EPA-650/2-73-037
November 1973
Environmental Protection Technology Series
,•.•.•.*
Wx'
ll^
-------
EPA-650/2-73-037
FEASIBILITY
OF A CW LIDAR TECHNIQUE
FOR MEASUREMENT
OF PLUME OPACITY
by
Richard A. Ferguson
Stanford Research Institute
Menlo Park, California 94025
Contract No. 68-02-0543
Program Element No. 1A1010
EPA Project Officer: William D. Conner
Chemistry and Physics Laboratory
National Environmental Research Center
Research Triangle Park, N. C. 27711
Prepared for
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
WASHINGTON, D . C. 20460
November 1973
-------
This report has been reviewed by the Environmental Protection Agency and
approved for publication. Approval does not signify that the contents
necessarily reflect the views and policies of the Agency, nor does
mention of trade names or commercial products constitute endorsement
or recommendation for use.
11
-------
SA in
FRONTISPIECE
MODULATED ARGON LASER BEAM TRANSMITTED FROM REMC
OUT OF VIEW AT LEFT BEING USED TO DETERMINE THE TARGET
CHARACTERISTICS OF THE THIN STEAM PLUME RISING FROM BUILDING
AT RIGHT
iii
-------
ABSTRACT
This report describes the work performed for the Environmental
Protection Agency during the initial proof-of-concept phase of a program
to develop an eyesafe CW lidar for remote measurement of the opacity of
smoke plumes from industrial smoke stacks. The analysis, design, con-
struction, and evaluation of n laboratory model CT lidar were performed
under SRI Project 1979 from 30 May 1972 to 30 May 1973 under KPA Contract
68-02-0543 to determine the limitations and potential of the technique.
The proof-of-principlc experiments combine what is called an FM-CW radar
technique with an argon laser. The technique involves modulating the
intensity of the laser beam at a frequency that changes rapidly and
linearly with time. A portion of the transmitted signal is mixed elec-
tronically with the light reflected from the targets in a device similar
to a radio receiver. Each target appears at a particular frequency.
By tuning the radar's receiver to these target frequencies, the re-
searchers were able to measure both the range and the opacity of semi-
transparent targets over distances of 100 to 200 meters.
-------
CONTENTS
FRONTISPIECE iii
ABSTRACT v
LIST OF ILLUSTRATIONS ix
ACKNOWLEDGMENTS xi
I INTRODUCTION 1
II SUMMARY AND RECOMMENDATIONS 3
III BACKGROUND 7
IV DESIGN OF THE CW LIUAR 9
A. An Idealized FM-CW Radar 10
B. An Idealized CW Lidar for Opacity Measurement ... 16
C. A Practical CW Laser Radar 22
1. Laser-Beam Modulation 23
a. Modulation Waveform 23
b. Modulation Techniques 39
c. Modulation Equipment 44
2. Laser Transmitter 47
3. Optical Receiver 51
4. Atmospheric Effects 57
5. Signal Processing 61
V EXPERIMENTAL PROGRAM 65
A. Design and Construction of a Laboratory-Model CW
Lidar 65
B. Target Display Characteristics 69
1. Discrete Targets 71
2. Steam-Plume Characteristics 74
3. Clear-Air Targets 77
vii
-------
C. Remote Opacity Measurements 81
D. Discussion of Experimental Results 86
REFERENCES 89
BIBLIOGRAPHIC DATA SHEET
viii
-------
ILLUSTRATIONS
1 Basic FM-CW Radar Performance 11
2 Multiple-Target Performance of an Ideal FM-CW Radar. . . 13
3 Effect of a Moving Target on FM-CW Radar Performance . . 14
4 Distributed Targets and Changes in Their Attenuation
Characteristic 15
5 Doppler Spreading of Distributed-Target Spectrum .... 17
6 Typical Range Profile on Spectrum Display 19
7 Signal-Processing Block Diagram 20
8 Result of Finite Limitations on Linear Sweep 24
9 Mixer Output Spectrum 28
10 Sideband Envelopes of Target Between Node Ranges .... 33
11 Inaccuracies in Measuring Relative Power Due to
Spurious Sidebands 36
12 FM Sweep Using Triangle Waveform 37
13 Improvement in Target Range Resolution by Increasing
Deviation, A.f 40
14 Optical FM-CW Radar Modulation 41
15 FM-AM-CW Radar Modulation . -43
16 Electrooptical Modulator Performance 45
17 Spatial Integration of Beam Intensity Function by
Receiver 54
18 Laboratory Model FM-AM-CW Lidar Equipment Block
Diagram 66
19 Measured Modulator Transmission Characteristic 67
20 Data Displays Used for Remote Target Measurements. ... 70
21 Lidar Target Range Used for Remote Range and
Opacity Measurements 72
ix
-------
22 Discrete-Target Spectra Showing Spurious Sideband
Changes with Range 73
23 Rooftop Lidar Site 75
24 Distributed Steam-Plume Target Spectrum 76
25 C\V Lidar Clear-Air Returns 78
26 Effects of Spatial Integration on Signal-to-Noise
Ratio of Clear-Air Return. , 80
27 Changes in Relative Target.Power Levels Caused by
Various Target Opacities 82
28 Results of Remote Opacity Measurements Using a CW
Lidar 84
-------
ACKNOWLEDGMENTS
The author would like to thank his associates, D. W. Jackson, D. G.
Saraf, J.W.O. Knotts, P. Chancellor, D. E. Arnold, R. Gumming, and L. L.
Guild for their valuable assistance in the conception, design, fabrica-
tion, testing, and documentation of the instrument described herein. In
addition, the support and guidance of Mr. William Conner of the Chemistry
and Physics Laboratory at NERC/EPA is greatly appreciated.
xi
-------
I INTRODUCTION
This report presents the results of an experimental study to demon-
strate the feasibility of a continuous-wave laser radar (CW lidar) for
the remote measurement of smoke-plume opacity. The C\V radar technique
, *
employed was proposed by D. W. Jackson" in January 1970 as an alterna-
tive to the pulsed laser technique reported by W. E. Evans3 in 1967. The
work discussed in this report included the analysis, design, fabrication,
and limited field evaluation of a laboratory-model CW lidar. The tech-
nical work was conducted under Contract 68-02-0543 with the Environmental
Protection Agency (EPA) for the twelve months ending 30 May 1973.
The main thrust of this effort was to demonstrate the CW lidar con-
cept with a working hardware model in order to facilitate development of
a low-cost, field-tested, portable smoke-plume opacimeter within three
years. A distinct lack of published experimental data on C\V laser radar
performance in the atmosphere had frustrated earlier efforts to analyti-
cally determine the proposed system performance. The experimental work
reported here provides the results of investigations into actual system
performance, operating characteristics, and remote transmission measure-
ment capability using components available at reasonable cost within the
state of the art. These new data demonstrate the feasibility of the CW
lidar technique and indicate promising directions for continued development.
*
References are listed at the end of the report.
-------
II SUMMARY AND RECOMMENDATIONS
A theoretical model was developed for a CW (as opposed to pulsed)
laser radar capable of measuring laser radiation scattered back from
atmospheric particles at ranges between 150 ?^d 500 meters. The model
was then examined in light of practical constraints on implementing the
CW lidar technique using a low-cost gas laser and inexpensive high-
frequency (HF) electronics.
After a survey of this analytical groundwork, actual lidar components
wore specified and purchased, and a laboratory-model CW lidar was con-
structed. Transmitting a continuous green beam of light at 5145 ang-
stroms, the lidar was used in a series of proof-of-principle experiments
to demonstrate the CW lidar concept,
The technique involves modulating the intensity of the laser beam at
a frequency that changes rapidly and linearly with time. A portion of
the transmitted signal is mixed electronically with the light reflected
from the targets in a device similar to a radio receiver. Each target
appears at a particular frequency. By tuning the radar's receiver to
these target frequencies, the operator can measure both the range and
the opacity of semi transparent targets over distances of 100 to 200 meters,
The blue-green argon laser beam is safer than the powerful, pulsed
beams used in earlier lidar monitoring devices that send out megawatts
of peak power in pulses faster than the eye's protective blink reaction.
The CW laser radar transmits less than 400 milliwatts of power continu-
ously, and is no more harmful to the eye than headlights or reflected
sunlight.
-------
Since the argon laser beam is sale and continuously visible, it is
easier to align and calibrate the CW system than a pulsed laser radar
system. Unlike the data output from a pulsed system, the CW radar
measurements can be read directly and continuously at the receiver's
output, promising easier operation and more reliable data.
A laser radar offers an accurate method for measuring targets that
are more than 80 percent transparent. As smokestack emission control
laws grow more stringent, enforcement agencies will need accurate measure-
ments of thinner, semitransparent smoke plumes that may be in violation
of the law.
In the effort to facilitate development of a field-portable, low-
cost lidar opacimeter within three years, the laboratory model was used
successfully in actual outdoor experiments to:
(1) Detect and find range to discrete targets.
(2) Detect and find range to distributed scattering targets
such as a steam plume and clear-air volumes.
(.'•i) Make remote measurements of target opacities between 10
percent and 50 percent (using a plate glass reference).
It is concluded that there are no fundamental obstacles to employing
the CW lidar technique in the remote measurement of smoke-plume opacity
using clear-air volumes as detectable reference targets. The engineering
requirements for designing an engineering model of the system have been
identified.
It is recommended that in order to provide a useful field-portable
instrument for remote stack monitoring, research, and enforcement pur-
poses in the near future, a refined research model should be designed,
constructed, and field-tested. This model should employ an improved,
low-noise receiver detector and dedicated signal-processing electronics
based on the three-AM-radio-receiver concept described herein. These
-------
improvements will facilitate night and even daytime use of laser light
backscattered from clear-air volumes as references against which plume
opacity may be measured. The design of the research model should include
investigations of innovative concepts for providing low-cost system-
component functions. The completed model should be transportable for
remote field measurements of actual plumes and could be used as a mobile
research tool to complement existing pulsed lidars. The component and
field-performance data provided by such a model would also serve as a
basis for constructing a preproduction prototype of a simple^ low-cost
remote monitoring device for routine enforcement purposes.
-------
111 BACKGROUND
Remote measurement of the density of a smoke plume is an important
capability in the field of air pollution research and control.2~s The
single-beam lidar transmission technique has been judged a useful and
appropriate means for remote measurement of stack-effluent opacity. The
technique's comparison of optical-wavelength power backscattered from
"clear-air" regions surrounding a plume is similar to other existing
methods using the transmission of visible light to infer particulate
concentrations.2' 7~9
However, lidar developments to date have employed pulsed lasers in
systems similar to that described by Evans.2 These pulsed lidars have
certain drawbacks resulting from the transient and high-peak-power charac-
teristics generic to pulsed radars. A relatively expensive pulsed laser
with high transmitted power is needed to obtain the necessary signal
levels from the backscatter phenomenon. 'The powerful pulsed beam also
presents a potential eye hazard, and precautions must be taken in the
design of the lidar to prevent operators or bystanders from looking into
the beam. Anyone looking into the beam risks receiving a damaging dose
of light energy to the retina since the laser pulse is much faster than
the eye's protective blink reaction.
An additional limitation on the pulsed lidar technique is that, be-
cause of the transient nature of the low-repetition-rate pulsed lidar
echo, data must be recorded on videotape, videodisk, or Polaroid film.
Each of these media makes for a somewhat tedious and tardy extraction
of opacity values or Ringelmann ratings from the data. While it is now
possible to record the transient echo digitally and to analyze the data
-------
nearly in real time with a minicomputer, the expense of fielding a prac-
tical network of these digital pulsed systems is probably prohibitive.
Drawing on the pulsed-lidar experience of associates at SRI, D. W.
Jackson proposed eliminating the transient, high-peak-power drawbacks of
pulsed lidars by employing a CW radar technique.1 Since radar theory
shows that the information obtained from a target is a function of average
radiated power, a one-watt CW laser could in theory make measurements
equivalent to those made by a 50-megawatt pulsed laser transmitting one
20-nanosecond pulse per second. A CW lidar holds promise of greater eye
safety, lower cost, and easier operation. A dearth of applicable pub-
lished data on CW lidar experiments led to the conclusion that a laboratory
model should be built to determine the requirements and potential for de-
signing a field-portable CW lidar for remote measurement of smoke-plume
opacity. The EPA made funds available for this initial work beginning
in May 1972.
-------
IV DESIGN OF THE CW LIDAR
Designing a C\V laser radar requires an understanding of CW-radar
principles, a practical technique for applying those principles to a
laser's optical carrier frequency, and knowledge of the effects on the
system due to actual target characteristics and practical limitations
on component performance. Section IV-A below satisfies the first re-
quirement by briefly analyzing the idealized frequency-modulated
continuous-wave (FM-CW) radar models useful in understanding the lidar
design presented later. The adequacy of these models has been amply
demonstrated in practice with existing VHF and UHF C\V radars. Next, an
idealized lidar for measurement of smoke-plume opacity is described in
Section IV-B; with this lidar, HF CW radar signals are used to modulate
the intensity of the optical carrier. This modulation technique neatly
sidesteps the significant problems of optical-frequency modulation,
coherent mixing, and holographic signal processing, while making atmo-
spheric aerosol particles visible to a radar employing only existing,
low-cost HF electronics technology. Finally, without resorting to
rigorous and cumbersome mathematics, the effects of particulate targets
and practical component limitations are illustrated with diagrams, simple
formulas, and examples drawn from the working model built during this
study. The elements of Section IV-C form a basis for designing a low-
cost, portable, eye-safe lidar for the remote measurement of smoke-plume
opacity. The actual design and performance of a "proof-of-concept" C\V
lidar opacimeter are discussed in later sections.
-------
A. An Idealized FM-CW Radar
The FM-CW radar employs a frequency-modulated continuous carrier
*
signal to provide range measurement. If the transmitted frequency f
o
is a linear function of time, the range (transit time) to a target is
proportional to the difference (or beat) frequency, f , between the echo
signal and the transmitter signal. A block diagram of the simplest FM-CW
radar is shown in Figure l(a), while l(b) and l(c) indicate the idealized
relationships between the transmitted frequency, the echo frequency, and
(lie difference frequency as functions of time.
Figure l(d) shows the frequency spectrum F of the output of the
R
mixer used to obtain the difference frequency f representing the range
Rl to a discrete stationary target placed in the radar's field of view.
Just as a single target appears as a "blip" on the amplitude-vs.- time
display of a pulsed radar, the FM-CW radar target appears as a blip on
a spectrum analyzer's amplitude-vs.-frequency display of F . The target
R
appears as a frequency-component spike at f . The reader should note
Rl
here that if the FM sweep function in Figure l(b) is not perfectly linear,
the difference function f does not have a constant value during the
sweep. Should the target cross section increase, the delayed FM echo
signal received at the mixer input increases, and so does the amplitude
of the output of this idealized mixer. The open-bar component in Figure
l(d) represents this increase in the mixer output signal appearing in
the spectrum display of F .
R
*
In principle, the carrier can be either the optical wave itself or a
radio-frequency sinusoid that modulates the intensity of the optical
wave. For reasons that will be developed later, the latter method was
used in our experimental work. However, the present discussion is
general and applies to either method.
The capital letter FR is used to refer to the entire beat-frequency
spectrum of any signals appearing at the mixer output. The lower-case
f refers to one component in that spectrum.
10
-------
FREQUENCY
MODULATION
-—
CONTINUOUS
CARRIER
SOURCE
1. . MIXER
1
*-
BEAT
FREQUENCY
SPECTRUM
ANALYZER
(•) BASIC TM-CW RAOAR
z
LU
b
LU'
LT
u_
Q
LU
I-
10
Z
TRANSMITTED
SIGNAL
DELAYED
ECHO SIGNAL
TIME
(b)
*
>
o
TIME
o
a.
o
a.
en
INCREASE IN SPECTRAL POWER
DUE TO INCREASE IN POWER
:LECTED FROM T
SWEEP RATE
R - TARGET RANGE
2RS
BEAT FREQUENCY, fR
(d) TARGET SPECTRUM DISPLAY
SA-1979-1
FIGURE 1 BASIC FM-CW RADAR PERFORMANCE
11
-------
Where two or more targets of varying cross section lie at different
ranges in the radar's field of view, superposition theory demands that
their echo signal amplitudes add linearly at the mixer input. As a re-
sult, the mixer output spectrum will contain two or more Fourier com-
ponents representing the ranges to the various targets. Figure 2(a)
shows the relationship of the multiple delayed echo signals to the trans-
mit signal, while 2(b) presents the resulting mixer output spectrum.
Only because the KM sweep function is linear can each of the multiple
stationary targets be distinguished as a separate component during the
sweep.10
Moving targets present still another situation for analysis with
the FM-CW radar model. The well known Doppler effect shifts by f the
frequency of a signal reflected by a target moving at some velocity with
a component v along the radar-beam axis. Figure 3(a) shows how an
approaching target shifts the echo-signal frequency upward (reducing f ),
while a receding target causes a downward shift in the echo frequency
(increasing i ). The Doppler effects on a spectrum display of f are
presented in Figure 3(b). Random velocities in a multiple-target situa-
tion can cause incorrect range interpretation, as can be seen by comparing
the Doppler version of a moving target at R2 [Figure 3(b)] with its
stationary range equivalent [Figure 2(b)].
Before leaving the discussion of a simple, idealized FM-CW radar,
one more important target situation needs to be considered—multiple
scatterers distributed in a volume. Consider first a line of closely
spaced stationary targets each of which scatters back some of the transmit-
beam energy, causing a characteristic attenuation of the signal with in-
creasing range. Figure 4(a) shows this multiple-target situation on an
idealized FM-CW radar-spectrum display of F . As the target spacings
R
AR, grow smaller, the corresponding spectral line spacings, Af ., grow
12
-------
o
z
HI
h-
I
<
oc
I-
a.
UJ
O
o.
cc
o
Ul
CL
TRANSMITTED
SIGNAL
— *RlH
'H2
R3
i-'i-i-'-* > I-L^ tV^riw OIUIUML.O
FROM THREE TARGETS
i ^B-
TIME
(a) MULTIPLE-TARGET ECHO RELATIONSHIPS TO fo
f=
R2
BEAT FREQUENCY, fR
I I
'R3
I
R1
R2
R3
RANGE, R
(b) SPECTRUM OF MULTIPLE-TARGET BEAT FREQUENCIES
SA-1979-2
FIGURE 2 MULTIPLE-TARGET PERFORMANCE OF AN IDEAL FM-CW RADAR
13
-------
o
z
UJ
O
o
in
<
OC
h-
TRANSMITTED
SIGNAL
DOPPLER-SHIFTED
ECHO FROM
APPROACHING
TARGET AT R2
-,HO FROM
STATIONARY
TARGET AT R2
DOPPLER-SHIFTED
ECHO FROM
RECEDING
TARGET AT R2
TIME
(a) DOPPLER-SHIFTED ECHOES OF TRANSMITTED SIGNAL
a
DC
LLJ
5
o
a.
_j
<
OC
fc
UJ
CL
c/>
FIXED
TARGET AT R2
rf ^
TD — f~
APPROACHING
TARGET AT R2
•
— «_ ^l
"^ 'D ^n
1
1
1
1
I
fR2
BEAT FREQUENCY, fR
I I 1
0 R1 R2 R3
RECEDING
TARGET AT R2
RANGE, R
(b) SPECTRUM SHOWING DOPPLER SHIFT FROM CORRECT fR
SA-1979-3
FIGURE 3 EFFECT OF A MOVING TARGET ON FM-CW RADAR PERFORMANCE
14
-------
a.
K
LU
I
(a)
(b)
O
0.
O
a.
-------
Where the distributed target consists of moving elements, Doppler
spreading of spectral lines becomes an important consideration. If the
various random target velocities are small, with some approaching while
others recede, the lines will not shift far and their attenuation charac-
teristic will still be apparent as shown in Figure 5(a). However, if the
velocities are high enough, the random shifts f up and down can spread
the f component lines far enough to seriously distort the targets'
attenuation characteristic as a function of range [Figure 5(b)J. Since
tho determination of the attenuation characteristic of multiple scatterers
distributed in a volume is the key to adapting an FM-CW laser radar for
opacity measurements, the magnitude of Doppler spreading of spectral
lines is an important consideration.
B. An Idealized CW Lidar for Opacity Measurement
The FM-CW radar technique idealized above can be extended for use
with a CW laser to model an FM-CW lidar capable of making remote trans-
mission measurements on smoke plumes. The use of HF-CW radar modulation
signals to modulate the intensity of the optical carrier, and a photo-
diode device to detect the HF envelope on the received echo, permits
existing HF radar electronic signal-processing designs to measure optical
power reflected from various ranges. The approach lends itself well to
the opacity-measurement problem. Since the optical carrier wavelength
is on the order of the particle size or shorter, the amount of radar
power backscattered from either "clear-air" particles or plume particles
can be great enough for detection. Measurement of the backscatter from
particles in volumes at particular ranges can be made simply by examining
the spectral power (or amplitude) of the HF mixer output components at
the appropriate range frequencies f . Transmission can then be deter-
R
mined by comparing laser power backscattered by volumes in front of and
behind a semi transparent target.
16
-------
IT
UJ
I
o
LU
IT
III
I
o.
_i
<
1C
b
LU
BEAT FREQUENCY. fR
la) RANDOMLY DISTRIBUTED MULTIPLE TARGETS
IfTTrr-
BEAT FREQUENCY, fR
(b) SAME COMPONENTS AS IN (a) SHOWN RANDOMLY DOPPLER-SHIFTED.
Note loss of attenuation characteristic Information.
SA-1979-5
FIGURE 5 DOPPLER SPREADING OF DISTRIBUTED-TARGET SPECTRUM
17
-------
If the particle velocities v. are low enough to make their Doppler
shifts f « Af ., the distributed target modeled in Figures 4(a)
through 4(d) can be used to explain the CW-lidar-opacity-measurement
technique. Figure 6 shows a typical measurement situation with a plume
at R = 300 m. Notice that the attenuation characteristic of the clear-
air particles suffers a discontinuity and change in value at the range
of the plume and beyond. When the plume particles lie in both the trans-
mitter and receiver lidar-beam paths, the amount of transmitted power
reaching the far clear-air region is diminished twice by the transmission
01 the plume along the round-trip echo path. After the received light
power backscattered from the far region. P . and the power received
far
from a near region in front of the plume. P , are both measured with
near
and then without plume attenuation, the change in the atmosphere attenua-
tion characteristic can be determined. The advantage of the transmission
technique is that the logarithms of the relative power P = P /P
rel far near
for the plume and no-plume cases can give a direct and uniform measure
of plume opacity or Ringelmann number that is independent of sky color
or backlighting conditions.
Given a CW lidar with the signal-processing functions shown in
Figure 7, an operator would measure smoke-plume opacity according to the
following scenario. Training the laser beam on the plume, the operator
finds the range to the plume by tuning a simple AM radio receiver through
the K spectrum until he finds the maximum signal at f , . Two other
R 6 plume
identical receivers ganged with this receiver then tune automatically to
frequencies f and f just off the skirts on either side of the
near far
plume frequency. The rms voltage outputs of these two receivers corre-
spond to the light power scattered back by "clear-air" particles at
ranges R and R . The operator does not have to measure both of
near far
the scattered power levels. Instead, envelope-detecting log amplifiers
take the logarithms of the P and P receiver outputs, which are
near far
18
-------
1 —
I
f
10
10~
-1
10
,-3 _
plume
CLEAR AIR
STRONG RETURN FROM
SMOKE PLUME
0
III
CC
Ul
t 10-
<
i
z
$ io-5
a
in
§ 10-6
ui
DC
?
10
io-8
~^te
^^
VOLUME
1 1
•W^ 1 CLEAR AIR
•K "far
si ^
T-I
1 1
I 1
1
1
1
(1
1
1 1
1
r-4- ,—
r
^^ ^^^
-• FAR BACKSCATTER
VOLUME
1
100
200
300
400
500
R — meters
rr -iz
f»r
TA-655582-14
FIGURE 6 TYPICAL RANGE PROFILE ON SPECTRUM DISPLAY
-------
CO
o
SPECTRUM OF »^ from miner
BEAT FREQUENCIES S
+f-
DETECTOR
•w-
DETECTOR
ATMOSPHERIC
TRANSMITTANCE
INPUT CALIBRATION
-W-
DETECTOR
TA-655582-17R
FIGURE 7 SIGNAL-PROCESSING BLOCK DIAGRAM
-------
then subtracted in a difference amplifier. It is the low-bandwidth out-
put voltage of this difference amplifier that can be read at a meter or
chart recorder by the operator. But first the operator needs to cali-
brate the system to compensate for such biasing factors as atmospheric
attenuation due to clear-air particles in the plume volume.
The calibration is relatively simple. Once his scattering-volume
ranges are selected around the plume, the operator directs the lidar beam
through clear air alone slightly to one side of the plume. The output
voltage that he reads is simply:
= log P - log P (l)
near far
L— —' nr» Til limfs
V
'no plume
= |-log P^
'no plume
By adjusting a dc bias control on the differential amplifier, the operator
can return the amplifier output to zero. This effectively sets a constant
dc baseline shift in the signal-processing device, biasing the output by
the amount [+log P ] no plume.
rel
Now the operator points his lidar back through the plume and ob-
serves the output of the amplifier again:
v = log P - log P + dc bias calibration factor
o near far! .
•- -"plume
or
v = -logfp ) + log/P )
o \ rel/ °\ rel/
plume no plume
v = log/P ) - log/P ) (2)
o e\ rel/ \ rel/
no plume plume
21
-------
Because light signal power scattered from the far volume behind the plume
is reduced twice by the transmission T of the plume,
/P }
\ far/
} = T2/P (3)
far/ , p\ far/
plume no plume
Thus,
v = - - log T (4)
o 2 p
This output is directly proportional to the plume's optical density,
where optical density is defined for transmission T between zero and 1.0
as (-log T). In short, the idealized CW lidar allows an operator to
continuously read the opacity of the plume after a simple calibration of
the clear-air return.
C. A Practical CW Laser Radar
The laser radar designer must hurdle a number of practical obstacles
to building a working model of the basic CW laser radar discussed in
Section 1V-A, above. Perhaps the easiest obstacles to define and sur-
mount are those posed by equipment performance limitations such as non-
linearity and bandwidth. More difficult to examine at the design stage
are problems for which inadequate design data exist, as is the case in
estimating the severity of Doppler effects from the various atmospheric
particle velocities. Also, investigations into actual operating charac-
teristics such as ease of data interpretation cannot be made until hard-
ware is built for evaluation.
22
-------
In this section the reduction of the basic CW laser radar to prac-
tice will be discussed under several subheadings, as follows: laser-beam
modulation, laser transmitter, atmospheric effects, laser receiver, and
signal processing. The intent here is to characterize the practical de-
sign problems of a CW lidar. Detailed design data are best included in
a later section.
1. Laser-Beam Modulation
The choice of modulation parameters, technique, and equipment
is crucial to the performance of a CW lidar. It will be helpful to
examine first some important and generally applicable effects of realistic
modulation waveform parameters. The discussion of practical ways to put
these waveforms on an optical carrier can then more profitably follow.
a. Modulation Waveform
The linear frequency sweep shown in Figure Kb)—funda-
mental to the multiple/distributed-target detection capability of the
FM-CW radar—cannot be sustained forever. Practical modulators have both
upper and lower limits on f , as well as maximum sweep rates S and im-
o
perfections in sweep linearity. Since the modulated frequency f cannot
o
increase indefinitely, at some time it will have to return in some fashion
to the desired minimum "start" frequency. The fact that the sweep lasts
only for a period T and covers a finite total deviation Af throws an un-
certainty into the range measurement.
As is obvious from Figure 8(a), if the target's range delay
time t is greater than the duration of the sweep T, f (t) will be the
R R
difference between the transmitted FM signal and whatever frequency con-
tent exists at the receiver input prior to the echo return. An equally
useless difference function results when the FM echo mixes with the
23
-------
0
LU
D
O
LU
(L
Q
LU
I-
K
U)
Lt
1-
ECHO SIGNAL
TRANSMITTED SIGNAL
TIME
(•I
u
a
LU
Lt
LU
ID
-tR-
-( T-i
TIME
">' SA-1979-6
FIGURE 8 RESULT OF FINITE LIMITATIONS ON LINEAR SWEEP
24
-------
unspecified signal following the transmitted sweep. Sweep duration T
must exceed t by a margin oi a few periods t in order to generate a
R R
constant value of f long enough for signal processing to be completed.
R
For example, a radar useful up to a 500-m range, where t =3 i±s, will
R
need a linear sweep lasting 10 us in order to generate 7 ^s of constant
beat frequency f . But if f does not complete at least one cycle in the
R R
7 us, no signal processor can measure the value of f during (T - t ).
R R
Moreover, if the signal processor is to discriminate between a change in
f and an arbitrary uncertainty Af in measuring f , at least one cycle
R R R
of Af must also be completed during (T - t ) as well. Thus three prac-
R R
tical design constraints result from sweep-duration considerations:
T » t (where t is the maximum value expected) (5a)
R R
T > 1/f + t (5b)
K R
1/Af + tn (5c)
R R
In order to apply these constraints to achieve desired range resolution,
the sweep extent, or total deviation Af , must be considered.
o
We can define the difference frequency f from the
R
parameters listed in Figure 8(b). For a linear sweep rate S lasting T
seconds, where Af = ST, f is equal to the frequency change that occurs
o R
during transit time t , or
R
fR = StR (6)
Since
tR = 2R/c
25
-------
then
f - (7a)
R c
or
2RAf
~
Solving for range R:
cf T
(8)
The uncertainty in range (AR) is related to the uncertainty in difference
frequency (Af ) by
R
cAf T
and recalling that
it follows that
AR , - (10)
o
and
Af * TT (ID
o 2AR
26
-------
Equation (10) shows that range uncertainty is a function only of the
maximum deviation Af achievable in a sweep time T.
o
The foregoing analysis applies to measurement consisting
of a single sweep through Af in time T. As in the example above, T
needs to be only a few tens of microseconds in order to meet the con-
straint posed by the maximum 3-us transit time to a target at 500 meters.
Laser modulators such as that used in this study can readily achieve a
linear sweep through a Af of 5 MHz in 100 \±s. From the above expression
o
for uncertainty AR in terms of Af it follows that one sweep through a
5-MHz deviation gives a 30-meter range uncertainty. However, if a large
number, n, of 5-MHz sweeps can be made, a signal processor can make an
improved measurement of beat frequency f during the extended period nT,
R
and thereby reduce the uncertainty in R.
Generating repeated sweeps where Af cannot extend above
o
some practical limit f requires that f be reset at the end of one
o,max o
sweep. Figure 9(a) shows one obvious way of obtaining a series of n
single-sweep measurements of f . The adjoining plot of f versus time
R R
for a single-target situation {Figure 9(b)J shows that this sawtooth FN
waveform repeated at some modulation frequency f gives a relatively
m
constant output value f , with a short burst of error frequency f (t)
R e
occurring every T = 1/f . Since there is a finite reset or turnaround
m m
time t at the end of the sweep, this periodic deviation from the constant
f value lasts for a period t = t + t . Just as sweep duration T was
R d R t
constrained to be longer than t by at least 1/Af [Eq. (5c)], so T must
R R
be reconstrained in the repetitive sweep situation as follows:
T ;> t + 1/Af (12a)
d K
T 2 t + l/fo (12b)
d R
27
-------
Q
D
_!
O.
TIME
(a! SAWTOOTH MODULATION WAVEFORM
TIME
(b) MIXER OUTPUT FREQUENCY
s PHASE SHIFT INTRODUCED
/ DURING DEVIATION
A
•• I \ I \ •
%./ vx %./
TIME
(c) MIXER OUTPUT SIGNAL
of
UJ '
5
2
t
'ECTRAI
Q_
w (
3
fR1
•«-«m-*-
Mil
BEAT FREQUENCY. fR
(d) MIXER OUTPUT SPECTRUM
FIGURE 9 MIXER OUTPUT SPECTRUM
28
SA-1979-7
-------
The periodic deviation in f has a much more profound
R
effect on the design of a practical FM-CW lidar than this time constraint,
however. Notice in Figure 9(b) that f (t) is now a frequency-modulated
function, deviating some Af at a repetition frequency f around an
e m
average frequency that can be shown to be f = 0. Fourier analysis of
an FM signal waveform shows that the spectrum contains an infinite number
of sideband components separated by f and clustered about the average
m
frequency. Figure 9(d) shows such a spectrum for the function f corre-
sponding to the single-target situation modeled in Figures 9(a) and 9(b).
It is no longer obvious from the spectrum display that a single target
is located exactly at Rl. Unless it is possible to narrow the envelope
of the sidebands around the "true" range frequency f , or otherwise
Rl
eliminate spurious spectral peaks, useful single-target detection would
be difficult and multiple-target discrimination impossible.
In designing a practical solution to the sideband problem,
it will be helpful to develop an understanding of the sideband spectrum
and its causes. A rigorous mathematical analysis is lengthy and is not
needed here to define and solve the problem. However, a few key points
will be developed mathematically here. Figure 9(c) shows the amplitude-
vs.-time plot of the waveform whose frequency f (t) changes as in Figure
R
9(b). Notice the amplitude and frequency irregularity during t . Notice
d
that the constant-frequency signal after the deviation does not neces-
sarily pick up in phase with its last value before the discontinuity.
It is these amplitude and phase irregularities caused by the periodic
deviation at rate f that generate the numerous Fourier sideband
m
components.
The spectrum of this frequency-modulated mixer output
signal is presented in Figure 9(d). Notice that the spectral power levels
of the sidebands have an envelope characteristic. As the. nature of the
29
-------
f deviation changes, the modulated signal waveform and its irregularity
R
also change, requiring a new sideband power distribution to generate the
new Fourier sum. During this study we found that the factors determining
the deviation in f (t) can indeed be altered to "improve" the resulting
R
FM spectrum.
The deviation of f (t) shown in Figure 9(b) can be charac-
R
terized by:
• Its duration, t = t -I- t
d R t
• Its waveform, f (t)
' e
• Its repetition rate, f
m
• The initial phase value cp at the start of t .
The practical factors controlling t are the target range R = ct /2, and
d R
the modulator's bandwidth (or maximum slewing rate), which determines
reset time t . Factors controlling the form of f (t) during the devia-
t e
tion are the extent and linearity of Af achieved by the sweep rate S
during deviation time t . The repetition rate of the frequency devia-
d
tions is equal to the FM sawtooth repetition rate f . Finally, the factor
m
determining cp is the duration of the constant frequency f up to the
1 Rl
start of t .
d
In order to single out a target on the spectrum display
of F , spurious sideband power must be reduced in such a way as to single
out the real target frequency as calculated ideally by Eq. (7):
2RlAf
f
Rl cT
At first glance this appears to be a discouraging task since many con-
trolling parameters do not appear in the idealized equation. Fortunately,
several parameters can be combined to give a simpler, approximate but very
30
-------
useful result showing that Af can be adjusted to minimize sidebands.
o
The dependence of target resolution on Af is borne out in practice as
o
well.
It can be shown that the starting phases of all segments
of the mixer output f are equal to each other. Hence, for phase co-
Rl
herence and the consequent reduction in spurious sideband levels, it is
necessary and sufficient that
fRl
an integer.
Combining this with Eq. (7) gives
2RlAf T
Under these conditions, there must be Just enough time for an integral
number of cycles of f to occur during the full period T . in order for
Rl m
phase to match at the endpoints of each T . In other words, for minimum
m
sidebands:
T /T = n (15)
m Rl
where n is an integer and T = 1/f . Since n is also a function of Af
Rl Rl o
[Eq. (14)], a handy design rule is that a single target will best stand
out from surrounding sideband components when
Af ~ = n . (16)
o c
31
-------
But this also means that f_,/f = n is a condition for minimum sidebands,
Rl m
so that only when Af is adjusted to make f a multiple of f will a
o Rl m
target stand out of its own sideband "noise." A radar operator with a
single knob to adjust Af could therefore isolate his target by "detuning"
the spurious sidebands with the knob.
Conversely, if Af is fixed along with f , there will be
o m
several target ranges R such that for
n
- <»)
o
clean radar target displays will result with no adjustment of Af. These
minimum-sideband ranges can aptly be termed "node ranges," and the spacing
between two adjacent node ranges is:
AR = R - R
n n+1 n
° (n + 1 - n)
AR = — - • (18)
n
For a single stationary target between two adjacent nodo
ranges, the sidebands at nf are largest for the values of n that make
m
nf close to the value of f calculated ideally from Eq. (7). The
m R
apparent "average" location of these larger sidebands is f , and if
R
their envelope characteristic is drawn as in Figure 10(b) the envelope's
peak lies at f . although no spectral component is found there. As the
Rl
target moves between node ranges, the larger sideband amplitudes change
as if to roll their envelope characteristic along the spectrum with the
changing target range (ignoring Doppler effects, of course). When the
target reaches a node range where fR = nf . the envelope narrows around
"•n. m
32
-------
«r
LU
Q
D
0.
5
_j
cc
I—
o
LU
CL
VI
(a)
, 1
i i 1 E
/MINIMUM SPURIOUS
/ SIDEBAND AMPLITUDE
/ FOR TARGET AT
NODE RANGE
I , n
1 1 1 i
BEAT FREQUENCY. fR
1
1
<"• 1 ^.CALCULATED TARGET FREQUENCY
i
LU
Q
h-
tt.
5
-------
the correct node-range component. This effect is depicted in Figures
10(a) through 10(d). It is instructive to examine range uncertainty in
this situation.
Since the sawtooth sweep is repeated and sideband com-
ponents f exist continuously, a sufficiently sophisticated signal
processor could measure n of the larger sideband amplitudes, A , calcu-
late the product f A , and determine an average range frequency
k+n
which would, for increasing n, get arbitrarily close to the ideal value
of f resulting from a single continuous sweep during the long measurement
R
time. But to the eye, or to a simple data processor that merely cate-
gorizes the target frequency as lying between the two largest sidebands,
the uncertainty in range is just the nods range spacing AR = c/2£f.
n
The reader has probably recalled by now the result derived earlier for
range uncertainty AR in a measurement made during single-sweep T through
Af [Eq. (10)]:
AR = c/2Af
o
which is just the node range spacing for the repetitive-sweep case. It
is apparent that, even though the repeated sweep sawtooth makes data
available for calculating accurate f , an unsophisticated signal processor
R
may not allow range resolution any better than that possible from a
single-sweep measurement.
34
-------
Another kind of accuracy limitation is also created by
the sideband noise. Recall the intensity-modulation technique assumed
here. The amount of power reflected back by a single target is propor-
tional to laser-beam intensity, and is therefore proportional to the
amplitude (not power) of the HF radar signals that modulate the optical-
carrier intensity and are returning from the target range. Thus it is a
display of spectral-component amplitude that gives a direct measure of
received laser power reflected by targets at various ranges. Now, if
two targets lie in the lidar's field of view, the sidebands of one will
add to the main component of the other, making it difficult to determine
accurately the ratio between the two target power levels. Figure 11
shows how such an inaccuracy develops. It should also be obvious that
the inaccuracy is minimized when spurious sideband levels are lowest--
namely, when the discrete targets are located at node ranges.
It was at these node ranges, for example, that discrete,
semi transparent targets were placed in this experimental study to demon-
strate the remote opacity-measurement capability of an FM-CW lidar, as
discussed in Section V. In practice, spurious-component amplitudes
could be reduced to between 5 and 10 percent of the main node-range
component. In an effort to improve upon this by further reducing side-
band levels, another kind of sweep waveform was investigated—the triangle
wave.
The desired effect of the triangle-wave sweep is to reduce
the duration of the error-frequency burst f (t) during t (Figure 12).
e d
Moreover, since the sweep function during turnaround time t is symmetric
in every respect with the function during the basic sweep time T, the
value of f (t) during t is just
e t
f (t ) = -f (20)
e\ t/ R
35
-------
Q UJ*
0
cc o
uj 1-4
2 J
£ IE °-
0 ° <
RECEIVE!
SPECTRAL
3 M
P,-4
~
1 1
1 1 1 1
SIDEBAND AT (R2
,HAS HEIGHT =• »
1 1 1 l_
jr 6 r
UJ
Q
<6
BEAT FREQUENCY, fH
(a) TARGET AT R1 REFLECTS FOUR UNITS OF OPTICAL POWER
t 4
_J
a.
SPECTRAL
3 10
-
-
|
|
/SIDEBAND AT fR1
HAS HEIGHT - 1
, I
|
I
V3
1
L
BEAT FREQUENCY, fn
-------
LU
D
O
UJ
Q
LU
t
5
TIME
(i)
1
cc
*- +f~
2UENCY,
> 3
a u
a:
a.
< ,
LU 'R
CD
1
/ /
/ 180° PHASE SHIFT IN^-
{_ fR DURING DEVIATION
•*"*£)•*•
\
-\
Notic*
ld • IR
/
!/ TIME
(b)
SA-1979-10
FIGURE 12 FM SWEEP USING TRIANGLE WAVEFORM
37
-------
In this case f (t) is either a constant +f or -f except during the re-
R R R
maining deviation time t . A sweep rate during t equal to the sweep
R t
rate during T apparently reduces the undesired effects of bandwidth or
slewing limitations that occur in a practical sawtooth modulator due to
its rapid reset time. A linear sweep produces the desired constant value
|f | for a larger fraction of the repetition period T in a triangle
wave than in a sawtooth waveform. Indeed, in practice the spurious
sidebands generated in a real mixer output by triangle-wave modulation
could often be reduced to as little as 3 percent of the amplitude of the
!n;iin node-range frequency component.
It can be shown that the triangle-wave sweep unfortunately
introduces a strong dependence on f , the center frequency of the trans-
mitted frequency sweep (Figure 12). Up to now. f did not figure sig-
oc
nificantly in the analysis or performance. Suffice it to say here that
because triangle-wave measurement of If I is made during both its posi-
K
tive and its negative phases, the 180° phase shift it undergoes during
the deviation period must enter into the calculation of phase-coherence
conditions leading to Eq. (13). These values in turn define the condi-
tions for minimum sideband power. A designer may choose to accept the
f dependence in order to gain a few percent in the sideband reduction.
oc
He should, however, be prepared to deal with the f sensitivity in cali-
oc
brating and operating the lidar. The sensitivity to this parameter greatly
reduces the number of practical combinations of &f and f that can be
m
"tuned in" by the operator to put a node range at the actual target range.
Both sawtooth and triangle sweep waveforms are suitable
for multiple-target situations where each discrete target is spaced several
node ranges away from its nearest neighbor. However, closer target
spacings would cause one target's sideband envelope to overlap the other
target's envelope, making it difficult both to resolve the two targets
38
-------
and to measure their relative reflected power levels accurately. Where
one target reflects 20 to 100 times the power of another, its sidebands
alone could dwarf the main component of the weaker target. Given a
minimum target spacing, AR , one might attempt to improve target reso-
min
lution by increasing &f until range node spacing, AR = c/2Af , is
o no
narrowed to a fraction of AR . Eventually the practical limit on Af
mm o
will be reached and no further sideband envelope narrowing can be achieved
(Figure 13). If at this point, AR > AR . , the multiple-target range
n mm
profile becomes inaccurate unless other steps can be taken to reduce the
spurious sideband component levels.
b. Modulation Techniques
Once an appropriate sweep modulation waveform has been
chosen, the techniques for applying the waveform to the optical-wavelength
carrier must be specified. There are two alternative design routes:
(1) modulate the optical frequency (color) of the transmitted laser beam
itself, or (2) modulate a lower subcarrier center frequency that in turn
is used to modulate the intensity of the laser beam.
The necessary prerequisite for optical-frequency modula-
tion is a tunable laser transmitter (see Figure 14). In the visible
region, such sources are now available, offering Af greater than tens of
gigahertz with single-frequency stabilities of 100 MHz. Modulation is
usually accomplished by changing the laser's cavity length, etalon spacing,
or end-reflector wavelength. The laser outputs can be swept through
large Af at audio rates by using acoustic or piezoelectric transducers.
Such lasers are uniformly expensive, and only very recently have they
emerged from the laboratory development stage.11 The linearity and re-
peatability of the sweep available from these lasers is difficult to
measure and is seldom if ever specified or measured by the manufacturers.
39
-------
c
1 R1 RANGE
(a) APPEARANCE OF TARGET LOCATION ON IDEAL FM-CW RADAR DISPLAY
CROSS SECTION. O
SPACING ARn • -^— •
n 2Af,
I I I I j
I I
3 R1 RANGE
K (b) APPARENT LOCATIONS OF SPURIOUS SIDEBAND TARGETS AT NODE RANGES
Q
or
TARGET
APPARENT
(
(
c
2A.2 .
. , •*, ,•*- 1 1 1
I ,
3 R1 RANGE
(c) ADJUSTMENT OF Af CHANGES ARn SO THAT R1 BECOMES A RANGE NODE
c
IjS 2Afg
Inii
3 R1 RANGE
(d) GREAT INCREASE IN Af NARROWS SPURIOUS SIDEBAND ENVELOPE
SA-1979-11
FIGURE 13 IMPROVEMENT IN TARGET RANGE RESOLUTION BY INCREASING
DEVIATION, Af
40
-------
TO
TARGET"
FROM
TARGET'
CONSTANT INTENSITY
WITH OPTICAL FM
Illlllllllllll
BEAMSPLITTER
Illillliilll!
ACCURATELY ALIGNED
MIRROR AND COHERENT
MIXER OPTICS
TUNING
TRANSDUCER
FM SWEEP
WAVEFORM
11 til
COHERENT
MIXER
AND
DETECTOR
».
BEAT-FREQUENCY-SPECTRUM
SIGNAL PROCESSOR
SA-1979-12
FIGURE 14 OPTICAL FM-CW RADAR MODULATION
-------
Certainly the practicality of these lasers is questionable in a low-cost,
field-portable lidar to be developed within three years.
Moreover, optical FM requires optical-frequency mixing in
order to generate the beat-frequency spectrum showing the target range-
vs.-power profile. Such coherent mixing requires very carefully aligned
receiver optics and a photodetector capable of extracting a beat-frequency
envelope that could have a microwave bandwidth in its own right. These
requirements can all be met with advanced components and techniques
;iv;i liable today, but few if any lend themselves well to the low-cost,
durable, field-portable instrument design that is the primary goal in
this effort.
Finally, since the optical-carrier and FM-signal wavelengths
are on the order of target particle size and spacing, the coherent effects
of multiple scattering and optical Doppler shifts could combine to put
complicated and possibly deleterious noise components into the signal
spectrum. The kind of holographic signal processing required to extract
optical-frequency signals from such noise must come from a very new and
still unexplored technology. The design problems posed by a coherent,
optical FM system are very interesting, but their solution is signifi-
cantly more difficult than an overall change in technique.
Given this perspective on optical FM, another modulation
technique appears to be the more promising design approach. By using
the optical carrier only to obtain high lidar target cross section, and
employing much lower VHF or HF intensity modulation to obtain the range-
profile information, the best of existing CW laser and HF electronics
technology can be exploited. This hybrid technique, used successfully
in ranging experiments during this study, might more properly be called
FM-AM-CW radar modulation. Figure 15 shows a convenient way of envisioning
this modulation technique. A carrier signal centered, for example, at
42
-------
FM
FREQUENCY MODULATION'
IN HF REGION
TO
TARGET
~
00
CONTINUOUS CARRIER
AT OPTICAL FREQUENCY f
INTENSITY OF
TRANSMITTED BEAM
FROM
TARGET
OPTICAL POWER
REFLECTED
FROM TARGET
HF
FM SIGNAL
GENERATOR
LASER
INTENSITY
MODULATOR
OPTICAL LASER
OUTPUT FREQUENCY fL
LASER
TRANSMITTER
TRANSMITTED FM ENVELOPE
OPTICAL
POWER
SENSOR
(photodiode)
A/VWVA/>
RECEIVED
HF
MIXER
fe-
BEAT
FREQUENCY
SPECTRUM
ANALYZER
ENVELOPE
SA-1979-13
FIGURE 15 FM-AM-CW RADAR MODULATION
-------
10 MHz, is frequency-modulated through a deviation Af of 6 MHz—i.e.,
from 7 MHz to 13 MHz. This FM signal is then used to modulate the in-
tensity (AM) of a laser-beam carrier operating continuously (CW) at one
optical wavelength. As Figure 15 indicates, the FM signal envelope is
detected by the photosensor and processed by electronics consisting only
of low-cost, medium-frequency components. The radar system remains
linear in the sense that a target whose optical cross section, o, in-
creases by a factor X will reflect an HF signal envelope that is also
larger by X. When heterodyned with the transmitted FM signal, a target
signal at beat frequency, f , will result that is also X times greater
R
in amplitude than before. The equipment needed to implement the hybrid
FM-AM-CW lidar modulation is readily available from well developed com-
mercial technology.
c. Modulation Equipment
From examination of the block-diagram elements in the
modulation scheme shown in Figure 15, it is apparent that the key ele-
ment is the FM-signal generator used both to modulate the laser-beam in-
tensity and to beat down the received signal in the mixer. The requisite
sawtooth (or triangle) sweep can be generated as a voltage waveform by
any of a wide variety of function generators. When this waveform is
used to drive a voltage-controlled oscillator (VOO), the oscillator's
output frequency will have the appropriate FM sweep waveform [Figure
16(a)]. Simple variations in voltage offset and sawtooth amplitude can
be made as adjustments, causing corresponding variations in FM-sweep
center frequency and total deviation, respectively. Likewise, errors
such as drift or nonlinearity in the voltage-function generator will
cause frequency shifts and nonlinear-sweep errors in the FM signal.
Moreover, the VCO will have an input bandwidth that must be considered
in designing the input-voltage function generator. A triangle wave at
44
-------
SAWTOOTH
FUNCTION
GENERATOR
LASER
/MXI
BEAM INTENSITY
'in
*
VOLTAGE-
CONTROLLED
OSCILLATOR
1
WlM/l/VUl POWER
*" AMPLIFIER
CABLE
ELECTRO-OPTIC
MODULATOR
MODULATED BEAM INTENSITY I
(a) A PRACTICAL CW LIDAR MODULATOR DESIGN
100%
O
2
-------
repetition frequency f will be accurately reproduced by a system whose
m
bandwidth accommodates the triangle's primary Fourier components at f ,
m
3f . and 5f . For example, a VCO with an input bandwidth of 100 kHz will
m m
not pass the third and fifth harmonics of a 50-kHz triangle wave, and its
resulting FM sweep will resemble more a sinusoidal sweep at 50 kHz. As
discussed earlier, a nonlinear sweep like this tends to add spurious
sideband components to the target display of an FM-CW radar. Therefore,
it is imperative that the sweep repetition rate fall well within the VCO
input bandwidth.
Once a clean FM signal has been generated, it needs to be
sent to both the laser-beam intensity modulator and the mixer. If the
connecting cables do not have a flat response over the sweep frequency
range, an additional AM component will be introduced on the FM signal.
This FM-to-AM conversion can be a significant cause of spurious sidebands
in the target-spectrum display. If buffer amplifiers or power amplifiers
are used to interface the FM signal with other system components, the
amplifier response must also be uniform over the sweep region Af.
The next key modulation component to consider is the laser-
beam intensity modulator. The most well developed, off-the-shelf com-
ponent is a type of electrooptic crystal such as KDP whose optical trans-
mission changes with increasing voltage applied to opposite faces of the
crystal. The intensity, I, of the output changes with applied voltage,
V , according to the following equation:
S
out J
= sin
in
A judicious choice of the bias-voltage value V and the range of the argu-
b
ment [V + V ] allows operation along the nearly linear portion of the
2
sin X curve, so that the useful modulator function becomes:
46
-------
I
out
1— = KVs + Tc (22)
in
where K and T are constants related to the shape and center point of the
optical transmission characteristic of the modulator assembly (see Figure
16).
The applied voltage changes the transmission of the modu-
lator in a linear, analog fashion. KDP modulators typically offer excur-
sions of 100:1 in transmission at bandwidths up to 10 MHz, and even up to
25 MHz with special driving circuitry. Once again, the nonlinearities
2
introduced by the modulator's sin X characteristic and bandwidth can de-
grade the FM signal appearing on the optical carrier in such a way as to
introduce spurious sidebands in the final target display. The better the
reproduction of the original sawtooth frequency sweep, the closer the
lidar performance is to that modeled in Section IV-A. It should be clear
from the present discussion that the design and specification of modulator
components is very important in reducing system sideband noise. It is
also important to note here that since the transmission of the modulator
assembly is varied about an average value T , the average transmitted
v/
beam intensity is nonzero. The resulting dc component in a square-law-
detected signal envelope has an important effect on system performance
and must be carefully considered during the design of the lidar receiver.
2. Laser Transmitter
The incentive for using an optical radar wavelength lies in the
nature of the intended target. Since a radar designer's goal is to maxi-
mize the amount of signal power to be detected after reflection or back-
scattering by the target, the designer should choose an operating
wavelength smaller than the target dimensions. For particulate targets
such as the aerosols and molecular constituents of the lower atmosphere,
47
-------
relatively high backscattering and cross sections can be obtained for
wavelengths less than 10 ^im, at which a mix of Mie and Rayleigh scattering
phenomena occur. As the radar wavelength shortens, the bulk-scattering
cross section of the particles rises, but the corresponding path attenua-
tion increases until in the far ultraviolet regions the useful lidar
range may be too short to be practicable. On the other hand, in the in-
frared region the backscatterlng cross section is lower by an order of
magnitude or more. A convenient middle ground is the visible region,
particularly because presently acceptable techniques for judging smoke-
plume opacity are visual. A visual technique is one involving the human
eye, whose response to light peaks at green wavelengths. Because these
visual methods typically involve comparison of plume darkness against
the scattered blue light from the sky, a laser wavelength in the green
is a reasonable choice for a CW lidar opacimeter that is to complement
or even supplant existing remote-measurement techniques. Also important,
most photodetectors are more efficient and sensitive to green light than
to longer, red wavelengths.
Power output levels from 10 mW to 1 watt can be supplied by such
commercially available CW gas lasers as the HeCd and argon types operating
at various wavelengths between 4416 and 5145 A. The stability, durability,
ease of handling, and lower cost of HeCd and argon lasers also make them
likely prospects for use in a mobile or field-portable instrument. CW
tunable dye lasers have recently appeared on the market and may prove
suitable in the present application should their prices drop significantly.
They might also provide a means to identify the chemistry of the effluents.
However, for basic plume transmission measurements at a single wavelength,
a tunable CW laser would be more sophisticated than is required.
The power level actually required from the laser transmitter
is constrained at the upper extreme by laser-performance and eye-safety
48
-------
considerations and at the lower end by receiver sensitivity and back-
ground noise.
HeCd lasers are now available that can transmit up to 50 mW at
4416 A. Argon lasers with rated power outputs from 25 mW to 10 watts
can transmit at any of a number of wavelengths: 5145, 4880, 4765, and
4579 A. The trend toward standardizing transmissometer sources to the
o
green region favors use of the argon 5145-A line. Such factors as tube
life, cooling requirements, line stability, and electrical power drain
can limit the choice of such an argon laser to one of the lower power
models for mobile or field-portable applications. Also, the ability of
the electrooptical modulator crystal to handle the rejected beam power
during its attenuating half-cycles limits the amount of power that can
be transmitted by an FM-CW lidar.
Finally, eye-safety considerations indicate that the trans-
2
mitted beam intensity should not exceed approximately 1 m\V/cm . This
number is derived from data showing the interrelations between continuous
source radiance, retinal irradiance, retinal image size, and thermal
damage (if any) done before the eye's natural blink reaction occurs.12
For example, if the lidar receiver design shows that one watt of power
must be transmitted to a target in order to get a detectable return sig-
3 2
nal, the transmitted beam s cross-sectional area should be 10 cm or
greater. This in turn means that the transmit lens must be some 40 cm
in diameter in order for the one-watt beam to be eye-safe even for per-
sonnel looking directly into the transmit aperture. Given the receiver's
sensitivity and the target's lidar volume backscatter cross section, eye-
safety considerations then lead to the first step in the design of the
transmitter optics: choosing the output-aperture dimension. The re-
mainder of the transmitter's optical design is straightforward and,
since it has no peculiar dependence on the CW technique, will not be
pursued here.
49
-------
Internal laser noise is one additional characteristic that can
have a serious effect on CW lidar performance. Since the CW lidar target
information is contained in the frequency spectrum, the spectral charac-
teristics of noise become important. Appearing as either spectral base
line "grass" or as spectral line-broadening frequency instability, the
noise components can interfere with the desired target signal spectrum.
Laser amplitude noise can be caused by the 60-Hz harmonics
riding on top of the dc gas ionization or exciter voltage. Amplitude
noise can also appear in the VHF region of the spectrum at a beat fre-
quency between two laser modes. Since both of these AM noise sources
steal carrier power and introduce it as spurious sideband components,
the designer must either specify components that reduce the noise level
or choose modulation parameters so that the spurious noise components
will not fall in the target signal spectrum, or both.
In a similar fashion laser frequency instability can cause a
kind of FM noise that is converted into AM noise in the laser itself.
The severity of this internal FM-to-AM noise conversion is difficult to
estimate and even more difficult to measure in the presence of other
FM-CW system noise sources. Laser manufacturers claim that laser FM
noise is generated by acoustical disturbances to the laser cavity caused
by the flow of cooling air or water, loud noises, or random motions of
the laser mount. Since all these phenomena occur in the audiofrequency
range, such internal noise could show up as low-level audiofrequency
components on the detected high-frequency modulation envelope. On the
other hand, atmospheric effects can also convert the VHF frequency in-
stabilities directly into VHF AM noise, as discussed later.
50
-------
3. Optical Receiver
The design of an optical receiver for incoherent, square-law
detection of backscattered laser signals has been amply discussed by
others.2'13 A few receiver characteristics, however, have effects on
performance that are peculiar to an FM-AM Ctt' lidar and should be examined
here.
The modulated laser-beam intensity v^s expressed in Eq. (22) as
I
o
= I. JKV + T (22a)
in|_ S CJ
Since V = V sin[2nf t'j , where f is the FM sweep function, the average
S o o o
value of I KV is zero. However, the quiescent modulator transmission
in S
value T is nonzero, and therefore an average nonzero intensity equal to
O
I T is continuously transmitted. Since the laser modulator input in-
in C
tensity I is very nearly constant, consider I instead as
in o
-------
which can be redefined as two power components
P = P + P (24)
R S C
The most sensitive, low-cost light detectors for use in this
application are square-law detectors such as a photodiode or photomulti-
plier tube (PMT). The current output of the detector is proportional to
incident power, so that in the case above, output current is:
(25)
where i is the FM sinusoidal function and i is a constant. The rela-
S C
tive value of i /i is just the ratio of modulator parameters K T /K V .
C S C C o o
This value typically averages 1.5 and higher for practical electrooptic-
crystal laser-beam modulators. As a rule, the better the modulator's
reproduction of the input signal, the higher the value of K T /K V . The
C C S S
resulting effect is that the average dc component, i } could exceed the
\j
continuous current rating of the detector for a bright lidar target.
Another result for a slightly smaller target return would be that i
\s
could drive the detector nearly into saturation, thereby reducing the
effective gain for the smaller, information-bearing current, i .
O
Even when both i and i are in the dynamic range of the de-
S C
tector, i creates a steady, additional shot-noise component above that
Vx
of i . In short, the relatively larger reflected power component P
i> C
contains no information, yet it can significantly reduce the lidar system's
signal-to-noise ratio (SNR) in the case of a strong target echo. The dc
component can mask the returning, information-bearing FM sinusoidal sig-
nal, P . The design of the lidar receiver should therefore take into
b
account both the modulator parameters K T and K V and the maximum ex-
C C S S
pected target reflectivity. As a practical matter, since the lidar
52
-------
operates continuously, it lends itself well to incorporation of an auto-
matic gain control (AGO to reduce PMT gain for bright-target conditions.
An AGC could keep both i and i in a range where i can be separated,
s c s
undistorted, from i in signal-processing circuits.
\^r
Considering the distributed nature of the atmospheric target
particles, an experienced lidar designer would point out here that cer-
tainly the backscattered signal from ranges very close to the receiver
would be decades stronger than those from distant targets, and would
therefore mask remote-target information whether there is an AGC or not.
The solution for a pulsed lidar is, of course, to make sure that the re-
ceiver's field of view does not intersect the transmitted beam for some
appropriate distance from the receiver. This is not a complete solution
for a CW lidar, however. The continuous average signal effects discussed
in the foregoing paragraphs combine in the receiver's field of view to
introduce a spatial integration effect on the CW lidar's distributed-
target performance at any range.
The modulated laser beam, frozen in an instant of time, would
appear as a column of light with alternating light and dark sections
spaced at the modulation wavelengths (Figure 17). For a constant modula-
tion frequency f , the peaks are spaced by \ = c/f . A receiver viewing
o o o
one thin slice of one section will see the full maximum and minimum value
of the intensity as the column is allowed to proceed across the field of
view at the speed of light. As the field of view is increased to include
one-half wavelength, the receiver will capture greater total power, but
the relative ratio of the total's maximum to its minimum will fall be-
cause, just as a dark region is leaving the field of view, a light region
is entering. When the receiver views a full wavelength, again the total
average received power increases, but the minimum value of intensity for
the dark section is exactly offset by the maximum values being received
from the light section at any instant of time.
53
-------
INTENSITY I (X) = K sin |27T — I + K T
O I -v / c C
LASER
TRANSMITTER
WAVELENGTH OF INTENSITY MODULATION
RECEIVER
(at SIMPLIFIED ANALYTICAL MODEL
(b) A TYPICAL EXAMPLE OF SPATIAL INTEGRATION GEOMETRY
ia
FIGURE 17 SPATIAL INTEGRATION OF BEAM INTENSITY FUNCTION BY RECEIVER
54
-------
In the expression for intensity of the frozen column of light
as a function of length, it is clear that every time length x increases
by one wavelength \ , one cycle is completed:
(26)
Of course, when the beam is allowed to move again, the length x traversed
by a modulation peak at the speed of light is just x = ct, which, when
substituted in the sine argument above, gives Eq. (23) for I (t).
o
The lidar receiver field of view, Q, encompasses at any instant
some extent of beam length Ax where Ax = 0R at range R. The power re-
ceived is just the integral of the intensity function over the spatial
extent AX that is being viewed:
.x+Ax
Ki J
P = K J K V sin|2TT —I dx + K T dx (27)
R If SO I \ I C C
P = K K V
\>
R 1 S O 2n
\
/2n
| —
\
cos| — xl - cos
K T Ax . (28)
The power received at an instant of time is clearly a function of field
of view, because Ax is the extent of the beam subtended by the field of
view at range R. Also, since the value of x is arbitrary, depending on
where the beam is frozen in space and phase, Eq. (28) can be simplified
by letting x = X /4 so that:
o
*o /2n \
DR = KiVo to cos(r Axj+ KiKc
Ax . (29)
As the field of view is widened, the term K K T Ax increases linearly.
J. kj v>
However, the cosine term periodically peaks at K K V \ /2n every time
ISO o
AX reaches an integral multiple of wavelength \ . The peak value of the
o
55
-------
information-bearing cosine signal can be expressed as a fraction P of
the nonzero average power:
rel K K T <0R) rel i
ICC c
It is clear that as receiver field of view increases, the detected sinu-
soidal signal current drops in relation to the average dc current. Since
the average dc component is a noise source and a limiting factor as dis-
cussed earlier, an increasing field of view acts to degrade system SNR
independent of optical background noise.
Therefore, in a CW lidar system, the receiver ideally must view
less than one-half wavelength of the modulating frequency f . As the
field of view widens, the stronger de-induced noise swallows up the fre-
quency component f .
o
Similarly, when f is a sweep function of time, f (t), in order
for the receiver's response to a distributed target to be uniform for
each frequency in the sweep range Af, the receiver must view at most a
fraction of the wavelength of the highest modulating sweep frequency,
f
o,max
As an example, assume that the maximum sweep frequency is 10
MHz, so that Ax < 15 m is the spatial integration constraint. Assume
also that the target is at 200 m, and that the receiver can be no further
than 5 m from the transmitter [Figure 17(b)]. If the receiver is to view
only Ax, its beam diameter at the target must be
and therefore field of view Q should be smaller than
56
-------
D 0.4
R =
which is not an unreasonable practical constraint.
4. Atmospheric Effects
Just as pulsed-lidar performance is affected by atmospheric
turbulence and nonuniform aerosol or particulate concentration, the C\V
lidar beam undergoes these atmospheric effects. Unlike the pulsed sys-
tem, however, the CW lidar continuously reads out the magnitude of error
introduced by these effects. For example, if, in the idealized lidar
of Section IV-B, the backscattering coefficient of the near clear-air
volume changes while that of the far volume does not, the continuous
readout of [log P - log P ] would vary accordingly about an average
near far
value. A chart recording of this relative reading would show the magni-
tude of excursions from the average value, giving a quick and reliable
measure of system accuracy under changing atmospheric conditions. An
operator would know in a minute or two whether conditions were too turbu-
lent for adequate plume-opacity measurements. Indeed, microwave FM-CVV
radar has been used recently specifically to measure air turbulence,1 "*• lp
Additional atmospheric phenomena may have deleterious effects
on the beam's modulation envelope. In this experimental study, noise
from these peculiar effects did not appear significantly above other sys-
tem noise levels. Nevertheless the design of an FM-CW lidar should take
them into consideration.
Particles and molecular constituents of the atmosphere are not
at all stationary when they scatter light. For example, molecular oxygen
4
at 273° K consists primarily of molecules moving at speeds between 10
5
and 10 cm/s. Recall that a moving target or scatterer shifts the fre-
quency of any signal that it reflects along its velocity vector. For a
57
-------
radar signal frequency f , velocity v, and speed of light c, the Doppler-
shift frequency f is:
14
For a laser beam at f = 6 X 10 Hz, where f is the optical carrier
o o
frequency, oncoming oxygen molecules could shift the received signal fre-
quency by as much as 2 GHz. Since the various scatterers in a volume of
clear air are moving with a continuum of velocity components along the
lidar-beam axis, the Doppler shift in the color of the reflected laser
light would be a continuous spectrum from 0 to 2 GHz around the optical
carrier frequency. In an optical FM-CW radar such as that shown in Figure
14, the random Doppler shift in the optical FM sweep signal would be on
the order of the sweep deviation Af obtainable from the tunable laser.
Therefore a change in return-signal frequency (color) due to Doppler
shift would be similar in magnitude to the intended frequency changes
due to target-echo delay time. As a result, Doppler effects from the
distributed atmospheric scattering target could cause serious inaccuracies
in the target range profile in a coherent, optical heterodyne system.
In an FM-AM-CW lidar, however, the range information is con-
tained not in optical FM but in the VHP or HF modulation envelope around
the optical carrier's intensity (Figure 15). The photodiode receiver
throws away the optical-carrier frequency information; changes of 20 GHz
14
in a 6 x 10 Hz laser frequency would not be detected in such a receiver.
Only the HF envelope signal is processed. Its highest frequency during
the sweep is on the order of 10 MHz, and it too undergoes a Doppler shift
due to atmospheric particle velocity:
f = — — — - X 107 = 33 Hz
3 X 10
58
-------
This shift is small compared to the value of Af = 5 MHz used in a prac-
tical FM-AM-CW lidar, and is also small compared to the resultant beat
frequencies, typically 100 kHz or so for a target range of 200 m.
A second atmospheric effect, coupled with Doppler shift, could
conceivably add amplitude noise to the Doppler-frequency noise just dis-
cussed. The severity of this new effect is difficult to calculate;
measuring it was one of the main goals in<\hls recent effort. It is
another manifestation of the FM-to-AM noise conversion mentioned earlier.
The theory behind this second effect is cumbersome to express mathemati-
cally. A qualitative explanation will help here.
Those who have observed an object illuminated by a laser have
undoubtedly noticed a kind of speckle pattern overlaid on the object.
The pattern is caused by constructive and destructive interference of
wavefronts emanating from various point reflectors on the object. Whether
a particular point will appear to be a bright or dark speckle is deter-
mined for a given beam diameter by (1) spacing between points on the ob-
ject, (2) distance to the observer's eye, (3) observer's pupil size, and
(4) the wavelength of the laser light. When the observer moves his head
he sees the speckle pattern change, as if the object were glistening.
He is changing the angle between lines drawn from two given object points
to his eye. As the two path lengths change, the apparent source bright-
ness goes through bright and dark interference cycles. In effect,
changing the determinants of the speckle pattern causes intensity varia-
tions, or "AM noise" at the observer's eye.
In a similar fashion, a stationary lidar receiver viewing
moving point scatterers will see a glistening image. Its intensity
appears to vary as any two scatterers move into the field of view in such
a way as to cause interference fringes to move across the receiver aper-
ture.18 As the aperture is enlarged, the effect decreases because more
59
-------
fringes are being spatially integrated over the aperture dimension, and
any new fringe moving in amounts to a smaller percentage of the whole
pattern. Similarly, two scatterers in the field of view may have Doppler-
shifting wavelengths as well as shifting positions that can cause speckle
variations. The glistening speckle effect caused by billions of indivi-
dual scatterers at various speeds, ranges, and angles-off-axis might be
complicated and serious enough to cause significant AM noise in the re-
ceived signal. Theoretical analysis of this point becomes quite diffi-
cult, and in the absence of experimental data one could speculate that
the AM noise might be dominated by the optical Doppler shifts, which
vvould add noise with bandwidth up to 2 GHz on the return signal. Or the
AM noise might be caused only over relatively small spatial regions, so
that only a very small receiver aperture or narrow field of view would be
adversely affected. The effect might be serious for an optical FM-CW
lidar with its need for coherent heterodyning, while not showing up at
all in an FM-AM-CW system.
As a practical matter, during this experimental study, no
serious Doppler or FM-to-AM conversion noise occurred to degrade system
performance. If any of these noise phenomena occurred, their contribution
was below that of other readily identified and correctable system sources.
It is not certain, however, that an optical design for an FM-AM-CW system
somewhat different from the one used here would also be free of this
noise. Nor is it clear that a purely optical FM-CW coherent radar would
escape these effects. Future CW lidar designs should not ignore the
potential for such Doppler spreading and FM-to-AM noise conversion to
degrade system performance. However, actual results show that at least
one FM-AM-CW lidar design is largely free of these adverse atmospheric
effects.
60
-------
5. Signal Processing
The foregoing sections have discussed the design of an FM-AM-CW
lidar without specifying the manner in which the received signal frequency
is subtracted from the transmitted signal, nor the method for displaying
and measuring target range and returned signal power. Some theoretical
considerations of measurement uncertainty were discussed in Section IV-
C-l-a, but. an explanation of the signal-processing functions required in
an FM radar still needs to be presented.
In Section IV-C-1 the advantages of VHP or HF frequency mixing
as opposed to optical, coherent heterodyning were pointed out. Indeed,
in the FM-AM-CW technique any of a number of inexpensive, commercially
available HF mixers can be used. Such mixers contain nonlinear elements
that multiply the transmitted FM signal envelope with the delayed, re-
ceived envelope to generate sinusoidal signals at the sum and difference
of the input frequencies. Since only the difference component is of
interest the mixer must be followed by a low-pass filter that eliminates
the higher sum-frequency components.
Practical mixers have limits to the dynamic range over which
their conversion of input-signal to ou,tput-signal amplitude is linear.
Care must be taken to guarantee that the photodetector output does not
exceed the mixer's maximum linear input level. Otherwise, comparisons
of various frequency-component amplitudes will be inaccurate. The mixer
will also introduce internal conversion noise around the signal during
the. heterodyning process. If the input signal is small enough, mixer
noise can overwhelm it. There is also a conversion loss through the
mixer that typically attenuates the signal power by a factor of five or
ten. While mixers are functionally simple to use, their detailed input-
output characteristics can be quite important in the design of a CW lidar,
where low-level and wide-dynamic-range signals are expected.
61
-------
Once the mixer output has been obtained and sum-frequency com-
ponents filtered out, it remains to analyze the spectrum of the difference-
frequency signal and noise. As discussed in the foregoing sections, this
spectrum contains the target range and reflected-power information that
is of primary interest to the operator. Signal processing at this point
is required to perform the following functions:
(1) Determine the frequencies of interest.
(2) Extract the target signal spectrum from the noise
spectrum.
(3) Modify and format those signals in such a way as to
provide a convenient measurement of target charac-
teristics.
Function 1 is simple when the system's modulation parameters
and target characteristics are known. A narrow bandpass filter tuned to
the known target frequency will eliminate all components except the
target-frequency component. If there are three targets, three such
narrow-band filters can be used, as in the idealized CW lidar opacimeter
discussed in Section IV-B. A portable AM radio does much the same thing:
it can be tuned to a predetermined broadcast frequency to process the
voice or music signals found there.
When the target spectrum is not known in advance, as in this
research and development effort, it is necessary to scan through the
frequency spectrum until an overall picture of the target range has been
developed in order to find any interesting spectral characteristics. A
spectrum analyzer is a useful instrument for this purpose. It sweeps a
narrow-band filter of some bandwidth B from zero frequency to some maxi-
mum frequency, displaying component amplitude Y versus sweep position X.
While it is certainly beyond the scope of this report to explain spec-
trum analyzer theory, it can be said that the instrument is of great
utility as a development tool. But since it cannot view three separate
62
-------
frequencies simultaneously, one cannot always use the analyzer to make
instantaneous comparisons of the reflected power levels of FM-CW lidar
targets. The analyzer takes some time to sweep from the frequency of
the near target to that of the far target. During that sweep time, the
target characteristics may have changed, and therefore the analyzer sweep
rate introduces uncertainties in measurements referenced to time-varying
signals. For slowly changing signals, however, the scan rate may be fast
enough to reduce the uncertainty to an acceptable level.
Functions 2 and 3 are straightforward electronic design prob-
lems and are not peculiar to the FM-AM-CW lidar technique. The idealized
CW lidar for measuring smoke-plume opacity shown in Figure 7 includes
functional signal-processing blocks that are readily achieved with inex-
pensive commercial components. In the CW lidar design the need for compli-
cated signal processing can be reduced by careful design elsewhere in the
system. A good choice of modulation and frequency parameters would allow
the target spectrum components to fall away from, or in between, noise-
spectrum components. For example, the so-called semiconductor 1/f noise
spectrum diminishes in strength at frequencies above 10 kHz, while the
modulator's power-amplifier signals may leak into the processor circuits
as noise in the frequency range above '1 MHz. By choosing modulation rate
f and deviation Af carefully, the designer can cause the desired beat-
m
frequency spectrum for targets between 100 and 300 m to lie between 50 and
500 kHz, thereby avoiding the given 1/f and leakage-noise spectra.
63
-------
V EXPERIMENTAL PROGRAM
The objective of the experimental program undertaken here was to
determine the engineering requirements necessary to design, fabricate,
and operate a CW lidar with range and amplitude resolution sufficient
to obtain transmittance measurements on semitransparent targets at ranges
of 500 to 1500 feet.
After surveying some of the analytical groundwork described in the
foregoing sections, the intensity-modulation, or FM-AM-C\V, approach was
selected over the less easily implemented coherent optical FM-CW tech-
nique. The next task included the design and construction of a laboratory-
model CW lidar with flexibility for doing research over a wide range of
parameters. With the completion of the lidar, preliminary system-noise
and performance-limitation measurements were made. These measurements
allowed the specification of optimum operating regions, and in these
regions the first ranging and remote opacity measurements were made using
the FM-AM-CW lidar. The details of these experimental program activities
and their results are included in the following pages.
A. Design and Construction of a Laboratory-Model CW Lidar
The FM-AM-CW lidar technique chosen is depicted in Figure 15, and
the modulator technique is detailed in Figure 16. For convenient refer-
ence this system has been redrawn in Figure 18, where each component is
specified by function, manufacturer, and model number.
Item 1, the sweep function generator, is needed to provide very
linear voltage ramps because the required linear frequency sweep f (t),
is patterned after this voltage waveform by the voltage-controlled
65
-------
POWER
AMPLIFIER
©
TO BEAM-
EXPANDING OPTICS
ARGON LASER
©
BUFFER AMPLIFIER
VOLTAGE (?)
CONTROLLED
OSCILLATOR
FROM
RECEIVER OPTICS
©
PHOTOMULTIPLIER
TUBE
SWEEP
FUNCTION
GENERATOR
©
LOW-PASS ©
FILTER
ITEM
Alfred Model 325A tor Wavetek 114) Function Generator
Exact Electronics Model 7030 Voltage Controlled Generator
Coherent Associates Model 30 Modulator Driver Unit
Spectra-Physics Model 164 00 2-Watt Argon Ion Laser
Coherent Associates Model 27 Electro-optic Modulator
RCA Type 7265 (or 6810 A) Photomultiplier Tube
Relcom Model M1 Mixer. 0.2-500 MHz
2-MHz Five-pole Low-pass Filter
Hewlett-Packard Model 8552/85538 Scanning Spectrum Analyzer
Tektronix 453 Dual-Trace Oscilloscope
SCANNING (»)
SPECTRUM
ANALYZER
DUAL-TRACE
OSCILLOSCOPE
SA-1979-16
FIGURE 18 LABORATORY MODEL FM-AM-CW LIDAR EQUIPMENT BLOCK DIAGRAM
-------
oscillator, Item 2. The specified VCO had a remarkably linear voltage-
to-frequency transfer characteristic. Its primary limitation was its
input bandwidth, which tended to round off the sharp turnaround points on
a triangular or sawtooth sweep waveform. This input bandwidth of 200 kHz
limited the sweep repetition rate to f < 20 kHz in order to maintain
m
sweep linearity and small relative turnaround time t . The majority of
measurements were made using a 10-kHz modulation rate. Buffer amplifiers
were used to reduce the load at the VCO output and thereby preserve sig-
nal purity.
The electrooptical modulator crystal (Item 5) has the measured
voltage-to-transmission curve shown in Figure 19. The crystal presents
a 65-pF capacitive load, requiring a rather powerful driver stage to pro-
vide the ±50-volt excursion around a bias point at 75 volts that makes
maximum use of the linear portion of the modulator transfer curve. The
CO
cr
o
it
o
Q
O
z
o
in
Z
<
a:
v-
20 30 40
50 60 70 80 90 100 110 120 130
APPLIED MODULATOR VOLTAGE — V
140 150 160 170
SA-1979-17
FIGURE 19 MEASURED MODULATOR TRANSMISSION CHARACTERISTIC
67
-------
specified power amplifier, Item 3, is designed to drive the crystal at
frequencies between dc and 10 MHz, a convenient range since the maximum
frequency available from the VOO is 11 MHz. Moreover, spatial-integration
effects in the receiver field of view dictated that the half-wavelength
of the highest intensity-modulation frequency be no shorter than approxi-
mately 30 m to encompass reasonable smoke-plume targets. This factor
also supports the choice of 10 MHz as an upper limit to the modulating
frequency.
To achieve research flexibility, a durable, tunable, 2-watt Argon
io>% laser transmitter was specified (Item 4). The unit provided plenty
of power for most experiments at 5145 A. No internal laser noise could
be measured save for 60- and 180-Hz ripple variations in the laser output
intensity. These variations were well below the 60-Hz noise introduced
by other components in the system such as the photomultiplier-tube power
supply (Item 6). Only by purposely detuning the laser could one generate
a beat frequency between two laser modes, but even then these beat fre-
quencies were greater than 25 MHz, well out of the range of the overall
10-MHz system bandwidth.
The photomultiplier tube (Item' 6) originally chosen for use as the
laser receiver had an S-ll photocathode that introduced too much dark-
current noise for adequate measurement of low-level backscatter from
clear air. After switching to an S-20 phototube, it was possible to de-
tect light backscattered from clear-air volumes at 100-m ranges.
Photomultiplier output current signals were buffered and mixed with
the FM sweep drive signal supplied by the VCO. The output frequencies
of the mixer consist of the sum and difference of the two input-signal
frequencies. Since only the difference frequency is desired, a low-pass
filter (Item 8) was used to eliminate the higher sum frequencies and pass
68
-------
the lower beat-frequencies containing target-range information on to the
display devices (Items 9 and 10).
Figure 20(a) shows the two kinds of data display employed in the
operation of the CW lidar. The dual-trace oscilloscope display shows
the voltage waveform (bottom trace) that drives the VOO; therefore this
trace is an indicator of the FM sweep waveform, deviation Af> and repeti-
tion rate f . The top trace is the output of the mixer showing the
m
Fourier sum of all target-frequency components as well as the error-
frequency components arising during sweep nonlinearities (note the blur
in the trace corresponding to the turnaround points of the triangular
sweep waveform). The second photo in Figure 20(b) is the spectrum-
analyzer display (Item 9) of the frequency components comprising the top
trace of Figure 20(a). Three main target frequencies appear, each of
which is 10 dB above its surrounding sidebands. Indeed, when these photos
were taken the lidar was viewing three targets downrange, and Af had been
adjusted to reduce each target's sideband envelope so that the three tar-
gets could be resolved on a display of the beat-frequency spectrum.
Notice that the differences in signal power returned by various targets
show up as height variations from component to component. By measuring
Af, f , and f from these displays, one can calculate target range from
m R
Eq. (7). By measurement of the heights of frequency components on the
spectrum display, the relative power scattered back from the target ranges
can be determined directly.
B. Target Display Characteristics
The goal of the overall effort is to develop a CW lidar for the re-
mote measurement of smoke-plume opacity. The measurement of the distance
to the plume is only incidental to this goal. Nevertheless, since
analysis of FM radar performance shows that single, discrete targets will
69
-------
OSCILLOSCOPE TRACE SHOWING MIXER OUTPUT SIGNAL (top) AND VOLTAGE
WAVEFORM THAT DRIVES THE FM GENERATOR.
—50 kHz
BEAT FREQUENCY — kHz
(b) SPECTRUM ANALYZER DISPLAY SHOWING THREE MAIN TARGET FREQUENCIES
COMPRISING TOPMOST SIGNAL TRACE ABOVE. NOTICE SIDEBANDS BELOW
EACH TARGET SIGNAL. LARGE SPIKE IS THE ZERO FREQUENCY MARKER.
SA-1979-18
FIGURE 20 DATA DISPLAYS USED FOR REMOTE TARGET MEASUREMENTS
70
-------
generate spurious target components in the beat-frequency spectrum, and
since these spurious targets can introduce inaccuracies in the measure-
ment of the relative power reflected by two targets (Figure 11), it was
important to experiment with discrete targets at various ranges before
tackling the distributed targets presented by plume and clear-air scat-
tering volumes.
1. Discrete Targets
The CW lidar equipment shown in Figure 18 was used for these
initial target measurements. The modulated argon laser beam was expanded
to a three-inch diameter and transmitted via periscope through the roof
of the laboratory and out along the target range shown in Figure 21.
Also on the roof and next to the steerable output periscope mirror was
the tripod-mounted lidar receiver, consisting of a 6-inch-diameter front
lens, a collimating lens, and a photomultiplier tube. The mirror was
adjusted to direct the beams down the parking-lot range, and the receiver
was then aligned with the transmitted beam. With the laser power output
reduced well below eye-safe levels, an assistant standing 90 m from the
transmitter held a white card in the beam and began walking away from the
lidar site. A series of photographs were taken of the spectrum analyzer's
display of the difference-frequency spectrum as the target moved slowly
out in range.
The spurious sideband components generated around each of the
single-target returns are shown in Figure 22. The CW lidar's mixer
output-frequency spectrum is displayed on a spectrum analyzer for nine
target positions at successively greater ranges. Each photo presents
the frequency-component amplitude on a logarithmic vertical scale versus
range (beat frequency) on the horizontal. This scale tends to over-
emphasize the sideband strength on first impression. The large signal at
71
-------
BACK
TARGET
REFERENCE
SEMI-TRANSPARENT
TARGET
STEAM
PLUME
,•:••::.-. . BLDG. 1 •.:};'•••»:.•::'.
GLASS
PLATE
REFERENCE
'BLDG.
OBJECT
TARGET
GLASS PLATE
REFERENCE
SA-1979-19
FIGURE 21
LIDAR TARGET RANGE USED FOR REMOTE RANGE AND OPACITY
MEASUREMENTS
left is the zero-frequency marker generated by the spectrum analyzer.
Since the sawtooth sweep repetition frequency employed is f =10 kHz,
m
the sidebands appear at multiples of 10 kHz across the display. The
apparent sideband width is an artifact of the display device.
Notice in Figure 22(a) that the sidebands are approximately 20
dB below the main target component at 50 kHz. This is a result of the
fact that a target located at a range node meets the minimum-sideband
condition discussed in Section IV-C. In Figure 22(b) the target has
moved to a range halfway between range nodes. Notice the increase in
72
-------
(a)
(b)
— 20 kHz
SA-1979-20
FIGURE 22 DISCRETE-TARGET SPECTRA
SHOWING SPURIOUS SIDE-
BAND CHANGES WITH RANGE
73
-------
sideband levels. As the target moves out in range, alternate photos
show the changes in the sideband power envelope for range-node locations
and nonrange-node locations.
An interesting feature shown in Figure 22(h) is that just as
the target moves slightly off the range-node position in 22(g), the first
sideband to the left of the main target component drops below the under-
lying system noise level. That a particular sideband can be eliminated
in this way may be a useful feature in future lidar opacimeter designs.
When the clear-air signal is brought up above the system noise and the
spurious sideband induced by the target is eliminated in this way, any
remaining sideband amplitude will be contributed solely by energy back-
scattered from clear air near the target. In a real smoke-plume measure-
ment, however, the target would not be a white card but a volume of dis-
tributed particle scatterers. A tenuous steam plume near the experimental
target range provided an excellent opportunity for examining such a dis-
tributed target's characteristics when viewed by a CW lidar.
2. Steam-Plume Characteristics
Figure 21 shows the location of this steam plume near the ex-
perimental target range at SRJ. Figure 23 is a photograph of the rooftop
lidar site after the transmit mirror and receiver telescope had been
steered over to view the steam plume. Notice that the beam's path is
clearly visible due to clear-air scattering all the way through and beyond
the plume. The plume itself is very reflective, and its bright front-
surface reflection causes the plume target spectrum shown in Figure 24(a).
Here the two plume sidebands appear approximately 10 dB below the main
plume return. Additional weaker sidebands can be seen at higher fre-
quencies, just above the system and background noise spectrum.
74
-------
SA-1979-21
FIGURE 23 ROOFTOP LIDAR SITE
-------
•50 kHz
BEAT FREQUENCY — kHz
STEAM PLUME TARGET RETURN RISES ABOVE NOISE SPECTRUM BY NEARLY
20 dB. TWO LARGEST SIDEBANDS ARE 10 dB LOWER. LARGEST SPIKE IS
ZERO-FREQUENCY MARKER.
-------
The steam plume return is much like a discrete target return.
The moving water droplets in the plume do not appear to have caused any
severe Doppler spreading in the target spectrum. Figure 24(b) is the
spectrum of a discrete target—a wall—near the steam plume. The two
spectra are quite similar, leading to a preliminary conclusion that
multiple-scattering or Doppler effects on the beam as it passes through
a plume do not make the plume target spectrum significantly different
from a discrete-target display.
The detection of light scattered from clear-air regions in
front of and behind the seemingly discrete plume target proved more dif-
ficult with the elementary CW lidar system built for these initial efforts.
Unlike the bright return from the front surface of the plume, the back-
scatter from extended volumes of clear air posed a problem unique to the
FM-AM-CW technique.
3. Clear-Air Targets
The primary problem in detecting clear air in the presence of
a bright steam plume return was that the plume return masked the clear-
air return. Since target energy appears only as spectrally narrow side-
bands at multiples of f in the beat-frequency spectrum, the strong
m
sideband-energy contributions from the main plume return are much greater
in amplitude than the contributions by clear-air scattering. It was
necessary to study the clear-air backscatter first by eliminating any
bright targets other than the clear-air scatterers in the lidar's field
of view.
This was done by steering the beam shown in Figure 23 off to
the left of the plume. The receiver was then trained on a segment of
the beam at about the range of the plume. The target spectrum of this
return is shown in Figure 25(a) on a linear scale. The main target
77
-------
(a) MAIN-TARGET COMPONENT
FROM CLEAR-AIR RETURN.
LARGE SPIKE IS ZERO-
FREQUENCY MARKER.
NOTICE CHANGE IN COMPONENT
AMPLITUDE AS DISPLAY SWEPT
FROM -fR TO +fR
(b) WHEN THE LASER BEAM IS BLOCKED,
THE NOISE LEVEL IS SEEN TO DROP.
c) CLEAR-AIR RETURN SIGNAL SHOWN ON
LOGARITHMIC SCALE APPEARS 10 dB
ABOVE NOISE. LARGE SPIKE IS ZERO-
FREQUENCY MARKER.
(d) CLEAR-AIR RETURN ON EXPANDED
SCALE, SHOWING THE MAIN COM-
PONENT IS SPECTRALLY NARROW.
1.0 kHz
SA-1979-23
FIGURE 25 CW LIDAR CLEAR-AIR RETURNS
78
-------
spectral component is several times greater than the surrounding noise
spectral amplitude. When the beam is blocked, the spectrum shown in
Figure 25(b) results. The noise level is lower by the measure of white
noise generated by PMT detection of the nonzero average intensity of the
reflected beam.
On an expanded frequency scale (1 kHz/div) it can be seen that
the clear-air component is spectrally narrow [Figure 25(d)]. This photo
shows a spectrum scan made with a 30-Hz filter bandwidth. One can con-
clude that use of a similarly narrow filter to monitor the clear-air
sideband center frequency will give a similar SNR better than 6:1,
allowing adequate resolution of the 10-to-50 percent attenuation in the
clear-air return that would be caused by an intervening smoke plume.
In order to get these clear-air returns it was necessary to
narrow the receiver field of view to reduce spatial-integration effects.
As discussed in Section IV-C-3, a CW lidar receiver viewing more than a
fraction of the intensity modulation wavelength will suffer a loss in
SNR. The series of spectrum photos in Figure 26 shows the spatial-
integration effect on the SNR of the PMT output for constant field of
view and increasing frequency—i.e., increasing number of modulation
wavelengths in the field of view. Each photo shows the spectrum of a
single modulation frequency component as detected by the PMT. As this
frequency was increased, the component was recentered on the display.
A similar series would result if the frequency were left constant while
the receiver field of view was progressively widened.
The top photo shows a 25-dB SNR obtained from clear air when
the receiver views a length of beam equal to only about 3 percent of the
modulation wavelength. The following photos show how the SNR drops below
10 dB for a receiver viewing more than a full wavelength. One can reason-
ably conclude that spatial-integration considerations are indeed important
79
-------
EACH PHOTO SHOWS THE SPECTRUM OF
THE FREQUENCY COMPONENT USED TO
MODULATE THE INTENSITY OF THE
LASER BEAM, AS DETECTED BY THE PMT.
AS THIS FREQUENCY IS INCREASED, IT IS
RECENTERED ON THE SPECTRUM DIS-
PLAY. SPATIAL INTEGRATION OF THE
SIGNAL OVER THE CLEAR-AIR VOLUME
BEING VIEWED RESULTS IN LOWER
SIGNAL-TO-NOISE RATIO FOR HIGHER
FREQUENCIES. FREQUENCY INCREASES
FROM PHOTO TO PHOTO, MOVING DOWN
THE SERIES.
SA-1979-24
FIGURE 26 EFFECTS OF SPATIAL INTEGRATION ON S!GNAL-TO-NOISE RATIO
OF CLEAR-AIR RETURN
80
-------
in the design of a CW lidar and in its subsequent performance in
distributed-target situations.
C. Remote Opacity Measurements
Having established that a discrete target, a thin-plume target, and
a clear-air volume of limited linear extent all have similar target
spectra when viewed by a CW lidar, it was possible to proceed to the
proof-of-principle experiments in remote measurement of lidar target
opaci ty.
The lidar beam was again pointed downrange as shown in Figure 21.
An 18-inch glass plate was placed in the beam at the fifth range-node
position from zero so as to minimize its sideband-amplitude contributions
at other range nodes. The plate was oriented so that its specular re-
flection did not return directly to the lidar receiver; instead, the re-
ceiver viewed only the diffuse scatter from the glass and dust particles
on the glass. Since the glass was 90 percent transparent, most of the
beam continued downrange to the ninth range-node position, where a ply-
wood board was placed in the beam. The receiver viewed the back target
through the front glass. With both targets at node ranges, their com-
bined sideband levels were relatively low, as shown in the linear display
of spectral-component amplitude versus beat frequency, Figure 27(a). It
can be seen that the plywood reference target is about 1.26 times stronger
than the return from the front glass reference target.
When a second glass is placed at the seventh range-node position,
it adds a third component between the two reference-target components,
as can be seen in Figure 27(b). The important change from 27(a) to 27(b),
however, is the change in relative amplitude between the front and back
reference targets. Now their ratio has been reduced to 1.02 from 1.26,
meaning that the light from the back reference—passing twice through
81
-------
(a) FRONT GLASS REFERENCE AT FIFTH NODE
RANGE AND BACK REFERENCE TARGET
(plywood) AT NINTH NODE RANGE SHOWN
WITH NO INTERMEDIATE SEMITRANSPARENT
TARGET.
GRANGE
(b)
(c)
(d)
INTERMEDIATE GLASS PLATE PLACED BETWEEN
FIRST TWO TARGETS AT NODE-RANGE SEVEN.
NOTICE THAT BACK-TARGET AMPLITUDE
DROPS BY 20%, MEANING THAT ONE-WAY
ATTENUATION OF GLASS PLATE IS 10%.
ALUMINUM SCREEN REPLACES CENTER GLASS
PLATE. BACK-TARGET LEVEL DROPS 43%
BELOW ORIGINAL, MEANING THAT SCREEN
IS 69% TRANSPARENT. THE SCREEN ITSELF
DOES NOT APPEAR AS A SPIKE BECAUSE
ITS REFLECTIVITY WAS MORE THAN 10
TIMES LOWER THAN THAT OF THE GLASS
AND PLYWOOD.
FINER ALUMINUM SCREEN MESH REPLACES
CENTER TARGET. FURTHER REDUCTION
MEANS ONE-WAY TRANSMISSION OF 58%
FOR THE FINE MESH.
SA-1979-25
FIGURE 27 CHANGES IN RELATIVE TARGET POWER LEVELS CAUSED BY VARIOUS
TARGET OPACITIES
82
-------
the new center target—has been attenuated to 81 percent of its intensity
when only clear air lay between the references.
From this one can surmise that the two-way transmission through the
center glass is 81 percent of the transmission through only the clear-air
path of Figure 27(a). This corresponds to a one-way transmission of
/0.81, or 90 percent, which is remarkably close to the 90 percent trans-
mission of the plate-glass target when measured in the laboratory. Such
a remote measurement of the opacity of a semitransparent target at a
range of approximately 160 m was first made with this CW lidar in April
1973.
Other semitransparent center targets were measured as well. Figure
27(c) shows the spectrum for a 1/4-inch aluminum wire mesh sprayed black
to reduce its reflectivity. Indeed, the return from the center target is
not visible on this scale but the telltale drop in the amplitude of the
return from the back reference target allows calculation of a 69 percent
one-way mesh transmission. Figure 27(d) shows still further reduction
when a 1/16-inch mesh is substituted in the center target position. Its
transmission, calculated from the remote spectrum display, appears to be
58 percent.
More than 50 different measurements were made of the semitransparent
glass and mesh targets, separately and in combination. The results of
tests on these thin targets are given in Figure 28. The chart is labeled
in percent transmission (one-way) and equivalent Ringelmann scale on the
right. The error bars represent the range of opacity values measured for
each type of target, the small-dots represent the actual remotely measured
value, and the large dots represent the value measured in the laboratory.
The accuracy of measurements made with this first CW lidar trans-
mi ssometer was limited by a number of factors; the effects of each factor
can be mitigated or eliminated by suitable refinements in technique and
83
-------
100
80
z
w
<
-------
subsequent engineering development. The most important factor was the
nonlinearity of the so-called "linear" spectrum-analyzer display scale
used when making three-target opacity measurements. The bottommost divi-
sion on this scale can represent anywhere from 0.1 division to 1.5 divi-
sions, depending on the signal gain taken in the spectrum analyzer's
front end. Some of the extreme data points in Figure 28 result from
photos taken before this equipment limitation was discovered, for which
scale-correction factors were not known.
Another limiting factor was the physical nature of the reference-
target setup. High winds during two experiments tended to jostle the
front glass reference plate, steering its specular reflection near the
receiver aperture and causing a corresponding increase in the amplitude
of the front reference target. Interestingly, an operator watching the
continuously refreshed target spectrum display could easily monitor this
changing reference level and write off any measurements as being subject
to error.
A third measurement error resulted as predicted in Figure 11. The
spurious sideband generated at Node 9 by the target at Node 5 added to
the amplitude of the main target component at Node 9 (and vice versa) to
give an inaccurate measurement of the height of one component relative
to the other.
Finally, there was the inaccuracy inherent in making centimeter-
stick measurements of traces in polaroid oscilloscope photos. Earlier
project plans had called for construction of a bank of three narrow-band
intermediate-frequency amplifiers, each tunable to a target frequency,
for continuous electronic measurement of signal amplitudes. Unfortunately,
other equipment considerations required that the beat-frequency spectrum
be scaled down from the 0-to-l-MHz range to the O-to-100-kHz range. Since
85
-------
the breadboard IF amplifier could not be scaled down accordingly, use of
this otherwise promising signalvprocessing technique had to be postponed.
D. Discussion of Experimental Results
The successive milestones in this experimental study can be listed
as follows:
(1) Preliminary analysis of FM-fW lidar.
(2) Design of laboratory-modei lidar.
(3) Detection and ranging to discrete targets,
(4) Detection and ranging to distributed scattering targets
including clear air.
(5) Remote opacity measurement using discrete object and refer-
ence targets.
(6) Proof-of-concept analysis and report.
In reaching each of these milestones, no fundamental obstacle to the use
of the CW lidar was found, Every significant source of error or per-
formance limitation was traced to an engineering technique or component
that is amenable to improvement, given spfficient attention in a subse-
quent engineering development effort. In reviewing the results, it is
unfortunate that remaining time and funds were not sufficient to accom-
plish a seventh milestone—namely, remote opacity measurement using the
detected clear-air volumes as reference targets, Jn order to accomplish
this milestone, an improved, low-noise PMT receiver 4s required, as well
as a more refined signal-spectrum processor based on the three-AM-radio
concept depicted in Figure 7. More realistic smoke pj,umes than the highly
reflective steam plume used here need to be employed i.n testing a refined
engineering model.
These are all straightforward tasks for the fixture, That this first
hardware assembly was not quite up to reaching th$ ideal seventh milestone
-------
is hardly an indictment of the technique, which passed the first six
milestones in a most positive manner. It is encouraging to have learned
that the atmospheric Doppler and FM-to-AM noise-conversion effects feared
by theorists do not degrade system performance in practice. It was sur-
prisingly easy to align and operate the eye-safe CW laser system, as
opposed to its pulsed-laser predecessor at an equivalently early stage
in its development. And it was gratifying-*o be able to make reasonably
accurate remote measurements of target opacity between Ringelmann 0 and 1
even with the rudimentary system implementation possible within the con-
straints of this first-phase effort.
The principal objective of the program was successfully accomplished:
the concept of CW lidar for remote opacity measurement was demonstrated
along with its capacity for detecting backscatter from remote volumes of
clear air. The engineering requirements for design of a subsequent engi-
neering research model have been determined and reported.
It is concluded that the FM-AM-CW lidar technique for detection,
ranging, and transmissometry through atmospheric particulates is feasible,
and that straightforward improvements, primarily in receiver detector and
postdetection signal processing, will allow accurate measurements of smoke
plumes thinner than Ringelmann 1.
87
-------
REFERENCES
1. D. W. Jackson, "Development of a CW Lidar for the Remote Measurement
of Smoke-Plume Opacity," Proposal No. ELD 69-165, Stanford Research
Institute, Menlo Park, California (22 January 1970).
2. W. E. Evans, "Development of Lidar Stack Effluent Opacity Measuring
System," Final Report, SRI Project 6529, Stanford Research Institute,
Menlo Park, California (1967).
3. W. B. Johnson, "Lidar Applications in Air Pollution Research and
Control," J. Air Poll. Contr. Assoc., Vol. 19, pp. 176-180 (1969).
4. W. B. Johnson and E. E. Uthe, "Lidar Study of Stack Plumes," Final
Report, Contract No. PH22-68-33 (NAPCA), SRI Project 7289, Stanford
Research Institute, Menlo Park, California (1969).
5. P. M. Hamilton, "The Use of Lidar in Air Pollution Studies," Intern.
J. Air and Water Poll., Vol. 10, pp. 427-434 (1966).
6. E. W. Barrett and O. Ben-Dov, "Applications of the Lidar to Air
Pollution Measurements," J. Appl. Meteorol., Vol. 6, pp. 500-515
(1967).
7. W. D. Conner and J. R. Hodkinson, "A Study of the Optical Properties
and Visual Effects of Smoke-Stack Plumes," Cooperative Study Project,
Edison Electric Institute and U.S. Public Health Service.
8. J. E. Yocom, "Problems in Judging Plume Opacity—A Simple Device for
Measuring Opacity of Wet Plumes," J. Air Poll, Control Assoc., Vol.
13, No. 1 (January 1963).
9. C. S. Cook, G. W. Bethke, and W. D. Conner, "Remote Measurement of
Smoke Plume Transmittance Using Lidar," Applied Optics, Vol. 11,
No. 8 (August 1972).
10. M. I. Skolnik, Introduction to Radar Systems, Section 3.3 (McGraw-
Hill Book Company, San Francisco, California, 1962).
11. T. E. Honeycutt and W. F. Otto, "FM-CW Radar Range Measurement with
a COg Laser," IEEE J. Quantum Electronics (February 1972).
89
-------
12. D. H. Sliney and B. C. Freasier, "Evaluation of Optical Radiation
Hazards," Appl, Optics, Vol. 12, No. 1 (January 1973).
13. M. Ross, Laser Receivers (John Wiley and Sons, Inc., New York,
New York, 1966).
14. E. T. Ebersol, "Radar Sees Gnats and StuffJ" Microwaves, p. 9
(February 1973).
15. "FM-CW Radar Detects Atmospheric Pressure," Electronic Design, Vol.
21, No. 15, p. 24 (19 July 1973).
16. A. Siegmann, Proc. IEE, pp. 1350-1356 (October 1966).
90
-------
BIBLIOGRAPHIC DATA 1- Report No. Z
SHEET EPA-650/2-73-037
4. Title and Subtitle
FEASIBILITY OF A CW LIDAR TECHNIQUE FOR MEASUREMENT OF PLUME
OPACITY
7. Author(s)
Richard A. Ferguson
9. Performing Organization Name and Address
Stanford Research Institute
333 Ravenswood Avenue
Menlo Park, California .94025
12. Sponsoring Organization Name and Address
EPA, Chemistry and Physics Laboratory
National Environmental Research Center
Research Triangle Park, North Carolina 27711
3. Recipient's Accession No.
5. Report Date
November 1973
6.
S. Performing Organization Kept.
No.
10. Project/T«sk/Worlc Unit No.
1 1. Contract/Grant No.
EPA 68-02-0543
13. Type ol Report & Period
Covered
Final Report
14.
15. Supplementary Notes
16. Abstracts
This report describes the work performed for the Environmental Protection Agency during the initial proof-of-
concept phase of a program to develop an eyesafe CW lidar for remote measurement of the opacity of smoke plumes
from industrial smoke stacks. The analysis, design, construction, and evaluation of a laboratory model CW lidar were
performed under SRI Project 1979 from 30 May 1972 to 30 May 1973 under EPA Contract 68-02-0543 to deter-
mine the limitations and potential of the technique. The proof-of-principle experiments combine what is called an
FM-CW radar technique with an argon laser. The technique involves modulating the intensity of the laser beam at
a frequency that changes rapidly and linearly with time. A portion of the transmitted signal is mixed electronically
with the light reflected from the targets in a device similar to a radio receiver. Each taiget appears at a particular
frequency. By tuning the radar's receiver to these target frequencies, the researchers were able to measure both the
range and the opacity of semi-transparent targets over distances of 100 to 200 meters.
17. Key Words and Document Analysis. 17o. Descriptors
Remote Measurement
Laser
Lidar
Optical Detection
Smoke-Plume Opacity
Measuring Instruments
Air Pollution Monitoring
Smoke Abatement Enforcement
17b. Identifiers/Open-Ended Terms
Air Pollution
Lidar
Clear-Air Scattering
Measuring Instruments
I7c. COSATI Field/Group 14/02, 20/05, 4/01, 13/02
18. Availability Statement
Release unlimited.
19.. Security Class (This
Report)
UNCLASSIFIED
20. Security Class (This
Page
UNCLASSIFIED
21. No. of Pages
102
22. Price
FORM NTIS-SS (1O-70)
USCOMM-DC «OS20-r»7!
91
-------
INSTRUCTIONS FOR COMPLETING FORM NTIS-35 (10-70) (Bibliographic Date Sheet based on COSATI
Guidelines to Format Standards for Scientific and Technical Reports Prepared by or for the Federal Government,
PB-180 600).
1. Report Numb«r. Each report shall carry a unique alphanumeric designation. Select one of the following types: (a) alpha-
numeric designation provided by the sponsoring agency, e.g., FAA-RD-68-09; or, if none has been assigned, (b) alphanu-
meric designation established by the performing organization e.g., FASEB-NS-87; or, if none has been established, (c)
alphanumeric designation derived from contract or grant number, e.g., PH-4J-64-932-4.
2. Leave blank.
3. Recipient's Accession Number. . Reserved for use by each report recipient.
4- Title and SubHtla. Title should indicate clearly and briefly the subject coverage of the report, and be displayed promi-
nently. Sec subtitle, if used, in smaller type or otherwise subordinate it to main title. When a report is prepared in more
than one volume, repeac the primary title, add volume number and include subtitle for the specific volume.
5. Report Dote. Eoch report shall carry a date indicating at least month and year. Indicate the basi* on which it waa selected
(e.g., date of issue, date of approval, date 01 preparation.
6. Performing Organization Code. Leave blank.
7. Aut! V,..). Give name(s) in conventional order (e.g., John R. Doe, or J.Robert Doe). List author's affiliation if'it differs
from the performing organization.
8- Performing Organization Report Number- Insert if performing organization wishes to assign this number.
9- Performing Organization Name and Address. Give name, street, city, state, and ?ip code. List no more than two levels of
an organizational hierarchy. Display the name of the organization exactly as it should appear in Government indexes such
as USGRDR-I.
10. Project/Task/Work Unit Number. Use the project, task and work unit number* under which the report was prepared.
11. Contract/Grant Number. Insert contract or grant number under which report was prepared.
12. Sponsoring Agoncy Nome and Address. Include zip code.
13. Type of Report and Period Covered. Indicate interim, final, etc., and, if applicable, dates covered.
14. Sponsoring Agency Cod*. Leave blank.
IS. Supplementary Notes. Enter information not included elsewhere but useful, such as: Prepared in cooperation with . . .
Translation of ... Presented at conference of ... To be published in ... Supersedes . . . Supplements . . .
16. Abstract. Include a brief (200 words or less) factual summary of the most significant information contained in the report.
If the 'eport contains a significant bibliography or literature survey, mention it here.
17. Key Words and Document Analysis, (a). Descriptors. Select from the Thesaurus of Engineering and Scientific Terms the
proper authorized terms that identify the major concept of the research and are sufficiently specific and precise to be used
as index entries for cataloging.
(b). Identifiers and Open-Ended Terms. Use identifiers for project names, code names, equipment designators, etc. Use
open-ended terms -written in descriptor form for those subjects for which no descriptor exists.
(c). COSATI Field/Group. Field and Group assignments are to be taken from the-1965 COSATI Subject Category List.
Since the majority of documents are multidisciplinary in nature, the primary Field/Group assignment(s) will be the specific
discipline, area of human endeavor, or type of physical object. The application(s) will be cross-referenced with secondary
Field/Group assignments that will follow the primary posting(s).
18. Distribution Statement. Denote releasability to the public or limitation for reasons other than security for example "Re-
lease unlimited". Cite any availability to the public, with address and price.
19 & 20. Security Clossificotion. Do not submit classified reports to the National Technical Information Service.
21. Number of Pages. Insert the total number of pages, including this one and unnumbered pages, but excluding distribution
list, if any.
22. Pric«. Insert the price set by the National Technical Information Service or the Government Printing Office, if known.
FORM NTt*-SB (10-70) USCOMM-OC 4O3ZO-P71
92
-------
ENVIRONMENTAL PROTECTION AGENCY
Technical Publications Branch
Office of Administration
Research Triangle Park. N.C. 27711
POSTAGE AND FEES PAID
ENVIRONMENTAL PROTECTION AGENCY
EPA - 335
OFFICIAL BUSINESS
AN EQUAL OPPORTUNITY EMPLOYER
Return this sheet if you do NOT wish to receive this material I I,
or if change of address is needed f~~l. (Indicate change, including
ZIP code.)
------- |