Remote Sensing For Wind Energy
Remote Sensing For Wind Energy
29, 2022
Peña, Alfredo; Hasager, Charlotte Bay; Lange, Julia; Anger, Jan; Badger, Merete; Bingöl, Ferhat;
Bischoff, Oliver; Cariou, Jean-Pierre ; Dunne, Fiona ; Emeis, Stefan
Total number of authors:
27
Publication date:
2013
Document Version
Publisher's PDF, also known as Version of record
Citation (APA):
Peña, A., Hasager, C. B., Lange, J., Anger, J., Badger, M., Bingöl, F., Bischoff, O., Cariou, J-P., Dunne, F.,
Emeis, S., Harris, M., Hofsäss, M., Karagali, I., Laks, J., Larsen, S. E., Mann, J., Mikkelsen, T., Pao, L. Y., Pitter,
M., ... Würth, I. (2013). Remote Sensing for Wind Energy. DTU Wind Energy. DTU Wind Energy E No. 0029(EN)
General rights
Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright
owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.
Users may download and print one copy of any publication from the public portal for the purpose of private study or research.
You may not further distribute the material or use it for any profit-making activity or commercial gain
You may freely distribute the URL identifying the publication in the public portal
If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately
and investigate your claim.
Remote Sensing for Wind Energy
DTU Wind Energy
E-Report
Technical University
of Denmark
Frederiksborgvej 399
4000 Roskilde
Denmark
Tel. +4546775024
[email protected]
www.vindenergi.dk
Contents
Page
3 Atmospheric turbulence 52
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
3.2 Definitions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
3.3 Rapid distortion theory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
3.3.1 RDT and surface layer turbulence . . . . . . . . . . . . . . . . . . . 56
3.3.2 Eddy lifetimes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
3.3.3 The uniform shear model . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4 Fitting spectra to observations . . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4.1 Uncertainties on spectra . . . . . . . . . . . . . . . . . . . . . . . . 58
3.4.2 Spectral fitting and prediction of coherences . . . . . . . . . . . . . . 59
3.5 Model spectra over the ocean and flat land . . . . . . . . . . . . . . . . . . 61
3.5.1 Code and textbook spectra . . . . . . . . . . . . . . . . . . . . . . . 62
3.6 Comparison with the spectral tensor model . . . . . . . . . . . . . . . . . . . 65
3.7 Wind field simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
Notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70
6 Remote sensing for the derivation of the mixing-layer height and detection of
low-level jets 122
6.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
6.2 Mixing-layer height . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 122
6.2.1 Acoustic detection methods (Sodar) . . . . . . . . . . . . . . . . . . 123
6.2.2 Optical detection methods . . . . . . . . . . . . . . . . . . . . . . . 127
6.2.3 RASS . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 130
6.2.4 Further techniques . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
6.2.5 Comparison of acoustic and optical MLH detection algorithms . . . . 134
6.3 Boundary-layer height . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 134
6.4 Low-level jets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
6.4.1 Formation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 135
6.4.2 Frequency and properties of low-level jets . . . . . . . . . . . . . . . 136
6.5 Summary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 138
Notation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 139
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 140
shear due to the proximity of the Earth’s surface. The wind speed at the ground is always
zero, both on and offshore.
Sodar is a RS methodology for measurements of the wind speed and direction aloft at
various heights in the atmosphere. Sodars are ground-based instruments that transmit a
sequence of short bursts of sound waves at audible frequencies (2000–4000 Hz) upward in
three different inclined directions into the atmosphere.
The sodar measurement technology was well established and in operational use for decades
by now, starting in the 1980’s where they served environmental protection issues and has
been extensively applied to atmospheric research for environmental protection air pollution
prediction measures well before the present burst in wind energy research and application.
In Germany for example, sodars have been commissioned on several nuclear installations to
replaced tall meteorological towers and serve now as operational monitoring devises of the
local wind speed, direction and atmospheric stability.
As the sound waves from a sodar propagate forward a small fraction of the transmitted
sound energy is scattered and reflected in all directions from temperature differences and
turbulence in the atmosphere. A very small fraction of this scattered energy reaches back into
the sodar’s detector, which in principle is a directional-sensitive microphone.
The height at which the wind speed is measured is usually determined by the time delay
in the backscatter from the transmitted pulse. Under standard atmospheric conditions with
sound propagation speed of about 340 m s−1 backscatter from a sodar measurement at 170
m height above the ground will reach back into the detector after 1 s delay time.
The wind speed component in the transmitted beam direction is subsequently determined
from the Doppler shift observed as frequency difference between the transmitted frequency
and the frequency of the received backscattered sound wave. By combining the measured
wind speed components obtained in this way from three differently inclined sound path direc-
tions, e.g. from one vertical and two inclined sound paths, the three-dimensional wind vector
including wind speed and direction and tilt can be measured by sodar from preset heights from
the ground and up to the limit determined by the sodar’s lowest acceptable Carrier-to-Noise
(C/N) ratios.
The above description is for a mono-static system, where transmitter and receivers are
co-located on the ground. But alternative configurations, e.g. in the form of so-called bi-
static sodar configurations exist as well, where the transmitter and receivers are separated
e.g. 100–200 meters on the ground.
Bi-static configurations have significant C/N-ratio advantages over mono-static configu-
rations for wind energy applications. Received backscatter in a bi-static configuration is not
limited to direct (180◦ ) backscatter from temperature (density) fluctuations only, but enables
also backscatter contributions from the atmospheric turbulence. And the higher the wind
speed the more turbulence.
As a consequence significant improvements of the C/N- ratios can be obtained from a so-
called “bi-static configuration”, in which the transmitter and the receiver are separated from
one another on the ground. This becomes in particular relevant during strong wind situations,
where the background noise level increases with the wind speed.
A particular configuration considered for wind energy applications is therefore the bi-static
continuous wave (CW) sodar configuration. Alternatively to the range gating in a pulsed
system, the range to the wind speed measurement in a CW system can be determined by
well-defined overlapping transmission and receiving antenna functions. At Risø DTU we have
build and investigated such a sodar system for wind energy applications.
Pros Cons
Portable Low duty cycle (1 pulse transmitted every 3 s,
and up to 6–10 s lapse times before all
three wind components have been sampled)
Build on well developed and well-proven Limited by low S/N- ratio at:
audio-frequency “low tech” technology 1) high wind speed conditions
2) during neutrally stratified conditions
Sodars are relatively cheap (priced down Prone to solid reflections from the
to some 25% of a corresponding wind lidar) surroundings (including wind turbines)
Low power consumption (one solar powered Prone to high background environmental noise
version uses less than 10 W)
Sound backscatter: Relatively high yield Low wavelength/aperture ratio (1:10)
(backscattered power at the detector of the results in undefined broad antenna beams
order of 10−10 W) Prone to beam bending with wind speed of
the order of 5% or higher of the speed of sound
Huge measurement and sample volumes
Signal processing limited by pulsed sodars
relative long data acquisition times
(sampling time per pulse of the order of 1 s)
ment standards for the replacement of in-situ reference met masts. Work within the IEC
is at the moment aiming at establishment of a new international IEC-standard for remote
sensed wind measurements, as e.g. obtained by lidars, for power curve measurements.
2. Wind energy resource measurements – The global wind resources are now being mapped
globally on shore, off shore, over hilly and in mountainous terrain, etc. Here also, high
accuracy is of uttermost importance for accurate site and resource assessments. Mea-
surement errors in excess of 1% are unacceptable by project developers and investment
banks.
3. Wind turbine control – RS lidar instruments that are directly integrated into the wind
turbines hub or spinner or even into the blades are also seen as a forthcoming RS mea-
surement technology that can help improve the wind turbines power performance and
possibly also diminish fatigue wear from extreme gusts and wind shear via active steering
the wind turbines individual blade pitch or, to come one day maybe, its trailing edge
flaps.
Researchers at Risø DTU have during decades now followed and contributed to the devel-
opment of improved instrumentation for RS of wind. Starting out already in the 60’s with
more general boundary-layer meteorological investigations of flow and diffusion our present
research and experimental developments within the meteorology and test and measurement
programs at Risø DTU has recently become more and more directed towards applications
within wind energy. Wind lidars and lidar-based wind profilers, their measurement princi-
ples, their measurement performances, and also their possible future integration within wind
turbines themselves are here addressed.
A pulsed lidar on the other hand transmits a sequence of many short pulses, typical 30 m
in effective length, and then it detects the Doppler shift in the backscattered light from each
pulse as they propagate with the speed of light. While a CW lidar measures from one height at
a time a pulsed lidar measures wind speeds from several range-gated distances simultaneously,
typically at up to 10 range gates at a time.
The pulsed lidar’s spatial resolution, in contrast to the CW lidar, is independent of the
measurement range. The pulse width and the distance the pulse travels while the lidar samples
the detected backscatter control its resolution. The spatial resolution in the beam direction
obtainable with the 1.5 µm wavelength pulsed lidar in the market today are of the order of
30–40 m.
In addition, while a CW lidar’s upper measurement distance is limited progressing uncon-
fined measurement volume at long distances, a pulsed lidar’s maximum measurement range
is limited by deteriorating C/N-ratios in measurements from far distances (height).
Moreover, while a CW lidar equipped with a 1 W 1.5 µm eye-safe laser has been tested
able to sample and process up to 500 wind speed measurements per second, a corresponding
powered pulsed lidar can handle only 2–4 wind speed samples per second. Each of these
samples, however, then on the contrary contain wind speeds from up to 10 range gates
(ranges) measured simultaneously.
CW vs pulsed lidars Overall, CW lidar features high spatial resolution in the near range and
very fast data acquisition rates, features that are well suitable for turbulence measurements.
Today’s commercial available CW lidar profilers measure radial wind speeds at ranges up to
∼ 200 m and wind vectors at heights up to 150 m.
The pulsed lidar configuration on the other hand features lower but always constant spatial
resolution properties (30–40 m) at all ranges. They are also inherently slower in their data
acquisition rate, but then they measure wind speeds at multiple heights simultaneously, and
they hold also potential for reaching longer ranges (heights) than corresponding powered CW
lidars. At the test site in Høvsøre Risø DTU, commercial available pulsed wind lidar profiles
have regularly measured the wind vector profiles up to 300 m height.
Figure 7: Pulsed wind lidars (six WLS7 WindCubes) and one Galion (far back) during testing
at Høvsøre, Risø DTU
Notation
CW continuous wave
C/N carrier-to-noise
lidar light detection and ranging
RS remote sensing
SAR synthetic aperture radar
sodar sound detection and ranging
References
Bingöl F., Mann J., and Foussekis D. (2008) Lidar error estimation with WAsP engineering. IOP Conf. Series:
Earth and Environ. Sci. 1:012058
Mikkelsen T., Jørgensen H. E., and Kristiansen L. (2007) The Bi-static sodar “Heimdall” – you blow, I listen.
Risø-R-1424(EN), Roskilde
Mikkelsen T., Courtney M., Antoniou I., and Mann J. (2008) Wind scanner: A full-scale laser facility for
wind and turbulence measurements around large wind turbines. Proc. of the European Wind Energy Conf.,
Brussels
Wagner R. and Courtney M. (2009) Multi-MW wind turbine power curve measurements using remote sensing
instruments – the first Høvsøre campaign. Risø-R-1679(EN), Roskilde
Wagner R., Jørgensen H. E., Poulsen U., Madsen H. A., Larsen T., Antoniou I., and Thesberg L. (2008) Wind
characteristic measurements for large wind turbines power curve. Proc. of the European Wind Energy Conf.,
Brussels
Kristiansen L. (2010) My own perception of basic statistics. Available on request from Torben Mikkelsen,
Roskilde
2.1 Introduction
The atmospheric boundary layer, ABL, is the lower part of the atmosphere, where the atmo-
spheric variables change from their free atmosphere characteristics to the surface values. This
means that wind speed goes from the free wind aloft to zero at the ground, while scalars,
like temperature and humidity approach their surface value. An illustration of the profiles is
shown on fig. 8.
Characteristics of the atmospheric boundary layer, ABL, are of direct importance for much
human activity and well being, because humans basically live within the ABL, and most of
our activities take place here. The importance stems as well from the atmospheric energy
and water cycles. Because the fluxes of momentum, heat, and water vapour between the
atmosphere and the surfaces of the earth all pass through the ABL, being carried and modified
by mixing processes here. Since these mixing processes mostly owe their efficiency to the
mechanisms of boundary layer turbulence, a proper quantitative description of the turbulence
processes becomes essential for a satisfying description of the fluxes and the associated profiles
between the surfaces and the atmosphere.
Description of the structure of the flow, relevant scalar fields, turbulence and flux through
the atmospheric boundary layers necessitates that almost all types of the flows, that occur
there, must be considered. For these objectives, there are very few combinations of character-
istic boundary layer conditions that are not of significant importance, at least for some parts
of the globe.
Important processes for ABL produced turbulence is the production of variability from the
average velocity shear that has to exist in the ABL, as illustrated in fig. 8, because a fluid like
the atmosphere gases cannot be dynamical stable with a mean shear as shown in fig. 8, and
will start producing swirling motion, called eddies, see fig. 9. The characteristic spatial scale
is the height where it happens, and the boundary layer height h. The thus created variability
is called boundary layer turbulence, and is essential for the ABL mixing processes mentioned
in the introduction.
For the purpose of mathematical treatment, one separates the variable in mean values and
fluctuations like:
• For velocity components: ui =< ui > +u′i , < ui > in a mean value, u′i is fluctuating
turbulence., i = 1, 2, 3
• For scalars, T , and q, temperature and humidity: T =< T > +T ′ . q =< q > +q ′
• co-variances and turbulent transport: < u′i u′j >: Transport of ui in the j direction (and
vice versa).
• < ui T ′ >, and < ui q ′ > :Transport of temperature and water vapor in the ui direction.
Multiplying the velocity co-variances by the air density, ρ ,we can say that the velocity
transport can be considered a momentum transport, similarly multiplying the temperature
transport by ρCp is converted to a heat transport, and multiplying the water vapor transport
by ρL is converted to transport of latent heat. Here Cp is the heat capacity of the air at
constant pressure, L is the heat of evaporation for water. Indeed these terms are often used
in the description, since they reflect the physical importance better than the statistical term
correlation.
Alternatively <> can be denoted with capital letters or over-bars. The coordinate system
can be described at xi , i = 1, 2, 3 or with x, y, z , with the corresponding ui or u, v, w where
u now is along the mean wind direction, w is vertical and v lateral (the second horizontal
component).
A typical behavior of 600 seconds of ABL velocity components and temperature are shown
on fig. 10.
Figure 11: The ideal simplest atmospheric boundary layer, ABL, which still provides realistic
results. The conditions are statistically stationary and horizontally homogeneous. The wind
speed is forced to zero at the surface, and attain the free wind value Ua in the free atmosphere
above the ABL height, h. The values of wind, temperature and humidity are constant at the
surface and above h, giving rise to a vertical flux of momentum, heat and water vapor between
the surface and the h-level
The flow is assumed statistically stationary, meaning there will be variations. But these
will be statistically horizontally homogeneous and stationary. However, as seen from equation
1 we must allow a pressure gradient that, again somewhat unrealistic, is taken as constant.
Without a pressure gradient, there will be no wind.
g
θ = T + Γ · z, with Γ = (2)
Cp
With Γ being about 0.01 K/m. It is noted that the only difference between T and θ is the
linear height variation.
Equation 1 shows that the vertical fluxes of scalars are constant with height, while the
momentum fluxes are slightly more complicated. Focusing on the two first equations for
the horizontal velocity components, we define the geostrophic wind, G, from the pressure
gradient, and perpendicular to the direction of this gradient:
1 ∂p 1 ∂p 1 ∂p 1 ∂p
G = (U1G , U2G ) = (− , ) = (− , ) (3)
f ρ ∂x2 f ρ ∂x1 f ρ ∂y f ρ ∂x
The two first equations in eq. 1 now take the form:
∂ ′ ′
0 = f (u2 − U2G ) − ∂x3 (u1 u3 )
∂ ′ ′
0 = −f (u1 − U1G ) − ∂x3 (u2 u3 ) (4)
This equation shows that the wind velocity approaches the geostrophic wind at the top of
the boundary layer, where the turbulence disappears. Down through the boundary layer the
wind solution depends somewhat on the additional assumptions, but generally it undergoes a
spiraling motion as it reduces to zero at the ground, see fig. 12.
The detailed behavior of the wind speed through the ABL, is simplest for atmospheric
neutral stability, meaning that there is no heat flux and water vapor flux between the surface
and the free atmosphere (fig. 11). This can be assured by keeping the potential temperature,
∆Θh , difference at zero, meaning that T a−T s = Γh, compare eq. 2, and as well qs = qa . For
such situations one can derive the equations in eq. 7, using so called scaling laws, where the
momentum transport, is a new scale introduced from the co-variance, the so called friction
velocity u∗ as :
21
κG u∗ 2 2
ABL : = (ln( − A) + B
u∗ f z0
−B
ABL : α = tan−1 ( ) (6)
(ln( fuz∗0 − A))
u∗ z
SBL : U (z) = ln
κ z0
The first equations relate the conditions at the top of the ABL to the wind at lower heights,
the so called Surface Boundary Layer, SBL. The first of these called, the resistance laws of the
ABL, relating the friction velocity, u∗ , to the geostrophic wind, G, and a so called roughness
length, z0 , plus two three additional “universal” constants, the von Kármán, κ, and the two
constants A and B. The term fuz∗0 is denoted the surface Rossby number. Since f is of the
order of 10−4 A is about 2 and B about 5 , z0 is small. The Rossby number term dominates
the first of the resistance laws. The second of the resistance laws estimates the angle between
the geostrophic wind and the wind in the lower part of the ABL, the ASL. In this part the
wind profile is described by the last of the equations, the so called logarithmic profile.
The set of equations allow us to estimate the wind profile in the ASL for a given geostrophic
wind speed for varying roughness length. Therefore the roughness length is a very important
parameter, see fig. 14. In the principle it is a characteristic length, where the velocity extrap-
olates to its surface value, which is zero, but since it is a measure of the “roughness” of a
given landscape, much work has been done to establish consensus about the roughness of
characteristic real landscapes.
As seen from fig. 13, the roughness generally follows the intuitive images of what that
roughness is associated with. The larger and the sharper the protruding elements, the larger
the roughness. Although this image is simplistic, it still summarizes the main aspects of the
roughness characteristics of landscapes, including season variation of some landscapes. To
emphasize the importance of the roughness length for wind energy, we compare in fig. 14 the
winds in the ASL for for different roughness values and a given geostrophic wind.
Until now we have considered only situations, where the potential temperature differences
between the surface and the free atmosphere is unimportant for the dynamics, this state is
called thermally neutral. Changing T a and T s to make ∆θh different from zero, and sim-
ilarly for qs and qa , a heat flux and water vapor must flow between the top of the ABL
and the surface, and there will be a density gradient between top and bottom. For such a
situation the variation of wind speed, temperature and humidity can be described by an ex-
tension of the simple scaling used for neutral conditions in eq. 7. These scaling formulations
are normally denoted the Monin-Obukhov formulation. The set of scales is summarized below.
′ ′
Water vapor scale:q∗ = − −w
u∗
q
T u3∗
Monin-Obukhov stability length:L = κg(T∗ +0.61q∗ ) (7)
The water vapor concentration, q, enters into the stability measure, because both q and
Θ influence the density and thereby the stability based on density fluctuations. This will
be repeated throughout this text, but not always, because temperature is typically more
important than humidity for the stability.
The Monin-Obukhov scales are widely used, not only in the atmospheric surface layer (ASL)
that was introduced in fig. 12, but also for extension into the full ABL. In eq. 7, g is the
acceleration due to gravity, and the stability length, L, is a scale that is derived from the
energy budget as a measure of the importance of the heat and water vapor fluxes relative to
the momentum flux. As in eq. 7, κ is the von Kármán constant (∼ 0.4). The scales defined in
eq. 7, change relatively little within the ASL, above this one typically uses the ASL measured
parameters.
Corresponding to the neutral wind profile in eq. 7, we can now formulate the set of stability
influenced profiles in eq. 9, as:
u∗ z z
u(z) = (ln − ψ( ))
k z0 L
θ0 z z
θ(z) − θ0 = (ln( ) − ψT ( )) (8)
k z0T L
q∗ z z
q(z) − q0 == (ln( ) − ψq ( ))
k z0q L
Stability influenced profiles of < u(z) >, < Θ(z) > and < q(z) > . Notice, that separate
“roughness length” parameters are introduced for the scalars. These are of the order of z0
for smooth surfaces, but typically of the order of 10 times less for rough surfaces.
From the definition of L, it is seen that Lz ∼ 0 for neutral conditions, with no scalar fluxes.
ψ(0) = 0 for all three ψ-functions, and hence the logarithmic profiles for all variable are
recovered for neutral. Close to the ground Lz ∼ 0 simply because z is close to zero. Hence
all profiles start as being logarithmic close to the ground. For stable conditions, we have:
u(z)
= (z/10)α (9)
u(z = 10)
where the relevant α will be function of height, stability and roughness, as seen by compar-
ing with eq. 9. The power law profile has the advance of being simple and that α is a direct
measure of the relative shear, as can be seen by differentiation of eq. 9.
For neutral conditions we have in eq. 7 not only a wind profile in the ASL, but also for the
whole ABL a resistance law relating the G to u∗ and the geostrophic wind angle. This can in
the principle be extended also to different stabilities and to the scalar variable, letting the A
and B constants, and corresponding parameters for scalars, be function of stability, typically
in terms of h/L (Zilitinkevich, 1972),(Zilitinkevich, 1975). However, the quality of the data
fit for these extensions are worse than for the simple neutral expression in eq. 7, and therefor
they are not much used.
For our ideal boundary layer discussion, we have simply fixed the boundary layer height,
just as we could fix the temperature difference between the surface and the air on the top of
the ABL. In the real world the ABL height is determined as the outer range of the boundary
layer turbulence. For neutral conditions one find that it must be proportional to u∗ /f , since
these two parameters are the only parameters available to characterize the ABL turbulence.
Indeed one finds that this fits the data moderately well with a coefficient equal to about 0.3.
For moderate stabilities, one can use an expression like h ∼ 0.7(u∗ L/f )0.5 , and for strongly
stable conditions h ∼ L, the Monin-Obukhov stability length.
For unstable and many stable situations, it is common to use an independent rate equation,
to determine h, at least over land, that h = h(t) is now derived from an equation like
dh/dt ∼ F (··). Especially for unstable situations, a very simple and successful equation
has been developed. Assuming an unstable ABL growing, in response to sunrise, against a
background stable potential temperature gradient, γ, using the in-stationary temperature
equation in eq. 1. The result is:
0.5 Zt
Q(t)
h(t) = 2 with Q(t) = w′ θ′ |0 (1 + 2A)dt′ (10)
γ
0
σx z h
≈ F ( , ) with x = ui , θ, q (11)
x∗ L L
Where it has been found necessary to include h/L even within the ASL. In general the
function F increase with increasing instability, and become constant for neutral-stable condi-
tion. For velocity one often considers eq. 11 in term of the Turbulence Intensity (T I) derived
from eq. 11 by division with the local wind speed. Thereby the T I becomes an expression
for the likely relative deviation from the mean speed one will encounter for given situations,
Also distributions of the vertical instantaneous shear has become an important parameter,
e.g. distributions of α in eq. 9. Inserting the respective expressions we get, for the two:
σu κF (z/L,n/L) 1 z
TI = u = ln(z/z0 )−ψ(z/L) ≈ ln(z/z0 ) for L →0
z du φ(z/L) 1 z
α = u dz = ln(z/z0 )−ψ(z/L) ≈ ln(z/z0 ) for L →0 (12)
where ϕ(z/L) is the derivative of ψ(z/L). The behavior of data on the two functions are
shown in fig. 14.
Additional information about the fluctuations can be seen in the correlations or the spectral
structure of the signals. Atmospheric signals vary as function of both time and space. For
boundary layer turbulence, one can mostly assume that the fluctuations vary with space only,
and that measured time variation at a stationary measuring station is due to advection of a
spatial variation of the signal.
This, surprisingly simple assumption, is denoted the Taylor hypothesis. Take a measured
u(t) as an example:
Z
SXY (ω) = X ′ (t)Y ′ (t + τ ) exp(iωτ )dτ
Z
1
X ′ (t)Y ′ (t + τ ) = SXY (ω) exp(−iωτ )dω
2π
(14)
Where X and Y are two turbulent variables as function of time, with a correlation as
function of a lag time, τ , with a corresponding cross spectrum SXY (ω) of frequency ω.
Similar expression could be formulated for the spatial correlation of < X ′ Y ′ (xi + δi ) > and
its corresponding wave number cross spectrum, SXY (ki ). Here the wave number analysis is
different from the frequency analysis, in that the spatial lags, δl , and the wave numbers, ki ,
are vectors. Taylors hypothesis is relating frequency to k1 only, the wave number one along
the mean velocity direction.
The spectra are in fig. 16 are scaled with the same scales as the profiles < u′ w′ >,
< w′ T ′ >, etc. Thet are plotted versus a normalized frequency n = f z/ < u >= k1 z, with
f in Hz, and, where we have again used Taylors hypothesis for the frequency-wave number
relation, the z appears from the Monin-Obukhov similarity (Kaimal et al., 1972).
The spectra in fig. 16 are empirical and other spectral expressions exist, but there is a
general consensus on their form and intensity, such that the different forms agree broadly
on the behavior of the spectra, although there are low frequency differences, that can be
important in connection with some load modeling on structures (Andersen et al., 2010). The
spectra vary systematically with height through the ABL and with stability in widely accepted
ways. Olesen et al. (1984), Hojstrup (1982), Mann (1998), and the Mann-lecture at this
course.
At present, there is a strong activity of extending profiles and turbulence expressions all the
way through the ABL, and even further up. This is both because the growing wind turbines
makes the information important and also because the breakthrough in the remote sensing
technology, has made systematic data gathering in these heights possible. Some progress
has been made, as indicated by the illustration in fig. 17. However, it is still a speciality in
development, because the profiles aloft become sensitive to many new important scales, like
the ABL height, but also on the detailed characteristics of the entrainment zone above the
ABL.
An often used scale for strongly unstable conditions, is a velocity scale constructed from the
turbulent heat flux, to be used when the heat flux is more important than the momentum flux,
and often in the upper part of the unstable ABL: w∗ = (h < w′ θ′ > g/T ) ∼ u∗ (−h/κL)1/3 ,
u∗ will typically at most be a few tens of cm per second, while w∗ can reach several meter per
seconds. Additionally, to the new scales entering the problem,the demands to homogeneity
becomes more severe, and several aspects of remote instationarity and inhomogeneity, not
really influencing the profile in the ASL, will influence wind speed and wind turning aloft,
like e.g baroclinity and remote changes in surface characteristics. We shall consider these
aspects later. Here we present experimental and theoretical effort to describe the wind profile
to greater heights.
We finally return to the concept of thermal stability, which as we have seen is a measure of
the importance of the atmosphere’s thermal stratification relative to the wind shear dynamics
for the flow structure. In the SBL it can be described by the stability length L (Eqs. 7 and
9). Commonly used are the so called Richardson numbers, that measures the ratio between
the potential and the kinetic energy across a layer. In the ABL it is mostly defined in terms
of gradients of the mean temperature and the mean wind as shown in the following equation
eq. 17.
p ∂η ∂u
−w′ u′ z = u2∗z = ρ ∂x +ν ∂z
η η
∂θ
−w′ θ′ z = νθ ∂z (16)
η
The real irregular surface as indicated on the figure and denoted, η, while ν and ν0 are the
molecular diffusivity of momentum and heat respectively. On the left hand side of eq. 16 the
transport taking place in the turbulent ASL have been indexed with a z. In the viscous layer
ASL turbulence cannot exist due to the nearness of the surface, and the transport is taken
over by molecular gradient diffusivity and for momentum also by the pressure perturbation
around the irregularities. The correlation obviously becomes larger with steeper and sharper
irregularities, roughness elements. If there is no roughness element the transport is limited by
the molecular diffusion, and a scale analysis indicates the z0 ∼ δ ∼ ν/u∗ . With roughness
elements present the roughness becomes larger. On the other hand the heat transport in eq. 16
includes no pressure term, since pressure appears in the momentum equation only. Therefore
the transport of heat in the viscous sublayer is basically due to the molecular diffusion only,
and z0T , and correspondingly z0q , must be expected to be smaller than z0 , and more so, the
more rougher the surface. As mentioned earlier typically z0T /z0 ∼ 0.1 for rough surfaces.
Next we turn to the surface temperature θ0 appearing in eq. 9 and implicit in all our
arguments about the effects of stability. The surface temperature is a boundary condition in
eq. 9 but it is as well a result of the energy balance at the surface and controlled by processes
above, on and below the surface. At the surface the energy balance per unit area can be
written as (Bodyko, 1974):
dθ0
C = NR + H + E + G (17)
dt
Where, C · dθ0 /dt the heat accumulation in a thin layer at the surface, with C being the
heat capacity of this surface layer.
N R is the net radiation at the surface, consisting mainly of the incoming short wave radia-
tion from the sun and outgoing long wave radiation from the surface itself. N R is positive day
time, most positive around noon, and it is negative during night, when the surface undergoes
additional cooling. Cloud cover, season and latitude is obviously important here.
H is the turbulent heat flux = ρCp < w′ θ′ > , cooling the surface when it is warmer
than the air, and heating it when it is colder, typically due to radiation at night and day.
Obviously it can reflect also changes in the air temperature due to advection, from eq. 9 we
have H == ρCp < w′ θ′ >∼ −ρCp (θ(z) − θ0 )κu∗ / ln(z/z0T ).
E is the latent heat flux=ρL < wq ′ >, reflecting that the surface can regulate heat by
evaporating water or condensing water vapor. Just as for temperature these processes can be
driven by radiation, but also by advection.
G is the ground layer heat flux, which as all the other fluxes in eq. 17 can be both negative
and positive. But it will typically be directed downwards during day time, and upwards during
night.
Figure 20: Measurement and modeling of seasonal variation of roughness for different types
of vegetated surfaces (Hasager et al., 2003).
heated the surface water evaporates, it will start sinking; now being heavier, because it retains
the salt from the evaporated water. If the surface water cools, it also becomes heavier due to
the cooling and sinks. All this give rise to an intense mixing in typically the upper 10 meter
of the ocean. In the heat exchange with the atmosphere the water therefore constitute a very
large heat reservoir that only can change its temperature slowly, and additionally has its own
heating and cooling from the ocean currents. Indeed when an air mass moves over an ocean it
always ends up at the temperature of the ocean. For these reasons the homogeneous marine
ABL is always close to neutral. The diurnal radiation cycle shows very little influence on the
water surface temperature, although it can be measured, but typical amplitudes are less than
a few tenths of a degree (Peña et al., 2008). The annual radiation cycle on the other hand
has significant influence on the sea temperature, because they involve enough heat to change
both the temperature and the depth of the mixed layer. However, stable and unstable condi-
tions happens over the ocean as well on shorter timescales, but they are mostly transitional,
associated with air masses moving across water surfaces with a different temperature, either
coming from a nearby land or associated with moving weather systems. We shall return to
these phenomena when coming to the inhomogenous and instationary ABLs.
The sea is also an obvious source of water vapor evaporation, indeed over the ocean since
the q0 is derived from the saturated pressure at the surface temperature. The ocean is also
a source of liquid water in the form of sea spray converting to marine aerosols. In wintertime
ν u2
z0 = 0.11 + β(c/u∗ , −) ∗ (19)
u∗ g
Where the first term reflects the molecular diffusion limit, discussed above, when only few
roughness elements are present. The coefficient, β, denotes the Charnock’s constant. It is a
function of the phase speed c of the dominant waves and u∗ . Since the roughness elements
will be mowing with the phase speed of the waves in the direction of the wind. The term
c/u∗ is denoted wave age, because c increases with the duration of the acting wind. β is
varying between 0.01 and 0.07, being smallest for mid-ocean mature waves with large phase
speed. A “typical” value is 0.015. β can be a function of other parameters as well: e.g bottom
topography, swells, and very high winds (e.g. hurricanes) resulting in foam covered waters
that reduces β further (Makin, 2005).
Figure 21: The wind profile close to the water surface, with the wave induced vorticity and
the small scale roughness element riding on the larger scale waves, with a phase speed c.
In spite of the functions shown in eq. 19 the roughness of the sea surface still remains one
of the smallest, one can encounter in nature. This means that high wind speeds will be less
efficient in forcing the stability towards neutral over water than over land, although over water
frequency of neutral stability increases with wind speed. Still high winds can be encountered,
associated with strongly stable flows over water, again reflecting an inhomogeneous situation
where warm air is advected over cold water. Here the friction almost disappears. Again we
shall return to this issue. Just as winter snow can modify the roughness of a land surface
strongly, the winter will some part of the world cover the water with ice, see again fig. 13,
and the roughness now will depend on the characteristics of the ice surface, ranging from
extremely low for smooth solid ice, to quite rough for pack ice. The small z0 also means that
the turbulence typically is lower over the water than over land, reflected also in a lower ABL
height over water than over land. The small z0 also means that the z0T and z0q are close
to z0 for low wind speeds, and start deviating only for rough pack ice or larger wind speeds,
with rough sea.
Starting with an abrupt roughness change as depicted in fig. 23. Here the wind blows from
a surface with one roughness over a surface with another roughness. The stability is neutral.
The turbulence associated with the new roughness growths into the upstream boundary layer
by diffusion. This new boundary layer is called and Internal Boundary Layer, IBL, because it
growths within the already existing boundary layer.
Figure 23: Description of the structure of the internal boundary, of height h(x), due to a step
change in terrain roughness. Rule of thumb h ∼ 0.1x
dh ∂h ∂h κA
= u(h) = Au∗0 → = → Cx h ln(h/z0 ) (20)
dt ∂x ∂x ln(h/z0 )
Where we have used that u(h) follows the logarithmic profile. C is found to be about 1.
Notice, we measure u∗ as u∗ at the surface, because u∗2 must be expect to vary with height
in the IBL. The system works for both a smooth- to -rough transition, and the opposite, as
long as one uses the largest roughness of the two in the IBL growth equation eq. 18. With
h(x) determined in eq. 18, we can find the new surface, the new u∗2 , matching the upstream
and downstream profile at h to yield:
Equation 21 provides a very successful estimate of the ratio between the upstream and
downstream u∗ values with fetch, x, since h = h(x). Additionally it provides us with a way
of characterizing magnitude of roughness changes, through the factor M (in the last term).
From eq. 13 we get that the roughness change between a water surface and a smooth land
surface with a roughness of 1 cm has the same roughness change magnitude as the roughness
change between the same land surface and a city or a major forest. Both roughness changes
are associated with a roughness ratio of about 1000. We may also use the results to estimate,
how wide a homogenous area that is needed for the assumption about homogeneity, need to
be. This is height dependent. The growth rate of the IBL given by eq. 20 is slightly slower than
h = 0.1x, meaning the measuring height should be less than 10 times the fetch, x, to feel
the new IBL. It can show that the height below, which the flow is in approximate equilibrium
with the new surface, is about 0.01x, meaning that we will need a homogeneous new fetch
of 100 times the measuring height to consider the upstream conditions homogeneous. Under
all circumstances the demands to a homogenous fetch, depends on the measuring height.
This is one explanation of why it has been found to be more difficult to obtain consistent
estimates of the profiles aloft. For example the profiles obtained up to 200 m in fig. 14,
demands homogeneous fetches of 20 km to be in equilibrium with the underlaying surface.
Uncritically, we could extend the theory to a full ABL, find as a rough estimate we should
have a fetch of 100h to have a new equilibrium situation.
From eq. 21 one can estimate also the ratio between the upstream and down stream wind
speeds. Unfortunately, the formulation is wrong for very large fetches. Where it predicts the
u∗ -ratio, it becomes unity, which is not in accordance with the resistance law for a new ABL
with the new z0 . Also stability effects are not included. Both deficiencies can be repaired, but
the results miss the appealing simplicity of eq. 21, and will not be treated here.
Instead we turn towards the flow over low hills: As the flow approach the hill a pressure
perturbation develops, reaching both up-stream and down-stream and breaking the air on the
front slope and back slope of the hill and and developing an acceleration on the top of the
hill.
Figure 24: The principles for flow over a low hill according to (Jackson et al., 1975) from
(Hunt et al., 1982), with the width, L, the height, H, defined.
The flow over the hill is described for three regions in the vertical. At a height L the pressure
perturbation has disappeared, and the flow is undistorted. In the inner region the maximum
flow perturbation takes place, while the wind is still forced to zero at the ground. Therefore,
this is the region with the largest wind shear. The maximum wind perturbation takes place at
the top of the inner layer, and gradually decreases to zero at the height, L. The hight of the
∆u H
≈2 (22)
u0 L
Troen & Petersen (1989) have developed the theory into a Fourier form, allowing one to
build an arbitrary landscape as a superposition of low hills, which when combined with the
roughness changes, becomes the model in the WAsP system Troen & Petersen (1989). If the
slope becomes too large, of the order of 20%, flow separation develops on both the upstream
and the downstream slope.
We next turn towards the situations which are driven by a changing heat flux. Here we
start with the diurnal cycle for a horizontally homogeneous area at midlatitude, where the
boundary layer show obvious diurnal changes. As we have discussed the changes, they will
be smaller if the wind is high, and the atmosphere is cloudy, than when the wind is low and
the sky is clear. But it will always be present, and limiting the degree of stationarity one can
expect for such a surface. A similar changing rate can be found in frontal passages, which
take a time of the order of one day or longer.
Figure 25: A heat flux driven diurnal change of the ABL for a homogeneous midlatitude land
surface (Stull, 1988).
We next turn towards an IBL controlled by the surface heat flux, where the IBL, here
denoted TIBL, growths against a stable upstream ABL. Here, we can use the growth of the
unstable boundary layer equation in eq. 10, also depicted in fig. 25. Just assuming that the
upward heat flux now happens as function of fetch, x, instead of time, inserting t = u/x
in eq. 10. A more comprehensive description of these types of models are found in Gryning
et al. (1990). The situation can happen for both land and water surfaces, being dependent
on differences in surface heat flux and/or surface temperature of the two surfaces involved.
However, it is quite easy to imagine, that it is developing at a coast line. Here the surface
temperature of the land and the sea is quite often different. Due to the two different surfaces
response in the surface energy budget (eq. 17) on a diurnal scale. As the time and space scales
for such coastal system increase, the differential heating between two surfaces, will influence
the atmospheric dynamics, because air will tend to rise over the relatively warm surface and
sink over the colder area. The resulting circulation is called a land-sea breeze system. At
night the land cools relatively to the sea, and the flow reverses. The land-sea breeze is the
best known one of the breeze systems. However, breeze systems occur also as breeze systems
around major cities, because of the “urban heat island” effect. When the spatial scale of a
land-breeze system increase to continental scale, one talk about Monsoons, which will typical
be seasonal rather than diurnal. However breeze systems of many scales can exist at the same
location, and will often interact.
A more complex situation appears for a transition from an unstable ABL passing over a
colder surface, where the more intensive unstable ABL has to decay before the stable IBL can
establish itself. Also, these situations occur quite often in some coastal areas.
Especially the transition from relatively warm unstable/neutral land ABL into a stable
IBL over the water (Melas, 1998). Several things happens. The air gradually looses heat to
the water and reach the water surface temperature at the interface, where a small neutral
layer establishes itself. This neutral layer slowly growths up through the stable IBL, until the
original temperature gradient between the two air masses is concentrated at the top of a new
boundary layer that however can take many hundred kilometers to establish itself.
Figure 27: Warm neutral air flowing over cold water, and develops a stable IBL, SIBL, gradually
transitioning to a neutral ABL with an inversion on the top. Over the water for off-winds.
(Lange et al., 2004). The broken triangular line indicates the possibility for a low level jet at
the transition.
A low level jet may happen in the transition between land and water, because the sur-
face friction in the boundary layer suddenly disappears, for the air crossing the coastal line.
Considering eq. 1, we see that the wind above the boundary layer is unchanged,friction is
unimportant here. The velocity within the ABL must accelerate because the friction is re-
duced, finally the wind has to go to zero at the surface. It is a transitional phenomenon that
gradually disappears as the full profile and stress profile reasserts themselves downstream from
the coast line. These low level jets can appears, wherever the surface stress reduces because of
increasing stability or decreasing roughness or both. It is quite well known over the nighttime
Plains of US, where they are a result of the growing night time stability.
The heat/temperature controlled IBLs may happen also especially over water, when syn-
optic air masses are advected over relatively cold or warm water. Here especially, the second
one, will give rise to a thermally induced convective IBL, which over the water typically will
1 ∂p gz ∂Tl 1 ∂p gz ∂Tl
G = (U1G , U2G ) = (− − , − ) (23)
f ρ ∂y f ρ ∂y f ρ ∂x f ρ ∂x
In eq. 23 the original G is present at the surface only, and the deviations are seen to increase
with z. Also they will influence both the magnitude and the direction of the geostrophic wind
for baroclinic situations. Indeed a low level jet can result from this type of baroclinicity as well
(Garreaud et al., 2005). The subscript l refers to that the temperature must be considered a
layer average. The thermal contribution is often called “The Thermal Wind”.
We finally mention coast line changes that firstly happens in the water, but may still
influence the air above the water: For deep water coasts, currents and wind may carry the
surface water outwards, and upwelling of deeper colder water replaces it closer to the shore.
This means that the flow approaching the land is in the process of cooling and developing a
stable internal boundary layer over the coastal water, as in fig. 27, before making landfall. A
well known example is the Californian west coast upwelling sea- breeze system. For shallow
coasts, on the other hand the coastal water reflects the land surface temperature more,
because the radiation balance reaches all the way to the bottom. Additional wave breaking
may also enhance the mixing through the water of both momentum and heat (Johnson et
al., 1999).
Figure 29: Different reanalysis sets, with start time, and horizontal grid size and time resolu-
tion. The vertical resolution corresponds to 20-30 levels.
Stability (temperature and humidity) and other possible data are here again “Nice to know”,
for the same reasons as for the wind resource estimation.
Finally we mention that it is quite normal that wind farm operators keep a meteorology
mast running when the farm has started to operate, both to support the monitoring of the
farm performance, and to support the service and maintenance function.
2.10 Summary
We have all most entirely avoided discussions of the different measuring systems in this pre-
sentation and concentrated on meteorological climatological aspects.
We have summarized knowledge about atmospheric boundary layers, starting with the
simplest ideal form, being statistically stationary and horizontally homogenous. We have
thereafter introduced aspect of reality into the picture, considering stability and unstability,
typical terrestrial boundary layers and marine boundary layers, and the important points of
instationarity and inhomogeneity.
Over land we have seen that stability changes have a clear diurnal cycle, the importance
of which becomes the smaller for the boundary layer structure, the higher the wind speed
and the larger the cloud cover. On the other hand the thermal properties become more and
more important the greater the height within the boundary layer. Over water stability is more
associated with advection of air masses than daily variation of insolation. On the other hand
high wind does not force the stability towards neutral to the same degree as for the terrestrial
boundary layers.
When the surface change, the boundary layer reacts by forming an internal boundary layer.
Orography changes, induce a pressure perturbation that reaches out in all directions within
a range of the same scale as the width the terrain change. The pressure perturbation also
perturbs the oncoming flow. Changes due to roughness and thermal changes on the other
hand diffuse into the flow as it moves across the surface. The transition zone for these changes
will typically be between 10 and 50 kilometers before a new boundary layer establishes itself
as a homogeneous boundary layer. The highest levels in the boundary layer reach equilibrium
the latest. The stable internal boundary layer over water for offshore flow takes the longest
fetch to reach equilibrium up to more than 100 km, meaning for example that an enclosed
sea as the Baltic Sea, with relatively cold water, must be considered entirely coastal.
We find that for complex terrain some general methodologies are possible, but one must
expect to have into invest more specific studies for projects within such areas, before having
Notation
A Coefficient in the resistance laws (6)
an entrainments coefficient in (10), area (15)
coefficient of proportionality in (17)
ABL, ASL Atmospheric Boundary Layer, Atmospheric Surface Layer
α Power in power law wind profile (9), Geostrophic angle (6)
B Another coefficient in the resistance laws (6)
β Charnock’s coefficient as function of wave age in (16)
C Surface heat capacity in (14), coefficient of proportionality in (17)
Cp The heat capacity at constant pressure for air
c Phase velocity of surface water waves (16)
δ Length scale lag-length, height of interfacial layer fig. 17
E Latent heat flux (14)
es : es = es (T ) saturated water vapor pressure in the atmosphere
η detailed surface height (13)
F unspecified function used in (11) and in deriving (10)
F Coriolis paramert (1)
G, U1G , U2G Geostrophic wind (3)
Soil energy flux (14)
g acceleration due to gravity (1)
H Height of hill (Fig 8), turbulent heat flux (14)
h height of ABL, and height of IBL (Figs. 11 and 20)
height of roughness element (15)
IBL Internal Boundary Layer
κ von Kármán constant in Monin-Obukhov turbulence description (6)
< psi(z/L) Stability function for profiles (8)
L Monin-Obukhov stability length scale (7)
width of the hill (fig. 21)
heat of evaporation
Γ The dry adiabatic lapse rate (2)
NR Net radiation (14)
ν, νθ Molecular viscosity and heat conductivity (13)
ρ density of air (1)
p air pressure (1)
q water vapor mixing ratio defined as density of water vapour/air density,
surface value of q (1)
q∗ Turbulence scale for q (7)
S Cross wind area of roughness element (15)
σ Standard deviation
T Temperature in K (1)
T∗ Turbulence scale for T or θ (7)
t time (1)
θ, θ0 Potential temperature in K, surface value of θ , (2)
Ui = u1 , u2 , u3 = u, v, w Three components of the wind velocity
Xi = x1 , x2 , x3 , = x, y, z Three spatial coordinates
z0 , z0T , z0q Roughness length for u, T and q
Bodyko, M.I. (1974) Climate and Life. International Geophysical Series, Vol. 18, Academic Press, New York
18:508
Bredmose, H., S.E.Larsen, D. Matha, A. Rettenmeier, E. Marino, L. Saettran (2012) D2.04: Collation of
offshore Wind Wave Dynamics. MARINET Report, EU-FP7 Grant no. 262552, 50p
Brutsaert, W.H. (1982) Exchange processes at the earth atmosphere interface Engineering Meteorology, Ed.
E. Platte, Elsevier, 319-369
Charnock H. (1955) Wind stress over a water surface. Quart. J. Roy. Meteorol. Soc. 81:639–640
Design Standards, (2007-2011) IEC 61400-3: Wind turbines- Part 3: Design of offshore wind turbines,
2009 DNV-OJ-J101. OFF SHORE STANDARD: Design of offshore wind turbine structures. DNV, 2007.
ABS(American Bureau of Shipping), Offshore Wind Turbine Installation, American Bureau of Shipping,
2010 ABS(American Bureau of Shipping), Design Standards for offshore wind farms, American Bureau of
Shipping, 2011 IEC 61400-1, Design requirements, 2008. IEC 61400-1, Design requirements-amendment1,
2010
Garreaud, R. D. and R. Munoz, 2005. The low level jet off the subtropicalwest coast of South America,
Structure and Variability. Monthly Weather Review 133:2246-2261
Gryning, S.E. and E. Batchvarova, (1990) Analytical model for the growth of the coastal internal boundary
layer during onshore flow. Quart. J. Roy. Met. Soc. 116:187-203
Gryning S.-E., Batchvarova E., Brümmer B., Jørgensen H., and Larsen S. (2007) On the extension of the
wind profile over homogeneous terrain beyond the surface layer. Bound.-Layer Meteorol. 124:251–268
Hahmann, A. N., J. Lange, A. Peña and C. Hasager The NORSEWInD numerical wind atlas for the South
Baltic DTU Wind Energy-E-Report-0011(EN)
Hasager, C.B.; Nielsen, N.,W.; Jensen, N.O.; Bøgh, E.; Christensen, J.H.; Dellwik, E.; Søgaard, G.. (2003)
Effective roughness calculated from satellite-derived land cover maps and hedge-information used in a
weather forecasting model Boundary-Layer Meteorology 109:227-254
Hunt, J.C.R. and J.E.Simpson, (1982) Atmospheric Boundary Layers over Non-Homogeneous terrain Engi-
neering meteorology (E.: E. Plate) Elsevier, Amsterdam 269-318
Håkanson, L. (1991) The Baltic Sea Environment Session 1: Physical Oceanography of the Baltics. The Baltic
University Programme, Uppsala University, Swedem, 35 p
Højstrup, J., (1982) Velocity spectra in the unstable planetary boundary layer J. Atmos. Sci. 39:2239-2248
Jackson, P.S. and J.R.C. Hunt, (1975) Turbulent wind flow over a low hill J. Roy. Met. Soc. 101:929–955
Jones, I S F, Y Volkov, Y Toba, S Larsen, N Huang, M Donelan (2001) Overview. Chapter 1 Wind stress over
the ocean ( Ed. I.S.F.Jones and Y.Toba). Cambridge University Press, 1 - 31.
Johnson, H K, H.J Vested, H. Hersbach, J Højstrup and S E Larsen (1999) The coupling between wind and
waves in the WAM model Journ. Atmos. Ocean. Tech 16(11):1780 - 1790
Kaimal, J.C., J.C. Wyngaard, Y. Izumi, and O.R. Cote, (1972) Spectral characteristics of surface-layer turbu-
lence Quart. J. Roy. Met. Soc., 563–589
Karagali, I,, A. Peña, M. Badger and C.B. Hasager (2012) Wind characteristics of the North Sea and the
Baltic Seas from QuickScatt satellite Wind Energy
Lange, B., S.E.Larsen, J. Højstrup and R. Barthelmie, (2004) The influence of thermal effects on the wind
speed profile of the coastal marine boundary layer Boundary-Lay. Meteorol. 112:587-617
Larsen, S.E. and Jensen, N.O. (1983) Summary and Interpretation of Some Danish Climate Statistics Risø-
R-399, 76 pp
Larsen, S.E., N. E. Tarp-Johansen, S. Frandsen and E. R. Jørgensen (1990) Södra Midsjöbanken Environmental
Data-Initial analysis Risø-I-2505. FP6 Integrated project DOWNVInd
Lettau H. H. (1962) Theoretical wind spirals in the boundary layer of a barotropic atmosphere. Beitr. Phys.
Atmos. 35:195–212
Makin, V. K. (2005) A note on the drag of the sea surface at hurricane winds Boundary Layer Meteorology
115:169-175
Melas, D., (1998) The depth of the stably stratified internal boundary layer over the sea Geophysical Research
Letters 25(13):2261-2264
Mortensen, N. G. and E.L. Petersen (1997) Influence of topographical input data on the accuracy of wind
flow modeling in complex terrain European Wind Energy Conference, Dublin, Ireland
Olesen, H.R., S.E. Larsen, and J. Højstrup, (1984) Modelling velocity spectra in the lower part of the planetary
boundary layer Boundary-Layer Meteorol. 29:285–312
Peña A., Gryning S.-E., and Hasager C. B. (2008) Measurements and modelling of the wind speed profile in
the marine atmospheric boundary layer. Bound.-Layer Meteorol. 129:479–495
Peña A. and Gryning S.-E. (2008) Charnock’s roughness length model and non-dimensional wind profiles over
the sea. Bound.-Layer Meteorol. 128:191–203
Peña A., Mikkelsen T., Gryning S.-E., Hasager C. B., Hahmann A. N., Badger M., Karagali I., and Courtney
M. (2012) Offshore vertical wind shear. Final report on NORSEWInD’s work task 3.1 DTU Wind Energy-
E-Report-005(EN), DTU Wind Energy, Risø¸ Campus, Roskilde
Smedman, A-S., H. Bergström and B. Grisogono, (1997) Evolution of stable internal boundary layers over a
cold sea Journ. Geophys. Research 102:1091-1099
Stull R. B. (1988) An introduction to boundary layer meteorology, Kluwer Academic Publishers, 666 pp
Tennekes, H. and J.L. Lumley (1972) A first course in turbulence MIT, Cambridge, Mass. 300p.
Troen, I. and E.L. Petersen (1989) European Wind Atlas Risø National Laboratory, Roskilde, Denmark
Wang, H., R. Barthelmie, S.G. pryor, and H-G. Kim (2012) A new Turbulence model for offshore wind
standards Subm. Geophysical Researh
Zilitinkevich, S.S. (1972) On the determination of the height of the Ekman boundary layer Boundary-Layer
Meteorol. 3:141–145
Zilitinkevich, S.S. (1975) Resistance laws and prediction equations for the depth of the planetary boundary
layer J. Atmos. Sci. 32:741–752
Zilitinkevich, S. , I. Mamarella, A. Baklanov and S. Joffre (2009) The effect of Stratefication on the aerody-
namic roughness length Meteorological and Air quality models for urban areas. Springer Berlin
3.1 Introduction
For many civil engineering structures, including wind turbines, dynamic wind loading caused by
the atmospheric turbulence is a serious concern for the designer. Gust loading on streamlines
bridge decks requires knowledge of the vertical wind field fluctuations not only in one point, but
also how the fluctuations are correlated in space Simiu and Scanlan (1996); Larose and Mann
(1998). Also the horizontal components may be of importance in bridge aerodynamics. For
dynamical load calculations on a wind turbine, for example at an off-shore location knowledge
of all three wind components and their spatial correlations are needed because the gusts are
‘sampled’ in a complicated way by the sweeping blades. Yet other structures such as tension
leg platforms used for extracting oil on deep waters are sensitive to slow variation in the
direction of the wind. Thus various engineering structures are sensitive to various components
of wind fluctuations at a wide range of frequencies and also to the spatial correlations of these
fluctuations.
The spatial structure of turbulence is also important in order to understand how remote
sensing instruments such as lidars measure in a turbulent flow fields. That is because the lidar’s
sampling volume is rather extended and thus very far from the almost point-like measurements
of a ultra-sonic anemometer. The description of how lidars measure turbulence may be found
in Mann et al. (2009) for a pulsed lidar, or in Sjöholm et al. (2009) for a continuous wave
(cw) lidar.
The purpose of this contribution is to model the spectral tensor of neutral atmospheric
surface layer turbulence. The spectral tensor contains all information on spectra, cross-spectra
and coherences, which usually are the input requested by wind engineers. We also want to
devise a general algorithm to simulate three-dimensional fields of all three components of the
wind velocity fluctuations. Such simulations are particular useful for time domain simulations
of gust loading of wind turbines and other structures.
In section 3.3 rapid distortion theory (RDT) is used to estimate the tendency of shear to
make turbulence anisotropic. RDT is a linearization of the Navier–Stokes equations and has
as such limited applicability. The influence of the non-linearity is modeled by postulating some
limit as to how much shear is allowed to make the turbulence anisotropic. This modelling uses
the concept of eddy lifetime. Despite the various assumptions and postulates the tensor model
only contains three adjustable parameters: a length scale describing the size of the energy
containing eddies, a non-dimensional number used in the parametrization of eddy lifetime,
and the third parameter is a measure of the energy dissipation.
These three parameters are estimated by comparing the model to measurements over the
sea in section 3.4. In section 3.5 the model is compared to various widely used wind engineering
spectral formulations. Finally, in section 3.7 the spectral tensor is used in a numerical algorithm
to simulate three-dimensional fields of all three components of the wind vector. This is done
by recasting the Fourier representation of the wind field in the discrete wave-vector space,
i.e. as a trigonometric series, where the statistics of the random coefficients are determined
by the spectral tensor. The method is considerably simpler, faster and in some aspects more
physical than many other currently used simulation algorithms. The method is now used in
bridge aerodynamics and in load calculations on wind turbines.
Much of the material presented here has previously been reported in Mann (1994, 1998),
and more details on many aspects may be found in these papers. Newer comparison with
neutral atmospheric data taken from Risø’s test station Høvsøre may be found in Peña et al.
3.2 Definitions
The atmospheric turbulent velocity field is denoted by ũ(x ), where x = (x, y, z) is a right-
handed coordinate system with the x-axis in the direction of the mean wind field and z as
the vertical axis. The fluctuations around the mean wind, u (x ) = (u1 , u2 , u3 ) = (u, v, w) =
ũ(x ) − (U (z), 0, 0), are assumed to be homogeneous in space, which is often the case in the
horizontal directions but is only a crude approximation in the vertical. Since turbulence over
the sea at high wind speeds is primarily shear-generated, the mean wind field is allowed to
vary as a function of z. Because of homogeneity, the covariance tensor
Rij (r ) = hui (x )uj (x + r )i (24)
is only a function of the separation vector r (h i denotes ensemble averaging).
We shall use Taylor’s frozen turbulence hypothesis to interpret time series as ‘space series’
and to serve as a ‘dispersion relation’ between frequency and wave number (Panofsky and
Dutton, 1984). Therefore, we can suppress the time argument in u .
We only aim at modelling the second order statistics of turbulence, such as variances, cross-
spectra, etc. For simulation purposes the velocity field is otherwise assumed to be Gaussian
(see section 3.7). It is still not clear how much influence the statistics of third order, such as
skewness, has on load calculations.
All second order statistics can be derived from the covariance tensor or its Fourier transform,
the spectral tensor: Z
1
Φij (k ) = Rij (r ) exp(−ik · r )dr , (25)
(2π)3
R R∞ R∞ R∞
where dr ≡ −∞ −∞ −∞ dr1 dr2 dr3 . The spectral tensor is the basis of the Fourier
simulation in section 3.7.
The stochastic velocity field can be represented in terms of a generalized stochastic Fourier-
Stieltjes integral: Z
u (x ) = eik ·x dZ (k ), (26)
where the integration is over all wave number space. The orthogonal process Z is connected
to the spectral tensor by
hdZi∗ (k )dZj (k )i = Φij (k )dk1 dk2 dk3 , (27)
∗
which is valid for infinitely small dki and where denotes complex conjugation (Batchelor,
1953).
Is it very difficult to measure the spectral tensor directly. Instead cross-spectra, defined as
Z ∞
1
χij (k1 , ∆y, ∆z) = Rij (x, ∆y, ∆z)e−ik1 x dx (28)
2π −∞
are often measured, say by two instruments separated by ∆y in the horizontal direction
perpendicular to the wind and ∆z in the vertical, and are used in practical applications. The
connection between the components of the spectral tensor and the cross-spectra is
Z ∞Z ∞
χij (k1 , ∆y, ∆z) = Φij (k )ei(k2 ∆y+k3 ∆z) dk2 dk3 . (29)
−∞ −∞
When the two indices i and j are the same and ∆y = ∆z = 0 Eq. (29) becomes the one-point
spectrum Fi (k1 ) = χii (k1 , 0, 0). This definition implies that spectra are two-sided, i.e. we get
the variance by integrating from −∞ to ∞. This convention is used throughout this chapter.
To distinguish between spectra as functions of wave number k1 (= 2πf /U ) and frequency
f we use F for the former and S for the latter, i.e. Si (f )df = Fi (k)dk . The coherence is
defined as
|χij (k1 , ∆y, ∆z)|2
cohij (k1 , ∆y, ∆x) = , (30)
Fi (k1 )Fj (k1 )
which can be interpreted as a normalized cross-spectrum.
Figure 30: Interpretation of the interplay of shear and turbulence: Two differently oriented
eddies are followed over three successive times. Shear stretches (along the axis of rotation)
and speeds up the upper eddy while the lower eddy is compressed and slowed down.
Physically, the last term on the right hand side of Eq. (32) may be interpreted as the
stretching of vorticity by the mean shear (see Figure 30). The first term is a distortion of the
mean vorticity by velocity fluctuations.
In order to solve Eq. (32) we have to Fourier transform the equation. In order to do so, it
is important to notice that wave fronts are advected by the mean flow i.e.
dk
= −(∇ U )k . (33)
dt
The solution to this wave front advection equation is
k (t) = exp(−∇ U t)k0 (34)
where exp means the matrix exponential.
For a general linear U Eq. (32) does not have analytic solution. However, for many simple
situations such as unidirectional shear, non-rotational stretching or compression, etc. such
solutions exists (Townsend, 1980).
and " #
1
k2 k02 βk1 (k12 + k22 ) 2
C2 = 2 3 arctan . (46)
(k1 + k22 ) 2 k02 − k30 k1 β
Eqs. (38) and (43) give the temporal evolution of individual Fourier modes.
1 10
kL
Figure 31: Eddy lifetimes as functions of the magnitude of the wave vector. The lifetimes
given by Eq. (51) give the most realistic results.
Comte-Bellot and Corrsin (1971) give another lifetime model which has the right asymptotic
behaviour for k → ∞, the ‘coherence-destroying diffusion time’ :
Z ∞ − 21
τD (k) ∝ k −2 p−2 E(p)dp
k
− 21 2
− 32 4 17 7 k− 3 for k → ∞
∝ k 2 F1 , ; ; −(kL)−2 ∝ (52)
3 6 3 k −2 for k → 0
which was constructed as the square of the eddy size divided by a k-dependent ‘turbulent
viscosity’.
Further, the inverse ‘eddy-damping rate’
( 2
3
− 12 k− 3 for k → ∞
τE (k) ∝ k E(k) ∝ 7 (53)
k− 2 for k → 0
4 × 10−4
k1 Fw (k1 )
U2
2 × 10−4
Figure 32: Spectra of w from the Great Belt Coherence Experiment. Mean wind speeds are
between 16 and 20 m s−1 and directions are in a narrow interval around the South. Dashed
spectra have slightly unstable stratification, gray have stable, and the thin have neutral.
Fv
0.2
Fw
ℑ(χuw )
ℜ(χuw )
Figure 33: Average u-, v-, w-, and cross-spectra of all the neutral runs present in Figure 32.
The ragged curves are measurements while the smooth are the model spectra. The model
has zero imaginary part of the cross-spectrum (quadrature spectrum).
Figure 34: The mast array on Sprogø viewed from SSE. The tiny dots at the top of the masts
are the omni-directional sonic anemometers.
0.5
1.0
∆y = 32.5 m
coh(k1 )
0.5
1.0
∆y = 47.5 m
0.5
1 2 1 2 1 2
k1 ∆y
Figure 35: The dots are measured coherences from the same set of data as used for Figure
33 for various horizontal separations ∆y and for all three velocity components. The lines are
the coherences predicted by the model.
These parameters are in turn used to predict the coherences as shown in Figure 35. As seen
from this figure the predictions agree well with the measurements except for the w coherence,
especially at the largest separation.
1000CDN
1
0
10 20 30
U (10 m) [m s−1 ]
Figure 36: The neutral drag coefficient CDN as a function of mean wind speed at z = 10 m.
The broad line is from Charnock’s relation Eqs. (64) and (62). The thin lines are empirical
relations from Geernaert (1987) and the dotted line is from NDP (1998), see Eq. (73).
30 m Eq. (62) is a good approximation to Eq. (63). Throughout this comparison we use
f = 10−4 s−1 .
Charnock (1955) argued that over the sea the roughness length is related to g = 9.8 m s−2
the acceleration due to gravity and the friction velocity by
u2∗
z0 = A (64)
g
where A, the Charnock constant, must be determined experimentally. On basis of an extensive
literature study of ocean data Garratt (1977) found that the best fit of Eq. (64) is A = 0.0144.
A slightly newer value is given by ESDU International (1982):
A = 0.0167, (65)
which will be used here. Over the ocean the neutral drag coefficient
2
u∗
CDN = (66)
U (10 m)
increases monotonicly with U as can be seen by solving Eqs. (64) and (62). This is shown
in Figure 36 as a broad line together with several recent empirical relations. The figure gives
a good impression of the uncertainty in estimates of drag coefficients. Among the various
reasons for this variability are atmospheric stability, surface currents, ‘wave age’, length of
the fetch over water, and water depth (Garratt, 1977; Geernaert, 1987; Brown and Swail,
1991). The spectral density of velocity fluctuations is in general proportional to the drag
coefficient so the uncertainty of the former is probably of the same order of the latter.
5
ESDU
u Kaimal
Simiu
Højstrup
NPD
1
0.5
0.1
5
v ESDU
Kaimal
k1 F (k1 ) [m2 s−2 ]
Simiu
0.5
0.1
5
w ESDU
Kaimal
Simiu
0.5
0.1
0.0001 0.001 0.01 0.1 1
−1
Wave number k1 [m ]
Figure 37: Comparison of spectral models. For the comparison z = 40 m and U = 40 m s−1
(over the sea) is chosen. For u (ESDU International, 1985), Eqs. (67), (70), (80), (76) are
used. For v and w (ESDU International, 1985), Eqs. (68) and (71), and (ESDU International,
1985), Eqs. (69) and (72), respectively. Eq. (63) together with Eq. (64) gives u∗ = 1.78 m
s−1 and z0 = 0.0054 m.
The spectra of Simiu and Scanlan (1996) have the same functional shapes as Kaimal’s but
the numerical constants are different:
f Su (f ) 100n
= , (70)
u2∗ (1 + 50n)5/3
f Sv (f ) 7.5n
2
= , (71)
u∗ (1 + 9.5n)5/3
and
f Sw (f ) 1.68n
= . (72)
u2∗ 1 + 10n5/3
Deviations from surface layer scaling are found in the model spectra from ESDU Interna-
tional (1985). Also the spectra of Norwegian Petroleum Directorate NDP (1998) and Højstrup,
u
0.5
v
w
0.2
k1 F (k1 )/u2∗
0.1
0.05
0.02
0.1 1 10 100
k1 z
Figure 38: The ‘sheared spectral tensor’ of section 3.3.1 (curves with dots) fitted to the
models by Simiu and Scanlan Eqs. (70)–(72). The result is given by Eq. (83).
The u-spectrum of NDP (1998) applies to winds over oceans and assumes the drag coef-
ficient to be
CDN = 0.525 × 10−3 (1 + 0.15U10), (73)
√
see Figure 36. Integrating dU /dz = u∗ /(κz) = CDN U10 /(κz) Eq. (73) implies that
z
U (z) = U10 1 + C ln (74)
10 m
with
C = 0.0573(1 + 0.15U10 )1/2 (75)
where U10 has to be measured in meters per second. While discussing the NPD spectrum we
also assume the unit of z to be meter, f is Hz and Su is m2 s−2 Hz−1 . The spectral density
of the longitudinal wind component is
2
U10 z 0.45
160
10 10
Su (f ) = 3n
5 (76)
1 + f˜n
with
z 2/3 U −3/4
10
f˜ = 172f (77)
10 10
where the ‘neutral length scale’ A = 3000 m and nt = f A/U . The second term in the
parenthesis is the Kaimal spectrum Eq. (67).
All spectral models are compared in Figure 37 for a specific choice of U and z. Generally,
ESDU has larger length scales compared to those by Kaimal and by Simiu & Scanlan, which
are similar. NPD and Højstrup support ESDU’s large u-scale. ESDU, though, has the most
peaked spectra and, at high wave numbers, slightly lower spectral densities. All spectra agree
fairly well at high wave numbers but have substantial scatter at low wave numbers.
u
2
v
1 w
k1 F (k1 ) [m2 s−2 ]
0.2
0.1
Figure 39: Example with z = 40 m and U = 40 m s−1 of the fit of the spectral tensor model
(curves with dots) to the ESDU models.
Here we fit the spectral tensor of section 3.3.1 to models that describe all three component
spectra, namely the ones by Kaimal, Simiu & Scanlan and ESDU.
We obtain the parameters Γ , L and αε2/3 by making a simultaneous least squares fit to
the u-, v- and w-model spectra for wave numbers in the range 0.05 < k1 L < 100. For the
3 Højstrup, Larsen and Madsen (1990) also gives a model for the v spectrum, but it was never compared
100
40
z [m]
25
50
L [m]
30 15
20
10
15 20
25
10 30
40
10 50
U [m/s] 60
70
5
80
0.01 0.1 1 10
αε2/3 [m4/3 s−2 ]
5
U
60 70
80
50
4 40
30
3 25
Γ
20
2
10 m/s
1 15
Figure 40: The parameters of the spectral tensor model derived from fits to the ESDU model
spectra for turbulence over the sea. Given U and z, all three parameters can be extracted
from these plots.
It is more complicated to get the parameters from the ESDU models because the spectra no
longer depend on U and z in a simple way. For each set {U, z}, a fit to the tensor model has
Table 4: Parameters of the spectral tensor derived from different sources for U (40 m) =
40 m s−1 at sea.
Another example is shown in Table 4 where the Great Belt data from Mann (1994) are
extrapolated using neutral surface layer scaling to U (40 m) = 40 m s−1 . The spectral fit for
these values of U and z is shown in Figure 39.
Literature coherences and coherences derived from the spectral tensor by Eqs. (29) and
(30) are compared in Mann (1998). Generally, the agreement is good.
where sinc x ≡ (sin x)/x. For Ll ≫ L, the sinc2 -function is ‘delta-function-like’, in the sense
that it vanishes away from kl much faster than any change in Φij , and the area beneath the
sinc2 -curve is 2π/Ll . Therefore, we get
∗ (2π)3
Cik (k )Cjk (k ) = Φij (k ). (89)
V (B)
The solution to Eq. (89) is
(2π)3/2 1/2
Cij (k ) = Aij (k ) = (∆k1 ∆k2 ∆k3 ) Aij (k ) (90)
V (B)1/2
with A∗ik Ajk = Φij and ∆kl = 2π/Ll . This result should be expected when comparing Eq.
(26) to (83).
Two problems occur by simulating a field by the Fourier series Eq. (83) with the coefficients
Eq. (90). The first is that for many applications the dimensions of the simulated box of
turbulence need not to be much larger than the length scale of the turbulence model L.
Therefore Eq. (89) may not be a good approximation to Eq. (88). The second problem is
that the simulated velocity field Eq. (83) is periodic in all three directions. Both problems
have been addressed in Mann (1998).
The algorithms above simulate a three-dimensional vector field on a three-dimensional
domain, but it can easily be modified to simulate one- or two-dimensional vectors in a 2- or
3-D domain (Mann, 1998). The algorithms are not needed for a one-dimensional domain, i.e.
simulation of wind fluctuations in one point as a function of time.
The implementation of the model includes three steps:
1. Evaluate the coefficients Cij (k ), either by Eq. (90) or a modification of this (Mann,
1998).
2. Simulate the Gaussian variable nj (k ) and multiply.
3. Calculate ui (x ) from Eq. (83) by FFT.
The time consumption in the first step is proportional to the total number of points N =
N1 N2 N3 in the simulation. The required time to perform the FFT is O(N log2 N ) (Press
et al., 1992). In practice, simulating a three-dimensional field, used for load calculations on
wind turbines, with millions of velocity vectors takes of the order of a few minutes on a
modern pc.
References
Batchelor G. K. (1953) The theory of homogeneous turbulence, Cambridge University
Bendat J. S. and Piersol A. G. (1986) Random data: Analysis and Measurement Procedures, 2 edn, Wiley-
Interscience
Brown R. D. and Swail, V. R. (1991) Over-water gust factors. Ocean Engng. 18(4):363–394
Charnock H. (1955) Wind stress on a water surface. Q. J. Royal Meteorol. Soc. 81:639–640
Comte-Bellot G. and Corrsin S. (1971) Simple Eulerian time correlation of full- and narrow-band velocity
signals in grid generated, ‘isotropic’ turbulence. J. Fluid Mech. 48:273–337
Derbyshire S. H. and Hunt J. C. R. (1993) Structure of turbulence in stably stratified atmospheric boundary
layers; comparison of large eddy simulations and theoretical models. in S. D. Moobs and J. C. King (eds),
Waves and Turbulence in Stably Stratified Flows, Clarendon, Oxford, 23–59 pp
ESDU International (1982) Characteristics of wind speed in the lower layers of the atmosphere near the
ground: strong winds (neutral atmosphere), ESDU International, London
ESDU International (1985) Characteristics of atmospheric turbulence near the ground. Part II: single point
data for strong winds (neutral atmosphere), ESDU International, London
Garratt J. R. (1977) Review of drag coefficients over oceans and continents, Mon. Wea. Rev. 105:915–929
Geernaert G. L. (1987) On the importance of the drag coefficient in air-sea interactions. Dynamics Atmos.
Oceans 11:19–38
Grigoriu M. (1993) On the spectral representation method in simulation. Prob. Eng. Mech. 8:75–90
Højstrup J., Larsen S. E. and Madsen P. H. (1990) Power spectra of horizontal wind components in the neutral
atmospheric boundary layer. in N. O. Jensen, L. Kristensen and S. E. Larsen (eds), Ninth Symposium on
Turbulence and Diffusion, American Meteorological Society, 305–308 pp
Hunt J. C. R. and Carruthers D. J. (1990) Rapid distortion theory and the ‘problems’ of turbulence. J. Fluid
Mech. 212:497–532
Kaimal J. C. and Finnigan J. J. (1994) Atmospheric Boundary Layer Flows, Their Structure and Measurement,
Oxford University Press, New York
Kaimal J. C., Wyngaard J. C., Izumi Y. and Coté O. R. (1972) Spectral charcteristics of surface-layer turbu-
lence, Q. J. R. Meteorol. Soc. 98:563–589
Kristensen L. and Kirkegaard P. (1987) Puff kinematics, Technical Report R–548, Risø National Laboratory
Larose G. L. and Mann J. (1998) Gust loading on streamlined bridge decks, J. Fluids Structures 12(5):511–536
Mann J. (1992) Investigation of atmospheric low-frequency turbulence over the ocean, Technical Report
I–634(EN), Risø National Laboratoty, Roskilde
Mann J. (1994) The spatial structure of neutral atmospheric surface-layer turbulence. J. Fluid Mech. 273:141–
168
Mann J., Cariou J.-P., Courtney M. S., Parmentier R., Mikkelsen T., Wagner R., Lindelöw P., Sjöholm M.
and Enevoldsen K. (2009) Comparison of 3D turbulence measurements using three staring wind lidars and
a sonic anemometer. Meteorol. Z. 18(2):135–140
Mann J., Kristensen L. and Courtney M. S. (1991) The great belt coherence experiment – a study of atmo-
spheric turbulence over water, Technical Report R–596, Risø National Laboratory, Roskilde
NDP (1998) Acts, regulations and provisions for the petroleum activities, Vol. 2,, Norwegian Petroleum
Directorate
Panofsky H. A. and Dutton J. A. (1984) Atmospheric Turbulence, John Wiley & Sons, New York
Panofsky H. A., Larko D., Lipschutz R., Stone G., Bradley E. F., Bowen A. J. and Højstrup J. (1982) Spectra
of velocity components over complex terrain. Q. J. R. Meteorol. Soc. 108:215–230
Peña A., Gryning S.-E., Mann J. and Hasager C. B. (2010) Length scales of the netural wind profile over
homogeneous terrain J. Appl. Meteor. Climat. 49:792–806
Press W. H., Flannery B. P., Teukolsky S. A. and Vetterling W. T. (1992) Numerical Recipes, 2nd edn,
Cambridge University Press
Shinozuka M. and Deodatis G. (1991) Simulation of stochastic processes by spectral representation, Appl.
Mech. Rev. 44(4):191–203
Shinozuka M. and Deodatis G. (1996) Simulation of multi-dimensional gaussian stochastic fields by spectral
representation. Appl. Mech. Rev. 49(1):29–53
Shinozuka M. and Jan C.-M. (1972) Digital simulation of random processes and its applications. J. Sound
Vibration 25(1):111–128
Simiu E. and Scanlan R. H. (1996) Wind Effects on Structures, 3. ed., John Wiley & Sons
Sjöholm M., Mikkelsen T., Mann J., Enevoldsen K. and Courtney M. (2009) Time series analysis of continuous-
wave coherent Doppler lidar wind measurements. Meteorol. Z 18(3):281–287
Townsend A. A. (1976) The Structure of Turbulent Shear Flow, 2nd edn, Cambridge University Press
Townsend, A. A. (1980) The response of sheared turbulence to additional distortion. J. Fluid Mech. 98:171–
191
4.1 Introduction
Remote sensing offers the wind industry an attractive alternative or complement to the tra-
ditional methods for obtaining accurate wind measurements that involve the siting of tall
meteorological masts. Laser anemometry (lidar) has now demonstrated its capacity for re-
source assessment, wind turbine power curve measurement, and turbine mounted deployment
for advance wind speed detection. Widespread acceptance of lidar by the industry has re-
quired that this technique be extensively validated, aiming towards a certifiable and traceable
measurement standard and formal accreditation of lidar methods for different applications in
a range of terrain types. This chapter outlines the lidar measurement process and capabilities
specifically for the case of continuous wave (CW) systems.
Wind lidar systems were first demonstrated in the 1970’s (Jelalian, 1992) and have since
been applied to a wide variety of applications including aviation and meteorology. Although
applications to wind energy were explored in the 1980’s (Hardesty and Weber, 1987; Vaughan
and Forrester, 1989), the lidar systems that existed at that time were too large, complex
and expensive to achieve serious acceptance in the industry. The situation has now changed
dramatically, with rapid growth of the wind generation industry coinciding with development
of a new generation of lidars based on optical fibre and other components that emerged from
the telecommunications boom of the 1990’s. The first all-fibre lidars were demonstrated in
the late 1990’s, and a commercial prototype unit (ZephIR) was mounted on a turbine to
demonstrate wind speed detection in front of the rotor plane in early 2003. A demonstration
of ground-based wind profiling followed shortly afterwards. ZephIR is a CW coherent lidar
system, and this approach was selected as a means to combine simplicity with high sensitivity
at measurement ranges relevant to wind energy resource assessment, and hence achieve a
robust, reliable system at relatively low cost. Following this pioneering work, the pace of
development has accelerated, and lidars have increasingly become an established tool in the
wind industry.
Section 4.2 of this chapter provides an overview of lidar techniques and technology. Dif-
ferent types of lidar system are surveyed, and the generic physical principles underlying their
operation are reviewed. The specific case examined in detail here is that of wind profiling by a
ground-based conically-scanned continuous-wave lidar, which has rapidly become established
as a powerful tool in the wind energy industry, and is exemplified by the ZephIR lidar, initially
developed by QinetiQ, the Natural Power, and now by Zephir Limited.
A number of assumptions must be made in order to extract values of wind speed from raw
lidar data; these are reviewed in section 4.3. The necessary steps that are required during
the end-to-end measurement process in order to arrive at a value of wind speed are detailed
in section 4.4. It is important to understand the potential sources of error and uncertainty,
and these are reviewed and analysed in section 4.6. Section 4.7 examines the important
requirement for lidar calibration and traceability and 4.8 discusses the emerging area of turbine
mounted lidar. Finally, section 4.9 draws together some conclusions and a summary of the
future outlook for lidar in wind energy.
to allow very precise velocity determination; as a consequence the Doppler shift measurement
by a CLR system is inherently calibrated. Finally, the LO amplifies the signal via the beating
process to allow operation at a sensitivity that approaches the shot-noise (or quantum) limit.
This very high sensitivity permits the operation of CLR systems in an unseeded atmosphere,
relying only on detection of weak backscattering from natural aerosols.
CW systems are the simplest form of lidar possessing the advantage of reduced complexity,
and their performance can be tailored closely to the requirements of the wind industry. How-
ever, the overall trade-off between the pulsed and CW options for each specific application
must take into account a number of factors including sensitivity, cost, velocity resolution,
and maximum and minimum ranges. Unlike pulsed lidar systems, which use time of flight
to discriminate between returns from different ranges, a CW lidar achieves operation at a
given range by beam focusing. Wind profiling is achieved by continuously scanning the beam,
focusing in turn at a number of predetermined ranges. For each range, a circular scan is
typically used. The rapid sampling rate inherent to CW lidar permits 1-s “snapshots” of the
flow across the scan disk at each measurement range. Focusing of the lidar beam results in a
Lorentzian spatial weighting function along the beam axis, with the sensitivity peak located
at the beam waist (Sonnenschein and Horrigan, 1971; Karlsson et al., 2000). This function
has a half-width given by the Rayleigh range (the distance from the waist at which the beam
area has doubled).
The beam diameter at the waist increases linearly with range while the Rayleigh range
increases roughly as the square. Hence the effective probe volume varies as the 4th power of
the focus range, and this strong dependence has some implications for the signal statistics
at shorter ranges (Harris et al., 2001b). The minimum range for a CW lidar is very short
with detection possible in principle at zero range, whereas a pulsed system is blinded while
the pulse is leaving the transmitter leading to a minimum range of 10’s of metres, typically
around 40–50 m.
Figure 42: Conical scan pattern as used for lidar wind profiling. Left: ground based vertical
scanning. The cone half-angle (θ) is typically of order 30◦ . The lidar can operate successfully
even when part of its scan is obscured, e.g. by an adjacent met mast. To build up a wind
profile, the lidar operates in a repeating sequence during which all the heights are interrogated
in series. Right: one of several wind turbine mounted configurations, where the lidar scans
around a horizontal axis, usually pointing into the wind.
Figure 43: Stages of evolution of the ZephIR lidar. The upper left-hand picture (a) shows
the lidar head mounted on the nacelle of a Nordex N-90 wind turbine (March 2003). The
top central picture (b) shows a prototype ground-based wind profiler at the Risø wind energy
test site at Høvsøre, Denmark. The top right picture (c) shows an early ZephIR production
model deployed in the field. The lower pictures show a more recent dual mode ZephIR DM300
deployed on an offshore platform (d) and on a turbine nacelle (e).
The system achieved several weeks of successful operation. It was then converted into a
ground-based scanning unit for wind profiling (Figure 43-upper middle frame). The ground
based system was first trialled in December 2003, and subsequently deployed in numerous
campaigns in the UK, Europe, and other parts of the world. The experience gained through
these trials has built confidence in the robustness and reliability of the core ZephIR design. In
late 2004, work started on a production instrument (Figure 43-upper right frame), designed
to perform autonomous wind profiling measurements at heights up to 200 m (Smith et al.,
2006), primarily for site surveys at proposed wind farm sites. The technology was transferred
to Natural Power in 2007, and subsequent development resulted in the more integrated ZephIR
Z300 system (Figure 43-lower frames), and the dual mode DM300 which can be both turbine
and ground mounted. ZephIRs have logged more than 3.6 million hours of deployment (May
2013 figures) around the world and Zephir Limited was eventually established as a standalone
member of the Fred Olsen Limited group in 2012.
Figure 44: Theoretical lidar sensitivity curves at focus heights 25, 50, 75 and 100 m for the
two cases listed above with A = 20 mm and 28 mm, corresponding to respectively the original
(red curve) and current (blue curve) ZephIR design. The peak is normalised to unity in each
case; the absolute peak value decreases as the inverse of height squared, so that the area under
each curve (representing the overall sensitivity) is always the same. This illustrates a useful
feature of focused CW coherent lidar that in uniform scattering, the signal-to-noise ratio is
independent of focus range. Data obtained in calibration measurements (black squares) at a
calibration range R = 68 m are in close agreement with the corresponding theoretical values
(dashed curve) at the equivalent height 58 m (=68 m× cos 30◦ )
4.4.6 Photodetection
The beat signal is detected by directing the optically-mixed beam onto a photodetector
which measures fluctuations in the light’s intensity. In the telecommunications wavelength
band around 1.55 µm, reliable photodiodes are readily available that are well suited to this
purpose. The photodiode converts the incident photons into photoelectrons, which generate
a measurable current (or voltage) that is normally amplified before digitisation. There are
generally four contributions to the output of the photodetector module:
• Dark noise – the intrinsic wideband noise floor generated by the detector and amplifier
combination in the absence of any incident light. Dark noise is due to the random
generation of electrons and holes within the depletion region of the photodetector device
that are then swept by the photodetector’s electric field.
• Photon shot noise (Bleaney and Bleaney, 1976) (sometimes called quantum noise) –
the random generation of photoelectrons by the incident LO beam leads to a wideband,
spectrally flat (white) Gaussian noise source. The shot noise power spectral density
increases in proportion to the optical power of the LO beam.
• Laser relative intensity noise (RIN) – intensity fluctuations that are in excess of shot
noise, caused for example by relaxation oscillation of the laser output (Siegman, 1986).
For a RIN-dominated noise floor, the power spectral density increases as the square of
LO power. Such oscillation is typically at relatively low frequency, peaking below 1 MHz,
and hence only affects the sensitivity of the lidar at low line-of-sight wind speeds around
1 m s−1 . In some systems it is possible to cancel the RIN by use of a dual-channel
balanced detector.
• Beat term resulting from the wind signal – this is the contribution that contains the
information on Doppler shifts from which the wind speed is derived. Its power spectral
density increases in proportion both to the LO power and the signal power.
The requirements for the detector are high quantum efficiency, sufficient bandwidth to
cope with the maximum Doppler frequencies of interest, and for the shot noise contribution
to significantly exceed that of dark noise. This latter requirement depends on a combination
of the detector’s intrinsic noise floor and the optical saturation threshold.
Figure 45: Stages in typical lidar signal processing: DFT analysis is carried out by a computer
integrated into the lidar system. As an example, 4000 individual spectra might be averaged
to achieve high sensitivity and measurable returns even in very clear air. This entire process
takes only 20 ms, giving ∼ 50 measurements of line-of-sight wind velocity per second
(CNR)
Figure 46: Wind lidar output for a ground-based, vertical scan ZephIR, illustrating many of
the features of a wind profile measurement. This example has been obtained at a height
150 m above ground level, one of several heights being probed in sequence. The lower trace
shows 147 individual line-of-sight wind speed values, obtained over a total period of 3 seconds
(plotted as white squares against azimuth scan angle). The same data, along with the least-
squares fit in red, are displayed above in polar coordinates on the figure-of-eight plot showing
the wind bearing to lie slightly to the East of North. The wind parameters, derived from the
fit, appear in the table on the right; the horizontal wind speed at this height is determined
to be 9.1 m s−1 . The plot on the upper left shows one of the spectra from which each point
on the other graphs is derived.
The high level of redundancy in the fitting process is advantageous and can be used to
identify non-uniform flow. The root mean square deviation of the points from the optimum
solution gives an indication of the quality of fit, and this can be related to the value of
turbulent kinetic energy (TKE; see Wagner et al. (2009)). More work is needed to establish
a full understanding of the turbulence information available from lidar signals (Banakh et al.,
1999). Note that information on turbulence is also available from the spectral widths of the
Figure 47: ZephIR lidar wind speed correlation with instrumented wind tunnel pitot tube.
Courtesy of LM Wind Power, DTU Wind Energy, and NKT Photonics.
A greater source of error arises from uncertainty about what provides the scattering from
which the Doppler shift is derived. The scattering is assumed to originate from atmospheric
particles moving at the same speed as the wind and positioned close to the focus of the
lidar beam (section 4.3). An obvious example where this breaks down is when the beam
intersects a solid object (e.g. a bird) that is moving at a different speed from the wind
giving a measurement which could be in error. However, in such a case the value of hVLOS i
so derived will stand out as clearly anomalous on the polar plot (Figure 46). The presence
of such points will be diluted by approximately 50 correct values of hVLOS i obtained from
uncontaminated parts of the atmosphere, and their inclusion should not introduce any bias.
A further safeguard against these erroneous points is provided by a simple “outlier removal”
algorithm. This identifies points that lie anomalously far from the best fit solution to Eq.
(100) and eliminates them. The least-squares routine is then rerun on this slightly reduced
set of hVLOS i, data pairs. The presence of precipitation within the probe volume leads to a
different source of uncertainty. The downward motion of rain and snow inevitably leads to
some error in the vertical component of wind speed. However, the presence of rain and snow is
normally easily identified from the measurement process (for example by detecting activation
of a rain sensor), and the resulting values of vertical wind eliminated from the data. Other
wind parameters are unaffected and can still be correctly inferred.
Table 5: Results of correlation analysis of 10-minute averaged horizontal wind speed for a
ZephIR 300 trial at Høvsøre, Denmark in March 2011. Gradients, m, (forced through the
origin) and coefficients of determination, R2, of unity would imply perfect agreement be-
tween the lidar and the mast-mounted cup anemometers. (It should be noted that the slopes
very close to unity are slightly fortuitous, since the cup anemometer measurements have
uncertainties at least of order ±1%, due to calibration and mounting/shadowing effects).
Table 6: Results of a comparative trial of a ZephIR lidar against a very tall mast, equipped
with two types of cups at each height. The data indicate that the extended probe length at
greater heights did not result in excessive bias or errors. [1]: Forced through the origin; [2]:
Only hourly average containing 6 valid 10 min measurements are compared
Figure 50: Graph showing results from routine re-calibration of four cups to be used at an
accredited met mast site. The vertical axis shows differences when the same cup was calibrated
at two independent wind tunnel standards facilities. The tunnel-to-tunnel calibrations of the
same cups show variability of the order of 1%.
be mounted on the nacelle roof, within the rotor spinner (Figure 51) or built into the rotor
blades.
Figure 51: Continuous wave dual mode ZephIR 300 lidars mounted on a nacelle roof (left)
and within a wind turbine spinner (right).
The high sensitivity and hence fast data rates of circular scanning CW lidar make it very
well suited to turbine mounted applications. By adopting different signal processing and data
analysis strategies; hub height wind speed, direction, shear exponent and turbulence can be
measured, or speed, direction and turbulence at discrete heights across the rotor can calculated
(or both). The latter technique can be used to generate rotor equivalent wind data which has
been shown to produce more precise and representative wind turbine power curves, especially
for large rotor diameters (Wagner et al., 2008; Wharton and Lundquist, 2012).By measuring
the inclination of the lidar in real-time, the circular scan CW approach allows the effects of
nacelle motion on both line-of-sight velocity and measurement height to be negated.
4.8.1 Least-squares fitting routine for horizontal scanning (turbine mounted) oper-
ation
The use of a CW lidar for turbine mounted applications is fundamentally a quite different
arrangement when compared to a ground based, vertical scanning configuration. Unlike the
latter, the scan axis is approximately horizontal, and the lidar is almost always predominantly
staring into the wind. A consequence of this is that the polar plot (of the measured line-of-
sight wind speeds as a function of scan angle) is no longer a figure-of-eight shape, but instead
takes on a more circular appearance (Figure 52).
Figure 52: An example of visualisation and analysis of data from a turbine-mounted ZephIR.
Left: polar plot of raw data, showing line-of-sight wind speeds with scan angle. The radial
axis is the LOS speed. The breadth and structure of plotted distribution gives an indication of
the spatial turbulence within the scan volume e.g. ground induced shear and turbulence can
be seen in the lower range of angles. Low level wind jets and wakes from other turbines can
also be detected in this manner. Centre: real time analysis of the received polar plot, showing
centroids of the received line-of-sight spectra (red dots) and fitted wind parameters (indicated
by the green curve). The central red dots are turbine blade returns and are automatically
filtered prior to fitting. Right: reference data and wind characteristics calculated from the fit.
An example of filtering that can be required is for the case of turbine blades. For a turbine
mounted CW lidar, situated on the roof of a turbine’s nacelle, and scanning upwind through
the turbine blades, the lidar must contend not only with quasi-periodic blocking of the beam,
but also strong Doppler returns from the blades themselves. Although the intensity of the
back reflected laser signals can be very high from these blades (typically 50 times higher than
the wind returns), this can help distinguish them from the line of sight Doppler returns from
the incoming wind. Additionally, the relatively slow, near perpendicular path of the blade
surfaces means that the Doppler shifts are relatively low frequency (giving Doppler returns
corresponding to typically < 2 m s−1 ). Hence efficient blade rejection filters, which remove
these signals from the wind field fitting process, are simple to implement. However, blade
effects do reduce the number of data points around the scan, and for this reason, hub (or
spinner) mounted CW lidars can have some advantages.
4.8.2 Turbine mounted CW lidar for wind turbine power curve measurement
Continuous wave lidar have now been deployed in a number of wind turbine power curve
measurement campaigns. Some results from two of these are reproduced in this section. The
first campaign discussed here was organised by ROMO Wind and carried out in flat terrain
Figure 53: (a) Measured turbine yaw misalignment before and after nacelle vane recalibration.
(b)Wind turbine power curves measured by nacelle mounted lidar before and after nacelle vane
recalibration at a measurement range of 180 m.
Measurements were also made at a series of other ranges between 10 m and 180 m and
are shown in Figure 54). While the measurements at close range ”flatter” the power curve
due to the rotor induction effect, they demonstrate that high quality power curves can be
obtained very close to the rotor. While these would have to be corrected to be quantitative,
they are very representative of the true wind field incident on the turbine rotor and can be
expected to be immune to say, the effects of complex terrain, as compared to more remote
measurements.
Figure 54: (a) time series of wind speed measurements taken at a series of measurement
ranges. The measurements are highly correlated, but rotor induction (blockage) effect is
clearly visible.(b) Series of wind turbine power curves generated from measurements at several
ranges. The upper curve is from the closest range (10 m), the lower curves were measured at
30, 50, 100 and 180 m respectively.
The second study reported here took place on a DTU test turbine situated at Roskilde,
Figure 55: Correlation plot of hub height wind speeds measured by a nacelle mounted CW lidar
and a pair of cup anenometers on a 500 kW Nordtank wind turbine at Roskilde, Denmark.
The lidar data was then reprocessed to measure wind speed and direction at the height of
each of the cups listed above, taking full account of any nacelle motion on both measurement
height and line-of-sight velocity. Turbine SCADA power output data was corrected for air
density and temperature and used to generate rotor equivalent power curves from both the
met mast and the lidar data (Figure 56). This preliminary result indicates that turbine mounted
CW lidars are capable of measuring power curves of comparable quality to met masts, but
with the advantage of easy redeployment. These are believed to be the first reported rotor
equivalent power curve measurements using a commercially available nacelle mounted lidar.
Subjects of current research in the turbine mounted lidar area include feed-forward control
for stress load reduction and short range (i.e. less than 2 rotor diameters) wind speed mea-
surement for turbines situated in more complex terrain where long range measurements are
less representative of the actual wind field incident on the turbine rotor.
requirements in the field of wind energy. Examination of the measurement process reveals
that the basic acquisition of line-of-sight Doppler spectra is a well-established method with
little scope for gross errors and miscalibration. The subsequent steps required to convert these
spectra into a profile of wind speed are more complex, however, and their validity relies on a
number of well-established assumptions. Much work has been performed to test the validity
of the assumptions outlined in section 4.3, and to understand the uncertainties and other
issues discussed in section 4.6.
Complex terrain remains a topic of great interest as it becomes increasingly necessary to
explore less ideal locations as potential wind farm sites. In such sites the horizontal wind
speed deduced by conically-scanned lidar can be subject to differences in comparison to that
measured by co-located cup anemometers when the flow is non-uniform across the lidar mea-
surement disk. A method has recently been developed in which the impact of inhomogeneous
flow at complex flow sites is examined using computational fluid dynamics (CFD modelling to
predict the bias that will be experienced by a lidar in comparison to a conventional met mast
equipped with cup anemometers. Similar percentage changes in wind speed as measured by
a mast are shown to occur if the mast were to be moved by ±50 m from its original location.
This suggests a methodology for resource assessment in complex terrain in which lidar is
used in combination with CFD modelling in order to (i) adjust the lidar data for the impact
of non-uniform flow and (ii) investigate the wind variations across the site that are a major
source of uncertainty for current techniques.
Lidar offers some potential advantages in turbine power curve measurement. The measure-
ment over an extended volume may give a more representative estimate of the wind energy
content of the air interacting with the blades, and the ability to re-position the lidar quickly
is clearly advantageous. A study reported by Wagner et al. (2008) has shown that exploiting
the lidar wind profile data can reduce the scatter of points in a measured power curve. In
another recent study (Cayla, 2010) a ZephIR lidar gave an almost identical power curve to
an IEC-instrumented power performance mast. The scatter of the points in the power curve
obtained using the ZephIR data at hub height was somewhat lower than that for the mast.
This result needs further investigation and possibly is a consequence of the more effective
sampling of the wind around the scan disk. It follows, interestingly, that remote sensing equip-
ment that agrees perfectly with the mast would therefore have provided higher scatter in the
power curve than ZephIR!
The extraction of turbulence data relevant to the wind industry from lidar signals is an area
that will benefit from further research and verification through field comparisons. Turbulence
Table 8: Correlation analysis from the first SeaZephIR trial in 2009: the table shows gradient
(m) and coefficient of determination R2 for plots of 10-minute wind speed for SeaZephIR on
a floating platform versus those measured by a second ZephIR unit positioned 800 m away
on land
Forward-looking turbine mounted lidar, either on the nacelle or in the hub, is another ex-
citing lidar development. Applications include turbine power curve measurement, energy yield
optimisation (e.g. by reducing turbine yaw misalignment) and gust and fatigue load reduction
allowing longer turbine lives and/or turbine build cost reduction. As already remarked, CW
lidar seems particularly well suited to this type of application, owing to its high sensitivity
(high average photon flux), high sample rate (50 Hz), and scan path that probes the wind
around the rotation path of the rotor. Another of its benefits is its flexibility in terms of
turbine mounting. In addition to nacelle roof mounting, it is, to date, the only class of lidar
that has been installed in a rotor hub (or spinner). Interest in the concept has increased
significantly since the world-first proof-of-principle demonstration of turbine-mounted lidar in
2003 (Harris et al., 2006, 2007), with several groups currently working towards evaluating
Figure 57: CW lidar experiments in a wind tunnel, prior to blade lidar deployment experiments.
Left: wind tunnel schematic. Centre: CW lidar twin telescopes. Right: External view of the
wind tunnel, showing the ZephIR 300 base unit.
In connection with turbine mounted lidars, significant recent efforts in the industry have
focussed on quantifying their potential benefits, as well as looking at the optimum lidar
configurations to use. For CW lidars, the cone scan angle, the number of ranges to scan over
(if indeed more than one is required), scan rates and the LOS processing algorithms are all
being investigated. Recent results in the literature have included:
1. Conical scan CW lidar was used to determine yaw alignment of a lidar (Kragh et al.,
2013) and demonstrated the ability to achieve a sub 4◦ yaw error over a 2 hour period,
even during periods of high turbulence.
2. Simulations examining the ability of turbine mounted lidar for accurate yaw alignment
(Kragh and Hansen, 2011) indicated yield, at below rated power, could be raised by 1%
to 5%.
3. A study reported in Schlipf et al. (2011) comparing conventional nacelle based wind vane
with lidar yaw alignment control, indicated that the yearly energy output of a 5 MW
turbine could be enhanced by ˜ 2% using the lidar.
4. Schlipf and Kuhn (2008) modelled the benefits of a nacelle mounted lidar for feed-forward
control, in particular turbine speed control. The study found reductions in standard
deviations of 91%, 90% and 71% for rotor speed, tower fore-aft moment and blade root
flap moment for gusts. For turbulent airflows, the reductions in standard deviations were
77%, 32% and 17% respectively.
5. Simley et al. (2011) simulated a conical scan CW lidar and showed accurate yaw align-
ment should be possible. Even in highly turbulent airflow, a precision of a few degrees
was achievable. The same paper also showed that RMS wind speed measurement errors
were lower for a CW system than a pulsed system for ranges ¡125m
6. Simulations using lidar feed-forward control (Laks et al, 2011) showed turbine fatigue
load reductions of approximately 20%.
Acknowledgements
The authors are grateful for the support and enthusiasm of their colleagues and collaborators.
Without them, wind power lidar technology could not have evolved to its current advanced
state.
Notation
a floating parameter for the fit of the line-of-sight velocity
A beam radius at the output lens
ADC analogue-to-digital converter
b floating parameter for the fit of the line-of-sight velocity
B wind bearing
c speed of light
floating parameter for the fit of the line-of-sight velocity
CFD computational fluid dynamics
CLR coherent laser radar
CNR carrier-to-noise ratio
CW continuous wave
DFT digital Fourier transform
D(ν) power spectral density from dark noise
ELO LO field
Es stable signal field
FFT fast Fourier transform
FPGA field-programmable gate array
h Planck constant
i fluctuating detector power output
IR infrared
LO local oscillator
m slope of the linear regression
Ps time-average optical signal power
PT transmitted laser power
R distance of the beam focus from the lidar output lens
R2 coefficient of determination
RIN laser relative intensive noise
D(ν power spectral density from RIN
SNR signal-to-noise ratio
t time variable
TI turbulence intensity
TKE turbulent kinetic energy
u wind speed component in the x-direction
v wind speed component in the y-direction
VAD velocity-azimuth-display
VH horizontal wind speed
VLOS line-of-sight wind speed
w wind speed component in the z-direction
References
Angelou N., Mikkelsen T., Hansen K. H., Sjöholm M., and Harris M. (2010) LIDAR Wind Speed Measurements
from a Rotating Spinner: SpinnerEx 2009. Project Upwind and Risø report 34: Risø-R-1741(EN)
Banakh V. A., Smalikho I. N., Köpp F., and Werner C. (1993) Representativeness of wind measurements with
a CW Doppler lidar in the atmospheric boundary layer. Appl. Opt. 34:2055–2067
Banakh V. A., Smalikho I. N., Köpp F., and Werner C. (1999) Measurements of turbulent energy dissipation
rate with a CW Doppler lidar in the atmospheric boundary layer. J. Atmos. Oceanic Tech. 16:1044–1061
Barker W. (2009) Analysis of ZephIR data from Mason City, Iowa. Natural Power internal report
Barker W., Pitter M., Burin de Roziers E., Harris M., and Scullion R. (2012) Can lidars measure turbulence.
EWEA conference PO74
Bingöl F., Mann J., and Foussekis D. (2008) Lidar error estimation with WAsP engineering. IOP Conf. Series:
Earth and Environ. Sci. 1:012058
Bingöl F., Mann J., and Foussekis D. (2009) Conically scanning lidar error in complex terrain. Meteorol. Z.
18:189–195
Bleaney B. I. and Bleaney B. (1976) Electricity and magnetism, Oxford University Press, Section 23.4
Cayla M. (2010) Comparison of ZephIR measurements against cup anemometry and power curve measure-
ments. Natural Power internal report (http://www.zephirlidar.com/resources/publications)
Chanin M. L., Gariner A., Hauchecorne A., and Portneuve J. (1989) A Doppler lidar for measuring winds in
the middle atmosphere. Geophys. Res. Lett. 16:1273–1276
Clifford S. F. and Wandzura S. (1981) Monostatic heterodyne lidar performance: the effect of the turbulent
atmosphere. Appl. Opt. 20:514–516
Courtney M. and Gottschall J. (2010) ZephIR 145 validation test. Risø report for Natural Power
Hardesty R. M. and Weber B. F. (1987) Lidar measurement of turbulence encountered by horizontal-axis wind
turbines. J. Atmos. Oceanic Tech. 67:191–203
Harris M., Pearson G. N., Hill C. A., and Vaughan J. M. (1994) Higher moments of scattered light fields by
heterodyne analysis. Appl. Opt. 33:7226–7230
Harris M., Constant G., and C Ward (2001) Continuous-wave bistatic laser Doppler wind sensor. Appl. Opt.
40:1501–1506
Harris M., Pearson G. N., Ridley K. D., Karlsson C. J., Olsson F. A., and Letalick D. (2001) Single-particle
laser Doppler anemometry at 1.55 µm. Appl. Opt. 40:969–973
Harris M., Hand M., and Wright A. (2006) Lidar for turbine control. Tech. Report NREL/TP-500-39154
Harris M., Bryce D. J., Coffey A. S., Smith D. A., Birkemeyer J., and Knopf U. (2007) Advance measurements
of gusts by laser anemometry. J. Wind Eng. Ind. Aerodyn. 95:1637–1647
Karlsson C. J., Olsson F. A., Letalick D., and Harris M. (2000) All-fiber multifunction CW 1.55 micron
coherent laser radar for range, speed, vibration and wind measurements. Appl. Opt. 39:3716–3726
Kragh K. A., and Hansen M. H (2011) Individual Pitch Control Based on Local and Upstream Inflow Measure-
ments. Proceedings of 49th AIAA Aerospace Sciences, American Institute of Aeronautics and Astronautics,
AIAA 2011–813, Orlando, Florida
Kragh K. A., Hansen M. H, and Mikkelsen T. (2013) Precision and shortcomings of yaw error estimation
using spinner-based light detection and ranging. Wind Energy 16:353–366
Lading L., Hanson S., and Skov Jensen A. (1984) Diffraction-limited lidars: the impact of refractive turbulence.
Appl. Opt. 23:2492-2497
Laks J., Pao L. Y., Simley E., Wright A., Kelley N., and Jonkman B. (2011) Model Predictive Control Using
Preview Measurements From LIDAR. Proceedings of 49th AIAA Aerospace Sciences, American Institute
of Aeronautics and Astronautics, AIAA 2011–813, Orlando, Florida.
Loudon R. (2000) The quantum theory of light, Oxford University Press, 3rd edition
Mikkelsen T., Hansen K., Angelou N., Sjöholm M., Harris M., Hadley P., Scullion R., Ellis G., and Vives G.
(2010) Lidar wind speed measurements from a rotating spinner. European Wind Energy Conf., Warsaw
Pearson G. N., Roberts P. J., Eacock J. R., and Harris M (2002) Analysis of the performance of a coherent
pulsed fiber lidar for aerosol backscatter applications. Appl. Opt. 41:6442-6450
Pedersen A. T., Montes B. F., Pedersen J. E., Harris M., and Mikkelsen T. (2012) Demonstration of short-
range wind lidar in a high-performance wind tunnel, EWEA conference PO78, Copenhagen
Peña A., Hasager C. B., Gryning S.-E., Courtney M., Antoniou I., Mikkelsen T. (2009) Offshore wind profiling
using light detection and ranging measurements. Wind Energy 12:105-124
Pitter M., Abiven C., Harris M., Barker W., and Brady O. (2012) Lidar and computational fluid dynamics for
resource assessment in complex terrain. EWEA conference PO80, Copenhagen.
Rogers T., Byme A., McCoy T., and Briggs K. (2012) Expected impacts on cost of energy through lidar based
wind turbine control. EWEA 2012.
A Rutherford A, M Harris M, W Barker W, E Burin de Roziers E., R Pitter M., R Scullion R., and C Slinger
C. (2012) Lidar calibration and performance validation process. AWEA 2012, Atlanta.
Schlipf D., and Kuhn M. (2008) Prospects of a collective pitch control by means of predictive disturbance
compensation assisted by wind speed measurements. German Wind Energy Conference (DEWEK 2008),
pp. 1–4
Schlipf D., Anger J., Kapp S., Bischoff O., Hofsäß M., Rettenmeier A., and Kühn M. (2011) Prospects of
optimization of energy production by lidar assisted control of wind turbines Proceedings of European Wind
Energy Conference and Exhibition, EWEA 2011.
Simley E., Pao L. Y., Frehlich R., Jonkman B., and Kelley N (2011) Analysis of Wind Speed Measurements
using Continuous Wave LIDAR for Wind Turbine Control. Proceedings of 49th AIAA Aerospace Sciences,
American Institute of Aeronautics and Astronautics, AIAA 2011-263, Orlando, Florida
Slinger C., Leak M., Pitter M., and Harris M. (2013) Relative Power Curve Measurements Using Turbine
Mounted, Continuous Wave Lidar. EWEA 2013.
Smith D. A., Harris M., Coffey A. S., Mikkelsen T., Jørgensen H. E., and Mann J. (2006) Wind lidar evaluation
at the Danish wind test site in Høvsøre. Wind Energy 9:87–93
Sonnenschein C. M. and Horrigan F. A. (1971) Signal-to-noise relationships for coaxial systems that heterodyne
backscatter from the atmosphere. Appl. Opt. 10:1600–1604
Vaughan J. M. and Forrester P. A. (1989) Laser Doppler velocimetry applied to the measurement of local
and global wind. Wind Eng. 13:1–15
Wagner R., Jørgensen H. E., Paulsen U. S., Larsen T. J., Antoniou I., and Thesbjerg L. (2008) Remote sensing
used for power curves. IOP Conf. Series: Earth and Environ. Sci. 1:012059
S Wharton S., and Lundquist J. K. (2012) Atmospheric stability affects wind turbine power collection. Envi-
ronmental Research Letters textbf7: 014005.
Zak J. A. (2003) Atmospheric boundary layer sensors for application in a wake vortex advisory system.
NASA/CR-2003-212175
This section complements the description of the measurement process for a Doppler lidar
for wind speed and direction determination, by focusing more on the pulsed lidar technology,
and particularly on the WINDCUBETM Doppler pulsed lidar, one of the most accurate remote
sensing devices available at the present time in the wind industry.
5.1 Introduction
There is a pressing need for good wind-speed measurements at greater heights and at more
site locations to assess the availability of the resource in terms of power production and to
identify any frequently occurring atmospheric structural characteristics that may impact the
operational reliability and lifetime of wind turbines and their components. To measure the
wind field up to the height of new generation wind turbines and at any location of interest,
remote sensors are needed to complement masts.
Different technologies are in use in this field, among them pulsed lidars (LIght Detection
And Ranging). The underlying principle of pulsed lidar measurement of wind and aerosols
is the use of optical heterodyne (coherent) detection, in which laser pulses are transmitted
into the atmosphere and scattered off of naturally-occurring small dust particles (aerosols)
entrained in the ambient flow field (Frehlich et al., 1994; Huffaker and Hardestry, 1996; Soreide
et al., 1997; Frehlich et al., 1998). Even though the measurement principle is well known and
similar to pulsed radars, pulsed lidars have only been used in wind energy site assessment
only since 2008. Their recent introduction is mainly due to laser revolution coming from fiber
telecommunication development in the late 1990s.
In this section, we describe the architecture of pulsed lidars, based on the WINDCUBETM
lidar developed by LEOSPHERE and ONERA, the French Aerospace Lab. We define the differ-
ent modules of a pulsed lidar, their specific functions and the individual level of uncertainties
they might bring to the lidar wind speed retrieval. We also show the differences between
pulsed and continuous wave (CW) lidars, used as well for wind sensing. We eventually give
the lidar equation in pulsed mode, giving the relationship between parameters, and we focus
on range and speed accuracy and resolution.
Laser source A pulsed lidar needs a continuous wave laser, called master oscillator (MO) to
generate the local oscillator (LO) beam and a pulsed laser to generate the powerful transmitted
pulse. The frequency offset between the two sources need to be stable with time to allow an
unbiased measurement of the Doppler shift.
The master oscillator provides the laser wavelength, the laser linewidth, the laser intensity
noise and the state of polarization. Each of these parameters has to be well known and stable
to guarantee the lidar performance. The required CW power is at least some milliwatts.
The pulsed laser delivers cyclic pulses of high energy. The pulse duration is some hundreds
of nanoseconds, that determines the length of the pulse in the atmosphere and so the spatial
resolution.
Pulse duration [ns] Pulse length [m] Minimum spatial resolution [m]
200 60 30
400 120 60
800 240 120
The pulse repetition frequency (PRF) is as high as possible, but cannot exceed a maximum
value PRFmax . To avoid ambiguity between return signals, the time between pulses (1/PRF)
must be longer than the round trip time of flight of the pulse to the greatest height to be
measured Zmax , PRFmax = c/(2Zmax), where c is the speed of light.
Circulator Figure 60 describes the circulator and the different elements which compose
it. The function of the circulator is to transmit the laser pulse from the laser (1) to the
telescope (2) and to direct the backscattered light from the telescope (2) to the receiver
(3). Power handling, transmission efficiency and isolation between (1) and (3) ports are
critical parameters. In coherent pulsed lidars, most circulators use polarization to perform
this function. Polarization is rotated in the telescope thanks to a quarter wave plate. The
transmitted polarization is circular.
Telescope The telescope magnifies the laser beam in order to reduce its divergence in
the far field, and focuses the beam at any distance. The larger the beam, the smaller the
divergence and the better the signal-to-noise ratio (SNR) at long range. Since the lidar SNR
is inversely proportional to the beam area, beam diameter must be minimized over the global
measurement range to ensure maximum efficiency.
The telescope can use reflective or refractive optics, perfectly corrected from geometrical
aberrations. In order to lose less than 3 dB on the detection efficiency, wavefront distortion
on the global roundtrip optical path has to be less than λ/4 RMS, including components and
atmospheric distortion.
Atmospheric turbulence as well creates wave distortion that degrades heterodyne efficiency.
To keep this distortion negligible, the telescope aperture has to be smaller than the coherence
diameter of the beam expressed by:
do = 2.4e10−8 λ6/5 Z −3/5 Cn−6/5 . (102)
Considering maximum values of index structure constant (Cn2
= 10 −13
m −2/3
) and propaga-
tion over about Z = 1 km, this limits the size of the telescope to about 10 cm at ground
level.
To measure the wind speed in any direction, and display 2D or 3D wind maps, a more
flexible scanner is needed. Double flat mirror scanners are mostly used (CTI, Halophotonics,
Leosphere), while double prism scanners, like those in Risø DTUs Windscanners (Mikkelsen,
2008), offer a more original and compact solution. A global hemispherical field of view can
be obtained with both solutions.
Tatm = exp(−2αZ) if the atmosphere is homogeneous and βpi (Z) is the backscattering
coefficient of the atmosphere at a distance Z, τ is the full-width-half-maximum (FWHM)
pulsed duration and Ω is the reception solid angle:
πσ 2
Ω= , (106)
Z2
where σ is the efficient telescope aperture radius and α and βπ are roughly proportional since
they both depend on aerosol concentration in the atmosphere.
However, because of the limited heterodyne efficiency, only a part of Pr (Z) is efficient for
the coherent detection. The lost part of Pr (Z) comes from phase and polarization mismatches.
An interesting way to estimate the detection antenna diagram is to propagate back in the
atmosphere the LO, and to compute the overlap integral of the signal and LO along the LOS.
Figure 63: Variations of I(Z) with focus distance Ft = 40, 80, 120, 160, 200, 240 m respec-
tively.
The current power on the detector can then be derived using the same equations as in CW
mode,
< i2het >= 2ηhet S 2 PLO Ps (Z), (109)
where PLO is the local oscillator power, ηhet the heterodyne efficiency, which depends on
phase, amplitude and polarization matching and S the detector sensitivity.
R R ∗
2
R R Ad Es (x, y)E
R LO (x, y)dxdy
R
ηhet = 2 2 , (110)
Ad Es (x, y)dxdy Ad ELO (x, y)dxdy
< i2het >= ηhet Tinst Tatm S 2 βπ (Z)PLO Ppeak τ cI(Z). (111)
The lidar equation shows that the signal power is proportional to pulse energy Ppeak τ and
proportional to LO power. LO amplifies the signal allowing it to be detected over the detector
noise.
Figure 64: Orthogonal frame of the WINDCUBE for retrieving the wind speed components
Figure 65: Variation of CNR vs distance (Altitude/cos(θ-dev)) for different beam radius of
curvature
Best focus corresponds to maximizing the data availability at all altitudes, from H = 40
m (Z = 46 m) to H = 200 m (Z = 230 m). Practically, it corresponds to balancing and
Figure 66: Pulse propagation and width of range gate represented on a time-distance plot
Velocity measurement at one point Vd (Z) depends on the velocity V r(R) at close points
and on the range weighting function RW F (R),
Z ∞
V d(Z) = RW F (R)V r(R)dR. (130)
−∞
For a pulsed lidar, Banakh and Smalikho (1994) gave an analytical equation for RWF, when
the pulse is Gaussian (FWHM = τ ) and range gate is flat (width τm ), as the convolution
between the pulse power profile and the range gate profile:
" √ √ ! √ √ !#
1 4 ln 2 ln 2 4 ln 2 ln 2
RW F (Z) = erf (Z − Zo ) + − erf (Z − Zo ) − .
τm c cτ τ cτ τ
(131)
Range resolution is defined as the FWHM of the function RW F which is roughly:
cτ
∆z1 = √ m . (132)
2erf ln 2τm /τ
Because of spectral extent of the signal and low frequency noise, this range is in fact a bit
smaller. For example, with λ = 1.55 µm, Fi = 68 MHz and Fs = 250 MHz, the practical
horizontal velocity range is [−50 m s−1 , +50 m s−1 ]. A passband filter cancels outband
Doppler shifts to avoid Doppler ambiguities.
Velocity resolution Velocity precision depends on both atmospheric parameters and lidar
parameters. Both are broadening the Doppler spectrum and hence limit the frequency es-
timation. Atmospheric parameters are wind gradient within the range gate and turbulence.
Lidar parameters are pulse duration CNR and number N of average spectra. The smaller the
pulse duration, the smaller the range gate and the velocity dispersion but larger the frequency
spectrum.
Eq. (135) gives the minimum velocity resolution as a function of relative parameters. It is
called Cramer Rao lower bound (CRLB),
√ √
2λ 1 + CNR
σvcrlb = √ . (135)
2τ N CNR
Eq. (135) assumes an infinite correlation time of the signal. In an actual lidar, σv is limited
by the finite correlation time τc , the smallest value between the pulse duration and the
Figure 68: Cramer Rao boundary for 100 averaged spectra (dark red) and 10000 averaged
spectra (light red). CNR is measured in narrow band (B = 1/τm )
For low CNR values, σv decreases as 1/CNR. Doubling the pulse energy divides by 2
the velocity resolution. For high CNR values, σv is constant. Spectral broadening comes from
speckle fluctuations in√the signal. The only way to reduce σv is to increase N . For intermediate
CNR, σv varies as 1/ CNR. It is therefore equivalent to increase pulse energy or number of
pulses. In this region, σv depends only on the average laser power. Using a low energy, high
pulse repetition rate (PRF) laser is then equivalent to using a high energy, low PRF laser,
assuming the average power is constant. This is one for the reasons of the recent raise of
high PRF fiber lasers for pulsed Doppler lidars.
Range ambiguity The rate at which pulses are transmitted, the PRF limits the range over
which heights can be unambiguously determined. To avoid ambiguity between return signals,
the inter pulse period (IPP=1/PRF) must be longer than the round trip time of flight of the
pulse to the greatest height. For example, to measure without ambiguity up to 5 km, the
PRF needs to be less than 30 KHz,
c
PRFmax = . (138)
2Zmax
Ando T., Kameyama S. and Hirano Y. (2008) All-fiber coherent Doppler lidar technologies at Mitsubishi
Electric Corporation. IOP Conf. Series: Earth and Environ. Science 1:012011
Banakh V. and Smalikho I. (1994) Estimation of the turbulence energy dissipation rate from the pulsed
Doppler lidar data. J. Atmos. Oceanic Technol. 10
Bingöl F., Mann J. and Foussekis D. (2009) Conically scanning lidar error in complex terrian. Meteorol. Z.
18:1–7
Boquet M., Parmentier R. and Cariou J. P. (2009) Analysis and optimization of pulsed Doppler lidar in complex
terrain. EWEC, Marseille
Boquet M., Parmentier R. and Cariou J. P. (2010) Correction of pulsed wind lidars bias in complex terrain.
ISARS
Boquet M., Callard P., Neve N. and Osler E. G. (2010) Return on investment of a lidar remote sensing device.
DEWI Magazin 37:56–61
Cariou J.-P., Augere B. and Valla M. (2006) Laser source requirements for coherent lidars based on fiber
technology. Comptes Rendus Physique 7:213–223
Courtney M., Wagner R. and Lindelöw P. (2009) Commercial lidar profilers for wind energy. A comparative
guide. EWEC, Marseille
Faghani D. et al. (2009) Remote sensing: practical use for wind power project. AWEA Wind Resource &
Project Energy Assessment Workshop, Portland
Frehlich R. G. and Kavaya M. J. (1991) Coherent laser radar performance for general atmospheric refractive
turbulence. Appl. Opt. 30
Frehlich R. G., Hannon S. M., Henderson S. W. (1993) Performance of a 2-coherent Doppler lidar for wind
measurements. J. Atmos. Oceanic. Technol. 11:1517–
Frehlich R. G., Hannon S. M. and Henderson S. W. (1998) Coherent Doppler lidar measurements of wind
field statistics. Bound.-Layer Meteorol. 86:233–256
Gottshall J., Courtney M. S., Lindelöw P. and Albers A. (2010) Classification of lidar profilers: How to introduce
lidars to power performance testing. EWEC, Brussels
Huffaker R. M. and Hardestry R. M. (1996) Remote sensing of Atmospheric wind velocities using solid state
and CO2 coherent laser systems. Proc. IEEE 84
Jaynes D. (2009) Lidar validation and recommendations for wind resource assessments. Wind Power Conf.,
Chicago
Lindelöw P. (2007) Effective sample volume fiber based coherent lidars for remote sensing of wind. PhD thesis
from Ørsted-DTU, Lyngby
Mikkelsen T. (2008) WindScanner: a facility for wind and turbulence measurements around large wind turbines.
In Nordic Conf. on: Global Challenges-Regional Opportunities: How can research infrastracture and Science
support Nordic competiveness?, Stockholm
Oldroyd A. et al. (2009) Testing and calibration of various lidar remote sensing devices for a 2 year offshore
measurement campaign. EWEC, Marseille
Siegman A. E. (1966) The antenna properties of optical heterodyne receivers. Appl. Optics 5:1588–1594
Soreide D., Boque R. K., Seidel J. and Ehernberger L. J. (1997) The use of a lidar forward looking turbulence
sensor for mixed-compression inlet unstart avoidance and gross weight reduction on a high speed civil
transport Nasa Technical memorandum
Westerhellweg A., Beeken A., Canadillas B., Neumann T. (2010) One year of lidar measurements at Fino1-
platform: comparison and verification to met-mast data. DEWEK
6.1 Introduction
This chapter gives an overview of the derivation of mixing-layer height (MLH) and the de-
tection of low-level jets (LLJs) by surface-based remote sensing instruments such as sodar,
lidar, ceilometer and RASS. The detection of vertical profiles indicating the structure of the
atmospheric boundary layer (ABL) is one of the principal tasks of experimental boundary-
layer research. MLH has become an important input parameter for the description of wind
profiles above the surface layer (Gryning et al., 2007; Peña et al., 2010). LLJs are secondary
wind maxima which occur at the top of a stable boundary layer. Over land, they are observed
at several hundred metres above ground at night-time. Over the sea, they are found with
offshore winds when warmer air flows over cooler water at the top of the shallow stable inter-
nal boundary layer at any time of the day. These shallow internal boundary layers are often
considerably less than 100 m deep.
The next section describes methods to detect the mixing-layer height, while Section 6.3
briefly mentions methods to capture the boundary-layer height. Section 6.4 gives more infor-
mation on low-level jets. Acoustic (sodar) and optical sounding techniques (lidar) have got a
broad coverage elsewhere in this volume. RASS techniques are still quite unusual in the assess-
ment of wind resources. Therefore, a subsection on technical details of this instrumentation
has been added in section 6.2. A more complete survey of remote sensing instrumentation
is given in Emeis (2010), an overview of applications of ground-based remote sensing is pre-
sented in Emeis (2011). The full scope of wind energy meteorology is presented in Emeis
(2012).
Figure 70: Sample time-height cross-section from acoustic sounding with a sodar. Left: acous-
tic backscatter intensity, right: horizontal wind direction. Thin black lines demark inversions.
valley air over a snow-covered valley floor. The multiple layering originated from an interlacing
of down-valley (wind direction around 190◦) and down-slope (wind direction around 230◦ )
flows. The layers are separated by temperature inversions and each higher layer is potentially
warmer than the next lower layer. They persisted nearly the whole day because no vertical
mixing took place in the stably stratified valley atmosphere.
Horizontal wind speed (HWS) method The HWS method requires the analysis of the
Doppler shift of the backscattered acoustic signals. The algorithm is based on the analysis
of the shape of hourly-averaged vertical wind speed profiles using the assumption that wind
speed and wind direction are almost constant within the mixing layer but approach gradually
towards the geostrophic values above the mixing layer. Beyrich (1997) listed this method in
his Tab. 2 but did not discuss it further. The applicability of the method is probably limited
to the well-developed convective boundary layers (CBL) due to the underlying assumptions.
Such CBLs are often higher than the maximum range of a sodar. Even if the CBL height is
within the range of the sodar the algorithm for the analysis of the Doppler shift often fails
above the inversion topping of the CBL due to too low signal-to-noise ratios. Today, small
Doppler wind lidars are available to derive wind speed and direction profiles through the whole
depth of the boundary layer. This facilitates the application of the HWS method.
Vertical wind variance (VWV) method The VWV method is also working only for CBLs.
It is based on the vertical profile of the variance of the vertical velocity component σw . In a
CBL σw reaches a maximum in a height azi . Typical values for a are between 0.35 and 0.40.
Thus, in principle, this is an extrapolation method. It has been tried for sodar measurements
because it permits a detection of MLH up to heights which are 2.5 times above the limited
maximum range (usually between 500 and 1000 m) of the sodar. Beyrich (1997) classified
this method as not reliable. A related method, which is based on power spectra of the vertical
velocity component, is integrated in the commercial evaluation software of certain sodars
(Contini et al., 2009). The application of the VWV method is now also been facilitated by
the easy availability of small Doppler wind lidars.
Enhanced acoustic received echo (EARE) method The EARE algorithm has been pro-
posed by Emeis and Türk (2004) and Emeis et al. (2007). The method is an enhancement of
the ARE method in two ways. Firstly, it includes further variables into the MLH algorithm that
are available from Doppler-sodars. The benefits of the additional usage of the variance of the
vertical velocity component have been demonstrated by Emeis and Türk (2004). Secondly, it
determines not only MLH from sodar measurements but also the heights of additional lifted
inversions. Especially in orographically complex terrain, the vertical structure of the ABL can
be very complicated. Emeis et al. (2007) have shown that several persistent inversions one
above the other which form in deep Alpine valleys can be detected from sodar measurements
(Fig. 70).
EARE determines three different types of heights based on acoustic backscatter intensity
and the variance of the vertical velocity component. Because the horizontal wind information
above the inversion is not regularly available from sodar measurements, horizontal wind data
have not been included into this scheme. In the following a letter “H” and an attached
number will denote certain derived heights which are related to inversions and the MLH;
while the variable z is used to denote the normal vertical coordinate. The EARE algorithm
detects:
• the height (H1) of a turbulent layer characterised by high acoustic backscatter intensities
R(z) due to thermal fluctuations (therefore having a high variance of the vertical velocity
component σw ),
H3 = z, if (R(z) > 105 dB and R(z + 1) < 105 dB and σw (z) < 0.3m s−1 ) or
if (σw (z) < 0.3 m s−1 and σw (z + 1) > 0.3 m s−1 and R(z) > 105 dB).(141)
The σw values used in Eqs. (140) and (141) have been determined by optimizing the
automatic application of the detection algorithm. In doing so it turned out that no lifted
inversions occurred with a variance σw higher than 0.7 m s−1 and that the variance σw in
nocturnal stable surface layers was always below 0.3 m s−1 . The first σw threshold made it
possible to distinguish between inversions and elevated layers of enhanced turbulence. The
latter σw threshold made it possible to differentiate between nocturnal stable surface layers
and daytime super-adiabatic surface layers although both types of surface layers yield more
or less the same level of backscatter intensity. Finally MLH from the acoustic remote sensing
is determined as the minimum of H1, H21 , and H3:
Gradient or derivative methods Hayden et al. (1997) and Flamant et al. (1997) proposed
to use the largest negative peak of the first derivative of the optical attenuated backscatter
intensity (B(z)) for the detection of H4 from LIDAR data (height of gradient minimum
H4GM ):
H4n = z, if ∂ 2 B/∂z 2|z−1 ≤ 0 and ∂ 2 B/∂z 2 |z > 0 and B(z − ∆h/2) ≥ Bmin
and ∂B/∂z|z ≤ ∂B/∂zmax for n = 1, ..., 5. (148)
The MLH from optical remote sensing is taken as the lowest height H4n :
Wavelet method A Wavelet method has been developed for the automatic determination
of mixing layer height from backscatter profiles of an LD-40 ceilometer by de Haij et al.
(2006). Before that wavelet transforms have been applied in recent studies for MLH deter-
mination from LIDAR observations (Cohn and Angevine, 2000; Davis et al., 2000; Brooks,
2003; Wulfmeyer and Janjić, 2005). The most important advantage of wavelet methods is the
decomposition of the signal in both altitude as well as vertical spatial scale of the structures
in the backscatter signal.
The Wavelet algorithm in de Haij et al. (2006) is applied to the 10 minute averaged range
and overlap corrected backscatter profile B(z) within a vertical domain of 90–3000 m. For
each averaged profile the top of two significant aerosol layers are detected in order to detect
MLH as well as the top of a secondary aerosol layer, like e.g. an advected aerosol layer or
the residual layer. This Wavelet MLH method uses the scale averaged power spectrum profile
WB (z) of the wavelet transform with 24 dilations between 15 and 360 m and step size 15 m.
The top of the first layer, H41 , is detected at the first range gate at which the scale averaged
power spectrum WB (z) shows a local maximum, exceeding a threshold value of 0.1. This
threshold value is empirically chosen, based on the analysis of several cases with both well
pronounced and less clearly pronounced mixing layer tops. H42 is optionally determined in
the height range between H41 and the upper boundary of detection. A valid H42 is detected
at the level with the strongest local maximum of WB (z) provided that this maximum is larger
than the WB (z) of H41 . MLH is set equal to H41 .
However, problems with this method arise e.g. in case of multiple (well defined) aerosol
layers, which renders the selection of the correct mixing layer top ambiguous. Furthermore,
in spring and summer the detection of the MLH for deep (convective) boundary layers often
fails. This is mostly due to the high variability of the aerosol backscatter signal with height
which limits the range for MLH estimation in those conditions (de Haij et al., 2006).
Variance method At the top of the CBL we have entrainment of clear air masses from
the free troposphere into the ABL. The entrainment process is temporarily variable and leads
locally to considerable fluctuations in the aerosol concentration. Therefore the maximum in
the vertical profile of the variance of the optical backscatter intensity can be an indicator for
an entrainment layer on top of a CBL (Hooper and Eloranta, 1986; Piironen and Eloranta,
1995). The method is called variance centroid method in Menut et al. (1999). The variance
method for the CBL height is also described in Lammert and Bösenberg (2006). Due to
the assumptions made this method is suitable for daytime convective boundary layers only.
An elucidating comparison between the gradient method and the variance method can be
found in Martucci et al. (2004) although they used a Nd:YAG LIDAR at 532 nm instead of
a ceilometer and thus suffered from a high lowest range gate in the order of 300 m.
6.2.3 RASS
The acoustic and optical methods for MLH determination, which have been described in the
sections above, are all indirect methods that try to infer the mixing-layer height from other
variables which usually adapt to the vertical structure of the ABL. The only direct and key
variable for the analysis of the presence of a mixing layer is the vertical profile of virtual
temperature. Temperature profiles can directly be measured with a radio-acoustic sounding
system (RASS). Fig. 72 shows an example. We start here with a short description of the
available RASS methods.
a Doppler-RASS.
ca = cs + w. (153)
Figure 75: Sodar-RASS. The acoustic antenna is in the middle, the electro-magnetic antennas
to the left and right.
The vertical air speed component w can be determined separately from the Doppler shift
of the backscattered electro-magnetic clear-air signal when operating a Bragg-RASS or from
the Doppler shift of the backscattered acoustic signal when operating a Doppler-RASS. Using
the definition of the acoustic temperature the height profile of cs can then be converted in a
height profile of the acoustic temperature Ta . For many purposes this acoustic temperature
is a sufficiently accurate approximation of the virtual air temperature. The maximum range
of a sodar-RASS is usually less than thousand metres so that such an instrument covers the
lower part of the boundary layer. Due to the high vertical resolution and the low minimum
range of about 30 to 40 m, which is comparable to the abilities of a sodar, a sodar-RASS is
well suited for the detection of shallow nocturnal boundary layers.
MLH detection algorithms MLH can be determined from the lowest height where the
vertical profile of potential temperature increases with height indicating stable thermal strat-
ification of the air. The great advantage of RASS measurements is that the magnitude of
6.4.1 Formation
The formation of low-level jets requires a temporal or spatial change in the thermal stability
of the atmospheric boundary layer which leads to a sudden change between two different
equilibria of forces. The flow must transit from an unstable or neutral condition where fric-
tion, pressure-gradient and Coriolis forces balance each other to a stable condition where only
pressure-gradient and Coriolis force balance each other (see Fig. 78). The sudden disappear-
ance of the retarding friction in the equilibrium of forces leads to an inertial oscillation of
the horizontal wind vector. Wind speed shoots to much higher values and the increased wind
speed leads to a stronger Coriolis force which provokes a turning of the wind vector as well.
The relevant equations for the wind components u and v at the moment of the disappearance
of the frictional force read:
∂u
= −f (v − vg ) (155)
∂t
∂v
= f (u − ug ) (156)
∂t
where f is the Coriolis parameter and ug and vg are the components of the geostrophic wind.
The terms on the left-hand side involve a dependence on time. Therefore, the analytical
solution of Eqs. 155 and 156 describes an oscillation with time, t:
u − ug = Dv sin f t + Du cos f t (157)
v − vg = Dv cos f t − Du sin f t (158)
where Du and Dv are the ageostrophic wind components at the beginning of the oscillation
in the moment when the retarding frictional force vanishes.
In the temporal domain this corresponds to a sudden change from an unstable daytime
convective boundary layer to a nocturnal stable boundary layer. This requires clear skies
in order to have rapid changes in thermal stratification but still non-vanishing horizontal
synoptic pressure gradients. Therefore, nocturnal low-level jets usually appear at the edges of
high-pressure systems (see shaded area in Fig. 79).
In the spatial domain this corresponds to a sudden transition of the flow from a surface
which is warmer than the air temperature to a smooth surface which is colder than the air
temperature. This may happen when the flow crosses the coast line from warm land to a
colder ocean surface or from bare land to snow or ice-covered surfaces.
Figure 80: Annual frequency of low-level jet events ordered by large-scale synoptic weather
patterns (Grosswetterlagen). From two years of sodar data for Hannover, Germany.
the frequency of occurrence of low-level jets over Northern Germany as function of these 29
large-scale weather types. The two most relevant types (the two left-most columns in Fig. 80)
are a high-pressure bridge over Central Europe (type “BM”) and a high-pressure area over
the British Isles (type “HB”). All in all a low-level jet appeared in 23% of all nights
Figure 80 showed the frequency of occurrence of a low-level jet as function of the weather
type. The relevance of a certain weather type for the formation of a low-level jet can be
assessed when comparing the frequency of low-level jet occurrence with the overall frequency
of occurrence of the respective weather type. Figure 81 has been produced by dividing the fre-
quencies shown in Fig. 80 by the occurrence frequency of the respective weather types during
the same observation period. There are two weather types where the occurrence frequency
is identical to the occurrence frequency of the low-level jets during this weather type. This
means that in every night when this weather type prevailed a low-level jet was observed. This
is indicated by a low-level jet efficiency of 1.0 in Fig. 81. Small deviations from unity are due
to the limited sample size evaluated for this purpose. These two weather types are “HNFA”
and “HFZ” which are both related to high-pressure systems to the North of the investigation
site
Such a high efficiency for forming a low-level jet allows for a quite certain forecast of the
occurrence of a low-level jet. Once such weather types are forecasted a low-level jet will form
with a very high probability. The values given in Fig. 81 can be used to give the low-level jet
Figure 82: Correlation between the height of the core of a low-level jet and its maximum wind
speed over Northern Germany (from two years of sodar measurements).
formation probability for Northern Germany for each of the weather types. For other areas
the investigation has to be repeated with local low-level jet data.
The height of the core of a low-level jet and its maximum wind speed seems to be correlated.
Figure 82 shows a correlation analysis from two years of sodar measurements in Northern
Germany. Weaker jet cores may appear at heights between 150 m and 200 m above ground
while stronger events (20 m s−1 core wind speed) usually appear at about 400 m with
considerable scatter. Typical wind shear values for the layer between the surface and 160 m
above ground vary between 0.04 and 0.10 s−1 . An analysis of daily wind speed amplitudes
at 160 m above ground from these sodar data shows that amplitudes of up to 14 m s−1 are
possible due to low-level jet events.
6.5 Summary
Wind resources depend on the large-scale weather conditions as well as on the local vertical
structure of the atmospheric boundary layer. Ground-based remote sensing is now a viable
technique to monitor the vertical structure of the atmospheric boundary layer. Three different
techniques are presently available: acoustic sounding (sodars), optical sounding (wind lidars
and ceilometers) and the combination of acoustic and electro-magnetic sounding (RASS).
Notation
a factor of proportionality
A1 (BM + BU ) /2
A2 (BM − BU ) /2
ABL atmospheric boundary layer
ARE acoustic received echo (method)
Bmin threshold value for B(z)
BM mixing-layer mean of B(z)
BU value of B(z) above the mixing layer
B(z) optical backscatter intensity
Bi (z) idealized optical backscatter intensity
∂B/∂zmax threshold value for the vertical derivate of B(z)
c speed of light
ca RASS speed of sound (= cs + w)
cs true speed of sound
CBL convective boundary layer
Du initial ageostrophic wind component
Dv initial ageostrophic wind component
DR1 threshold value for the vertical gradient of R(z)
DR2 threshold value for the vertical gradient of R(z)
EARE enhanced acoustic received echo (method)
f Coriolis parameter
fa acoustic frequency
fe electro-magnetic frequency
Hn analysed height (n is a one-digit number)
H4GM height of minimum of vertical gradient of B(z)
H4IPM height of minimum of second vertical derivate of B(z)
H4LGM height of minimum of logarithmic vertical gradient of B(z)
HWS horizontal wind speed (method)
MLH mixing-layer height
MLHac mixing-layer height from acoustic sounding
MLHop mixing-layer height from optical sounding
q specific humidity
Rc threshold value for R(z)
R(z) acoustic backscatter intensity
RASS radio-acoustic sounding system
t time
Ta acoustic temperature
u horizontal wind component (Eastwards)
ug horizontal component of geostrophic wind (Eastwards)
v horizontal wind component (Northwards)
vg horizontal component of geostrophic wind (Northwards)
VWV vertical wind variance (method)
w vertical wind component
WB (z) scale averaged power spectrum profile
z height above ground
zi boundary-layer height
∆fe Doppler shift of electro-magnetic frequency
∆h height interval
Θ potential temperature
Θa acoustic potential temperature
Θv virtual potential temperature
λa wavelength of sound wave
References
Angevine W., White A. B., and Avery S. K. (1994) Boundary layer depth and entrainment zone characterization
with a boundary layer profiler. Bound.-Layer Meteorol. 68: 375–385
Asimakopoulos D. N., Helmis C. G., and Michopoulos J (2004) Evaluation of SODAR methods for the
determination of the atmospheric boundary layer mixing height. Meteorol. Atmos. Phys. 85: 85–92
Beyrich F. (1995) Mixing height estimation in the convective boundary layer using sodar data. Bound.-Layer
Meteorol. 74: 1–18
Beyrich F. (1997) Mixing height estimation from sodar data - a critical discussion. Atmosph. Environ. 31:
3941–3954
Beyrich F. and Görsdorf U. (1995) Composing the diurnal cycle of mixing height from simultaneous SODAR
and Wind profiler measurements. Bound.-Layer Meteorol. 76: 387–394
Boers R., Spinhirne J. D., and Hart W. D. (1988) Lidar Observations of the Fine-Scale Variability of Marine
Stratocumulus Clouds. J. Appl. Meteorol. 27: 797–810
Bösenberg J. and Linné H. (2002) Laser remote sensing of the planetary boundary layer. Meteorol. Z. 11:
233–240
Brooks I. M. (2003) Finding boundary layer top: application of a wavelet covariance transform to lidar backscat-
ter profiles. J. Atmos. Oceanic Technol. 20: 1092–1105
Cohn S. A. and Angevine W. M. (2000) Boundary Layer Height and Entrainment Zone Thickness Measured
by Lidars and Wind-Profiling Radars. J. Appl. Meteorol. 39: 1233–1247
Contini D., Cava D., Martano P., Donateo A., and Grasso F. M. (2009) Comparison of indirect methods for
the estimation of Boundary Layer height over flat-terrain in a coastal site. Meteorol. Z. 18: 309–320
Davis K. J., Gamage N., Hagelberg C. R., Kiemle C., Lenschow D. H., and Sullivan P. P. (2000) An objective
method for deriving atmospheric structure from airborne lidar observations. J. Atmos. Oceanic Technol.
17: 1455–1468
de Haij M., Wauben W., and Klein Baltink H. (2006) Determination of mixing layer height from ceilometer
backscatter profiles. In: Slusser J. R., Schäfer K., Comerón A. (eds) Remote Sensing of Clouds and the
Atmosphere XI. Proc. SPIE Vol 6362: paper 63620R.
Devara P. C. S., Ernest Ray P., Murthy B. S., Pandithurai G., Sharma S., and Vernekar K. G. (1995)
Intercomparison of Nocturnal Lower-Atmospheric Structure Observed with LIDAR and SODAR Techniques
at Pune, India. J. Appl. Meteorol. 34: 1375–1383
Emeis S. (2010) Measurement methods in atmospheric sciences. In situ and remote. Vol. 1 of Quantifying
the Environment. Borntraeger, Stuttgart, XIV+257 pp., 103 Figs, 28 Tab. ISBN 978-3-443-01066-9.
Emeis S. (2011) Surface-Based Remote Sensing of the Atmospheric Boundary Layer. Series: Atmospheric and
Oceanographic Sciences Library, Vol. 40. Springer Heidelberg etc., X+174 pp. 114 illus., 57 in color., H/C.
ISBN: 978-90-481-9339-4, e-ISBN 978-90-481-9340-0, ISSN 1383-8601, DOI: 10.1007/978-90-481-9340-0
Emeis S. (2012) Wind Energy Meteorology - Atmospheric Physics for Wind Power Generation. Springer
Heidelberg etc., XIV-196 pp. ISBN 978-3-642-30522-1, e-book: ISBN 978-3-642-30523-8
Emeis S. and Schäfer K. (2006) Remote sensing methods to investigate boundary-layer structures relevant to
air pollution in cities. Bound.-Layer Meteorol. 121: 377–385
Emeis S. and Türk M. (2004) Frequency distributions of the mixing height over an urban area from SODAR
data. Meteorol. Z. 13: 361–367
Emeis S., Jahn C., Münkel C., Münsterer C., and Schäfer K (2007) Multiple atmospheric layering and mixing-
layer height in the Inn valley observed by remote sensing. Meteorol. Z. 16: 415–424
Emeis S., Schäfer K., and Münkel C (2008) Surface-based remote sensing of the mixing-layer height - a review.
Meteorol. Z. 17: 621–630
Emeis S., Schäfer K., and Münkel C (2009) Observation of the structure of the urban boundary layer with
different ceilometers and validation by RASS data. Meteorol. Z. 18: 149–154
Eresmaa N., Karppinen A., Joffre S. M., Räsänen J., and Talvitie H. (2006) Mixing height determination by
ceilometer. Atmos. Chem. Phys. 6: 1485–1493
Flamant C., Pelon J., Flamant P. H., and Durand P. (1997) Lidar determination of the entrainement zone
thickness at the top of the unstable marin atmospheric boundary-layer. Bound.-Layer Meteorol. 83: 247–284
Gerstengarbe F.-W., Werner P. C., and Rüge U. (Eds.) (1999) Katalog der Großwetterla-
gen Europas (1881 - 1998). Nach Paul Hess und Helmuth Brezowsky. 5th edition. Ger-
man Meteorological Service, Potsdam/Offenbach a. M. (available from: http://www.deutscher-
wetterdienst.de/lexikon/download.php?file=Grosswetterlage.pdf or http://www.pik-potsdam.de/ uw-
erner/gwl/gwl.pdf)
Grimsdell A. W. and Angevine W. M. (1998) Convective Boundary Layer Height Measurement with Wind
Profilers and Comparison to Cloud Base. J. Atmos. Oceanic Technol. 15: 1331–1338
Gryning S.-E., Batchvarova E., Brümmer B., Jørgensen H., and Larsen S. (2007) On the extension of the
wind profile over homogeneous terrain beyond the surface layer. Bound.-Layer Meteorol. 124:251–268
Hayden K. L., Anlauf K. G., Hoff R. M., Strapp J. W., Bottenheim J. W., Wiebe H. A., Froude F. A., Martin
J. B., Steyn D. G., and McKendry I. G. (1997) The Vertical Chemical and Meteorological Structure of the
Boundary Layer in the Lower Fraser Valley during Pacific ’93. J. Atmos. Environ. 31: 2089–2105
Hennemuth B. and Kirtzel H.-J. (2008) Towards operational determination of boundary layer height using
sodar/RASS soundings and surface heat flux data. Meteorol. Z. 17: 283–296
Hooper W. P. and Eloranta E. (1986) Lidar measurements of wind in the planetary boundary layer: the method,
accuracy and results from joint measurements with radiosonde and kytoon. J. Clim. Appl. Meteorol. 25:
990–1001
Lammert A. and Bösenberg J. (2006) Determination of the Convective Boundary-Layer Height with Laser
Remote Sensing. Bound.-Layer Meteorol. 119: 159–170
Marshall J. M., Peterson A. M., and Barnes A. A. (1972) Combined radar-acoustic sounding system. Appl.
Opt. 11: 108–112
Martucci G., Srivastava M. K., Mitev V., Matthey R., Frioud M., and Richner H. (2004) Comparison of lidar
methods to determine the Aerosol Mixed Layer top. In: Schäfer K., Comeron A., Carleer M., Picard R. H.
(eds.): Remote Sensing of Clouds and the Atmosphere VIII. Proc of SPIE 5235: 447–456
Melfi S. H., Spinhirne J. D., Chou S. H., and Palm S. P. (1985) Lidar Observation of the Vertically Organized
Convection in the Planetary Boundary Layer Over the Ocean. J. Clim. Appl. Meteorol. 24: 806–821
Menut L., Flamant C., Pelon J., and Flamant P. H. (1999) Urban Boundary-Layer Height Determination from
Lidar Measurements Over the Paris Area. Appl. Opt. 38: 945–954
Münkel C. and Räsänen J. (2004) New optical concept for commercial lidar ceilometers scanning the boundary
layer. Proc. SPIE Vol 5571: 364–374
Neff W. D. and Coulter R. L. (1986) Acoustic remote sensing. In: Lenschow D. H. (ed) Probing the Atmospheric
Boundary Layer. Amer. Meteorol. Soc., Boston, MA, 201–239
Peña A., Gryning S.-E., and Hasager C. B. (2010) Comparing mixing-length models of the diabatic wind
profile over homogeneous terrain. Theor. Appl. Climatol. 100:325–335
Piironen A. K. and Eloranta E. W. (1995) Convective boundary layer depths and cloud geometrical properties
obtained from volume imaging lidar data. J. Geophys. Res. 100: 25569–25576
Schäfer K., Emeis S. M., Rauch A., Münkel C., and Vogt S. (2004) Determination of mixing-layer heights
from ceilometer data. In: Schäfer K., Comeron A. T., Carleer M. R., Picard R. H., and Sifakis N. (eds.):
Remote Sensing of Clouds and the Atmosphere IX. Proc. SPIE Vol 5571: 248–259
Schäfer K., Emeis S., Junkermann W., and Münkel C. (2005) Evaluation of mixing layer height monitoring
by ceilometer with SODAR and microlight aircraft measurements. In: Schäfer K., Comeron A. T., Slusser
J. R., Picard R. H., Carleer M. R., and Sifakis N. (eds) Remote Sensing of Clouds and the Atmosphere X.
Proc. SPIE Vol 5979: 59791I-1–59791I-11
Seibert P., Beyrich F., Gryning S.-E., Joffre S., Rasmussen A., and Tercier P. (2000) Review and intercompar-
ison of operational methods for the determination of the mixing height. Atmosph. Environ. 34: 1001-1-027
Sicard M., Pérez C., Rocadenbosch F., Baldasano J. M., and Garcı́a-Vizcaino D. (2006) Mixed-Layer Depth
Determination in the Barcelona Coastal Area From Regular Lidar Measurements: Methods, Results and
Limitations. Bound.-Layer Meteorol. 119: 135–157
Steyn D. G., Baldi M., and Hoff R. M. (1999) The detection of mixed layer depth and entrainment zone
thickness from lidar backscatter profiles. J. Atmos. Ocean Technol. 16: 953–959
White A. B., Senff C. J., and Banta R. M. (1999) A Comparison of Mixing Depths Observed by Ground-Based
Wind Profilers and an Airborne Lidar. J. Atmos. Oceanic Technol. 16: 584–590
Wulfmeyer V. (1999) Investigation of turbulent processes in the lower troposphere with water-vapor DIAL
and Radar-RASS. J. Atmos. Sci. 56: 1055–1076
Wulfmeyer V. and Janjić T. (2005) 24-h observations of the marine boundary layer using ship-borne NOAA
High-Resolution Doppler Lidar. J. Appl. Meteorol. 44: 1723–1744
7.1 Introduction
Power performance measurement is central to the wind industry since it forms the basis for
the power production warranty of the wind turbine. The power curve measurement has to
be realised according to the IEC 61400-12-1 standard. The power curve is obtained with
10-min mean power output from the turbine plotted against simultaneous 10-min average
wind speeds. The standard requires the wind speed to be measured by a cup anemometer
mounted on top of a mast having the same height as the turbine hub and located at a distance
equivalent to 2.5 rotor diameters from the turbine.
Such a plot usually shows a significant spread of values and not a uniquely defined function.
The origin of the scatter can mainly be grouped into three categories: the wind turbine com-
ponents characteristics, sensor error and the wind characteristics. Within the last group, the
current standard only requires the wind speed at hub height and the air density measurement.
However, other wind characteristics can influence the power production like the variation of
the wind speed with height (i.e. wind speed shear). The influence of wind speed shear on
the power performance was shown in several studies: some based on aerodynamic simulations
(Antoniou, 2009; Wagner et al., 2009) others based on measurements (Elliot and Cadogan,
1990; Sumner an Masson, 2006).
A major issue is to experimentally evaluate the wind speed shear. The wind speed profile
is usually assumed to follow one of the standard models such as the logarithmic or power law
profiles. However, these models are valid for some particular meteorological conditions, and
therefore, cannot represent all the profiles experienced by a wind turbine. Measurements are
then a better option but are also challenging. Indeed characterising the wind speed profile
in front of the rotor of a multi-MW wind turbine requires measurements of wind speed at
several heights, including some above hub height, i.e. typically above 100 m. Remote sensing
instruments such as lidar or sodar then appear as a very attractive solution.
This chapter starts with a description of the influence of the wind speed shear on the power
performance of a multi-MW turbine. The challenge of describing the wind speed profile is
then discussed followed by a description of an experiment using a lidar for its characterisation.
This is followed by the introduction of the definition of an equivalent wind speed taking the
wind shear into account resulting in an improvement of the power performance measurement.
Finally, some recommendations about remote sensing instruments are given to successfully
apply this method.
height @mD
140
120
100
80
60
40
wind speed @msD
4 6 8 10
Figure 83: Wind profiles used as input for the wind speed shear aerodynamic investigation.
Black curve: no shear; grey curve: power law profile with shear exponent of 0.5. Both profiles
have the same wind speed at hub height
Figure 84 shows the free wind speed (i.e. the absolute wind speed in absence of a turbine)
seen by a point at a radius of 30 m from the rotor centre, rotating at the same speed as
the rotor as a function of time for the 2 inflow cases. Whereas in a uniform flow the blade is
subjected to a constant wind speed, in a sheared flow, the point is exposed to large variations
of wind speed even though the inflow is laminar. The variation of the wind speed seen by this
rotating point in time is only due to the fact that it is rotating within a non uniform flow
(wind speed varying with height).
Figure 85 shows the variations of the free wind speed seen by the same rotating point as
function of the azimuth position (0◦ = downwards). The point experiences the hub height
wind speed (same as uniform inflow) when it is horizontal (±90◦ ), lower wind speed when it
is downward (0◦ ) and higher wind speed when it is upward (180◦ ).
A rotating blade does not experience the free wind speed because of the induction from the
drag of the rotor. In reality, a rotating blade is directly subjected to the relative wind speed w
(i.e. the speed of the wind passing over the airfoil relative to the rotating blade) and the angle
of attack φ (i.e. the angle between the blade chord line and the relative wind speed) with the
effects of the induced speed included. The variations of these two parameters as function of θ
are shown in Figure 86. As these two parameters directly depend on the free wind speed, they
vary with the azimuth angle in a sheared inflow, whereas they remain constant in a uniform
inflow.
The relative speed and the angle of attack are derived from the rotor speed and the induced
velocity. Therefore, they depend on the way that the induction is modeled and it is difficult
9.0
8.5
8.0
7.5
7.0
time@sD
120 140 160 180 200
Figure 84: Time series of free wind speed seen from a rotating point, positioned at a radius
of 30 m, rotating at rotor rotational speed (no induced velocity). Black curve: no shear; grey
curve: power law profile with shear exponent of 0.5
u y @msD
9.0
8.5
8.0
7.5
7.0
Θ @°D
-150 -100 -50 50 100 150
Figure 85: Free wind speed seen from a rotating point, positioned at a radius of 3 0m, rotating
at rotor rotational speed, as function of the azimuth angle θ. Black curve: no shear; grey curve:
power law profile with shear exponent of 0.5
to evaluate their variations due to a non uniform flow in a simple way. However, some basic
considerations (ignoring the induction) can give a basic insight to the variation of the relative
speed and the angle of attack as the blade rotates. In case of uniform inflow, the free wind
speed is the same at any point of the swept rotor area. Therefore, the angle of attack and
the relative speed are the same at any azimuthal position (see Figure 87).
In case of sheared inflow, the free wind speed depends on the position of the blade. When
the blade is horizontal, the free wind speed is the speed at hub height and the speed triangle
is the same as in Figure 87. Below hub height, the wind speed is lower than the hub height
speed, see Figure 88 (left). Consequently w and φ are lower than those at hub height. Above
hub height, the wind speed is higher than the hub height wind speed. Consequently w and φ
are higher than those at hub height, see Figure 88 (right). The variations in φ and w cause a
variation of the local lift and drag as the blade rotates, which finally results in the variation of
the local tangential force contributing to the wind turbine power (see Figure 89). For a given
φ, local lift dFL and local drag dFD , the local tangential force dFT is given by (Manwell
4.0 44.5
3.5 44.0
3.0 43.5
2.5 43.0
Θ @°D Θ @°D
-150 -100 -50 50 100 150 -150 -100 -50 0 50 100 150
Figure 86: phi (left) and w (including induction) (right) as a function of θ, seen from a point
at radius r = 30 m on a rotating blade. Black curve: no shear; grey curve: power law profile
with shear exponent of 0.5
rotation
rW
Φ
u
w
Figure 87: Speed triangle for a blade element at radius r. rΩ is the blade speed and w
corresponds to the sum of the pitch angle, the twist angle and φ. As the twist angle is
constant for a given position on the blade and the pitch angle is 0◦ for wind speeds below
rated speed, φ represents the variation of the angle of attack
et al., 2002):
dFT = dFL sin φ − dFD cos φ. (160)
As the wind speed increases with height (e.g. in the case of the power law profile), the
amplitude of the variations of the free wind speed seen by a rotating point increases with the
radius (not shown here). The local tangential force consequently varies with the radius too.
As the torque results from the integral of the tangential force over the whole rotor, it thus
depends on the wind speed profile.
Figure 88: Speed triangles at the top and bottom of the rotor showing the effect of wind
speed shear. The speed triangle at hub height is shown with dashed arrows
dFT @kNmD
0.18
0.16
0.14
0.12
0.10
0.08
0.06
Θ @°D
-150 -100 -50 50 100 150
Figure 89: Local tangential force seen from a point r = 30 m on a rotating blade as function
of θ. Black curve: no shear; Gray curve: power law profile with shear exponent of 0.5
where ρ is the air density, R the rotor radius and c is the chord (of the circle defined by the
rotor swept area) as function of z which varies between the bottom and the top of the rotor:
p
c = 2 2Rz − z 2 . (162)
In order to compare to the power output variations, Figure 90 also shows the difference
between the kinetic energy flux for a power law profile and a constant profile, normalised with
the power obtained with a constant profile KEprofile − KEhub /KEhub (see gray dashed line)
5
. Figure 90 shows two other interesting results:
1. The kinetic energy flux varies with shear exponent.
2. The power output of the turbine does not follow the same trend as the kinetic energy
flux.
Despite the high uncertainty in the modeled power output for a sheared inflow, the results
highlight that the influence of the shear on the power performance of a turbine can be seen
as the combination of two effects:
• The variation in kinetic energy flux (power input).
• The ability of the turbine to extract the energy from the wind, which depends on the
details of the design and the control strategy of the turbine.
5 The kinetic energy flux is not an output of HAWC2Aero, it has been here estimated and for power law
profiles, the normalized difference does not depend on the wind speed
-10
Figure 90: Normalised difference in power output and kinetic energy flux between shear and
uniform case as function of the shear exponent, for various wind speed at hub height
These results clearly indicate that wind speed profiles encountered by the turbine during a
power performance measurement should be known and taken into account.
æ æ
1.4 1.4
model: Α= 0.33 model: Α= 0.4
mast measurement æ mast measurement æ
Height HeightHub
Height HeightHub
1.2 1.2
1.0 æ 1.0 æ
0.8 0.8
æ æ
0.6 0.6
æ æ
0.4 0.4
0.6 0.7 0.8 0.9 1.0 1.1 1.2 0.6 0.7 0.8 0.9 1.0 1.1 1.2
U Urated U Urated
Figure 91: Two examples of profiles and their fit to the power law model (using the wind
speed values at 40 and 80 m). The model fits very well the measured profile on the left, but
it does not work for the profile on the right
10
% Hprofiles numberL
8
6
4
2
0
-0.4 -0.2 0.0 0.2 0.4
wsp error at 1.45 hub height
Figure 92: Distribution of the error made on the wind speed at 116.5 m height when assuming
a power law profile with a shear exponent estimated with the wind speeds at 40 and 80 m
evaluated the goodness of fit with the residual sum of squares RSS, defined as:
X 2
RSS = (uf it (zi ) − ui ) , (164)
i
where uf it is the fit function defined in Eq. (163) and ui the wind speed measured by the
lidar at height zi (i = 1 to 9). Figure 93 shows two examples of measured profiles with their
shear exponents and RSSs.
The profiles were then divided into two groups according to the RSS:
• Group 1: RSS ≤ 0.1 – the profiles from this group have a shape close to a power law
profile and are referred to as the power law profiles.
• Group 2: RSS > 0.1 – the profiles from this group have a shape that cannot be well
represented by a power law profile and are referred to as the non power law profiles.