12.215 Modern Navigation: Thomas Herring

12.
215 Modern Navigation
Thomas Herring
Review of last class
• Atmospheric delays are one the limiting error sources in GPS
• In high precision applications the atmospheric delay are nearly
always estimated:
– At low elevation angles can be problems with mapping
functions
– Spatial inhomogenity of atmospheric delay still unsolved
problem even with gradient estimates.
– Estimated delays are being used for weather forecasting if
latency <2 hrs.
• Material covered:
– Atmospheric structure
– Refractive index
– Methods of incorporating atmospheric effects in GPS
11/29/2004 12.215 Modern Naviation L20 2

Today’s class
• Ionospheric delay effects in GPS
– Look at theoretical development from Maxwell’s
equations
– Refractive index of a low-density plasma such as
the Earth’s ionosphere.
– Most important part of today’s class: Dual frequency
ionospheric delay correction formula using
measurements at two different frequencies
– Examples of ionospheric delay effects

Microwave signal propagation
• Maxwell’s Equations describe the propagation of
electromagnetic waves (e.g. Jackson, Classical
Electrodynamics, Wiley, pp. 848, 1975)
4π 1 ∂D
∇ • D = 4 πρ ∇ × H = J+
c c ∂t
1 ∂B
∇ •B = 0 ∇×E+ =0
c ∂t

Maxwell’s equations
• In Maxwell’s equations:
– E = Electric field; ρ=charge density; J=current
density
– D = Electric displacement D=E+4πP where P is
electric polarization from dipole moments of
molecules.
– Assuming induced polarization is parallel to E then
we obtain D=εE, where ε is the dielectric constant of
the medium
– B=magnetic flux density (magnetic induction)
– H=magnetic field;B=μH; μ is the magnetic
permeability
Maxwell’s equations
• General solution to equations is difficult because a
propagating field induces currents in conducting
materials which effect the propagating field.
• Simplest solutions are for non-conducting media with
constant permeability and susceptibility and absence
of sources.

Maxwell’s equations in infinite
medium
• With the before mentioned assumptions Maxwell’s
equations become:
1 ∂B
∇ •E = 0 ∇ ×E + =0
c ∂t
με ∂E
∇ • Bcomponent
• Each cartesian = 0B satisfy the
= 0 ∇ × B −of E and
c ∂t
wave equation

Wave equation
• Denoting one component by u we have:
1 ∂ 2
u c
∇ u− 2 2 =0
2
v=
v ∂t με
• The solution to the wave equation is:
ik.x−iωt ω ω
u=e k= = με wave vector
v c
k×E
E = E 0e ik.x−iωt B = με
k

Simplified propagation in ionosphere
• For low density plasma, we have free electrons that
do not interact with each other.
• The equation of motion of one electron in the
presence of a harmonic electric field is given by:
m ⎡⎣ x + γ x + ω02 x ⎤⎦ = −eE(x, t )
• Where m and e are mass and charge of electron and γ

is a damping force. Magnetic forces are neglected.

Simplified model of ionosphere
• The dipole moment contributed by one electron is p=-
ex
• If the electrons can be considered free (ω0=0) then the
dielectric constant becomes (with f0 as fraction of free
electrons):
4 πNf 0e 2
ε(ω ) = ε0 + i
mω (γ 0 − iω )

High frequency limit (GPS case)
• When the EM wave has a high frequency, the
dielectric constant can be written as for NZ electrons
per unit volume:
ω 2p 4 πNZe 2
e(ω ) = 1 − 2 ω 2p = ⇒ plasma frequency
ω m
• For the ionosphere, NZ=104-106 electrons/cm3 and ωp
is 6-60 of MHz
• The wave-number is
k = ω 2 − ω 2p /c

Effects of magnetic field
• The original equations of motion of the electron
neglected the magnetic field. We can include it by
modifying the F=Ma equation to:
e
mx − B 0 × x = −eEe − iωt for B 0 transverse to propagation
c
e
x= E for E = (e1 ± ie 2 ) E
mω (ω ∓ ωB )
e B0
ωB = precession frequency
mc

Effects of magnetic field
• For relatively high frequencies; the previous equations
are valid for the component of the magnetic field
parallel to the magnetic field
• Notice that left and right circular polarizations
propagate differently: birefringent
• Basis for Faraday rotation of plane polarized waves

Refractive indices
• Results so far have shown behavior of single
frequency waves.
• For wave packet (ie., multiple frequencies), different
frequencies will propagate a different velocities:
Dispersive medium
• If the dispersion is small, then the packet maintains its
shape by propagates with a velocity given by dω/dk as
opposed to individual frequencies that propagate with
velocity ω/k

Group and Phase velocity
• The phase and group velocities are
1
v p = c / με v g =
ω
( )
d
με (ω ) + με(ω ) /c
dω c
• If ε is not dependent on ω, then vp=vg
• For the ionosphere, we have ε<1 and therefore vp>c.
Approximately vp=c+Δv and vg=c-Δv and Δv depends
of ω2

Dual Frequency Ionospheric
correction
• The frequency squared dependence of the phase and
group velocities is the basis of the dual frequency
ionospheric delay correction
R1 = Rc + I1 R2 = Rc + I1 ( f1 / f 2 ) 2
φ1λ1 = Rc − I1 φ 2 λ2 = Rc − I1 ( f1 / f 2 ) 2
• Rc is the ionospheric-corrected range and I1 is

ionospheric delay at the L1 frequency

Linear combinations
• From the previous equations, we have for range, two
observations (R1 and R2) and two unknowns Rc and I1
I1 = (R1 − R2 ) /(1 − ( f1 / f 2 ) 2 )
( f1 / f 2 ) 2 R1 − R2
Rc = ( f1 / f 2 ) 2 ≈ 1.647
( f1 / f 2 ) 2 − 1
• Notice that the closer the frequencies, the larger the

factor is in the denominator of the Rc equation. For
GPS frequencies, Rc=2.546R1-1.546R2

Approximations
• If you derive the dual-frequency expressions there are
lots of approximations that could effect results for
different (lower) frequencies
– Series expansions of square root of ε (f4
dependence)
– Neglect of magnetic field (f3). Largest error for GPS
could reach several centimeters in extreme cases.
– Effects of difference paths traveled by f1 and f2.
Depends on structure of plasma, probably f4
dependence.

Magnitudes
• The factors 2.546 and 1.546 which multiple the L1 and
L2 range measurements, mean that the noise in the
ionospheric free linear combination is large than for L1
and L2 separately.
• If the range noise at L1 and L2 is the same, then the
Rc range noise is 3-times larger.
• For GPS receivers separated by small distances, the
differential position estimates may be worse when
dual frequency processing is done.
• As a rough rule of thumb; the ionospheric delay is 1-
10 parts per million (ie. 1-10 mm over 1 km)

Variations in ionosphere
• 11-year Solar cycle
400
350 Approximate 11 year cycle Sun Spot Number

Smoothed + 11yrs
300
Sun Spot Number
250
200
150
100
50
0
1980 1985 1990 1995 2000 2005 2010
Year

Example of JPL in California
1.0
0.5
Ionospheric Phase delay (m)
0.0
-0.5
-1.0
-1.5
-2.0
-2.5
-3.0
-8 -6 -4 -2 0 2 4 6
PST (hrs)

PRN03 seen across Southern
California
0.6
0.4
Ionospheric Phase delay (m)
0.2
-0.2 CAT1 LBCH

CHIL PVER
HOLC USC1
JPLM
-0.4
-0.6
-2 -1 0 1 2 3 4
PST (hrs)

Effects on position (New York)
500
Kinematic 100 km baseline
250
(mm)
-250 L1 North L2 North LC North
RMS 50 mm L1; 81 mm L2; 10 mm LC (>5 satellites)

-500
0.0 0.5 1.0 1.5 2.0 2.5
500
250
(mm)
-250 L1 East L2 East LC East
RMS 42 mm L1; 68 mm L2; 10 mm LC (>5 satellites)

-500
0.0 0.5 1.0 1.5 2.0 2.5
Time (hrs)
Equatorial Electrojet (South America)
o
0 Site at -18 Latitude (South America)
-2
Ionospheric L1 delay (m)
North Looking
-4
-6
-8
-10
South Looking
-12
0 1 2 3 4 5 6 7 8
Hours

Summary
• Effects of ionospheric delay are large on GPS (10’s of
meters in point positioning); 1-10ppm for differential
positioning
• Largely eliminated with a dual frequency correction
(most important thing to remember from this class) at
the expense of additional noise (and multipath)
• Residual errors due to neglected terms are small but
can reach a few centimeters when ionospheric delay
is large.

12.215 Modern Navigation: Thomas Herring

Uploaded by

Copyright:

Available Formats

12.215 Modern Navigation: Thomas Herring

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

12.215 Modern Navigation: Thomas Herring

Uploaded by

Copyright:

Available Formats

12.

215 Modern Navigation

11/29/2004 12.215 Modern Naviation L20 2

11/29/2004 12.215 Modern Naviation L20 3

11/29/2004 12.215 Modern Naviation L20 4

11/29/2004 12.215 Modern Naviation L20 6

11/29/2004 12.215 Modern Naviation L20 7

11/29/2004 12.215 Modern Naviation L20 8

• Where m and e are mass and charge of electron and γ

11/29/2004 12.215 Modern Naviation L20 9

11/29/2004 12.215 Modern Naviation L20 10

11/29/2004 12.215 Modern Naviation L20 11

11/29/2004 12.215 Modern Naviation L20 12

11/29/2004 12.215 Modern Naviation L20 13

11/29/2004 12.215 Modern Naviation L20 14

11/29/2004 12.215 Modern Naviation L20 15

• Rc is the ionospheric-corrected range and I1 is

11/29/2004 12.215 Modern Naviation L20 16

• Notice that the closer the frequencies, the larger the

11/29/2004 12.215 Modern Naviation L20 17

11/29/2004 12.215 Modern Naviation L20 18

11/29/2004 12.215 Modern Naviation L20 19

350 Approximate 11 year cycle Sun Spot Number

11/29/2004 12.215 Modern Naviation L20 20

11/29/2004 12.215 Modern Naviation L20 21

-0.2 CAT1 LBCH

11/29/2004 12.215 Modern Naviation L20 22

-250 L1 North L2 North LC North

RMS 50 mm L1; 81 mm L2; 10 mm LC (>5 satellites)

-250 L1 East L2 East LC East

RMS 42 mm L1; 68 mm L2; 10 mm LC (>5 satellites)

11/29/2004 12.215 Modern Naviation L20 24

11/29/2004 12.215 Modern Naviation L20 25

You might also like