EITN90 Radar and Remote Sensing

Lecture 3: Propagation Effects and

Mechanisms
Daniel Sjöberg

Department of Electrical and Information Technology

Spring 2018
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

2 / 54
 Learning outcomes of this lecture
In this lecture we will
I Get an overview of propagation phenomena.
I See how they can be quantified using the propagation factor.
I Learn about the basic structure of the atmosphere and how it
affects electromagnetic waves.
I See the basics of diffraction phenomena and multipath
propagation.

3 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

4 / 54
 Propagation factor
Propagation effects in addition to free space are modeled by the
complex voltage propagation factor Fv
E00 = Fv E0 = (F ejφF )E0
where E00 is the one-way received electric field strength, and E0 is
the corresponding field strength when only free space effects are
considered. When N effects are considered, we have
Fv = Fv1 Fv2 · · · FvN = F1 F2 · · · FN · exp[j (φ1 + φ2 + · · · + φN )]
| {z } | {z }
= F , amplitude = φF , phase

In terms of power we need to consider

|Fv |2 = F 2 = F12 F22 · · · FN2
and the two-way propagation in the radar range equation implies
Pt G2 λ2 σ 4
Pr = F
(4π)3 R4
5 / 54
 Line-of-sight and shadow regions

A target is in the line-of-sight (LOS) region if a straight line can be

drawn to it from the transmitter without passing an obstacle.

There can still be interaction with targets in the shadow region,

due to refraction and diffraction effects. However, this is usually
significantly weaker than LOS.

6 / 54
 The atmosphere
The layered structure of the atmosphere can significantly affect
the propagation of electromagnetics waves.

I Troposphere: 4/5 of atmosphere mass, most weather

processes (and water vapor) occurs here.
I Stratosphere: little water, little weather, increasing
temperature.
I Mesosphere: decreasing temperature, strong winds.
I Thermosphere: high temperature region.
I Ionosphere: atoms and molecules ionized by radiation. 7 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

8 / 54
 Attenuation by scattering and absorption
When a wave interacts with a cloud of particles, each particle can
both scatter and absorb the wave, leading to attenuation.

Dense collections of particles, large compared to wavelength,

attenuate more. The attenuation is characterized by the (two-way)
attenuation coefficient α:
αL
F 2 = 10αL/2 , F 2 [dB] = , [α] = m−1
2
The two-way attenuation in dB is F 4 [dB] = αL.
9 / 54
 Heterogeneous atmosphere

For long range propagation, we need to consider combinations of

several propagation regions.

F 2 = 10α1 L1 /2 10α2 L2 /2 · · · 10αN LN /2

The attenuation can be due to atmospheric molecules as well as
rain, fog, dust, and can vary strongly with time.

10 / 54
Typical attenuation coefficients

11 / 54
 Atmospheric gases and water vapor

Peaks correspond to molecular resonances (rotational or

vibrational), and may help isolating short-range systems.
Long-range systems typically operate in regions of low attenuation.
12 / 54
 Rain
Attenuation due to rain drops depend strongly on rain rate and
frequency.

EV

EH

The small difference between polarizations is due to the flattening

of falling rain drops, becoming larger for horizontal polarization.
13 / 54
 Fog
Using the formula (α in dB/km, M in g/m3 , f in GHz, T in ◦ C)
 
11.152
α = M −1.347 + 0.66 f + − 0.022 T , f > 5 GHz
f
the fog attenuation α as function of water vapor concentration M
is as below. Data for different kinds of fog are found in Table 4-3.

g/m3 14 / 54
Fog parameters

15 / 54
 Snow and hail
Frozen water in crystalline particles changes electromagnetic
properties. Larger effects than for rain at higher frequencies.

Similar characterizations apply to attenuation due to dust and

smoke. 16 / 54
Snow parameters

17 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

18 / 54
 Standard atmosphere
There are several models of “standard” atmosphere, for instance
US Standard Atmosphere and International Standard Atmosphere.
The models give a baseline for predictions. Typically, the refractive
index decreases as height increases, due to the thinning of the
atmosphere.

When the real scenario leads to smaller or larger refraction than

the standard atmosphere, it is referred to as anomalous refraction,
further divided into subrefraction, superrefraction, and trapping or
ducting.
19 / 54
 Snell’s law

Snell’s law of refraction for a planar, layered structure, states that

n0 cos α0 = n1 cos α1 = n2 cos α2 = · · · = ni cos αi

If the refractive index is decreasing, n0 > n1 > n2 · · · , then cos αi

must be increasing, that is, αi must be decreasing. Hence, the
propagation direction of the wave is bent towards the earth
(standard refraction).
20 / 54
 Spherical form of Snell’s law

For a radial structure, Snell’s law must be weighted with the radius,

n0 r0 cos α0 = n1 r1 cos α1 = n2 r2 cos α2 = · · · = ni ri cos αi

This is a small correction, since the earth radius a = 6 371 km is

much larger than the thickness of the atmosphere (around 100 km).
21 / 54
Angle and range estimation errors due to refraction

Elevation angle measured by radar appears larger than true

elevation.
Range measured by radar appears longer than true range.

22 / 54
 Angle estimation errors

Unfortunately, the book does not state the position of the target
with respect to these curves. Largest error at high altitude and
small grazing angle.
23 / 54
 Range estimation errors

Unfortunately, the book does not state the position of the target
with respect to these curves. Largest error at high altitude and
small grazing angle. 24 / 54
 Effective earth model

The geometric radar horizon (assuming a spherical earth)

p
Rh = 2aht
The extended horizon Rh0 due to refraction is given by the same
formula if the effective earth radius ae is introduced (calculated
from dn −8 m−1 )
dh = −3.9 · 10

4
a → ae ≈ a
3
25 / 54
 Anomalous refraction
Propagation conditions differing from the standard model
( dn −8 m−1 )
dh ≈ −4 · 10
I Subrefraction: dn > 0, rays bend upward.
dh
I Superrefraction: dn more negative than standard atmosphere,
dh
rays bend more strongly downwards.
I Ducting and trapping: dn < −16 · 10−8 m−1 , rays may be
dh
trapped in regions 10–20 m (sometimes up to 200 m) in
height. This significantly extends the horizon.

The ducting phenomenon is often caused by temperature and

humidity effects, and can vary with time. Radar wave propagation
can be very different at different times of day. 26 / 54
 Effects of ducting on elevation coverage

The trapping of rays inside ducts depends on angle, and can

produce significant distortion to the intended free space coverage.
Requires good computer models to predict the behavior.
27 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

28 / 54
 Turbulence

Atmospheric fluctuations,
typically in clear, hot, humid
weather. Mostly a problem at
frequencies above 80 GHz,
fluctuations can be around
1–2 dB in amplitude and 300
microradians in angle of
arrival.
29 / 54
 Ionosphere

I D layer: exists only on daylight hours, bends and absorbs low

frequency (3–7 MHz).
I E layer: similar characteristics to D, but higher altitude and
exists at all hours.
I F1 layer: weaker than F2, blends into F2 at night.
I F2 layer: densest electron density, produced by UV. Bends
waves below 30–50 MHz, strongest by day.
30 / 54
 Ionosphere

The ionosphere is dispersive (frequency dependent) due to the

free-moving electrons (electron density Ne , in [electrons/m3 ])
s  2
fp p
n(f ) = 1− , fp ≈ 9 Ne ≈ 9 MHz
f

With Ne increasing with height, n is decreasing if f > fp . Using

reflections in the ionosphere, a radar can see over the horizon.
Requires low enough frequency, and low enough angle.

31 / 54
 Diffraction
Even though an obstacle is blocking the path, some power can be
diffracted into the shadow zone.

With a cylindrical object (long edge), the waves inside the shadow
region are typically cylindrical waves (power decay as 1/R):
F 2 (θ)
F2 =
kR
where k = 2π/λ is the wave number.
32 / 54
 Knife-edge and rounded tip, formulas

Exact solutions can be found for simple geometries, depending on

the radius of curvature b relative the wavelength λ, or kb = 2πb/λ:
    
2 1 θ+π θ+π
F (θ, kb = 0) = √ sec + csc b < λ/50
2 2π 2 2
C0 h i p
F 2 (θ, kb) = (kb)1/3 √ exp −τ0 (kb)θ/3 sin(π/3) 1/k b > λ/50
2
For a real dielectric surface, C0 = 0.910719 and
τ0 = 1.8557 exp(π/3) = 5.2881. Strong dependence on curvature!
33 / 54
Knife-edge and rounded tip, graphs

34 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

35 / 54
 Multipath

The electric field is propagated over four different paths:

I Edd : path RTR (direct-direct, or DD).
I Edi : path RTPR (direct-indirect, or DI).
I Eid : path RPTR (indirect-direct, or ID).
I Eii : path RPTPR (indirect-indirect, or II).
All the possibilities need to be added to find the total field.
36 / 54
 Multipath

Using the reflection coefficient Γ , the one-way and two-way

propagation factors are
F 2 = |1 + Γ e−jkδR |2
F 4 = |1 + 2Γ e−jkδR + (Γ e−2jkδR )2 |2
The book uses only cos(kδR), but ends up with correct results, for
2
instance Fmax 4
= 6 dB and Fmax = 12 dB. 37 / 54
 Coherent summation

Depth of oscillation depends on amplitude of reflection coefficient.

38 / 54
 Multipath signal lobing, one-way

Receiver height 10 m, frequency 10 GHz. Note that multipath

propagation can extend maximum range (signal above threshold).

39 / 54
 The reflecting surface

The roughness of the surface (compared to wavelength) needs to

be taken into account. Relative permittivity r , conductivity σ+ ,
2σh
rms roughness σh , rms slope β0 . Spatial correlation is T = tan β0 .

40 / 54
 Multipath reflection coefficient

The multipath reflection coefficient

Γ = Γ0 D(ρs + ρd )

has several different contributions:

I Γ0 : Fresnel reflection coefficient of smooth earth surface.
I D: spherical earth divergence factor.
I ρs : specular roughness modifier.
I ρd : diffuse roughness modifier.
These are described shortly in the following.

41 / 54
 Fresnel coefficients, smooth surface
At a flat interface between two materials, (rc1 , µrc1 ) and
(rc2 , µrc2 ), the reflection coefficient of a plane wave incident from
material 1 at grazing angle γ can be explicitly calculated,
q
rc2 rc2 µrc2
rc1 sin γ − rc1 µrc1 − cos2 γ
Γ0VV = q rc1 , µrc1
rc2 rc2 µrc2 γ γ
rc1 sin γ + rc1 µrc1 − cos2 γ
q
µrc2 rc2 µrc2
µrc1 sin γ − rc1 µrc1 − cos2 γ rc2 , µrc2
Γ0HH = q
µrc2 rc2 µrc2
µrc1 sin γ + rc1 µrc1 − cos2 γ

Most often, the materials are non-magnetic (µrc1 = µrc2 = 1), and
the complex permittivity can be written rc (ω) = r + σ+ /(jω).

As γ → 0, we have Γ0VV,HH → −1, that is, at small grazing angles

there is complete reflection.
42 / 54
 Fresnel coefficients, smooth surface
At a flat interface between two materials, (rc1 , µrc1 ) and
(rc2 , µrc2 ), the reflection coefficient of a plane wave incident from
material 1 at grazing angle γ can be explicitly calculated,
q
rc2 rc2
rc1 sin γ − rc1 − cos2 γ
Γ0VV = q rc1
rc2
sin γ + rc2 − cos2 γ γ γ
rc1 rc1
q
sin γ − rc2
rc1
− cos2 γ rc2
Γ0HH = q
sin γ + rc2
rc1
− cos2 γ µrc1 = µrc2 = 1

Most often, the materials are non-magnetic (µrc1 = µrc2 = 1), and
the complex permittivity can be written rc (ω) = r + σ+ /(jω).

As γ → 0, we have Γ0VV,HH → −1, that is, at small grazing angles

there is complete reflection.
42 / 54
 Fresnel coefficients, example

rc1 = µrc1 = µrc2 = 1, rc2 = 2.5 . 43 / 54

 Critical angle, Brewster angle

I For a wave coming from a denser medium to thinner

(| rc2
rc1
| < 1), there is total reflection (|Γ | = 1) if the grazing
p
angle γ < γc . The p critical angle is γ c = arccos( rc2 /rc1 )
(or θc = arcsin( rc2 /rc1 ) in terms of angle to the normal).
This explains how waves can be reflected in the ionosphere.
I For a vertically polarized wave,
p there is zero reflection at the
Brewster angle γB = arctan( rc1 /rc2 ), whereas the
horizontally polarized wave is significantly reflected.

44 / 54
 Divergence factor

−1/2
2r12 (r − r1 )

D ≈ 1+
ae r(h1 − (r12 /2ae ))
Important at long ranges, beyond the horizon. Otherwise, D ≈ 1.
45 / 54
Surface roughness

σh
σh0 = sin γ
λ
|ρs | = exp[−(4πσh0 )2 ]
p
|ρd,limit | = 0.5 1 − |ρs |

46 / 54
 A word of caution

The end of Section 4.9.3 becomes quite technical, with much new
terminology. Do not dive too deeply into this if you find it
confusing.

47 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

48 / 54
 Wave propagation in a material
An electromagnetic wave propagating in a material with the factor
ω √
e−jkc x = e−αx e−jβx , kc = β − jα = nc , nc = rc µrc
c
The complex permittivity c = 0 − j00 has a typical frequency
dependence as below:
 = 0 − j00
+
(0) vibration +
+ rotation
−
conduction + −
− transition
−
relaxation (∞)

103 106 109 1012 1015 f (Hz)

49 / 54
 Skin depth
Electromagnetic waves propagating in a lossy medium attenuate as
e−x/δ , where δ = 1/α is the skin depth. For good conductors
s s
2 1
δ= =
ωµσ+ πf µσ+

50 / 54
 Propagation delay
In addition to attenuation, material introduce delay due to reduced
phase velocity, vp = c/ Re(nc ). This is a source of errors in range
estimation.

51 / 54
 Outline

1 Propagation basics

2 Atmospheric attenuation and absorption

3 Atmospheric refraction

4 Turbulence, ionosphere, diffraction

5 Multipath

6 Penetration in materials

7 Conclusions

52 / 54
 Conclusions

Propagation in addition to free space is related to

I Atmospheric attenuation due to scattering and absorption in
molecules and particles; rain, fog, snow, dust etc.
I Refraction due to the layered structure of the atmosphere,
gradient of refractive index.
I Turbulence, dispersive effects in ionosphere, diffraction at
edges.
I Multipath propagation and interference can significantly
modify the received power.

53 / 54
 Some left-out topics

See the book for suggestions on further reading on the following

topics:
I Atmospheric emission
I Surface wave propagation
I Ground-penetrating radar
I Atmospheric turbulence sensing
I Trans-ionospheric propagation

54 / 54