Light Detection And Ranging

(LIDAR)

Aman Raj
Regd.No.: 1501289117

Department of Computer Science &

Engineering Trident Academy of
Technology Bhubaneswar-751024,
Odisha, India.

October 2018
 Seminar Report on

Light Detection And Ranging

Submitted in Partial Fulfilment of the
Requirement for the Award of the Degree of

Bachelor of
Technology in
Computer Science & Engineering

Submitted by

Aman Raj
Regd.No.: 1501289117

Under the Guidance of

Mr. Jyoti Ranjan Sahoo

Professor, Dept. of CSE

Department of Computer Science & Engineering

Trident Academy of Technology
Bhubaneswar-751024, Odisha, India.

Oct0ber 2018
 CERTIFICATE

This is to certify that this Seminar Report on the topic entitled ) Light
Detection And Ranging which is submitted by Aman Raj in partial
fulfillment of the requirement for the award of the of Bachelor of
Technology in Computer Science & Engineering of Biju Patnaik
University of Technology, Odisha, is a record of the candidate's own
work carried out by him under my supervision.

Supervisor Head of Deartment

(Mr. Jyoti Ranjan Sahoo) (Dr.Abhay Kumar Samal)

Professor, Dept. of CSE Asso.Prof., Dept. of CSE
Trident Academy of Technology Trident Academy of Technology
Bhubaneswar, India. Bhubaneswar-751024, Odisha, India.

i
 ABSTRACT

LIGHT DETECTION AND RANGING

LiDAR is an acronym for Light Detection And Ranging, sometimes also

referred to as Laser Altimetry or Airborne Laser Terrain Mapping (ALTM).
The LiDAR system basically consists of integration of three technologies,
namely, Inertial Navigation System (INS),LASER, and GPS. The Global
Positioning System (GPS) has been fully operational for over a decade, and
during this period, the technology has proved its potential in various
application areas. Some of the important applications of GPS are crustal
deformation studies, vehicle guidance systems, and more recently, in
LiDAR. Geo Spatial Information is an important input for all planning and
developmental activities especially in the present era of digital mapping
and decision support systems. LiDAR is much faster than conventional
photogrammetric technology and offers distinct advantage over
photogrammetry in some application areas. Its development goes back to
1970s and 1980s, with the introduction of the early NASA-LiDAR systems,
and other attempts in USA and Canada (Ackermann, 1999). The method
has successfully established itself as an important data collection
technique, within a few years, and quickly spread into practical
applications. Early 1980's, second generation LiDAR systems were in use
around the world but were expensive and had limited capability. With the
enhanced computer power available today, and with the latest positioning
and orientation systems, LiDAR systems have become a commercially
viable alternative for development of Digital Elevation Models (DEM) of earth
surface.

AMAN RAJ
ROLL NO:11
REGD NO:1501289117
BRANCH:CSE
 ACKNOWLEDGMENTS

I take this opportunity to express my gratitude to the people who have been instrumental
in the successful completion of this project. I am, in the first place, obliged and
grateful to my parents without whose support and care I could not have completed
this report. I express my deep gratitude towards my guide, Mr. Jyoti Ranjan Sahoo
,Professor, Dept. of CSE,Trident Academy of Technology, Bhubaneswar, for his
tremendous sup- port,encouragement and help.
I convey my sincere thanks to Dr. Abhay Kumar samal, Asso.Prof., Dept. of CSE, De-
partment of Computer Science & Engineering, Trident Academy of Technology, Bhubaneswar,
for his cooperation in the completion of the project.
I would like to extend my gratitute to the Department of Computer Science & Engineer-
ing, Trident Academy of Technology, Bhubaneswar, for their support and cooperation.
Finally, I extend my appreciation to all my friends, teaching and non teaching sta , who
directly or indirectly help me in this endeavour.

Place: Bhubaneswar Aman Raj

Date: October 20, 2018 1501289117

iii
 CONTENTS

i. Certificate

ii. Abstract

iii. Acknowledgements

Contents

1 Introduction

2 General History of Lidar

3 Lidar Overview

 Lidar System Components

 Lidar Fundamentals
 Lidar Equations

4 Various Lidar Systems

 Raman Lidar
 Differential-absorption Lidar (DIAL)
 Resonance Lidar
 Doppeler Lidar

5 The Qinetiq ZephIR CW Doppler Lidar

 ZephIR Overview
 ZephIR PhotoDetection
 ZephIR Fourier Analysis
 ZephIR Cloud-Correction Algorithm
 ZephIR wind Velocity Estimation
 ZephIR and Range-Resolved Measurements
6 Conclusion

References
 CHAPTER 1

Introduction
The Renewable Energy Research Laboratory’s (RERL) Light Detection And Ranging (LIDAR)
system is capable of remotely measuring wind speed and direction by the technique of laser
remote sensing. This system merges established laser technology with more affordable internal
components to make it available for commercial use. Manufactured by Qinetiq of England, the
instrument is specifically designed for wind energy resource assessment applications. The lidar
offers great promise in terms of its ability to provide wind speed data at the hub height of a
modern wind turbine. This instrument is also attractive because it is small and capable of
deployment by a team of only two people.

Another remote wind speed measurement alternative is Sound Detection And Ranging
(SODAR). This measurement technique uses the same basic measurement principle as lidar.
However, sodar wind speed measurement is based on the analysis of acoustic signals rather than
laser radiation backscatter. The sound waves that are emitted by the sodar come in the form of
pulsed “chirps” that are subject to echo interactions from nearby structures or trees that can
corrupt the data [1], [2]. Sodars must also estimate the true value of the speed of sound at a given
temperature in real time. Since local temperatures can vary appreciably with time and height,
sodar wind speed measurement is somewhat more complex than that of the lidar.

Before the lidar is dispatched for data collection, its ability to accurately measure wind speed and
direction is first explored in the lidar data validation measurement campaign. To evaluate the
validity of the measurements that are collected by the lidar, experimental data are compared to a
control. This report contains an analysis of the procedure by which the lidar is validated with
respect to cup anemometers that are installed on a nearby tower.

1
 CHAPTER 2

General History of Lidar

In 1930, E.H. Synge was the first to suggest that atmospheric

density measurements could be obtained by analyzing the light
return scatter obtained from searchlights that illuminate the sky
[4]. Early lidar systems, such as the type proposed by Synge,
operated in a biaxial configuration that allowed range-resolved
measurements. In a biaxial setup, the lidar detector is located some
distance (up to several kilometers) away from the point where light
is transmitted to the atmosphere. The receiver’s field-of-view
(FOV) can be scanned along the searchlight beam to obtain a
height profile of the scattered light’s intensity by applying simple
geometry. Figure 1 illustrates the biaxial setup as well as other
configurations that will be discussed later.

2
Chapter 2: General History of Lidar 3

Six years later, in 1936, Duclaux was able to acquire atmospheric

density measurements at an altitude of 3.4 km by applying the
method that Synge proposed [5]. Hulbert later extended this work
to obtain measurements at 28 km [6].
Further developments in lidar technology introduced the
monostatic configuration where the transmitter and receiver are
grouped together in a monostatic-coaxial or monostatic-biaxial
arrangement, shown above in Figure 1. This design improvement
allowed lidar systems to incorporate transmitters that pulse the
light source, thereby permitting the measurement of round-trip
time of flight of the scattered light pulse. In 1938 Bureau was the
first to use a pulsed, monostatic system to determine cloud base
heights which signaled the berth of range-resolved lidar
measurement techniques as we know them today [7].
With the flexibility that monostatic configurations gave to the
experimentalist in terms of obtaining vertically profiled
measurements, the next major advance came from the advent of the
modern laser in 1960 by Theodore H. Maiman [8]. Lidar
technology took yet another leap when the Q-switched, or giant
pulsed, laser appeared as a result of the work done by F.J.
McClung and R.W. Hellworth in 1962 [9]. The ability to transmit
electromagnetic radiation at specific wavelength and frequency
characteristics has inextricably linked the advance of the lidar to
laser technology developments. Smullins and Fiocco were the first
to incorporate modern laser technology in a lidar system when they
used a pulsed ruby laser to detect light that was scattered from the
surface of the moon and later, from the lower atmosphere [10],
[11].
 CHAPTER 3

Lidar Overview
In the study of lidar technology, it is important to understand the
concepts of basic lidar components in order to gain a better
understanding of the various types of lidar systems. This section
introduces the essential lidar system components that are included
in all forms of lidar instruments. Additionally, the fundamentals
of lidar wind speed measurements are presented as well as an
overview of the various types of lidar systems in use today.

3.1. Lidar System Components

In general, a lidar system consists of three main components: the
transmitter, the receiver and the detector. A simple block diagram
of these basic components is shown below in Figure 2.

Figure 2: Block Diagram of a Generic Lidar System [12]

4
Chapter 3: Lidar Overview 5

3.1.1. Lidar Transmitter

The transmitter includes the laser that is used to generate a
continuous or pulsed laser beam at a variety of wavelengths
ranging from the infrared through the visible and into the
ultraviolet [12]. The wide variety of wavelengths that are used in
lidar systems gives it a capacity to measure a number of
atmospheric variables. Many systems also incorporate a beam
expander in the transmitter module that can help reduce the beam
divergence and increase the beam diameter, which in turn
diminishes unwanted background return scatter that can add noise
to the return signal. This beam expander typically comes in the
form of a convex lens. Also, a portion of the laser beam is sampled
and used as a reference to which the backscatter signals can be
compared.
The laser wavelength is an important design characteristic that
must be considered for the application of lidar technology to wind
resource monitoring. Laser wavelengths longer than 1.4µm do not
penetrate the eye and cannot reach the retina. Thus, an important
design criterion for laser remote wind speed measurement is the
choice of wavelength that will not require special personnel or
equipment to safely install and operate [13]. Given the importance
of this requirement, available technology allows a choice of three
realistic wavebands centered around 1.5µm, 2µm and 10µm [14].
More delicate 10µm and 2µm systems require larger and more
expensive optical equipment. These systems are therefore
inappropriate for autonomous wind speed measurement purposes.
The
1.5µm systems have recently been developed to perform reliably
and accurately while incorporating fiber-optic components that are
inexpensive and widely available [15]. These traits in parallel with
the satisfaction of the eye safety requirement make the 1.5µm
waveband a suitable laser wavelength for wind power resource
monitoring applications.
 Chapter 3: Lidar Overview 6

3.1.2. Lidar Receiver

In the receiver, a telescope collects the photons that are scattered
by the body that is being measured and directs them to a
photodetector that converts the light into an electrical signal. The
size of the telescope plays an important role in the accuracy of the
lidar since the strength of the electrical signal depends on the
amount of light that can be collected by the telescope. Naturally,
the larger the optical telescope, the larger proportion of photons
that can be detected after scattering. The diameter of most lidar
telescopes range from approximately 10 cm to a few meters [12].
Smaller telescope diameters can be used when lower heights (less
than 150 m) are being probed because the intensity of light that is
returned at these heights is more substantial. However, the focal
ratio of the telescope can cause uncertainty in range-resolved
measurements1. Before the collected light is directed to the
detector, many lidars introduce spectral filtering based on
wavelength, polarization and/or range [12]. The simplest case of
spectral filtering involves an interference filter that transmits light
in a certain pass-band around the wavelength of interest while
discarding any signal that falls outside of this band.
Chapter 3: Lidar Overview 7

3.1.3. Lidar Detector

The detector is the system component that records the intensity of
the light that is collected by the receiver. The detector in various
lidar systems can record information about the return signal by
using either a photon counting method, analog signal detection or
coherent detection method. While each method has its advantages,
coherent detection is most commonly used for wind velocity
measurements. Coherent detection allows the frequency shift of
the return signal to be determined by a relatively straightforward
method that isolates the difference between the frequency of
emitted light and that of the backscattered light. Lidar systems that
employ coherent detection tend to be cheaper and more robust
because this method eliminates several sensitive and costly
components that are associated with photon counting.

3.2. Lidar Fundamentals

The fundamental principle that governs lidar operation is based on
various forms of light scattering. This section introduces and
describes each type of light scattering that pertains to lidar
measurement. The objective of this section is to introduce the
concepts needed to more completely understand lidar operation.
3.2.1. Rayleigh Scattering
Rayleigh scattering is one form of light scatter. It is defined as the
elastic scattering of light from particles that are very small
compared to the wavelength of the scattered radiation. Elastic
scattering occurs when there is no change in energy between the
incident light and the target molecule. In the context of lidar
operation, Rayleigh scattering is used as a synonym for molecular
scattering. The intensity of Rayleigh scattered light is proportional
to λ-4 and dominates the elastic-backscatter signals at short laser
wavelengths [16].
 Chapter 3: Lidar Overview 8

3.2.2. Raman Scattering

Conversely, Raman scattering is the inelastic scattering of light
where the energetic state of the molecule is changed and thus the
wavelength of scattered light is shifted as well.
1
The focal ratio is expressed as N = Df where f is the focal length of
the telescope lens and D is the
diameter of the entrance pupil; the ratio expresses the
diameter of the entrance pupil in terms of the effective focal
length of the lens.

In any instance of light scattering, the majority of the light

scatters elastically (Rayleigh scattering). However, a small
fraction of scattered light is scattered with optical frequencies
different from the frequency of the incident photons. The Raman
effect corresponds to the absorption and subsequent emission of a
photon via an intermediate electron state, having another energy
level.
3.2.3. Mie Scattering
Mie scattering is another form of light scatter. Described by
Gustav Mie [17], this form of light scattering is not limited to a
certain size of particles. Furthermore, Mie scattering is based on
the assumption that the behavior of scattered light is a result of
contact with spherical aerosols.
Chapter 3: Lidar Overview 9

3.2.4. Light Scatter for Wind Speed Measurement

Lidar-based wind speed measurements are often based on Raman,
or inelastic, light scattering although it is also common to find
lidar instruments that are based on Rayleigh scattering. Lidar
instruments that are designed for the wind energy sector are
typically based on the principle of Raman scattering. As such, the
lidar instrument analyzed in the lidar data validation experiment
(section 7) operates on the principle of Raman scattering to obtain
wind speed measurements by the detection of small changes in the
frequency of scattered light with respect to a reference beam with
the same frequency as the emitted light.

3.3. Lidar Equation

In its most general form, the detected lidar signal can be written as
[16]:
P(R) = KG(R)β(R)T(R)
Equation 1
Where the power P received from a distance R is made up of four
factors. The first factor, K, summarizes the performance of the
lidar system and the second, G(R), describes the range-dependent
measurement geometry. These two factors are completely
determined by the lidar setup and can thus be controlled by the
experimentalist. The information about the atmosphere, and thus all
of the measurable quantities, are contained in the last two factors of
the equation above. The term β(R) is the backscatter coefficient at
distance R. It stands for the ability of the atmosphere to scatter
light back into the direction from which it came. T(R) is the
transmission term and describes how much light is lost on the way
from the lidar to distance R and back again. Both β(R) and T(R)
are the subjects of investigation and are unknown to the
experimentalist. For further detail on the specific equations that
govern these five components, the reader is referred to [12] .
 CHAPTER 4

Various Lidar Systems

The lidar is a versatile instrument that can remotely measure a variety
of atmospheric properties. Accordingly, there are many different forms
of lidar systems that are available for a range of applications. This
section gives a brief overview of each of the most common types of
lidar systems.

4.1. Raman Lidar

A Raman lidar is a variation of a lidar system that is designed to detect
Raman scattering that results from the illumination of a target of interest
in the atmosphere. Today, Raman lidars are typically used to measure the
distribution of aerosols and other gaseous species in the atmosphere [18].
The Raman lidar technique harnesses the characteristics of the inelastic
scattering of light to measure data by detecting shifts in the wavelength
of scattered light. Because a much smaller proportion of light is scattered
inelastically, at shifted wavelengths, Raman lidars must be equipped with
very sensitive receivers. The receivers found in Raman lidar systems are
capable of detecting extremely small backscatter intensity levels. The
Raman lidar can be manipulated to measure the concentration of a wide
range of atmospheric molecules because the wavelength shift (caused by
Raman scattering) is different for distinct molecules.

10
Chapter 4: Various Lidar Systems 11

4.2. Differential-absorption Lidar (DIAL)

The DIAL technique is based on photon absorption by molecules in the
atmosphere. If a photon has exactly the right amount of energy to allow a
change in the energetic state of a molecule, then the photon is absorbed.
This characteristic can be applied to the detection of trace gases in the
atmosphere by selecting a transmission wavelength that corresponds to
the absorption line of the constituent of interest. A DIAL lidar transmits
two closely spaced wavelengths in tandem. One of these wavelengths is
known to correspond to the absorption line of the substance under
investigation while the other is emitted in the wing of the absorption line
where it will not be absorbed as strongly. During the transmission of these
two wavelengths into the atmosphere, the intensity of the light that
corresponds to the substance absorption line will be diminished. When the
lidar instrument detects the backscatter intensities, the concentration of
various substances can be determined.

4.3. Resonance Lidar

Resonant scattering is an elastic process that occurs when the energy of
the incident photon is equal to the energy of an allowed transition within
the atom of investigation. The process is elastic because when the atom
absorbs a photon, it also releases another photon with the same frequency
as the incident light. Because each atom has unique absorption
characteristics, the resonance lidar method can be applied to the
measurement of the concentration of a particular atom, ion or molecule
in the atmosphere.

4.4. Doppler Wind Lidar

The Doppler phenomenon relates to the frequency change of radiation as
perceived by an observer that is moving relative to the source of the
radiation. This effect, while most famously applied to sound waves, also
applies to electromagnetic waves. The measurement of a perceived
frequency shift is accomplished by the illumination of naturally occurring
aerosols that travel at approximately the same speed as the wind.
Examples of these aerosols include pollution particulates, pollen or dust.
 Chapter 4: Various Lidar Systems 11

The study of the wind velocity at approximately 0.1 to 60 m/s is

interesting in the field of wind power resource monitoring, but since the
frequency shift at these speeds relative to the speed of light corresponds
to a very small fraction, the measurements require extremely sensitive
equipment.

4.4.1.Doppler Wind Lidar: Continuous Wave or Pulsed Laser Operation

Doppler wind lidars typically employ the use of pulsed laser operation
because a larger amount of energy can be emitted in short pulses which
allow the time-of-flight to be measured, permitting operation at much
longer ranges. The duration of the pulses is typically on the order of a few
microseconds [19]. Because pulsed lidars emit powerful bursts of
radiation, they are subject to more stringent laser safety guidelines that
make the eye safety requirement more difficult to achieve [13].
With the introduction of fiber optics and telecommunication industry
components, the continuous wave (CW) coherent Doppler lidar has
become much more economical in recent years [15]. CW lidars are much
less complicated systems that can be used for wind measurements at
heights in the lower atmospheric boundary layer because scatter
intensities at shorter range are relatively strong. The CW Doppler lidar
emits a continuous beam of radiation that is sampled in discrete chunks
by a signal processor that is part of the detector. The measurement of
wind speeds at various heights is achieved by adjusting the laser focus
internally rather than by measuring the time of flight of the return signal
as is done in pulsed lidar systems. A disadvantage of CW systems is that
they are limited to sensing at a maximum range of approximately 200
meters because beam diffraction can cause measurement instability [14].

4.4.2 .Doppler Wind Lidar: Scanning Techniques

In order to measure horizontal wind speeds, the beam of the lidar must
be tilted from vertical by some angleθ . By tilting the lidar beam, the
horizontal and vertical wind speed components can be obtained by
purely geometric means. These wind speed components are extracted
from the line-of-sight (LOS) wind speed. The line-of-sight wind speed is
the vector component of the wind speed along the axis of the lidar laser
beam transmission path. It is common to refer to the
Chapter 4: Various Lidar Systems 11

line-of-sight velocity as vLOS or vR for radial velocity. Coherent Doppler

lidars are typically scanned via the Velocity Azimuth Display (VAD)
technique whereby the laser beam is swept in a circular pattern about
azimuth angle φ . When the beam of the lidar is swept in such a manner,
it intersects the wind at different angles. The resulting line-of-sight
velocity becomes a function of φ that behaves like a rectified sine wave
with maximum up-wind and down-wind speeds occurring at the peaks.
More details on the behavior of line-of-sight velocity measurements are
given in section 5.
 CHAPTER 5

The Qinetiq ZephIR CW Doppler Lidar

5.1. ZephIR Overview

The RERL’s lidar, manufactured by Qinetiq Ltd., is a fiber-optic based lidar
system. This particular instrument is marketed as the “ZephIR” lidar system.
It falls under the category of Doppler wind lidar systems. This instrument is
a 1.55µm CW coherent laser radar (CLR) that has a laser output power of 1-
Watt with a measurement range of 10m 150m (according to the manufacturer).
The ZephIR is a monostatic coaxial system where both the emitted and
backscattered light share common optics.

The system is specifically designed for autonomous wind resource

assessment purposes and includes a laser that emits at an eye-safe
wavelength. Since the eye safety requirement is satisfied, the assembly and
operation of this unit does not require the assistance of qualified laser
technicians. Figure 3 shows the ZephIR lidar instrument, which consists of
three “pods:” the battery pod (lowest pod), the electronics pod (middle pod)
and the optics pod (upper pod). The system is also equipped with a
meteorological mast that includes a thermometer and barometer as well as a
wind direction sensor.

14
 Figure 3: Qinetiq ZephIR Lidar System

The equation for the time-averaged optical signal power PS of a CW

CLR (such as the ZephIR) is shown below in Equation 2.
PS = πPT β (π )λ
Equation 2

In Equation 2, PT is the transmitted laser power, λ is the laser wavelength

at transmission and β (π ) is the atmospheric backscatter coefficient. It is
important to note that Equation 2 is independent of both focus range and
system aperture size [20].

Sections 5.2 - 5.6 outline the process by which the ZephIR obtains
a wind speed measurement after the system receiver has collected
the scattered light.

15
5.2. ZephIR Photodetection
When the lidar receiver collects scattered light, it is then optically
mixed with the reference, or local oscillator (LO), beam as shown in
Figure 4. While Figure 4 shows a generic lidar system in a bistatic
configuration (where the transmitter and detector are separated) the
same general principles apply to both monostatic and bistatic lidar
instruments.

Transmitted light
Laser
Laser Target
Local oscillator
(Reference beam) Target
Detector

Scattered &
Reflected light
(With Doppler frequency shift)
Figure 4: Generic Bistatic Lidar System

The detector creates an electric signal that is digitally sampled for the
purpose of determining the Doppler shifted frequency of the return light.
The conversion of incident photons to photoelectrons, which generate a
measurable current, is accomplished by a photodiode of the same type
that is commonly used in the telecommunications industry. The
photoelectrons can then be amplified and digitized for the subsequent
detection of the Doppler shifted return frequency. The output of the
photodetector is, however, comprised of many sources of noise. The main
noise components of the photodetector signal are:

16
 • Dark noise – The intrinsic wideband noise floor that is generated
by the detector and amplifier combination in the absence of any
incident light.
• Photon shot noise – Also known as quantum noise, this source of
noise is the random generation of photoelectrons by the incident LO
beam that leads to a wideband, and spectrally flat, noise source. The
shot noise power spectral density can be shown to increase in
proportion to the optical power of the LO beam [21].

Laser relative intensity noise (RIN) – The intensity fluctuations that are in
excess of photon shot noise. Such intensity fluctuations can be caused by (e.g.)
relaxation oscillations of the laser output where a small disturbance in the laser
power causes a damped oscillation of the laser output power before once again
returning to steady state [22]. Such oscillations typically occur at low frequency
levels and hence only affect the sensitivity of the lidar during low wind speed
events.
These sources of noise therefore require that the ZephIR photodetector have
a high level of quantum efficiency, sufficient bandwidth to measure
maximum Doppler frequencies of interest and a photon shot noise
contribution that sufficiently exceeds the dark noise intensity level. It is
desirable to have a dominant shot noise contribution because it is spectrally
flat and thus more predictable so that it can be treated as a “noise floor.” The
ZephIR’s InGaAs (indium gallium arsenide) photodiode is capable of
meeting these requirements for applications in wind resource monitoring
[15].

5.3. ZephIR Fourier Analysis

After the detector converts the backscattered light to an electric signal, the
signal is digitally sampled at a rate of 100 MHz by the data acquisition
system that is incorporated in the ZephIR. Next, the signal is sent through a
low-pass filter with a cut-off frequency of 50 MHz. A 512-point fast Fourier
transform (FFT) is then applied to the digitized signal to determine its
frequency content. The 512-point FFT yields 256 bins in the spectrum. In
order to increase the signal to noise ratio, 4,000 of these individual power
spectra are averaged to create each wind spectrum. After the averaging step,
a clear Doppler frequency peak appears in the wind spectrum as shown in
Figure 5.

17
 Frequency [Hz]

Figure 5: Doppler-Shifted Wind Spectrum [23]

5.4. ZephIR Cloud-Correction Algorithm

CW CLR systems such as the ZephIR do not create range -resolved
measurements by observing the “time-of- flight” of the emitted
radiation, as is the case with pulsed lidar systems. Instead, CLR systems
focus their beam at a specific height to obtain measurements. This
technique can lead to problems when the lidar’s beam intersects a cloud
base. When such an event occurs, the cloud’s frequency contribution to
the Doppler-shifted return signal can contaminate backscatter signal
from the aerosols at the desired height. If left unchecked, this
contamination can cause an overestimation of the true wind speed at the
height of interest. The severity of this phenomenon depends on a number
of factors. As such, the threat of measurement error increases for low
cloud height, high lidar altitude sensing, low aerosol scattering at the
desired height and high wind shear conditions.

In order to mitigate the risk of wind speed overestimation due to the presence
of clouds in the atmosphere, the frequency component that is associated with
the cloud base must be identified and isolated from th Doppler spectra. The
ZephIR employs an effective cloud-correction algorithm that is proven to
minimize measurement error [24].

18
The details of this algorithm are proprietary but the essence of its operation is
illustrated in Figure 6. The upper part of this figure shows the wind spectrum at
150 meters before the cloud-correction algorithm has been applied. A broad
aerosol return signal appears to the left of a narrow peak that is caused by the
presence of clouds at a slightly higher altitude. The middle plot shows the
corresponding spectrum that is obtained by focusing the lidar beam at 300
meters, where cloud density is assumed to be more intense. Notice that at 300
meters, the spectral peak from the cloud retains the same Doppler shift and its
peak is amplified. The lower plot of Figure 6 shows the resulting wind
spectrum after the cloud-induced frequency component (middle plot) has been
subtracted from the original spectrum (upper plot). The outcome of this process
is the elimination of the frequency component that was caused by the presence
of clouds.

Figure 6: ZephIR Cloud-Correction Algorithm

19
A detailed analysis performed by Albers in 2006 confirms that the
ZephIR cloud-correction algorithm is effective in dramatically
improving the quality of the wind velocity measurements at 65
meters and 124 meters [24].

5.5. ZephIR Wind Velocity Estimation

After each wind spectrum is produced, it is checked for corruption
resulting from the presence of clouds. Next, an algorithm that
determines the shifted frequency of the scattered light is applied to the
digitized signal that is generated by the photodetector. The dominant
frequency of the photodetector signal can be determined by a simple
procedure where the Doppler peak is chosen based on the location of
maximal power density in the spectrum. A better method is to employ
an algorithm that determines the first moment (centroid) of the spectra
around the largest peak [25]. The ZephIR system incorporates a similar
peak-picking algorithm, but the details of its operation are proprietary.
The frequency behavior of the scattered light is described by Equation
3 where vLOS is the line-of-sight velocity and λ0 is the wavelength of the
transmitted light.
fDoppler Shifted = 2vLOS
λ0
Equation 3
Equation 3 can be rearranged to show that the line-of-sight wind speed is
determined by multiplying the shifted Doppler frequency by a simple
conversion factor of 0.775 ms-1 per
MHz, or λ20 . A study performed by Frehlich contends that, for pulsed lidar
operation,
this calibration factor suffers negligible drift over long periods of time
[25]. Furthermore, Jorgensen et al. contend that the ZephIR (a CW
Doppler lidar) is capable of stable laser frequency transmission at
1.55µm with less than 0.2% drift over long periods of time [26]. Thus,
the ZephIR is an absolute instrument that does not require calibration

20
The ZephIR emits laser radiation in a circular pattern by reflecting the
laser beam off of a spinning optical wedge, via the VAD scanning
technique. The wedge is positioned such that the beam is transmitted at
an angle of 30 degrees from zenith, thereby creating an upside-down cone
shaped probe volume. The line-of-sight velocity data measurements
therefore become a function of scan angle, shown in Equation 4
vLOS = a cos(φ -b) + c ,
Equation 4
where angle φ is the azimuth scan angle. The parameters a, b and c in
Equation 4 are obtained by applying a non-linear least squares fit to
the line-of-sight data that are collected by the lidar. The wind speed
can then be determined by substitution in the following equations:

a
u=
sin(θ)
c
w=
cos(θ)
Bearing ± 180o = b

Equation 5

where u is the horizontal wind speed, w is the vertical wind speed and b is
the direction of approaching wind. The parameter b is directly obtained in
the curve-fitting operation. If the line-of-sight curve-fit is poor, then it is
possible that a 180° wind bearing ambiguity can occur. This potential
ambiguity is resolved by verification with the lidar met mast wind
direction sensor. When the lidar mast is unobstructed, wind direction
measurement errors are rare [24].

21
 In the curve-fitting algorithm, the wind speed data below a certain
threshold (approximately 0.5- 1m/s) are eliminated from consideration.
Low wind speeds are eliminated from the curve- fitting algorithm because
they are typically more variable and thus more likely to disrupt the
accuracy of the overall curve fit. If there are enough data for a valid fit,
then the curve-fitting algorithm proceeds. However, if there is excess
noise in the return signal, then a fit may not be possible until the noise
threshold is incremented and another fitting iteration is attempted. The
next step in the curve-fitting algorithm searches for large deviations from
the sine fit. Velocity data with large, non-Gaussian deviations are
separated and eliminated from consideration. After these steps, a
nonlinear least squares fit is performed on the filtered data. The result of
this process is illustrated in Figure 7 where the solid line represents the
best fit to the line-of-sight wind speed data.

VLOS [m/s]

Figure 7: Lidar Line-of-Sight Velocity as a Function of Scan Azimuth [27]

When the line-of-sight wind speed vs. azimuth angle is plotted on a polar
axis the result is shown in Figure 8 for a three-second measurement period.
When the atmospheric backscatter coefficient is large, the lidar calculates
up to a maximum of 150 line-of-sight data points for each three-second
measurement period.

22
However, when clear conditions are present, fewer data points are available
for the curve fitting process. In the event of extremely clear conditions, a
valid wind speed measurement may not be possible. The actual number of
data points in the curve fitting process is supplied in the ZephIR output data
file. More detail on ZephIR operation in clear conditions is given in section
10.3.6.
The relationship of the data in Figure 8 to the best-fit approximation (solid
line) suggests that the wind flow across the probe volume is uniform and
the slight asymmetry in the lobe sizes indicates that the presence of a
vertical wind speed component. Here, the atmospheric backscatter
coefficient is large because 147 data points were available for the curve
fitting process. The wind direction shown in Figure 8 is approaching from
the NNE direction.

Figure 8: ZephIR Polar Line of Sight Wind Speed Plot in m/s vs. Azimuth Angle

5.6. ZephIR and Range-Resolved Measurements

The ZephIR scans the wind at up to 5 user-programmable heights at a rate of
1 revolution per second for a period of three seconds at each height. Wind
speed and direction data can be measured at various heights by focusing the
laser beam at a preset range. In order to understand the way in which the
lidar measures at various heights, the hourglass-shaped beam geometry must
be considered.

23
 Conclusions
This report presented a complete analysis of the Qinetiq ZephIR lidar
system. In the lidar data validation experiment, the lidar accurately
measured the wind velocity at a test site in Hull, Massachusetts for
approximately 1.5 months with minor difficulties. The lidar data were
compared to simultaneous cup-anemometer measurements that were
collected at a nearby radio tower with sensors at 61 m, 87 m and 118 m.
The cup-anemometer data were found to exhibit temporary periods of
measurement irregularities, which were later removed by a rigorous data
filtering process. After filtering, the data from the two measurement
sources demonstrated strong correlations at all heights. Over the course
of the 1.5 month measurement campaign, the lidar was shown to
introduce an overall measurement bias of approximately +1% at heights
above 80 meters.
A head-to-head uncertainty analysis was also presented. The purpose of
this study was to present more information regarding the benefits and
limitations associated with cup anemometer and lidar-based wind speed
measurement. In this analysis, the various sources of measurement
uncertainty for each sensor were investigated. A generic CW lidar
instrument was shown to reduce the amount of measurement uncertainty
involved in wind resource assessment with traditional cup anemometers.
The lidar’s ability to reduce the error associated with tower-mounted cup
anemometer measurement (e.g. tower shadow effects, wind shear
extrapolation, etc.) makes lidar technology attractive for wind power
applications.

While the lidar is a promising new advance for the wind energy
industry, this report identifies the need for improvement in the
following areas:

• More robust battery backup system

• More restrictive data validity requirements during cold
temperature startup/operation

24
 • Expanded averaging options when data are directly retrieved
from the compact flash memory card

The advantages provided by the lidar are summarized as follows:

• Portability and ease of rapid deployment

• Small footprint that does not require permits to place at a given site
• Silent operation that is not subject to echo interactions
• Ability to provide accurate hub-height wind data
• Capacity to minimize measurement uncertainty
With these advantages, the resource assessment portion of any wind
energy project development process could be improved and streamlined
by the substitution of lidar-based measurement in place of cup
anemometers. The most important benefit associated with lidar wind
resource assessment is that the lidar is capable of measuring wind speeds
at the hub height of a modern wind turbine. When hub height data are
available for analysis, access to financial capital for wind project
development becomes less restrictive because wind shear extrapolation is
one of the most detrimental sources of measurement uncertainty in a
given wind resource assessment campaign .

Overall, the Qinetiq ZephIR lidar achieves a favorable review based on

the findings of this report. Thus, this system is approved to be added to
the RERL suite of wind speed measurement devices .

25
 REFERENCES

1. Antoniou, I., H.E. Jorgensen, T. Mikkelsen, S. Frandsen, R. Barthelmie, C.

Perstrup, and M. Hurtig, Offshore Wind Profile Measurements from Remote
Sensing Instruments, in European Wind Energy Conference. 2006: Athens,
Greece.
2. Hensen, W.L.W., A.L. Rogers, and J.F. Manwell, Wind Resource Assessment
Using SODAR at Cluttered Sites, in European Wind Energy Conference.
2006: Athens, Greece.
3. Lackner, M.A., A.L. Rogers, and J.F. Manwell, Wind Energy Site Assessment and
Uncertainty 2006, Renewable Energy Research Laboratory.
4. Synge, E.H., Philosophical Magazine, 1930. 52: p. 1014-1020.
5. Duclaux, J. Phys. Radiat., 1936. 7: p. 361.
6. Hulbert, E.O., J. Opt. Soc. Am., 1937. 27: p. 377-382.
7. Bureau, R., Meteorologie, 1946. 3: p. 292.
8. Maiman, T.H., Nature, 1960. 187: p. 493.
9. McClung, F.J. and R.W. Hellworth, Journal of Applied Physics, 1962. 33: p. 828-
829.
10. Smullins, L.D. and G. Fiocco, Nature, 1962. 194: p. 1267.
11. Fiocco, G. and L.D. Smullins, Nature, 1963. 199: p. 1275-1276.
12. Argall, P.S. and R.J. Sica, Lidar, in Encyclopedia of Imaging Science and
Technology. 2002, Wiley.
13. Safety of Laser Products - Part 1: Equipment classification, Requirements, and
User's Guide., IEC, Editor. 2001, IEC60825-1.
14. Harris, M., M. Hand, and A. Wright, Lidar for Turbine Control. 2005, National
Renewable Energy Laboratory.
15. Karlsson, C.J., F.A.A. Olsson, D. Letalick, and M. Harris, All-fiber Multifunction
Continuous-wave Coherent Laser Radar at 1.5 µm for Range, Speed, Vibration,
and Wind Measurements. Applied Optics, 2000. 39(21): p. 3716-3726.
16. Weitkamp, C., Lidar Range-Resolved Optical Remote Sensing of the
Atmosphere. 2005.
17. Mie, G., Articles on the optical characteristics of turbid tubes, especially
colloidal metal solutions. Annalen der Physik, 1908. 25: p. 377.
18. Guibert, S., V. Matthias, M. Schulz, J. Bosenberg, R. Eixmann, I. Mattis, G.
Pappalardo, M.R. Perrone, N. Spinelli, and G. Vaughan, The Vertical Distribution
 of Aerosol Over Europe - Synthesis of One Year of EARLINET Aerosol Lidar
Measurements and Aerosol Transport Modeling with LMDzT-INCA. Atmospheric
Environment, 2005. 39(16): p. 2933-2943.
19. Pearson, G.N., P.J. Roberts, J.R. Eacock, and M. Harris, Analysis of the
Performance of a Coherent Pulsed Fiber Lidar for Aerosol Backscatter
Applications. Applied Optics, 2002. 41: p. 6442-6450.
20. Sonnenschein, C.M. and F.A. Horrigan, Signal-to-noise relationships for coaxial
systems that heterodyne backscatter from the atmosphere. Applied Optics, 1971.
10: p. 1600-1604.
21. Bleaney, B.I. and B. Bleaney, Electricity and Magnetism. 1976:
Oxford University Press.

22. Siegman, A.E., Lasers. 1986, Mill Valley, California: University Science Books.
23. Smith, D.A. and e. al. Wind Lidar Evaluation at the Danish Wind Test Site
in Hovsore. in European Wind Energy Conference. 2004. London, UK.
24. Albers, A., Evaluation of ZephIR. 2006, Deutsche WindGuard Consulting GmbH.
25. Frehlich, R.G. and M.J. Yadlowsky, Performance of Mean-Frequency Estimators
for Doppler Radar and Lidar. Journal of Atmospheric and Oceanic Technology,
1994. 11: p. 1217-1230.
26. Jorgensen, H., T. Mikkelsen, J. Mann, D. Bryce, A. Coffey, M. Harris, and D.A.
Smith. Site Wind Field Determination Using a CW Doppler Lidar- Comparison
with Cup Anemometers at Riso. in The Science of Making Torque from Wind.
2004. Delft.
27. Danielian, R., H.E. Jorgensen, T. Mikkelsen, J. Mann, and M. Harris, Surface-
Layer Wind and Turbulence Profiling from LIDAR: Theory and Measurements,
in European Wind Energy Conference. 2006: Athens, Greece.
28. Banakh, V.A., I.N. Smalikho, F. Kopp, and C. Werner, Representativeness of
Wind Measurements with a CW Doppler Lidar in the Atmospheric Boundary
Layer. Applied Optics, 1995. 34(12): p. 2055-2067.
29. Fowles, G.R., Introduction to Modern Optics, Second Edition. 1989:
Dover Publications.
30. Harris, M., G. Constant, and C. Ward, Continuous-wave Bistatic Laser Doppler
Wind Sensor. Applied Optics, 2001. 40(9): p. 1501-1506.
31. NRG. NRG Maximum 40 Specifications. [cited; Available from:
http://www.nrgsystems.com/store/product_detail.php?cd=11&s=1899.