Ionosferic Radio Propagation - Article
Ionosferic Radio Propagation - Article
Ionosferic Radio Propagation - Article
CHAPTER 10
km
700
600
400
NIGHT
AM BAND
300
.25
2
DAY
6 7 8
9 10
20 MHz
Figure 10-1. Typical midlatitude day and nighttime ionograms, recorded by a C-4 ionosonde at Boulder, Colorado. The daytime ionogram shows reflections
from E, Es, F1 and F2 layers; the nighttime ionogram those from Es and F2 layers.
10-2
= 0.009
(10.4)
(10.5)
2X(1 - X)
(10.6)
where
with
(10.2)
rm f
(10.7)
and
e, Eo,
and m are natural constants, Ne is the electron density,
and f is the wave frequency. Below the ionosphere, Ne = 0,
and u = 1. Within the ionosphere, Nc > 0, and u < 1.
At a level where X = 1,
(10.3)
(10.3)
the refractive index u becomes zero. The wave cannot propagate any farther and is reflected. The quantity fN, which
relates the electron density to the frequency being reflected,
is called the plasma frequency. Inserting the natural constants into Equation (10.3) permits us to deduce the useful
relation between electrondensity and plasma frequency (which
is identical to the probing frequency being reflected)
eB
(10.8)
10-3
CHAPTER 10
plasma two and only two "characteristic" waves can propagate. These two characteristic waves are called the ordinary
or o-component and the extraordinary or x-component seen
in the ionogram shown in Figure 10-1. A radio wave with
arbitrary (often linear) polarization will split in the ionospheric medium into two characteristic or o-and x-components, which in general propagate independently.
The reflection condition u = 0 gives two solutions for
X; for the + sign (o-component)
X = 1
(10.10)
(10.16)
as in the no-field case, Equation (10.3); for the - sign (xcomponent)
X = I -
Y.
then
(10.11)
At the reflection level for the O-component the plasma frequency equals the probing frequency, fN = f. The x-component is reflected at a lower level that depends on the local
magnetic field strength. It can be shown that the critical
frequencies fo and fx, for fo > fH, are related by
h' > h.
(10.12)
(10.14)
f)
10-4
(10.17)
dz,
(10.18)
f] dz.
(10.19)
2I
Detectorand
I-f
Ii
Receiverm/)er
"
/T-Ampofer
V4?eo-Amrp/feer
,/ fi'
_CT
//morkersro
Figure
10-2. fo.Schematic
presentation
of the
components of
Ionospheric
major
Pulse
Sounder.
the
frequency
oscillator
The transmit
pulse
is amplified
inansynchronism
with
the transmission
for reception
ausedusing either a
is
matched with
bandwidth
in width
to of transmitted
width
the pulse
(overall
bandwitho
B=1/P,
P the
pulse).
the major
After
componeng
/Digital
Ionosonde
Whileand
verdetection
amplification,
thethe
video
signalof
modulates
the
amplified
incal
synchronismng
withtheofHybrid
transmissiontter
receiverng
theirm
in one orand
several
power stages
and transmitted,
using
a
slowly
in the direction
the X-axis inand
the focal plane
of
suitable wide-band antenna with a vertical radiation pattern.
an imaging optic results in an ionogram recording such as
The same antenna can be used for reception using either a
that shown in Figure 10-1. Since sounders based on the
polarization
or mode identification
[Biblup)
Reinisch,
The
are mixed
signals
down
to and
intermediate
frequencyreceived
in
an
and
IF(or
amplifier,
amplified
that
especially problematic
transistorized
for receivers. More
become commonplace.
This
therefore
less
costly
receiver
still
at many
ionospheric
observatories,
especially
theoperated
well-known
C3 and
C4
which
ionosondes,
were
devel-
permits
the use
antennas
phased
inof smaller
arrays and
for
angle-of-arrival
measurements
polarized
asprotectas
antennas
for
Transmit/Receive
or T/R during
switch, and
which
receiver
input
fromproblematic
overloading
transmission
ofthe
the
pulse,
especially
for transistorized
receivers.
More
re-
Deflecting the
for
a Y-axis
sawtooth
transmitter
with
sounders, in contrast
re-to the digital sounders
developed
in
oped
by NBS have
and which
distributed
were
as
the
respective
antennas
ollocated
generation,
synchr
reworldwide
eptioniz
of
cessing,
more
become
known
asation
analog
sounders,they
in contrast
to recently
the digital
sounders
developed
in
voltage
haveand reception
several places
relatively
easy, a much
sounders
moreanalog
are-
10-5
10-5
CHAPTER 10
manding task arose when investigators attempted to sound
the ionosphere over paths of varying distances to determine
the mode structure and the propagation conditions directly.
If the transmitted signal is to be received within the
receiver bandwidth, the systems must be started at a precise
time, and must have perfectly aligned frequency scans. This
was achieved using linear frequency scans and synchronous
motor drives, which derived their A/C voltage from crystal
oscillators IBibI, 1963]. A large step forward was the development of frequency stepping sounders such as the Granger Path Sounder IGowell and Whidden, 1968] which combined digital and analog techniques. Digital techniques
generated ionograms by stepping synthesizer/transmitter and
receiver through the desired frequency range, providing selectable frequency spacing (for example, 25, 50, or 100
kHz, linear or linear over octave bands). The frequency
synthesis itself and the data processing/recording however,
used the standard analog techniques. All digital and hybrid
pulse sounders currently available use these frequency stepping techniques.
10.1.1.5 Digital lonosonde. The rapid development of
integrated circuits, microprocessors and especially ReadOnly-Memories, and of inexpensive storage of large capacity, has led to the development of digital ionospheric
sounders. These systems have some analog components,
but use digital techniques for frequency synthesis, receiver
tuning, signal processing, recording, and displaying of the
ionograms. However, to the modern sounder, the digital
control of all sounder functions, the ability to digitally control the antenna configuration, and above all, the immense
power of digital real time processing of the data prior to
recording on magnetic tape or printing with digital printers
are of special importance.
A digital amplitude ionogram, recorded by a Digisonde
128 PS at the AFGL Goose Bay Ionospheric Observatory
is shown in Figure 10-3. This system developed at the University of Lowell [Bibl and Reinisch, 1978a,b] uses phase
coding, spectral integration, polarized receive antennas for
o/x component identification, and fixed angle beam steering
of the receive antenna array for coarse angle of arrival measurements to provide a rather complete description of the
properties and origin of the reflected echoes. Using a standard set of 128 range bins for each frequency, the sounder
integrates the sampled receiver output signals for a selectable number of integrations, improving the signal-to-noise
ratio and providing the samples for spectral analysis. Since
for each frequency-range-bin or FRB only one return is
recorded, a search algorithm determines from the set of
separate signals (o, x, several antenna directions, Doppler
lines) the signal with the largest amplitude and retains amplitude and STATUS, that is, special signal characteristics.
Using a special font [Patenaude et al., 1973], the resulting
digital amplitudes are printed out providing the analog presentation essential for the recognition of the detailed struc10-6
600
O-SIGNALS
VERTICAL
400
300
200300
600
500-
10
11
[MHz]
1
Goose Bay Ionospheric Observatory 16 June 1980 1720
Figure 10-3. Digital daytime amplitude ionogram recorded by a Digisonde 128PS at the AFGL
AST. Coarse angle of arrival and polarization information is used to separate the vertical ordinary trace shown in the upper part of the figure.
300
200
CHAPTER 10
NEGATIVE
-
POSITIVE
STATUS
DOPPLER
IONOGRAMS
AMPLITUDE AND
STATUS
ONOGRAMS
AMPLITUDE
DOPPLER
......
MHZ
Figure 10-4. Amplitude/status ionogram taken by the AFGL airborne ionospheric observatory with a Digisonde 128PS at Thule, Greenland 9 December
2231 (UT). The lower right panel shows the amplitude ionogram after removal of radio noise. The Doppler ionogram shown in the upper
right panel is produced by replacing each amplitude in the ionogram below with a number representing the measured Doppler shift. The
separation into positive and negative Doppler traces (approaching and receding reflection regions) is shown in the two panels on the left.
10-8
KM
600
600
500
500-
-400
00
100
MHZ
MHZ
-13.0
Frequency Characteiistics
13.0
12.0-
11.0
11.0
10
10.0
9.0
9.0
70
8.0
60
....
4.0 -
12.011.0
10.09.09
M13.0-
I
a
"i
'H tir~
,..
ate
Frn4.0
t... fifi"
a
-10.0
I H 1.
aill II~Ia
MH
....- .....
4.0
'rEiit
ic
Ant~~~~~~~~~~~~~~~~~iilllsilii
lill:-iiil HI !IIli
A"Iglflij
,;;i;?i~i~li,)i
!-i" -tHHE13
a
9'I~~~~~
'.0~
-1~~~~'1~
10.
9.0 -j
3.-8.0
6.0-
3.0
9
ji
1.0
7.0
l,
.ta
3.3.
6.0
5.0
4.0
5..
main
linear waveform
theFM/CW
of
interruption
12
12AST
instanthe very narrow
system is transmission.
the
advantageof
120The
419
Figure
10-5.
Ionospheric
131.0
7.0fl
i ..i
6.0..
:4.
il[ilIi
:!l:liiii'?$iiiiiil$'''
~-~.........
iiiii:!sbriiii$iiiiii'iiiLi~ii?
13.0
E
1.0 ~
12.0
.
~11.0
~
~
7.0
6.0- : ~:['
13.011.0,......
10D-
2.aIai
7.01.0 :-
-12.0
9.0
characteristics
spanning
three
days
produced
of 1Hz.
Although the
121digital integration and spectral analysis
digital
ionograms
recorded
at Goose
Bay
28-30
83The
U0AR
1980.
April
integrated
Figure 10-5. Ionosphericcharacteristicsspanningthree days producedfromdigitalionogramsrecordedat Goose Bay 28-30 April 1980.The integrated
E layers.
layers.
changesof
dynamicchanges
sporadic E
the dynamic
and sporadic
solar and
showsthe
the solar
of the
characteristicshows
appearance of
heightcharacteristic
the appearance
and the
F layer
layer and
the F
of the
height of
minimum height
of the
the minimum
height
evidenceof
as evidence
well as
as well
of these
these layers
variabilityof
events.
diurnal variability
the diurnal
auroral events.
some auroral
showthe
of some
characteristicsshow
E frequency
frequencycharacteristics
and E
layers as
The FF and
The
system
ceivers with
monostaticsystem
signal, aa monostatic
directsignal,
unwanteddirect
the unwanted
with the
ceivers
and
quasi-random
and a quasi-random
switch
T/R switch
using aa T/R
developed, using
was developed,
was
interruption of the linear waveform transmission. The main
advantage of the FM/CW system is the very narrow instantaneous bandwidth of the transmitted signal, allowing a
similarly narrow receiver bandwidth (nominal 100 Hz at
sweeprates
reduced
is further
furtherreduced
bandwidthis
This bandwidth
kHz/s). This
20 kHz/s).
of 20
sweeprates of
order
the order
of the
ofR
of 1 Hz.
Although the digital integration and spectral analysis
used in the modern digital pulsed ionosondes decreases the
effective bandwidth of a pulse receiver significantly (by a
bandwidth
effective bandwidth
an effective
to an
analysis to
spectrum analysis
by spectrum
by
10-9
CHAPTER 10
-I
inreceiver
Maine.
I--
-J
site
have FM/CW
been obtained
with transmit
as low
as I1Wfor
(CW).
The
system
definitely
is power
a good
solution
the
MacLean
[1969].
Since groundbased
ionosondes obtain ionospheric echoes
measuring geophysical
Re
rcie
site in Maine.
fB
fT
fN fxS 3fB
4fB
foF2
-...........
....
fxF2
I
200400-
1200i
it1200-
'
_u,,nt_ Si
T.e
PS.B,
'.
Ar
140005
Figure 10-7.
15 2025
45
55
65 70 75
85
95
FREQUENCY (Mc/s)
DAY 319 (15 NOVEMBER 1962) 0731/10 GMT (138 0 E, 30 0 S)
SATELLITE HEIGHT I011 Km
10.5 11.5
An Alouette I topside ionogram illustrating Z-, 0- and X-wave traces, cutoffs, resonance spikes, and earth echoes.
distributed over the whole globe, has produced a rather uniform analyzed data base which is archived at the World Data
Centers for Solar Terrestrial Research located at Boulder,
Colorado (WDC A), Izmiran, USSR (WDC B), Tokyo,
Japan (WDC C1) and Slough, UK (WDC C2). With some
exceptions, the individual world data centers store data originating in their respective regions. WDC A stores the data
from the western hemisphere and also data from France and
India.
To provide special instructions for the analysis of the
extremely complex ionograms from high latitude stations,
a High-Latitude Supplement to the URSI Handbook on Ionogram Interpretation and Reduction has been published
[UAG-50, 19751.
For special research efforts, it is often essential to go
back to the source data, the ionospheric films of a specific
station(s). For the western hemisphere, these films are stored
10-11
CHAPTER 10
at the World Data Center A for Solar Terrestrial Physics,
NOAA/NGSDC, Boulder, Colorado. A Catalogue of Ionosphere Soundings Data [UAG-85, 1982] provides access
to this data base, which spans the period from 1930 through
fmin
toE
foFI foF2fxF2
600
400
today. The longest and still continuing operation of an ionosonde station started at Slough, UK in January 1930.
Continuous operation starting before 1940 is still ongoing
at Canberra, Australia (1937); Heiss Island, USSR (1938);
Huancayo, Peru (1937); Leningrad, USSR (1939); Tomsk,
USSR (1937); and Tromso, Norway (1932).
To provide an overview of some of the more important
ionospheric parameters that can be derived from an ionogram and introduce their geophysical meaning, two ionograms are provided in the form of a sketch (Figure 10-8),
and the parameters are identified. Both ionograms depict
the same ionospheric conditions (taken from Figure 10-1)
with the exception of an Es layer that can suddenly appear,
possibly as the result of a windshear at E layer heights. This
Es layer can obscure parts of the trace from reflections at
higher regions of the ionosphere. A list of parameters and
their identification and interpretation is provided here as a
general reference and not as a guide for ionogram analysis.
300
Parameter
h' F2
200
00
o0
h'F
fbEs
fxEs
h'E
transmitted
-1-17l
- -
14
16 IBMH
600
soo400
selected
300ves
200
h'Es
6 7 8 9 10
15
20 MHz
Meaning/Comments
F2 layer extraordinary wave critical frequency. Can be used to infer foF2 using Equation (10.12)
if foF2 is obscured by interference.
foF1
F1 layer ordinary wave critical frequency. This layer is often smoothly merging with the F2 layer
resulting in the absence of a distinct cusp and in difficulties of determining the exact frequency
(L condition).
foE
fbEs
Es layer blanketing frequency. Returns from higher layers are obscured by the Es layer up to this
frequency. This frequency corresponds closely to the maximum plasma density in the (thin) Eslayer [Reddy and Rao, 1968].
fxEs
foEs
foEs can be inferred, applying Equation (10.12). If fbEs < foEs, the layer is semitransparent. Es
and higher layers are both observable. The determination of foEs and fxEs for all cases is subject
to a complex set of rules beyond the scope of this outline (see URSI Handbook on Ionogram
Interpretation). Modern Sounders, using polarized receive antennas, permit unambiguous foEs
determination.
10-12
b) Virtual heights
h'F
Meaning/Comments
Minimum frequency at which returns are observed on the ionogram. Since radio wave energy is
absorbed in the D region according to an inverse square law (Absorption - l/f2 ), the variation of
fmin is often used as a coarse indicator of the variation of D region ionization. fmin is not an
absolute value (as for example foF2), but depends directly on the transmitted power and the
antenna gain. Comparison between stations, therefore, can be only qualitative.
The minimum virtual height of the ordinary wave F trace taken as a whole. Due to the effects of
underlying ionization and profile shape on the travel time of the pulse, these minimum virtual
heights are only useful as coarse and "relative" height classifiers (high, average, or low layer,
compared to a reference day). True height analysis must be made to give more meaningful height
parameters, such as the height of the layer maximum (hmaxF2).
h'F2
The minimum virtual height of the ordinary wave F2 layer trace during the daytime presence of
the F1 layer. When an Fl layer is absent, the minimum virtual height of the F2 layer is h'F,
defined above.
h'E
h'Es
hpF2
The virtual height of the ordinary wave mode F trace at the frequency 0.834 x foF2. For a
single parabolic layer with no underlying ionization this is equal to the height of the maximum of
the layer, hmax.In practice hpF2 is usually higher than the true height of the layer maximum.
Useful as a rough estimator of hmaxbut strongly affected by a low foF2/foF1 ratio (< 1.3).
MUF(3000)F
A set of "transmission curves" [Davies, 1966 and 1969] developed for a selected propagation link
distance (the URSI standard is 3000 km) permits the determination of the Maximum Usable
Frequency, which the overhead ionosphere will permit to propagate over the selected distance.
The MUF is determined from the estimated transmission curve tangential to the F-trace. For this
ionogram MUF(3000)F would be 17.0 MHz.
10.1.1.10 Ionosonde Network. Even though the routinely operating ionosondes forming the worldwide network
are independent, generally operated as subchains or as individual stations by national or private organizations, their
operation is coordinated by the "Ionospheric Network Advisory Group (INAG)", working under the auspices of Commission G (On the Ionosphere), a Working Group of the
International Union of Radio Science (URSI). INAG publishes the "Ionospheric Station Information Bulletin" at varying intervals. The Bulletin provides a means of exchanging
experiences gained at the various ionospheric stations, discusses in detail difficult ionograms for the benefit of all
participants, and disperses information on new systems, new
techniques, special events (for example, eclipses), relevant
meetings, ahd general network news. URSI's International
Digital Ionosonde Group (IDIG), which provides a forum
for the discussion of standardization proposals, for the exchange of software, and for the general exchange of experiences with these rather new and still maturing systems
10-13
CHAPTER 10
180
150W 120W
90W
60W
30W
30E
60E
90E
120E
150E
180
690
-I ~~ <ad<>t -- a
_-
p30t
SASH
e- -X--- e
[;S4~
x60
90
Figure 10-9. Map of vertical incidence ionospheric sounder stations 1984.
10-14
D = (,
2 m
kTe/4iiN2e)l/
COMPLETESPECTRAFOR VARIOUS
VALUES
OF a
(10.22)
00
a10.003
10
'
-e
: 0
:
L 10-
0
DOPPLERSHIFT(
Figure 10-10. The variation of the overall spectrum for different values
of the ratio
. The ion has been assumed 0 +.
These curves assume that collisions are negligible and that
Te = Ti [Hagfors, 1961].
1/2
(10.23)
where mi is the mass of the dominant positive ion and Tj is
the ion temperature. Considering Te Ti, and the dominant
ion to be O+, Afj 2 x 10- 2 Afe. The echo energy is,
therefore, mainly concentrated in a relatively narrow spectral window rendering the radar investigation feasible with
apparatus of much lower sensitivity than initially envisaged.
Figure
how
10-10
the
spectral
Figure 10-10 shows
how
theshows
spectral shape
shape depends
depends on
on aa
parameter a - 47iiD/Afor the case Te = Ti. For a > 10,
the scattered energy is entirely due to the electronic component, whereas for very small values of a, the electronic
component decreases and the energy appears mainly in the
ionic component with a much smaller bandwidth. The electronic component now appears as a single line, known as
the plasma line, at a Doppler shift approximately equal to
the plasma frequency of the medium.
In the ionosphere, the electrons and the ions are at different temperatures and the spectrum of the ionic component
changes for different values of the ratio Te/Ti at a given
value of a. This is shown in Figure 10-11 for the case
a = 0.1 for O + ions which illustrates the double-humped
is encountered in the ionosphere, the total scattering crosssection (o) may approximately be given by
(10.24)
04
12
16
20
24
DOPPLERSHIFT (Af,)
Figure 10-11. Spectra of the ionic component for the case of
(=4iiD/A) = 0.1 [Evans. 19691. (Reprinted with per
mission from IEEE c 1969.)
10-15
CHAPTER 10
based radar system [Evans, 1969]. The most obvious measurement is the electron density (N) versus altitude (h) profile
made by recording the variation of echo power Ps as a
7000
4000
(10.25)
1550 EST
2500
2000
500
400
300
200
150
10
ELECTRON
DENSITY
N/m3)
Figure 10-12. An electron density profile obtained at Jicamarca that extends to almost one earth radius [Bowles, 1963].
Location
Frequency (MHz)
Jicamarca,
Peru
Arecibo,
Puerto Rico
St. Santin,
France
Millstone Hill,
USA
Sondrestrom,
Greenland
EISCAT*
Transmitter:
Tromso, Norway
Receiver:
Tromso, Norway
Kiruna, Sweden
Sodankyla, Finland
*European Incoherent SCATter facility
10-16
50
430
935
440
1300
1300
224
(monostatic)
933.5
(tristatic)
Antenna
6
Pulsed
2
Pulsed
0.15
Continuous
3
4
Pulsed
5
Pulsed
32 m parabola
71
30 m x 40 m
parabola cylinder
67
Power (MW)
300 m spherical
reflector
20 m x 100 m
reflector
68 m
25 m parabola
2
32 m parabola
32 m parabola
30
47
57
FLOATING
GUARD
GUARD
DRVEN
(DIA = 0.24m)
COLLECTCO
(DIA = 0.165cm)
(DIA = 0.058cm)
S/
probe is replaced with the entire exposure conducting surface of the rocket or satellite. If the exposed conducting
surface of the rocket or satellite is much greater than the
~2Tc
tC
2.3-
SPACECRFAFT
SURFACE
.89
DMSP SSIE
07
DMSP
SSIE
ELECTRON
SENSOR
ION SENSOR
GOLD PLATED
ALUMINUM
COLLECTOR
1.75" DIAMETER
COLLECTOR
art,//~~\
/
//
ok
\ GUARDRING
//
<~~~~\
YMOUNTING BOOM
\
\
XSWEPT
B~APERTURE
GRID
SUPPRESSOR
\\ ELECTRON
///
2.25"
AMPLIF
ER
// /
\\ ELECTRON
TO
I _LL
F I
AMP
GOLD PLATED
TUNGSTEN
TRANSPARENCY0.80
ON
E R
COLL:CTOR
COLD PLATED ALUMINUM.
AAMPL
A
CONDUCTING SURFACES GOLD
PLATED.
OND.
26
41'
15
DIAMETER
Figure 10-15.
10-17
CHAPTER 10
0 +
10-18
02
(10.26)
K2
O2 + e -- O + O('D)
O(1 D) -- O (3P) + hv(6300)
OP)
(10.27)
(10.28)
(1029)
Since K2
K1, and in regions where O+ is the dominant
ion (0+ = N) the 6300A volume emission rate is given
by
dh
IQ2N
(10.30)
where
Field of View
1550
Pass Band
25 A at f 1.4; 5 A at f 8
Resolution
Spectral Response
Picture Storage
Tube Gain
Threshold Sensitivity
20 R at 2 s exposure
I kR at 30 frames/s
Dynamic Range
Flatness of Field
Repetition Rate
Temporary Storage
Permanent Storage
Process Controller
Display Systems
Real-Time Display
Character Generators
Digital Encoding
K3
O+ + e -m O(CP)
5
0( P)-->
O(5S) + hv(7774)
(10.31)
(10.32)
700
600 -
1(7774):30R
\0(7774)
Ne
E 5001(6300} 4R
(10.33)
400
300
20
I0-
I
I0
PHOTONS/cm
- ec
' 0'
0'
'
EL
Figure 10-16. Electron density profile derived from true height analysis
of bottomside ionogram matched to a modified Chapman
function for the topside. Also shown are calculated 6300
and 7774 A volume emission rates and column intensities
in Rayleighs.
10-19
CHAPTER 10
decreasing O2 concentration, the 6300 A volume emission
rate is confined to the bottomside and reaches a maximum
value at 300 km, below hmax F2 (360 km). The 7774 A
volume emission rate is proportional to [Ne]2 and attains its
maximum value at hmax F2. Because of the broad altitude
exatent of the equatorial electron density profile the 7774
A volume emission rate displays a similar broad extent and
falls to 50% of the maximum value at 300 km and 450 km.
Thus the 7774 A airglow measurements provide information
over a broad altitude range, with approximately one half of
the emission produced above the F layer peak.
All sky imaging photometer measurements conducted
within a few degrees of the magnetic equator and near the
Appleton Anomaly region (--18 ML) have established the
two-dimensional horizontal extent of equatorial plasma depletions. These are also the regions of post-sunset equatorial
spread F, VHF radar backscatter plumes, and amplitude and
phase scintillation on transionospheric radio propagation.
Figure 10-17 shows an example of all sky images at 6300 A
and 7774 A near the equatorial edge of the Appleton Anomaly. The bright region over the southern two-thirds of the
image is airglow from the high-density anomaly region. The
North-South aligned dark band is a region of decreased
airglow emission. Comparison with simultaneous in situ
measurements from the Atmosphere Explorer satellite (AEE) shows this airglow depletion to the collocated with a
region of significantly decreased ion density. Having established the relation between airglow emission processes and
F layer density, the all sky images provide a two-dimensional map of these depleted regions. In addition, the dynamics of these regions can be monitored over extended
15 DECEMBER
1979
6300ooA
7774A
MLm
GNORTH
MAG NORTH
ACPOSITIONX
Eplitude
GLON 7W
0o~"aj ,
L,\
1o
. E BIMS
E
(..
434kin)
particles,
'
but are convected from a source region at or equatorward of the dayside cusp. These patches are also subject
to structuring processes that lead to scintillation.
......
98F~iljV
V
l
//
l?
rO',
*IA
/-
SOME ASPECTSOF
1031/10.2
102L~__,_
0035
129
03
0033
-138 -144
~k
oo
0034
36
hancement and structuring within these auroras leds to amand phase scintillation.
During Bz south conditions, large patches of F region
plasma are observed to drift across the polar cap in the antisunward direction. Coordinated satellite measurements show
70-0 14
0005
loi_--4
0031
0029
0027
002 5
147 -144 -139 -31
-122 -11 -100
27
36 17, 20o245-0624 23
2831
-7041 -8 73
2336
2306
2236
UT
MLAT
134
GLAI
3344GLON
(W
2208
L-(WTMV
It is convenient to refer to radio waves having frequencies below 3000 kHz as "long waves". These include Ex-
10-20
NORTHW
0800
0900 UT
Figure 10-18. Auroral images taken at 10 minute intervals from 0620 to 0950 UT on 29 January 1979. The upper row under each hour shows the 6300
A images; the lower row, the 4278 A images.
10-21
CHAPTER 10
Table 10-3. Long wave frequency bands.
Designation
Abbreviation
ELF
VLF
LF
MF
use, except for communications that require wave penetration beneath the surface of the ocean or earth. The VLF/LF
bands are used extensively for navigation and military communication. The standard AM broadcast systems utilize part
of the MF band (535-1606 kHz). Long radio waves are also
used in basic ionospheric research, meteorology and thunderstorm study and tracking, standard frequency and time
distribution, geological studies, and minerals exploration.
Long waves propagate by a number of different modes.
These include propagation over the surface of the earth by
diffraction modes, ELF propagation by transmission-line
type modes, propagation by ionospheric reflection (or earthionosphere waveguide modes) and propagation through the
ionosphere by so-called "whistler" modes. Each type of
mode requires a separate physical description and mathematical formulation.
(10.34)
providing
d Us2a cos '{a/(a+h)},
(10.35)
10-22
Frequency Range
0.003-3
3-30
30-300
300-3000
kHz
kHz
kHz
kHz
Wavelengths
105-102 km
102-101 km
10-1 km
1-0. 1 km
".
CURRENT
ELEMENT
10
1.0
0i0JC
EZ
IGI
2
E
/
-V
I
Figure 10-19. Vector field-components at a point P in a cylindrical coordinate system. The plane XOY represents the surface of
the earth.
H,= M(t')
i2dr 1A0.2
or G
0.5
argG. l=seconds
(10.36)
0.05
(10.37)
where
t' = (t -
d/c),
X is the
vacuum
of
space
is the
Eo
permittivity of free space
is the permittivity
108/f m),
m), Eo
(=3 xx 5108/f
(= 8.854 x 10-12 F/m), and u is the free space permeability (= 4ii x 10-7 H/m). The far-field components are
related by
(=3
(10.39)
Ez = -ZoHo
[Heckscher
and Tichovolsky,1981]
wavelength
free
2.0
1.0
Figure 10-20. Height variations in the amplitudes and phases of 100 kHz
(10.38)
Eo =0,
0.15
0.25
0.10
Zenneck wave, without any radiation field, requires an infinitely long source distribution [Hill and Wait, 1978].
For a finitely conducting earth, Equation (10.39) is still
r---
A-v-T-,
-r
1 -r
,.--
00
80s
60\
surface
in
accordance
with
the
Fresnel
formulas.
Norton
SPHERICAL EARTH
the Zenneck wave is important historically and conceptually, it is generally difficult to excite because of its rather
slow decay with height. In fact, the excitation of a pure
1000
100
DISTANCE,
10000
km
10-23
CHAPTER 10
Table 10-4. Ratio of groundwave radial and vertical electric fields for
SURFACE TYPE
(S/m)
a/o
EP/E,
Perfect Conductor
80 0.00118 144.99
Sea Water
Good Soil
10 2
20 0.0236 144.660
Fresh Water
10- 3
80 0.0713
132.890
143.25
10.2.1.2 Idealized Spherical-Earth Models. A mathematical treatment of the groundwave on a smooth spherical
earth of homogeneous, isotropic material was undertaken
5 0.403 114.760
2 x 10 5
Thick Ice
by Watson [1919] to determine if an atmospheric reflecting
layer (ionosphere) was required to explain the large fields
produced by distant transmitters. Such early theoretical analthe ground. The ratio of the radial and vertical electric fields
yses were handicapped by the poor convergence of the inat the surface is given by
finite series contained in the solutions. That difficulty was
largely overcome by Van der Pol and Bremmer's 11938]
(10.40)
EP/E,
(1 - /p)/p,
"residue series" solution, which has become a basis for
modern numerical analysis of the groundwave.
where
Figure 10-21 shows examples of field-strength vs dis(10.41)
tance curves for waves of selected frequencies propagating
p - e/E, - ia/Weo.
over "good soil," assuming a vertical source on the surface.
The earth-curvature causes the wave amplitude to decrease
Starting just below the surface of the assumed uniform
with distance more rapidly than it would on a flat earth of
earth, Ep and Ho decrease exponentially with depth. The
the same material, but near the source the flat- and roundfields are 1/e of their values on the surface at the "skin
earth models give essentially the same fields.
depth" 8, which may be estimated from
One way to present both amplitude and phase data is in
terms of the complex factor W, which is the ratio of the
(10.42) actual field component to an idealized one calculated as if
the earth were flat and perfectly conducting. The flat earth
10
Poor Soil
= 503/of
10 0.0745
-ii/2.
when arg p
(10.43)
o I
,FK
........
propagation distance is taken to be the same as the curvilinear distance on the sphere. Figure 10-23 shows curves
of W in the complex plane, for 100 kHz propagation over
sea water and good soil.
0.7
t100OOkHz
M 0.6
........
.''..'I
-_
go
-0.3
-0.
o;0,3
-0.1 O0S
0.1
0.2
0.3
0.4
.5
0.6
0.7
0.8
10
o1
o2
11
;1o
FREQUENCY,
H2
10-24
.0
10-25
CHAPTER 10
the guide's height decreases, the ratio of power leakage to
total power flow increases [Burrows, 1978].
ELF attenuation is low enough to support very longrange propagation, and a planar model of the earth is inadequate if the path length exceeds the earth's radius. The
most important effect of earth curvature is the closure of
the guide, which allows the field to return to the source
point after one complete encirclement. Therefore, the total
field is the sum of the field arising from propagation over
the shorter great-circle path from the source, and that arising
from propagation over the longer one. The local effect of
the curvature is small, however, and the wave propagates
in the curved guide with nearly the same parameters as it
would in a planar one.
The geomagnetic field interacts with the charged particles of the ionosphere to produce an anisotropic conducting
medium. However, the electrical mismatch between the atmosphere and the ionosphere is large at ELF, and the transition between them abrupt, so the ionosphere acts much
like a perfect reflector. The effect of the anisotropy is therefore small, and ELF attenuation and phase velocity depend
only slightly on the direction of propagation.
10.2.2.1 ELF Field-Strength Calculations. The expressions most often employed for calculating ELF fields in the
earth-ionosphere waveguide are based upon an idealized
model that assumes the earth and the ionosphereto be sharply
bounded and homogeneous. Experimental measurements of
the waveguide properties at ELF have shown consistently
that they can be represented accurately by formulas based
on such a simple model. Complicated calculations that account for vertical structure of the ionosphere [Field, 1969
and Pappert and Moler, 1978] also confirm that the simple
model is adequate for many purposes.
For the idealized model, the magnitude of the magnetic
field of the signal from a horizontal dipole is approximately
IH
IdA
2Xo
c/v - 0.985
C
cos
,
h1/ho
(10.46)
and
o.143f
(10.47)
(10.47)
(10.48)
where H is an (arbitrary) reference height. The corresponding values for ho and h1 are
2] ( p/a 1/2
J sin p/a
-7T
H -
Ln
2.5 x 10
2iTrf
(10.49)
and
(10.44)
where
2+ L (2.39 x
104)
(10.50
Xg(10.45)
riwpoXc/v
230.0
2-CALCULATED
* MEASURED
E 100
eIn
e MEA
/URE~
-
2.0
I
c/v
XFX,/h,,
h
(10.51)
or 2.0
DAYTIME
I-I
0.5
0.143f
0.2
'5 7 10
20
50
100 200
FREQUENCY,
Hz
Figure 10-24. Calculatedand measuredELF daytime attenuationrates
[Bannister,1982].
hE/ho
+
hrrkonEhE)
(10.52)
10.0
-CALCULATED
E
m 5.0
MEASURED
1.8
1.6_
CALCULATED
DAYTIME
MEASURED
z 2.0
C)
0
NIGHTTIME
0 l
5'
20
FREQUENCY,Hz
Figure 10-25. Calculated and measured daytime ELF phase velocities
[Bannister, 1982].
50
FREQUENCY,Hz
Figure 10-26. Calculated and measured ELF nighttime attenuation rates
[Bannister, 19821.
10-27
CHAPTER 10
10.2.2.3 Anomalous ELF Propagation Occasionally the
nighttime field strengths measured at 40 Hz to 80 Hz have
decreased by 4 dB to 8 dB in the northeastern United States
[Bannister, 1975, 1980]. Those relatively severe nighttime
fades sometimes occurred during the several days following
magnetic storms, when similar (but less pronounced) behavior was found to coincide with phase disturbances on
VLF paths across the northern United States. These shortpath (- 1.6 Mm) field strength reductions might have been
caused by enhanced ionospheric ionization due to precipitating electrons from the radiation belts. However, attempts
to correlate the fades with geomagnetic indices have met
with limited success. Simutaneous measurements taken in
Connecticut and the North Atlantic area during the magnetically quiet period of early March 1977 have indicated
some of these anomalies might have been caused by a moving nocturnal sporadic E layer.
Calculationsby Barr [1977] and Pappert and Moler [1978]
show that nocturnal sporadic E can produce marked maxima
and minima in the propagation characteristics of ELF radio
waves. One physical explanation for the effect is interference between waves reflected from the normal E region and
from the sporadic E region. Pappert [1980] showed that a
sporadic E patch one square megameter in extent could
account for the 6-8 dB fades that have been observed. Similarly, patches 0.5 Mm2 in extent could account for the more
commonly observed 3-4 dB fades.
Many other ionospheric disturbances can cause ELF
10.2.3
Very low and low frequency (VLF/LF) waves are reflected from the lowest regions of the ionosphere (the D
region during daylight and the lower E region at night), and
apart from the sunrise and sunset periods, exhibit propagation characteristics that are very stable in both phase and
amplitude. The LF band (30-300 kHz) is useful for communications to distances ranging from hundreds to several
thousand kilometers, shorter than the almost global ranges
achievable at VLF but much longer than the groundwave
distances normally associated with the MF band.
Beginning in 1911 with Austin, various empirical formulas, deduced from numerous measurements, have been
used to estimate the field strengths of these long waves. For
example, Pierce derived a semiempirical formula to describe
VLF propagation over water during the day. That formula,
which gives the vertical electric field strength at a distance
d from a transmitter having a radiated power P (kw) at
frequency f (kHz), is [Watt, 1967]
210 N/P
=
exp
da
(10.53)
TM -TRANSVERSE
MAGNETIC
ELTM
WAVEPOLARIZATIONS
ONOPHERE
'
'-
"M
EXCITED
Y
EXCITEDBY
VERTICAL
CURRENTS HORIZONTAL
CURRENTS
:.
TE-TRANSVERSE
ELECTRIC
FAMILIES
OF PROPAGATION
MODES
Figure10-27. Characteristicsof transversemagnetic(TM)and transverse
electric(TE)waves [Kosseyet al., 1982].
ionosphere, the TM and TE modes are not entirely independent, but are coupled. In general, when an electron acquires a velocity from the electric field of the wave the
magnetic forces cause it to have a component of motion
perpendicular to the electric field, thus causing polarization
conversion. The interaction is described by the AppletonHartree equation (for example, see Ratcliffe [1959]), and
the polarization conversion has been demonstrated experimentally (as reported by Bracewell, et al., [1951] Lewis,
et al., [1973] among others). The coupling effect is most
pronounced at night when the waves reflect higher in the
ionosphere, where the electron-neutral collision frequency
is smaller than the electron gyrofrequency.
10.2.3.2 Waveguide Propagation Equations. The detailed equations for VLF/LF waveguide propagation are
described, for example, by Wait [1970], Pappert and Bickel
[1970], Galejs [1972], and Field et al., [1976]. To illustrate the key dependencies and to define the commonly
used notations, they are given here for the case when
geomagnetic anisotrophy (that is, polarization conversion) effects can be neglected. That approximation is fairly
accurate for long-range VLF/LF propagation under normal daytime conditions.
Transverse Magnetic (TM) Modes
Conventional ground-based VLF/LF transmitters are
vertical and their fields are composed of a superposition of
TM waveguide modes. Following Field [1982c] the vertical
electric field is given by
e
exp
/ d/a
ir
IL cos 4
exp
sin d/a
exp
e
2id)G(hT)Gj(h)
Vj
d
8.7/
(10.54)
V/m,
10-29
CHAPTER 10
10-4
d/a
sin d/a
exp
(-
iTr)\
120(riILexp
sin ti
4EH=-l2TiXd exp
Sm 12
V-4
AI exp
(10.55)
8.7
-d)
v,"
G,(hr)G,,(hR) V/m.
10-30
11
10-9
'
E-
o
E 10-5-
lo-lo z
z
\
',
-\
m\\
\
E
106o
10
15
\\\
20
25
30
30
35
35
40
40
45,
45
FREQUENCY,kHz
\o~l
70
60
lo
',
,
\
10
kHz,
10
cr = -
et3
S/m
1976].
[Field
al
.,
1976
'%
=2
0:2
10
15
20
25
30
35
45
40
and, above a few kilometers, the TE mode height-gain factors increase sharply to values well in excess of one-hundred.
For elevated antennas, the large TE mode height-gain factor
mitigates the effects of the small excitation factor, and these
modes can be excited about as effectively as TM modes.
At frequencies above about 30 kHz the first TM mode
develops a broad maximum in its height-gain factor in the
40-60 km altitude range. Such "whispering gallery" behavior is not important for ground-based or airborne terminals, but may be significant for very high, balloon-borne,
agation analysis than is mode theory. In intense disturbances, however, higher order modes are much more severely
attenuated than lower, and the mode sum can be used well
into the LF regime. Figure 10-30 shows that the first TE
mode is slightly less attenuated than the first TM mode,
although that result depends on the specific normal daytime
conditions and ground conductivity assumed.
Figure 10-31 illustrates the ground conductivity depen-
10-31
CHAPTER 10
40
'
-TM
35
-- I
al0
\
TM MODE
i0~~~~~g~~
---
-- let TE MODE
\
~~~~\
a2
'
E0
. _\.MODERATE
--
w
5
104
1d73
o2
o
10-1
CONDUCTIVITY,
Siemens/m
1I o' /
,/30
/
30
t-
I< WEAK
MODERATE
INTENSE
SEVERITYOF DISTURBANCE
Figure 10-32. Attenuationratesof the lowest orderTM and TE modes,
for disturbedionosphericconditionsand a perfectlyconductingearth [Fieldet al., 1976].
10-32
INTENSE
DISTURBANCE
I
0I
1ll
DISTANCE,Mm
Figure 10-33. TE/TMsignal strengthsvs distancefor ambientday and
3
intenselydisturbedionospheric
conditions:
20kHz,o = 10S/m,
10atdegree
antenna
inclination
receiver
12.2km
[Field,
1982c].and transmitterand
l/./
2 lo
O"H
,/
-I
IOkHz/
E --TE
DAY
\Id\
0dB
al., 1976].
20
,the
,AMBIENT
-TM
A E
\DISTURBANCE
NORMALD6ATIhEI
TE
80
70
8C0
102897kHz
70
60
60
20
20
15 567 kHz
st00
0
l60
c
21 794 kHz
70
60
28021 kHz
DISTANCE,Mm
Figure 10-34. VLF signal calculations for a mostly seawater path from Hawaii to Sentinel, Arizona:-daytime,----nighttime
[Morfitt, 1977].
10-33
CHAPTER 10
\,
50
< 40
3 x~i~{'~
'E 20
i
k.
9\\
/
o0
-10,
0
.~---
\""<*~,
X-. ~
"o
250':.'
m0-
\
x( o
_propagation
,
1
DISTANCE , Mm
quencies as high as 300 kHz for daytime propagation [Pappert, 1981] and as high as 60 kHz for nighttime propagation
[Morfitt et al., 1982]. Figures 10-34 and 10-35 show representative TM signal strengths computed by NOSC for
nominal daytime and nighttime models of the ionosphere
appropriate for midlatitudes, and for propagationover water.
At VLF, many of the propagation predictions have been
validated by NOSC airborne measurements, such as those
described by Bickel et al. [1970].
10.2.3.5 Other VLF/LF Propagation Prediction
Techniques. In addition to the waveguide mode formulations a number of other mathematical techniques have been
developed for describing the propagation of VLF/LF waves.
These include the zonal harmonic or spherical wave analysis
method [Johler, 1964, 1966]and the wave-hopmethod [Berry,
1964, and Berry and Chrisman, 1965]. The spherical wave
technique has the attractive feature that it can model variations in height in the earth-ionosphere cavity, but it requires
the use of large-scale, very fast, digital computers for its
implementation. For the higher VLF and LF frequencies the
wave-hop method requires the least computing time, but it
has not been formulated in a way that lends itself to modeling
of discontinuities in the earth-ionosphere duct [Jones and
Mowforth, 1982].
In the wave-hop approach the field strength at any point
is the sum of the groundwave (see Section 10.2.1) and a
series of "hops," which represent waves that have been
reflected from the ionosphere and/or the ground. The hops
are numbered according to the number of times they have
been reflected from the ionosphere. Each reflection results
in a reduction in hop amplitude so that usually a relatively
small number of hops are needed to provide good fieldstrength estimates. For VLF propagation at distances less
than about 1000 km and for LF propagation, in general, the
10-34
the geomagnetic field are greatly diminished and the ionosphere can be assumed isotropic. Under such conditions the
mathematical inversion problem becomes somewhat simpler.
Field et al. [1983] have developed an inversion technique,
appropriate for isotropic propagation, which has been used
in conjunction with VLF/LF pulse reflection data to derive conductivity profiles of the severely disturbed polar
ionosphere.
A problem with profiles calculated by inversion is that
of nonuniqueness, which can be caused by either incompleteness of data or the nonlinear dependence of the reflected
signal on the propagation medium. In addition, the profiles
characterize narrow regions of the ionosphere, since the
propagation data contain information about only those altitudes where the ionosphere interacts appreciably with the
reflected wave.
The altitude constraints are even more severe if long
path propagation data are used, rather than steep-incidence
reflection data. Nevertheless, some effort has been devoted
to deducing the structure of the ionosphere from long path
data. As described by Crain [1970], the data in this case
are the attenuation rates and phase velocities of the propagating waveguide modes, and the analysis is a trial-anderror technique effectively to find an ionospheric conductivity profile which provides a waveguide mode or wave
hop structure that agrees with the observed distribution of
radio field strength.
In essence the long-path technique is similar to the steepincidence approach insomuch as ionospheric reflection coefficients are calculated as an intermediate step in obtaining
the mode constants. In order to synthesize the total field as
measured, one has to take care to add in as many modes as
contribute to the field. This can add a great deal of complexity to the application of the technique. Nevertheless,
the technique has been applied with much success to develop
phenomenological models of the lower ionosphere. Although such models may not be consistent in all respects
with those derived from detailed analyses of the aeronomy
of the upper atmosphere, they have found widespread application in long wave propagation prediction codes.
The results obtained by Bickel et al. [1970], Morfitt
[1977] and Ferguson [1980] are especially noteworthy in
that regard. They have performed detailed analyses of a
large volume and a wide variety of VLF/LF propagation
data and have derived analytic models of the lower ionosphere for propagation prediction. Those models are simply
exponential height-profiles of conductivity, which can be
specified by only two parameters, scale height and reference
height. Following Wait and Spies [1964] the conductivity
parameter, wr, depends on the ratio of electron density to
electron-neutral collison frequency, and is taken to be of
the form wr(z) = 2.5 x 105 exp (B(z-H')), where z (km)
is altitude, B is the inverse scale-height (km1) and H' (km)
is a reference height. The value of electron density N(z),
in electrons/cm3 , is calculated as a function of height by the
equation N(z) = 1.43 x 107 exp(B(z-H')-0.15z). The
10-35
CHAPTER 10
Table 10-5. Suggested exponential profiles for use in long wave propagation prediction codes. Frequencies, f, are in kHz [Morfitt et al., 1982].
Seasonal-Diurnal
Propagation
Condition
Summer day
Summer night
Winter day
Winter night
H' (km)
(km )
B
70
87
74
80
0.5
0.0077f + 0.31
0.3
0.035f - 0.025
(10 < f < 35)
Linear change between high and
middle latitude
0.0077f + 0.31
10.2.4 MF Propagation
At night medium frequency skywaves can propagate to
considerable distances with relatively little attenuation, but
during the day the skywaves are severely attenuated in passing through the ionospheric D-region, so that only the
groundwave provides usable signals. Thus, during the daytime MF signals are very stable, while at night they are
much less so owing to the variability of the lower E-region
of the ionosphere, and to interferences between the groundwave and skywaves.
10.2.4.1 MF Groundwave Propagation. The propagation of MF groundwaves can be described using the techniques discussed in Section 10.2.1. Because of the shorter
wavelengths, however, such factors as the earth's atmosphere (and hence, the effective radius of curvature of the
earth), terrain elevation, conductivity changes, and trees and
buildings along the propagation path usually influence MF
90-75
(high latitudes)
(transition latitudes)
<70
(middle latitudes)
INPHERIC
POLARIZATION
/
COUPLINGLOSS
COUPLINGLOSS
set of propa-
POLARIZATION
RUN
o
LOSS
TRASMITTER
RECEIVER
19821.
......... |
,,,,,
| W|
m
14
the Ionosphere
J0
0.7
i'
l.5
U),
Z
:t
4
'~T~i~-~
-es--------s---5
10
11
1213 15
'
20o
HOPLENGTH,kmx100
Figure 10-37. Computed MF ionospheric reflection losses:-east-west
propagation at all latitudes,---north-south propagation at
the magnetic equator [Knight, 1982].
10-37
CHAPTER 10
.I.l_l
90c
--- ABIENT
,I
DAY
-NIT/
60
/
I
4,) 40
V 40 -
i/
I
/
30
F
20
lo
_-
3
lo
O3 lo-2 lo 10 s1 12
14
O15 i6
~
07~108 109
1010
FREQUENCY, Hz
10-38
10-39
CHAPTER 10
of 3000 km between the (hypothetical) transmitter and the
receiver. The experimental data for the numerical factor
M(3000) (in terms of coefficients) come from 13 ionosonde
stations covering the geomagnetic latitude range from 7 S
to 880N.
The stability and predictability of the E layer results in
a 99% probability (highest) of supporting radio propagation
and communication via the E layer. The next highest probability is via the regular F layer. When neither of the above
modes is possible, the Es mode is considered for communication.
For computing the system performance the program allows a selection from 10 antenna patterns. The program
takes into account the ground losses, ionospheric losses,
free space losses, and the excess losses. The program computes the radio communication circuit reliability, service
reliability, and the multipath evaluations.
The ITS-78 model has several limitations. The results
from the model are useful only when the operating frequency
is below the maximum usable frequency. The model assumes that transmission will be by reflection from the ionosphere. For this the transmitter and the receiver must be
on the same side of the ionospheric layer (for example,
ground-to-ground communication. The model does not take
into account the daytime Fl layer which usually develops
between the E and F2 layers. The model does not adequately
account for the electron density above the altitude of hmF2.
Finally, the model does not take into account the dependence
of absorption on the operating frequency in considering the
D layer absorption.
10.3.1.2 The Bent Model. The Bent Model [Bent et al.,
1972; Llewellyn and Bent, 1973] is basically for groundto-satellite communication but can be adapted for groundto-ground or satellite-to-satellite communication. The main
purpose of the Bent model is to determine the total electron
content (TEC) of the ionosphere as accurately as possible
in order to obtain high precision values of the delay and
directional changes of a wave due to refraction. Ground-tosatellite communicationdemands operatingfrequencies which
are higher than the MUF. Thus the mode involves the transmision refraction characteristics of the F2 layer and the
electron density distribution above the height of the F2 peak
must be known.
The model provides (output) the vertical total electron
content above the transmitter, the profile of vertical electron
density with altitude, and the total electron content along
the path between the satellite and the ground. It also provides
the refraction corrections to the elevation angle, the range,
and the range rate.
The input parameters to the model are the date, Universal
Time, locations of the transmitter and receiver (ground and
satellite), rate of change in elevation and altitude of the
satellite, operating frequency, the solar flux (10.7 cm flux),
and sunspot number.
10-40
(10.56
aN N e' 3
h2 -
3;
N:
N=Ne
hi_
LtL__.
__ I
k202
ho.
h.L
-8i, L. 1-----
N N(
tI-
b. \2
-........
fib
2I
j\
Ym
I N Nm( N: I I;
I ~f,,F2
if
f. F2
I
N2
iI
II
N1
NoNm
are needed.
The schematic for the electron density distribution with
altitude
-____~for the IONCAP program is shown in Figure
10-40. The model has 3 parabolic layers, E, Fl, and F2.
2) The altitudes for the peak electron densities are hmE, hmF1,
and hmF2. The half thickness widths for the layers are YmE,
is assumed to be 4. The E layer has fixed altitudeshmE = 110
km and YmE = 20 km. IONCAP improves on the ITS-78
model by incorporating D and Fl layers. The D layer con-
ELECTRON
DENSITY,
N
Figure 10-39. Schematic for exponential and bi-parabolic profiles for the
electron density distributionwith altitudefor the Bent model.
HEIGHT,km
VIRTUAL
hmFz.. ....
T
ymF2
PARABOLIC
NOSE
mFI
predictions.
c
LINEARVALLEY
LINEAR
(or PARABOLIC)
LEDGE
ymF
-
hmE
YE_
70
PARABOLC
NOSE
-EXPONENTIAL
/70m I
TAIL
I
fOE
foFI
FREQUENCY
foF2
Figure 10-40. Schematic for the electron density profile and virtual height
10-41
CHAPTER 10
layers, the electron density is assumed to be linear for the
frequency range fv to fu where fv = xv x foE and f,, = xu
x foE. Typical values for xv and xu are 0.85 and 0.98
respectively. Thus the electron density decreases above the
parabolic nose 0.85 foE (<foE) and continues upwards up
to 0.98 foE (<foE) producing a linear valley in the transition
region. When the x,, = xv = 1 the valley is absent in the
transition region and the curve is a vertical line starting at
the tip of the parabolic nose of the E layer. The F1 layer
forms a linear or parabolic ledge depending on the magnitudes of hmF1,hmF2, foF1, and foF2. In the ITS-78 model,
the Fl layer is assumed to be absent. In the IONCAP model
the numerical coefficients for foE are functions of geographic latitude for both solar maximum and minimum from
the work of Leftin [1976]. The model uses foF1 maps of
Rosich and Jones [1973]. It also takes into account the
retardation below the F2 layer.
For the MUF computations the model uses the corrected
form of Martyn's theorem. As the absorption equations using the secant law do not work for lower frequencies at
altitudes below 90 km, these equations have been modified
in the IONCAP program. The IONCAP provides two programs (1) the ITS-78 short path geometry and (2) the long
path (>10 000 km) geometry. In addition to the ITS-78
model, the path computations now include the Fl mode,
the over-the-MUF mode, D and E region absorption losses,
and sporadic E losses. A correction to frequency dependence
is added for low frequencies reflected from altitudes below
<90 km.
The improvements over ITS-78 can be summarized (see
IONCAP) as follows:
1. The description of the ionosphere is now more complete.
2. The loss equations have been supplemented. This
includes E mode adjustments, sporadicE effects, overthe-MUF losses, and losses for low reflection heights.
3. The ray path geometry calculations have been revised. This was an empirical adjustment to Martyn's
theorem.
4. The loss statistics were revised to include the effects
of the sporadic E layer and of over-the-MUF modes.
5. A separate long path model was developed.
6. The antenna gain package was revised.
All models predict only the quiet ionosphere, which
shows a large systematic dependence on latitude, longitude,
season, time, and sunspot activity.
10.3.1.4 The Bradley Model. The Bradley model contains two modifications to the existing models: (1) the filling
of the valley between the E and F layers (Fl, or F2 if Fl
is absent), by parabolic layers [Bradley and Dudeney, 1973]
and (2) a simple formulation of the prediction of the probability of the high-frequency propagation [Bradley and Bedford, 1976].
The assumption that the electron density distribution in
the E and F layers is parabolic in shape, results in a valley-
10-42
301
(10.58)
(10.59)
NJ(e) =
akWk(0),
E
k =
CHAPTER 10
The programs are in FORTRAN-4 and ALGOL-60 computer codes. These programs are available from the World
Data Center, Boulder, Colorado.
Topside ( 1 )
HST
Intermed. (
--
the physical processes responsible for the observed ionospheric phenomena. The processes responsible for the ionosphere are production, maintenance, and decay of the ionsphere. As the approach deals directly with the physical
u,
IHST'"ting
-/ Xntermed.
(4)
/2[HEF
~B
-E - Valley 5 1 --
I HBR
HABR
t
HME/
_.
_HgOX
/D
HA
NME
log N
and Oran and Young [1977] also considered the minor constituents He and H. For the dissociation and ionization of
the neutral species, the incident solar EUV radiation in the
NMF2
F region
E region
Peak
Height Density
+ 15% 30%
5% 10%
Temperature
+ 30%
+ 10%
Table 10-6. Variations in the physical processes used in the theoretical models
Nisbet
[1971]
Processes
Stubbe
[1970]
120-1250
N2 , 02, 0
O +, NO ', 02 +
O +, NO +0, 02,
Hi
Neutral Winds
Horizontal Winds
Electric Fields
Magnetic Fields
Yes
Yes
Additional Features
10-44
5 Reactions
120-1500
0, He,
H
10 Reactions
N 2 , 02,
200 700
N2, 02, O, He
120-1200
N2 , 02, O, He, H
4 Reactions
24 Reactions
Horizontal Winds
O , NO +, 02,
H+
N + , He + , N2 +
Horizontal Winds
CHAPTER 10
Table 10-7.
Intercomparison
of the empirical-computer
based ionospheric
models.
Ionospheric Models
Parameter
ITS-78
D Region
E Region
foE
hmE
YmE
F1 Region
foF1
IONCAP
Bent
IRI 79
Non-deviative and
deviative absorptions only
Same as ITS-78
+ E Layer
exponential
extension down
to 70 km
Not modeled
Modeled
Modeled by Leftin
et al. [1968]
110 km fixed
20 km parabolic
shape
Same as ITS-78
+ exponential
down to 70 km
Leftin [1976]
coefficients
Not modeled
Modeled
Kouris and
Muggleton [1973a,b]
Not modeled
Modeled
Eyfrig 11955]
Ducharme et
al. [1973]
Bi-parabolic
Modeled
Rawer [1981]
Not modeled
Rosich & Jones
coefficients 11973]
hmF 1
YmF I
hmF1/ymF1 = 4
(fixed)
F2 Region
Bottomside
foF2
hmF2
YmF2
Haydon-Lucas
coefficients [1968]
Shimazaki eq [1955]
+ E layer
retardation
Same as ITS-78
Bent coefficients
Kelso [1964]
F2 Region Topside
Not modeled
Not modeled
Up to 1000 km
Modeled
Rawer [1981]
E-F Transition
Region
Not modeled
Modeled
Not modeled
Maeda [1971]
Electron-Density
Profile
Not computed
Available up
to hmF2
Available up
to 3500 km
Available up
to 1000 km
Electron, Ion
Temperatures
Not modeled
Not modeled
Not modeled
Modeled
Ion Composition
Not modeled
Not modeled
Not modeled
Modeled
Total Electron
Content (TEC)
Not computed
Not computed
Computed
Not computed
MUF
Not modeled
Not modeled
Function of Kp
Not modeled
Not modeled
Not modeled
Input Parameters
required
Sunspot number
Sunspot number
Sunspot number
and 10.7 cm
solar flux
Sunspot number
10-46
(Continued)
Ionospheric Models
Parameter
Noise Parameters
Galactic
Atmospheric
Manmade
Modeled
Same
Modified
Same
Modeled
Not modeled
Not modeled
Not modeled
Not modeled
System Performance
Modeled for
short path
(<3000 km)
Also has a
long path
option >
10,000 km
Not modeled
Not modeled
Antenna Patterns
Modified
ITS-78
package with
17 antenna
options
Not modeled
Not modeled
Sporadic E
Modeled in terms
of occurrence
frequency
Same as ITS-78
Not modeled
Not modeled
Circuit Reliability
Service Probability
Multipath Evaluation
Modeled
Modeled as
ITS-78
Not modeled
Not modeled
MUF 50%
FOT 90%
HPF 10%
10.4.1
Morphology
of the Ionosphere
SOLARMAXIMUM
CHAPTER 10
14
14
SUMMER 1958
12 -
oI
12 -
T= 188
O-
06
12
WINTER 1958
F2
18
24O
T= 159
06 8
12
6z
SUNSET
SUNRISE
O
0
06
12
18
HOUR (LT)
24
150E 0
SUNRISE
II
06
SUNSET
12
18
HOUR (LT)
24
SOLAR MAXIMUM
10
10
SUMMER 1964
8- ~~~~~N
lo "
4
WINTER 1964
I~~8
T=lO
'1OE/
~
06
10
cr
SUNRISE
Figure 10-42.
10-48
II[z
T=O
06
12
18
24
JSUNSET
SUNRISEI
SUNSET,
06
12
18
24
Diurnal variations of the critical frequencies of the E, F1, and F2 layers for solar maximum (1958) and solar minimum (1964) and for
summer (January) and winter (June), at a typical midlatitude station (Canberra). The parameter T is an ionospheric index related to the
level of solar activity.
90
60 E
180
120E
600 W
120W
J60
CD
0-
310090080070060*
a.
60070
300 40050
80
200
500400 30
~~~~~~~~~~~30,
-1
'~~~~~~~-J
:z ,,I
30
Crr
~~~5
\\\\\
SOLAR
60X/-ZENITH
e::~~~~~~~~~~0
ANGLE
700
::
-0.90SUN
"50
7~~~~~~~~~~~~00
80
901
0O
90
04
08
600 E
120*E
12
1800
16
20
1200 W
'"UW~~~~~~60
0.5
60 0 W
~.
1.
2.0
00
1.0
C,
303
3.0
3.0
C,,I
0 90104
90 o .
O'E
08
12
12O'E
180,
16
~~~~~2
.5
20
120OW
60'W
solar
(botton panel).
angle for the
W maximum
60 *
2.0 2.5The top
30panel shows the
E variation of the solar zenith
0~~~~~~~~~~~~~~
.5 same month.
~10-49
80
20E
W
1.5 60
3.0
O1.
it
3.5~~~~~~~~~~~.4.
o
0
a"30 - 9004
3C3
IzI"'
r-,
08
2~~
60
16
20
~0
2.5
30
90
12
~~~~~~.02.
-~~5
04
81
62
Figure
10-43The
Figure
10-43. The 00
geographical
frequency of the
solar minimum (center panel) and
gclvraino variation of the
h critical
rtclfeunyo
h E layer
ae as
a a function
ucino of local
oa time
iefrJnfor June attslrmnmm(etrpnl
n
J 90
60.5
O
B
60
12
soau
aiuhto
ae) 4 h oppnlsostevraino
h olrznt 16 nl
o h 20aemnh
"r'~~~~~~~~~~~~~~~~~~~04
CHAPTER 10
1. DIURNAL-variation throughout the day, which is
largely due to the variation of the solar zenith angle.
2. SEASONAL-throughout the year.
3. LOCATION-both geographic and geomagnetic.
:
4. SOLAR ACTIVITY-both long term and disturbances
5. HEIGHT-the different layers.
These variations have all been deduced experimentally,
by world-wide observations of the ionosphere over the past
few decades. The reader may refer to Davies 11966, Chapter
31and Hargreaves [ 1979, Chapter 5] for details. The diurnal,
seasonal, solar cycle and height variations of the ionosphere
may all be deduced by routine monitoring of the ionosphere
at one location. Figure 10-42 shows these four variations
for a typical midlatitude station. 1958 was a period of high
solar activity, as indicated by the high values of the ionospheric index, T (see Section 10.4.4). The figure also illustrates the mid-latitude seasonal anomaly, the name given to
the initially unexpected fact that foF2 is higher in the winter
than in the summer, in spite of the larger solar zenith angle.
This anomalyandothersare describedby Hargreaves[ 1979,
Chapter 5].
Once the diurnal, seasonal, solar cycle and height variations of the ionosphere at a given location have been deduced, the next step is to measure and understand the variation with location. This has been achieved through an
international effort of observations and data exchange, and
we now have reliable maps of the world-wide distribution
of the important ionospheric parameters. The accuracy of these mapsover the ocean areas, whereno observationsare
Z 40
w 20
/
4.2
-J 20
WO
60
so
8
I
00
I
04
I
08
I
12
16
80
60
10-50
20
24
4.5
5.5
/40/
4\0
04 0
z 40.5/
w 20
o
-J 2 0
6.8
6.4
/
/ GROUND
DAY-NIGHT
LINE
_
04
08
I I I I
12
20
16
24
LOCAL TIME
/DAY- NIGHT
D
LINE
4.0
3
X]11
MHz .
-=5
.___
50
,)
1
0_
-)
\(
9>
8
_~ ,o ~30
{ 50
'" '~6~-'
~
2_3
5 -
10 </ C)
"-
"
\\
/ /
<much
7
._
70
s90o )
0
500E180
L===step
1
I'
s
150-W 120W 90W 60W 30W 0c
GEOGRAPHIC LONGITUDE
N90
Lw=_
8
7
a
30
10I
8
_
98=
N5\\\\\(_=--
7c
.1
2.
30
,o50
C
a
>
70 -
'oo[
7c
,~~)'o~
sgo9,_|-
=-<
5
o
the
\\\4\\> lower end
/at
of the HF band. A review of the D region
\'4,.
30E
i~oo~
< 8'been
C
-,|I
60 L2~
'GEOGRAPHIC
ONGIT;;0UDE
10-51
CHAPTER
10
.. .:.'..: .. ::' :
';)~.-
2'
'1.: .::::...::::-... .
IONOSPHERE:
.tinuously
''-
'. '-IONO.P..
'-~
.':"'::*
.OFINCIENcE
ELETION
ANGLE
A7..... __.._./~
manycalculations
for HFpropagation.
E- MODE
F-MODE
I- hop 2-hop I-hop
9.0
20.5 28.1 48.4
4.4
2.6
1.9
1.3
CIRCUIT:FORA 1OOOkm
12-hop
elevation angle
K* sec ()
E- MODES
90
90
80 -
8.0
8_0-
70
7.0
70 -
aW 60
,,\
-J
F-MODES
9.0
8.0
--
7.0
ELEVATIONANGLE
6.0
\2-HOP
5.
50
5.0
W 60o
L
6.0
I-HOP
ELEVATION
ANGLE 2-HOP
50
\
(9
9.0
.x
4
4.0
o
v-Y
w"J~
,30V
20C
>
- 3,0
seca~~~
~--_
101
0
I 1 1111111
II1111111
I 1
No
/\2.0
~~~~~K*
PATH LENGTH(km)
%
o
o0
0.0
40-
4.0
1~i
>
I-HOP
I~I-HOP
I-HOP
, 30-I
ISA
20C
)~K (PHIO~
e
2.0
sec (PHI)
1.0
01
3.0
\,
1
11
1 1 1111111
0.0
U)
PATHLENGTH (km)
Figure 10-47. Plots showing the elevation angle and obliquity factor (k sec 0) for propagation via the E and F layers. For a given circuit length and
number of hops, these plots show the required elevation angle (necessary for selecting an appropriate antenna) and obliquity factor for the
two layers (useful, among other things, for consideration of possible E layer screening).
10-52
.......
LOWER.
..L . O
.. IONOSPHERE
o o.MHEIGHT
.-
OF MAX ION
OENSITY
. . . . .
*-
Figure 10-49. Sample ray path for a fixed frequency but varying elevation
angle. If the operating frequency
propagated ray.
Figure 10-48. Sample ray paths for a fixed distance but different fre-
quencies.As the frequencyincreases,the ray mustpenetrate further into the ionosphere.If the frequencyis too
high, that is, above the MUF for the circuit, the ray pen-
etrates.
10-53
CHAPTER 10
CHOICEOF MODE
::::::::.:.:...............................................................
Under
normal
conditions,
there
are threecomfactors
that must
be operating
considered
to achieve
successful
are
These
munications.
1. Choosing a suitable operating frequency.
2. Choosing a suitable antenna system.
3. Ensuring that the wanted signal is at a level above
that of the local radio noise at the operating frequency.
The choice of a suitable operating frequency is the subject of main concern in the present context since this is
where a knowledge of ionospheric physics and radio wave
m..,.....................\
.............................................
.. .......
................
. . ............
/
Tx
Rx
0-200
600-90
200-500
400-70
250-50
and 10-20
1000-2000
10-30
and low angles
2000-3000
5-15
and 20-30
500-1000
> 3000
10-54
Required Radiation
(Elevation) Angles
low angles
300
km
IOOO
Figure10-51.Sketchillustrating
thefactthattheantenna
shouldbechosen
tomatchthecircuitlength.Theantennapatternillustrated
is notappropriate
for shortcircuits(R - 300km).
propagation via the IF mode. The impact of the choice of
an antenna for a given circuit is illustrated in Figure 10-51
for hop lengths of 300 and 1000 km. In this case, the antenna
heavily favors the longer circuit and is quite inadequate for
short-haul circuits.
Given a frequency that the ionosphere will support and
an antenna which emits sufficient power in the direction
taken by the signal that arrives at the receiving site, the
third thing to ensure is that the signal strength is above the
strength of the local radio noise. This noise can be natural
or manmade. Below about 20 MHz, natural noise is caused
by either distant thunderstorms, which cause a general increase in background noise level, or local thunderstorms
which are usually much more obvious causes of poor signalto-noise ratio. Galactic cosmic noise becomes the dominating natural noise above -20 MHz when it penetrates the
ionospherefrom above at frequenciesabove foF2. The worldwide distribution and characteristics of atmospheric radio
noise can be obtained from the CCIR Report 332 [CCIR,
19631. Manmade noise includes such things as industrial
noise due to welders, diathermy machines, car ignition,
power lines and so on. Interference from other communicators using the same operating frequency can also be regarded as noise.
There are many techniques available to ensure an adequate signal-to-noise ratio. The more environmentally acceptable method is to aim at lowering the noise level. This
can be achieved by choosing as a receiver site some location
that is remote from the major sources of manmade noise,
which usually entails being away from major cities. Some
noise rejection is also possible using horizontally polarized
antennas which de-select the local noise that tends to be
vertically polarized. Careful attention to the azimuth of the
main lobe of the receiving antenna can also result in the
beam not being aimed at a nearby source of noise. In heavily
urbanized areas (Europe, East and West Coast of the United
States), the noise level is often set by the large number of
other communicators.
The first step in predicting HF communication conditions is to set up an appropriate model of the ionosphere.
To have practical application, this model must include all
five main variations of the electron density distribution of
the ionosphere (altitude, location, diurnal, seasonal, solar
cycle), must include some measure of the remaining variability of the ionosphere after these main variations have
been accounted for, must exist as a reasonably efficient and
fast computer code, and must possess some method for
projection forward in time. Ideally, it should also be capable
of a modelled response to short time scale events such as
shortwave fadeouts and ionospheric storms. (These are discussed in Section 10.4.6.)
Most of the ionospheric models used for communication
predictions are empirical models based on world-wide observations of the ionosphere over the past four decades.
Observations of the main parameters of the ionosphere, foE,
foF1, foF2, M(3000)F2, and Es have been used to produce
worldwide contour maps of monthly median values of these
parameters for each hour of each month and for low and
high levels of solar activity. An example of such maps is
given in Figure 10-45. To calculate the maximum usable
frequency (MUF) for a given circuit, hour and month, the
values of foF2 and M(3000)F2, for example, are determined
by interpolation in the appropriate world map for the expected point(s) of reflection of the signal. For a level of
solar activity other than the low or high levels for which
the maps are drawn, values of the ionospheric parameters
are obtained by interpolating in each of the low and high
activity level maps and then interpolating again in these
results to find the values appropriate to the given level of
activity.
In practice, the only prediction actually performed in
making predictions of suitable HF operating frequencies is
that of the general level of solar activity that can be expected
to pertain at the epoch for which the predictions have been
made. This prediction is usually made in the form of some
ionospheric or solar index, which is a single parameter describing the gross behavior of the ionosphere, and which is
used to drive the computer program. The most common
index used for prediction programs such as the U.S. programs IONCAP and ITS-78, and by the CCIR, is the average
sunspot number, R12. However, Wilkinson [1982] has shown
10-55
CHAPTER 10
A PLOTFORBRtSBANE
JUNEOOUT
200
______________________________
/r57
48
150
55
47
/59
a:
/79
49
69/70
X 60
w
Z
m
W
/67
46
71/
o
56
50
61
73
77/ 5
74years
6
43
0
53
54d756
6f/1
6Dudeney
10
12
14
FREQUENCY(MHz)
Figure 10-52. The variation of the monthly median value of foF2 for
June, 00 UT, at Brisbane, Australia, as a function of the
10-56
/l
'x,,
8_
'_-X
,-'
_ -,
N
I \
E
\,
>6 -
-.
.. '\-..
'
..
/-|/
z x
0o
.
x
4 -
x
x
....
UPPER
DECILE
MEDIAN
_-.-.
LOWERDECILE
XXXXXX MEASUREDVALUE
FOR A MAGNETICALLY
DISTURBED
PERIOD
00
00
04
08
Il
12
16
2
20
I
24
I50'EMT
10.4.5.2 SporadicE. Sporadic Eisatwo-facedphenomenon having both advantages and disadvantages. There is
no doubt that when a strong sporadic E layer is present, it
presents an ideal propagation mode suitable even for medium speed digital data transfer. However, a dense Es layer
can screen the F layer, preventing signals from ever reaching
that layer and causing Es-F combination propagation modes
that may have none of the desired characteristics-the signals may miss the target or may arrive at an angle not
specified in the original receiving antenna design. Many
prediction schemes for HF propagation include statistical
occurrence rates for Es propagation modes, but the Es layer
remains essentially unpredictable.
10.4.5.3 Problems in High Latitudes. After major solar
flares with particular characteristics, high energy protons
penetrate into the lower ionosphere at the poles and causes
wide-spread and long-lasting disruptions to HF communications. The polar cap absorption event is discussed briefly
in Section 10.4.6. The penetration of electrons with energies
in excess of 10 keV into the D region leads to increased
absorption in the auroral zone. This auroral absorption may
have severe consequences for circuits crossing the auroral
zone, but strong absorption is usually limited geographically
to patches a few hundred kilometers in extent and the duration is typically half an hour to a few hours.
An important feature of the undisturbed high latitude
ionosphere is the great variability in space and time. This
variability severely limits the usefulness of a median model
of the ionosphere. The F region trough, which marks a
transition between the midlatitude and high-latitude ionospheres, can have severe and detrimental effects on signals
propagating through it. The sharp gradients at the walls can
cause reflections and result in off-great-circle propagation.
Propagation in the auroral region may introduce rapid
and severe fading of HF signals. Diversity techniques or
some sort of real time channel evaluation technique (see
Section 10.4.9) therefore become almost mandatory.
The polar cap ionosphere in darkness, of importance to
long distance HF communication between higher latitude
stations such as in Canada and northern Europe, is essentially unpredictable for HF purposes. Enhanced F ionization
regions resulting from particle precipitation are randomly
distributed within a low density background ionization. Especially during slightly disturbed conditions polar cap plasma
convection moves high density particles anti-sunward at
10-57
CHAPTER 10
high velocities (Buchau et al, 1983]. These changes will
result in rapid MUF and mode variations and will in general
lead to poor channel performance. MUF variations of typically several MHz from hour to hour and day to day were
observed on an arctic link [Petrie and Warren, 1968]. In
summertime the Fl mode tends to dominate the mid-morning MUF, and during winter/solar minimum conditions Es
modes over transpolar circuits (for example, Andoya, Norway to College, Alaska) have a greater than 50% occurrence
[Hunsucker and Bates, 1969].
10.4.5.4 Problems in Low Latitudes. HF communication problems at low latitudes due to steep spatial and temporal gradients have been discussed by Lakshmi et al. [1980].
The very steep gradients in foF2 during sunrise hours give
rise to several difficulties [Lakshmi et al., 1980]:
1. HF link operators are expected to get their frequencies
cleared from the appropriate governmental authority well in
advance and it is usual practice to fix one frequency for the
daytime and another for the nighttime. The use of the night
frequency during sunrise will require much more power than
is normally permitted while the frequency allocated for the
daytime will be higher than the MUF during the transition
period.
2. Point-to-point links normally use inexpensive tuned
directional antennas, and frequent change of operational
frequency is deleterious from the point of view of antenna
efficiency.
3. In the case of long distance circuits in the east-west
direction involving multi-hop F region propagation, the
problem of the sunrise period will extend to a large number
of hours, because the different F region reflection points
will fall in the sunrise transition location at different periods.
The steep gradients associated with the equatorial anomaly cause problems with north-south circuits. For example,
if we consider the anomaly peak in the northern hemisphere
to be at 15N geomagnetic latitude, if a north-south
HF circuit is operating such that the reflection point is on
either of the sides of the peak, and if the frequency of the
link is very close to the MUF, a peculiar situation arises.
If the point of reflection is equatorward of this anomaly
peak, the radio waves incident on the ionosphere for the
northern circuit will continuously encounter increasing levels of electron density due both to the vertical gradient as
the radiowave penetrates higher into the ionosphere and the
horizontal gradient as the wave progresses in the direction
of increasing electron density. On the other hand for the
same link in the return direction, the horizontal gradient is
reversed. Thus the real MUF values for the two opposite
directions in the same circuit can vary by a large margin
depending on the angle of incidence and on the magnitude
of the horizontal gradient. In fact, rather frequently, especially when the operating frequency is close to the MUF
(calculated ignoring horizontal gradients), only one way
communication would be possible. This has been one of the
10-58
10.4.7 UnusualPropagationModes
The usual monthly-median HF predictions normally assume simple propagation modes such as the IF and 2F
modes and in general these predictions are quite successful.
Some prediction systems also include Es and such propagation modes as lEs and Es-F. However, it is found in
practice that other unusual propagation modes can also exist.
In general these unusual modes have one feature in common-they rely for their support on some particular feature
of the ionosphere that is restricted in latitude. We shall
consider here a few examples of such modes.
Possibly the most useful propagation mode not normally
predicted is the F2 super-mode encountered on transequatorial circuits during periods when the equatorial anomaly
is well developed, that is, during the afternoon and early
evenings, during equinoxes, at high levels of solar activity.
This mode involves two F region reflections without an
intervening ground reflection and is characterized by high
signal strength, low fading rates and an MUF 10-15 MHz
higher than for the normal (2F) mode. It is often described
as the afternoon-typeTEP (transequatorialpropagation)mode.
10-59
CHAPTER 10
A second mode associated with the equatorial ionosphere, and which can coexist with the afternoon-type TEP
mode in the early evening, is called the evening-type TEP
mode. It is characterized by strong signals, rapid flutter
fading and frequencies well above the normal 2F MUF. A
propagation frequency of 102 MHz has been regularly observed on a Japan-Australia circuit. The propagation mode
in this case is probably a ducted mode, the signals traveling
within the walls of equatorial "bubbles" [Heron and
McNamara, 1979 and Winkler, 1981]. This mode is also
most likely to occur during the equinoxes at solar maximum.
A review of TEP has been published by Nielson and Crochet
[1974].
Enhancements of the MUF on transequatorial circuits
have also been attributed to combination modes in which
one hop is via a reflection from the equatorial Es belt which
stretches a few degrees either side of the magnetic equator
[McNamara, 1974a]. In an Es-F combination mode, the
MUF is enhanced because the F layer hop is longer than
the usual hop length on a 2F mode. Similar MUF enhancements on nighttime circuits have been attributed to a combination mode in which one hop relies on scatter from F
layer irregularities [McNamara, 1974b]. Propagation above
the MUF due to scatter by small scale irregularities can
routinely be observed [Rawer, 1975].
The high latitude also produces its share of unusual
propagation modes, mostly because of the presence of the
midlatitude trough. Strictly speaking, the modes are not
different from those expected-they are just heavily affected
by the presence of the trough.
Another interesting propagation mode is the round-theworld (RTW) propagation mode. This mode relies on aspects of the vertical distribution of electron density rather
than on latitudinally localized features of the ionosphere.
Consideration of observed elevation angles, signal losses
and time delays have led to the conclusion that RTW signals
propagate via a chordal mode of propagation, rather than
by a uniform multihop mode.
10.4.8 Short-TermForecasting
of HF Conditions
Long-term predictions of monthly median parameters of
the ionosphere such as described in Section 10.4.4 are the
traditional approach to frequency management. Frequencies
are selected which should ensure at least a 90% success rate
for communications at all times of day, during any season,
and at all parts of the solar cycle. Engineering decisions
such as transmitter power and antenna configuration are also
made at this stage. Optimum working frequencies (OWF)
are derived for each hour of the day, usually for a month
at a time, using empirically derived frequencies and signalto-noise ratios. Use of predictions in this manner continues
to provide acceptable, though not high quality, communications for many purposes such as voice and telegraph con-
10-60
ditions. For a survey of the historical development of forecasting methods for HF propagation see Rawer [1975].
The effects on the ionosphere and communications of
solar disturbances which cause ionospheric storms have been
discussed in Section 10.4.6. However daily values of foF2
and MUF are known to vary by about 15% to 20% from
the monthly median values during quiet times as well as
during ionospheric/magnetic storms. See, for example, Figure 10-53. These variations may be superimposed on slower
upward or downward drifts in values over several days. It
is desirable to predict all of these variations for the purpose
of efficient radio communications.
There are several possible approaches to the short-term
forecasting of ionospheric parameters.
1. Using associated geophysical parameters.
2. Using the ionospheric parameters themselves.
3. Using ionospheric indices.
The day-to-day variability of foE and foF1 is so small that
monthly median values of these parameters can be used to
represent the daily variation in the E and Fl regions [Rush
and Gibbs, 1973] and therefore offer little difficulty.
Short-term forecasts of foF2 have been made by relating
changes in foF2 to corresponding changes in selected solar
and geophysical variables such as the 10.7 cm solar flux
and the geomagnetic index Kp [Bennett and Friedland, 1970;
Ichinose et al, 1980]. The disadvantage of this type of approach is that the independent variables upon which changes
in foF2 are assumed to depend must themselves first be
predicted. Even if this is done successfully, only limited
success is possible because the correlations of these parameters with parameters of the ionosphere are not very high.
A more successful approach to the short-term forecasting
of foF2, or alternatively of an operational MUF on a given
circuit, is a prediction scheme based directly upon immediate past observations of foF2 or an operational MUF. Such
prediction schemes for foF2 are described, for example, by
Rush and Gibbs [1973], Lyakhova and Kostina [1973],
McNamara [1976], and Wilkinson [1979].
Rush and Gibbs 11973] used a five-day weighted mean
value of foF2 to predict daily and hourly values of that
parameter. The method of Lyakhova and Kostina [1973] is
based on the observation that correlation coefficients between the deviations of foF2 from median values remain
greater than 0.5 for up to four hours. The high correlation
between hour-to-hour variations of foF2 has been discussed
by Lyakhova [1960], Radinov [1963], Gautier and Zacharisen [1965] and Rush [ 1972].
McNamara [1976] made predictions of foF2 at a particular location up to 3 hours ahead by projecting forward
the trend in the departures of the last few hours' observations
from a 15-day running median. Wilkinson [1979] on the
other hand, simply projected forward in time the deviation
of an observed foF2 value from the predicted monthly median value of foF2. He found this technique to be effective
for lead times of up to about 3 hours.
Similar techniques have been applied to oblique circuits
ditions requires the projection forward in time of the effective ionospheric index. In the limit of zero lead-time, the
"forecast" in fact becomes a real-time assessment of the
ionosphere.
A brief review of short-term forecasting for HF propagation has been given by CCIR [1981b]. The special case
of high latitude propagation is considered in CCIR [1983].
CHAPTER 10
referred to that report. One of the more recently developed
techniques of RTCE, which can serve as an example of the
technique, is that of pilot tone sounding.
In pilot tone sounding [Betts and Darnell, 1975], low
level CW tones are either inserted into the data spectrum
or transmitted in potentially available alternative channels.
At the receiver, simple measurements on the tones enable
the relative states of the channels to be specified in terms
of predicted (digital) data error rates. In a situation where
a multicomponent broadcast is being used to radiate identical
data simultaneously on each of several frequencies, this type
of RTCE system allows the best component to be selected
automatically at the receiver.
One of the main advantages of the pilot-tone RTCE
technique is the extreme simplicity of the concept and the
implementation when compared with other techniques for
RTCE, for example, oblique incidence ionosondes. A further advantage of the technique is that it permits RTCE and
data signals to be combined in a simple format, rather than
requiring a separate stimulus. The technique also lends itself
readily to automation. With an automated system, there will
be little difficulty in principle in including into the algorithm
for selecting the optimum working frequency any measurements of other pertinent properties of the pilot tone, for
example, amplitude, Doppler shift, and spectral width. The
noise level in the channel could also be included.
The pilot-tone sounder technique does not permit the
determination of the operational MUF or identification of
propagation modes. The former may cause difficulties when
the MUF is decreasing (for example, due to normal diurnal
variation) if the working frequency is just below the MUF.
Difficulties could be severe if the MUF changes significantly
during the RTCE time (typically 3 min per channel).
The signal optimization problem at HF is potentially
very complex and to date RTCE has been used almost exclusively as a frequency selection aid. Darnell [1975b] lists
the various parameters on an HF link which could be under
the control of the communicator and hence be made adaptive
in response to the RTCE data:
Transmission frequency
Transmitter power level
Bandwidth
Data rate
Modulation type
Time at which the transmission is made
Signal processing algorithm at the receiver
Elevation angles for antenna array beams
Diversity type
To this list may be added antenna polarization for near
vertical incidence propagation.
The potential advantages accruing from the use of RTCE
techniques have been summarized by Darnell [1975a]:
1. The effect of man made noise and interference can
be measured and specified quantitively [Darnell, 19781.
2. The facility for real-time, on line measurement of
propagation and interference allows the use of relatively
10-62
10.4.10 Conclusion
It has not been possible within the present limitations
of space to present more than a cursory overview of the
subject of HF propagation. However, it is hoped that most
readers will be able to follow up on topics of special interest
by going to the references cited. In this regard, the value
of the proceedings of the "Boulder Workshop" (Donnelly,
1979; 1980) cannot be stressed too highly. The CCIR documents and the AGARD Lecture Series No. 127 should also
prove very useful, but possibly somewhat harder to obtain.
of an over-the-horizon (OTH) radar designed for the surveillance of aircraft at ranges of approximately 1000 to 3500
km.
The radar reflections are the result of scattering from
electron density irregularities aligned along the lines of force
of the earth's magnetic field. The characteristics of the fieldaligned scatterers are such that the radar echoes originate
in a small range of angles about perpendicular incidence to
the magnetic field lines and that the amplitude of the auroral
echo is aspect-angle sensitive. In addition to the orthogonality condition, it is necessary that this geometric configuration take place at ionospheric altitudes, that is, 80 km
and above.
The probability of observing radar reflections from ionization irregularities is also dependent upon the frequency
of occurrence of E layer auroral or spread-F irregularity
activity in the region of interest. Reflections can be expected
to occur in regions where both the conditions of near-perpendicularityat ionsophericheightsand high auroral or spreadF irregularity activity are satisfied. The fulfillment of only
one requirement is not sufficient to warrant a radar reflection.
An appropriate model to use to define the condition of
auroral activity is the Feldstein-Starkov [1967] auroral belt
(oval) model. Since the location and extent of the auroral
oval are a function of time and geomagnetic activity, E layer
auroral echoes can be expected to appear over a wide area
during magnetically disturbed periods for an HF backscatter
radar located at high latitudes.
The F layer irregularities can be described in terms of
the probability of occurrence of spread-F derived by Penndorf [1962] or in terms of Aarons' [1973]irregularity scintillation region.
In this section, an estimate is made of the characteristics
of ionospheric clutter that could be observed by an HF
backscatter radar operating in the midlatitude with the antenna beam oriented towards the polar region. The topics
to be discussed are the amplitude, the backscatter crosssectional area, the angular extent, the Doppler frequency
variation, the frequency of occurrence, the diurnal and seasonal variation and the correlation with solar-geophysical
conditions.
10-63
CHAPTER 10
For aurora which fills the antenna beam, o is replaced
the volume scattering coefficient, ov,that is, radar cross
section per unit volume, such that, for a rectangular antenna,
30
AURORAL
20
x1
-by
~2o0~~~~
- a d
rREGIONOF
&
,B
200
300 40o0
500
FFREOQUENCY
(MHz)
o100
uJ
/
_
1> 15
x15
600
-
700
800
800-MHz ORTHOGONAL
POLARIZATION
---.
-10AURORA
800-MHz TRANSMIT.
POLARIZATION
-20
(10.61)
-30
PtGtGrG r
(47r)R KTB NF LsLp
(10.60)
3 4
frequencies in April 1960. The region of uncertainty between 30 and 401 MHz results from the fact that there was
no way to specify to what degree the auroral scatterers filled
the 30-MHz antenna beam. At 401 and 80 MHz, the antenna
beamwidth was 1.2 , while at 30 MHz it was 15 . Thus, it
was assumed that the scatterers completely filled the beam
at the two higher frequencies. However, it was doubtful
that the beam was completely filled at 30 MHz. The 1 x 10
aurora is assumed to be a point target for the 30-MHz antenna.
The data in Figure 10-54 are normalized with respect to
401 MHz. For a filled beam antenna, the relative signal
strength of the 30-MHz echoes was 6-dB greater than those
at 401 MHz, while, for the point target case, there was a
28-dB difference between the 30- and 401-MHz data.
The absolute magnitude of the 30-MHz auroral echoes
can be deduced by comparing the system sensitivities of the
30- and 401-MHz radars. According to Table 10-9, the 401MHz radar was 48.9-dB more sensitive than the 30-MHz
radar assuming that the aurora was a point target and 26.9dB more sensitive, assuming that the 30-MHz antenna beam
was completely filled with scatterers. Thus, it follows that
Table 10-9. Comparison of SRI International-Scotland 30-MHz and 401-MHz radar sensitivities utilized in the April 1960 radar-auroral observations.
Point Target
(401 MHz/30 MHz)
(dB)
Scatterers Fill
Antenna Beam
(401 MHz/30 MHz)
(dB)
30-MHz
Radar
401-MHz
Radar
x2 (m2 )
Pt (Mw)
G,Gr(dB)
(10.0)2
(0.748)2
-22.5
-22.5
0.0015
40
0.12
84
+ 19
+44
+ 19
+44
15
1.2
- 11
15
1.2
- 11
T(IJS)
300
B (kHz)
NF (dB)
L, (dB)
6
9
7
300
6
4.5
3.1
Parameter
(DEG)
3
(DEG)
10-64
0
0
+4.5
+3.9
+48.9
0
0
+4.5
+3.9
+26.9
50
45
SATURATION
1
Z
.-
40-
79
35
Z
LEVEL-2 30 -NOISE
c
-133 dBm
u25
u. 20
0
70
50
30 MHz-
NORMAL
-115
-105
-95
-85
AMPLITUDE
-75
1dBm)
-65
-55 .
0.1
-45
_
_
Figure 10-55. Amplitude histogram of SRI International 401-MHz Scotland radar auroral data [Leadabrand et al., 1965].
0.1
0.01
-120
the magnitude of the 30-MHz radar-auroral data was 20.9dB less than that of the 401-MHz data.
The amplitude distribution of the 401 -MHz radar-auroral
__
I , I
I I
-88
-72
-104
AMPLITUDE (dBm)
-56
30-MHzradar.
Table 10-10. Comparisonof the sensitivityof an HF backscatterradar to that of the SRI International-Scotland
Point Target
(HF Radar/SRI)
(dB)
Scatterers Fill
Antenna Beam
(HF Radar/SRI)
(dB)
HF
Backscatter
Radar
SRIScotland
Radar
X2 (m2 )
Pt (MW)
G, Gr
(dB)
(10.0)2
0.8
45
(10.0)2
0.0015
40
13 H
(DEG)
13v
(DEG)
T (pLS)
15
-4.8
15
15
400
10
2.5
100
12
3
300
Parameter
B (kHz)
NF (dB)
L, (dB)
6
9
7
0
+ 27.3
+5
+ 3.5
- 12.2
-3
+4
(r = 400 Ls) + 37.1 dB
(T = 10 ps) +21.1 dB
0
+ 27.3
+5
+ 1.2
- 14.8
+3.8
- 12.2
-3
+4
+ 33.5 dB
+ 1.5 dB
10-65
CHAPTER
10
Table 10-11. Estimate of radar-auroral clutter levels of a hypothetical HF backscatter radar at 30 MHz based on SRI International data.
Scatterers Fill
Antenna Beam
Point Target
= 4 0 0 ps
T = 10 ps
Parameter
(dBm)
(dBm)
(dBm)
T = 1]0RS
(dBm)
Upper Decile
Upper Quartile
Median
-42.0
- 51.0
-61.0
- 58.0
- 67.0
- 77.0
- 45.6
- 54.6
- 64.6
- 77.6
- 86.6
- 96.6
Lower Quartile
Lower Decile
Standard
Deviation
-71.0
- 80.0
14.8
- 87.0
- 96.0
14.8
-74.6
- 83.6
14.8
- 106.6
- 115.6
14.8
Statistical
f-",
T = 400ps
imately the 10- to 20-MHz range. From simultaneous auroral echo measurements at frequencies of 49.7, 143.5 and
226 MHz, Flood [1960] has deduced a value of n = 3.5
between 49.7 and 143.5 MHz and n = 6.5 between 143.5
and 226 MHz.
Measurements of E region radar aurora at frequencies
in the HF band and at VHF (143.8 MHz) from locations in
Finland and Germany indicated a frequency dependency of
n = 3 [Oksman et al., 1979].
Since there is a wide discrepancy in the experimental
measurements of the frequency dependence of auroral backscatter, which could be due to the characteristics of the
auroral ionization, that is, inhomogeneous distribution of
auroral electrons, varying scale sizes of ionization irregularities and different scattering altitudes, a value of 4 for
the exponent in Equation (10.62) is assumed in this analysis.
The estimated median auroral clutter amplitudes that
could confront the hypothetical HF backscatter radar in the
5- to 30-MHz band are shown in Figure 10-57. The external
(10.62)
-20-T = 400/ts
. -0
10-54, it can be shown that for this power law n = 2.5 for
a point target and n = 0.5 for the filled beam case. However, an examination of other radar-auroral data taken at
Scotland revealed a value of n = 6.8 for a point target and
n = 2.7 for a filled beam [Leadabrand et al., 1965]. An
analysis of auroral echoes from simultaneous multiple frequency observations in Alaska by Leadabrand et al., [1967]
revealed than n = 2 for frequencies between 50 and 850
MHz and n = 5 between
-60
80
E
10
i
'
POINTTARET
PFILLNTA
PULSE LENGTH
<-120
--
10-66
-160
100
6 8 101
2
FREQUENCY
(MHz)
,
4
8 102
100
--0
_100
s
90
70
g901
TARGET
o -POINT
--FILLED BEAM
= PULSE LENGTH
70
T= 10sS
.....
POINT TARGET
---FILLED BEAM
T= PULSE LENGTH
os
60
0- 6[
C-
--
05o
404
10
2 4
6
8102
100
8 101
102
FREQUENCY (MHz)
FREQUENCY (MHz)
cT
,2
(10.63)
where Ah is the thickness of the aurora in the vertical direction. Since Ah was an unknown factor in the Polar Fox
it was assumed to be unity (1m). Thus, Equaexperiment,
II
tion (10-63) can be written as
=
or =
0`AA
('AR[
(10.64)
(10.64)
10-67
CHAPTER 10
Table 10-12. Average area scattering coefficient of E layer and F layer
irregularities based on Polar Fox II data, January 1972.
Area Scattering
Coefficient (dBsm/m2 )
Frequency
(MHz)
Statistical
Parameter
E Layer
F Layer
8.125
Upper Quartile
Median
Lower Quartile
Upper Quartile
Median
Lower Quartile
- 34.2
-47.6
-58.6
-56.5
-63.9
- 74.0
- 27.2
-33.6
-42.0
-37.1
-47.2
- 54.3
14.875
Table 10-13. HF backscatter radar-estimated average cross-section of E layer and F layer irregularities based on Polar Fox II data, January 1972.
Cross Section
Ratio
(dB)
Average
Cross Section
(dBsm)
Frequency
(MHz)
Pulse Length
(Ais)
Statistical
Parameter
E Layer
F Layer
F Layer/E Layer
8.125
400
Upper Quartile
Median
Lower Quartile
Upper Quartile
Median
Lower Quartile
Upper Quartile
Median
Lower Quartile
Upper Quartile
Median
Lower Quartile
61.2
48.0
37.0
45.4
32.0
21.0
39.1
31.7
21.6
23.1
15.7
5.6
71.8
65.4
57.0
55.7
49.3
40.9
61.9
51.8
44.7
45.8
35.7
28.6
10.6
17.4
20.0
10.3
17.3
19.9
22.8
20.1
23.1
22.7
20.0
23.0
10
14.875
400
10
10-68
120
__-___-r_
DOPPLER
[T
I
I
SPREAD (LEFT SCALE)
60
100
0.6
o.2
<0.4
~::~-
='
,0.2
16
24
2 DE-AND F-LAYER-
g~~~~~~~~od
__I
32
40
48
56
64
72
80
80
40
60
<~c
88
nr
.j
-
/+20
40 .
a.
0
Figure 10-60. Cumulative distributionof the angular extent of ionospheric
clutter based on Polar Fox II data.
_+300
-/z-6 dB
ca
20
+10
0
FREQUENCY
(MHz)
100
-18 d
_+80
+60
-6 dB
+40
Statistical
Parameter
E Layer
Clutter
F Layer
Clutter
E and F Layer
Clutter
20
Upper Decile
Upper Quartile
Median
87.7
83.8
60.2
76.0
42.4
18.8
87.0
74.4
27.6
Lower Quartile
Lower Decile
19.0
11.6
10.0
2.8
13.6
3.4
10
15
20
25
30
FREQUENCY (MHz)
Figure 10-62. Estimate of E layer HF radar-auroral total doppler frequency
10-69
CHAPTER
10
10.6.6 Conclusions
The radar-auroral clutter-to-noise ratio in the 5-30MHz
frequency range, as predicted from the SRI International
radar data acquired in Scotland, is dependent on the frequencydependencelaw of the auroral backscatterecho power.
Assuming a frequency variation of f-4 , the clutter level maximizes at approximately 10 MHz. When a f 2 law is assumed, the maximum clutter level is shifted to a frequency
of 30 MHz.
For an assumed f-2 dependence, the radar-auroral clutterto-noise ratio at 5 MHz is 15.6-dB less than that derived
on the basis of an f-4 law. Assuming that the radar aurora
is a point target, the clutter-to-noise ratio for a 400 us pulse
is estimated to be 3.6-dB greater than that deduced for
auroral scatterers that completely fill the antenna beam. For
a 10-us pulse, the corresponding difference is 19.6 dB.
An analysis of the HF ionospheric backscatter data obtained in Maine under the Polar Fox II program reveals that
the F layer cross section could be 10 to 20 dB greater than
the E layer cross section at 8.125 MHz and 20 to 23 dB at
14.875 MHz.
The median value of the angular extent of E layer, F
layer and combined E and F layer ionospheric clutter appears
to be on the order of 60.2 , 18.8 and 27.6 , respectively.
Extrapolating from the SRI International-Scotland UHF
radar-auroral data, it is estimated that for 5-MHz radar transmissions reflected from E layer auroral ionization, the
1-sigma value of the total Doppler frequency deviation, that
is, Doppler frequency shift and spread, could be on the
order of 13 Hz at the -6 dB normalized signal level,
+ 17 Hz at -12 dB, and + 19 Hz at -18 dB. At 30 MHz,
the 1-sigma Doppler deviations should increase to approximately + 78 Hz, + 105 Hz, and + 114 Hz, respectively.
HF field-aligned backscatter echoes usually occurs most
often near local sunset. However, experimental observations
have shown, at times, peak activity near the midnight hours.
There is a strong correlation of backscatter echoes with
the solar cycle. That is, the percentage of days displaying
backscatter reflections decreases with decreasing sunspot
number.
The characteristics of sunset echoes are found to be
different from those occurring during the night. The sunset
echoes are generally of the diffuse type, and are independent
of magnetic activity while the nighttime echoes are discrete
and correlate with magnetic activity.
10-71
CHAPTER 10
oo0
LOCALTIME
-I
GEOMAGNETIC
\
y\
9/ 95
753
TO
) No
S
PERCEUDENT-DWELL
Figure
le
SUNSET
6IM
j.I
10-64.
S
ample
fading
SECONDS
MET
ORDINATE
SGNXCEEDED
LEEL
trbution
-'3
intensity
Of
FAOEOdRT1o
TIME
9 SS
produc
WA
SUNRISE
-6
20"
----- -----IY)~50
216
DENSITY
OF
LOW
CROSSHATCHING
Moderate
solar During
activityEXCED.D
70
90,95
763
99
FADE DORATION
Io
SECONDS
irregularities. Fade duration and cumulative probability distribution are also shown.
SI
(10.65)
mama
of continuoussolar visibility.
where Pmaxis the power level of the 3rd peak down from
indicatedbelow.
10-72
S4
dB
0.075
0.17
0.3
3
6
0.45
10
Scaling of the chart records is facilitated by simply measuring the decibel change between the Pmaxand Pminlevels.
2k sin K2 LA
K
k
Frequency
Dependence
K2
re
(DN(K
KYK,)
In the presence of a relative motion between the propagation path and the irregularitiesthe spatialspectrum0(Kx,Ky)
is convected past the observer and a temporal variation of
signal phase and amplitude is observed. If the irregularity
structure does not change during the convection ('frozenin' hypothesis) and the irregularities have a uniform velocity
u in the x-direction then the power spectra of log amplitude
and phase in the frequency domain are given by
1 +
.U
z directions
K:
(10.66)
10.7.4
\ 2k /
KL
DX.
s ((K
f
x cosF-(z
KY)dKY
i(
)2k
(10.67)
(PAN,(,KOdK
Lkk
2/J
where
w-angular frequency of phase and amplitude fluctuation
2
q2 =
+ K.
u
The term within the bracket with the upper sign is known
10-73
CHAPTER 10
as the amplitude or Fresnel filter function and that with the
lower sign as the phase filter function. The Fresnel filter
function oscillates with the variation of frequency f and attains its first maximum at the frequency fF = u/ 2Xz .
The behavior of the phase filter function is very different
from the amplitude filter function as it fails to attenuate the
low frequency regime. Equation (10.67) provides a relationship between the irregularity power spectrum in the ionosphere and the amplitude or phase scintillation spectrum
obtained on the ground. Since the irregularity power spectrum has a power law variation [Dyson et al., 1974; Phelps
and Sagalyn, 1976] of the form K P, the power spectrum of
amplitude scintillation shows a maximum at the frequency
fF due to the Fresnel filter function. On the other hand, the
phase fluctuations are dominated by the low frequency regime. At f > fF, both amplitudeand phase scintillationspectra
show an asymptotic variation f-p when the three-dimensional irregularity spectrum is of the form K P. Thus from
a study of weak scintillation spectra the spectral form of the
irregularities causing scintillations may be deduced.
Figures 10-65a and 10-65b show two samples of weak
phase and amplitude spectra obtained from 244 MHz scintillation observations made at Goose Bay, Labrador by the
use of the Fleetsat geostationary satellite [Basu et al., 1982].
The phase and amplitude scintillation data were detrended
by a high pass detrend filter with a cut-off frequency of
0.0067 Hz and the data sample yields an rms phase deviation
of 2.4 radians and amplitude scintillation index S4 = 0.51
conforming to weak scintillation criterion. Both spectra show
an asymptotic variation in the high frequency region, the
amplitude spectrum showing f-2.35 variation, and the phase
spectrum showing f 2 .6 variation. If we consider that the
scintillation spectra have an average f-2.5variation, the corresponding three-dimensional irregularity power spectrum
is expected to have a power law wavenumber spectrum of
the form K-3. 5 . The decrease of power spectral density at
the low frequency end of the phase spectrum is caused by
the detrend filter. It may be noted that in the amplitude
scintillation spectrum the high frequency roll-off starts at
the Fresnel frequency of about 40 mHz. For an irregularity
layer height of 350 km, the observed Fresnel frequency
yields the irregularity drift velocity as 37 m/sec.
When scintillations become intense, the theory outlined
above does not hold, and strong single scattering as well as
multiple scatter theories appropriate to such cases have been
developed [Yeh and Liu, 1982 and references therein]. In
such cases it becomes difficult to relate in a straightforward
manner, the scintillation spectra to the irregularity structures
in the ionosphere. Figures 10-65cand 10-65d show a sample
of intense amplitude scintillation data and its spectrum
(S4 = 0.88) analyzed from 257 Hz scintillation data obtained at Ascension Island near the crest of the equatorial
anomaly by the use of transmission from the geostationary
satellite, Marisat [Basu and Whitney, 1983]. In contrast to
the weak amplitude scintillation spectrum shown earlier,
these spectra show a flat low frequency portion and increased
10-74
bandwidth of the frequency spectra, implying the development of shorter scales in the diffraction pattern. It is found
that the correlation lengths get progressively smaller with
increased strength of scattering [Rino and Owen, 1981; Basu
and Whitney, 1983]. Under conditions of strong scattering,
the phase scintillation spectra are also believed to suffer
from the refractive scattering from very large scale irregularities [Booker and MajidiAhi, 1981]. Since wave propagation through a strong irregularity environment has considerable systems applications, intensive work in this area
is in progress.
ll
[ IV 11 Ill
11
1 1111[
1 1I 1ll~ll
I I I1
I1
-20.
-26.8
P=
-2.60
II 1
20
T - 276
P
-235
10
10.
0.
-~~~~~~10
-10
T
0.00lO
lllI
IO.1
0.01
ll
-50
rl]
II[
TI
I I
[lllmI
Tl
0001
2158UT
0.1
lII
10
FREQUENCY-HZ
MIslandR.
257 MHZ
I lo
IOTi
(SE-HC)
FRE
122079MAR.
S410.88
UHF
10.88
20
257 MHz
TheASC.
I
2158 UT
I0
122079
S4=
n-6002established
-20
-50
-60
-20o
00rence
spread
andF
1 99Rastogi,
980]
that equatorial
159.is9817997
range
TIME (SEC)
c
that
nighttime
ionospheric
equatorial
irregularity
FREQUENCY- Hz
regions
Figure 10-65. Spectra of (a) phase scintillation and (b) intensity scintillation under weak scatter conditions at 244 MHz observed at Goose Bay, Labrador
on 31 March, 1979. (c) Data of 257 MHz intensity scintillation and (d) its spectrum under strong scatter conditions observed at Ascension
Island on 20 December, 1979.
The evidence from the corelation of scintillation occurrence and spread F [Rastogi, 1980] is that equatorial range
10-75
,o
> 000
km --
20
THREE
*C ,NS IOND
,o N.D
..
'IASCENSION
km
/(
DIMENSIONALE
ONC
GITU
O
PATCH MODEL
Figure 10-67.
Map of equatorial
Figure 10-66.
HUANCAYO, PERU
80
0 ---60 -
NOV-JAN
20
NOV-JAN
20
12
18
24
12 LT
60m
12
24
18
12 LT
60
',
o
GUAM
GUAM
K=O-I
40
FEB-APR
FEB-APR
40-
20 -
20
I
LZ
40
12
18
I~
12
18
20 [Z -
12
(Kp
Z_7N
iI
_,_~
18
12 LT
~1_1
24
1'addtb60=3
-
z.LI
24
____I
12
40
12 LT
18
24
12
12 LT
MAY-JUL
18
24
12 LT
24
12 LT
9mantcodios60
20
.eI_.___I
24
and disturbed (Kp
6
3'
12 LT
- 9) magnetic conditions.U
12
18
PERCENTOCCURRENCE
GREATER
THAN 2d8 (S4:SUNSET
(350 km)
SUNRISE
JUL
AUG
SEP
0CT
NOV
TEE
-0C
,,
JAN
~0
FEB
FEB
,0\ ,o
MAR
APR
.MA
-
JUN
,
15
L
21
03
09
HUANCAYO
APRIL
76
15
LT
54 GH,
OCT 77
10-69
igurePercentage
occurrenceof 1.5 GHz scintillation 2 dB
by the occurrence of L band, 1500 MHz activity, at Huancayo, Peru. That evidence is shown in Figure 10-69 [Basu
et al., 1980b] L-band activity at Huancayo does not suffer
from strong scattering or from saturation (as does 136 MHz
and 250 MHz data on occasion); the data show clear equinoctial maxima.
From available data it appears as if geomagnetic control
of the occurrence of scintillation varies with longitude. The
generalization can be made that increased magnetic activity
inhibits scintillation activity before midnight--except duringthose monthswith very low scintillationactivity(May-July
for the region (00-70W) and November-January in the
Pacific longitudes (135-180E). After midnight the scintillation activity in general increases slightly with the presence of magnetic storms. The data shown in Figures 10-70
and 10-71 are for a year's observation in each case. The
complexities of the magnetic control of scintillation occurrence are illustrated by the variations in the curves of occurrence at each station in each season. For further details
see Mullen [1973].
10.7.8.3 In-Situ Data. Basu and Basu [1980] have developed a model from in situ, theoretical, and scintillation
studies. In their morphological model of scintillations, in
situ measurements of irregularity amplitude, AN/N, as computed from T sees of data are utilized in conjunction with
simultaneous measurement of electron density N. A combination of AN/N and N data provides the required AN
parameter as a function of position and time.
10.7.8.4 Sunspot Cycle Dependence. From the viewpoint of electron density variations, the equatorial region
around the magnetic equator displays a complex pattern.
During the day an increase in maximum electron density
occurs away from the equator. The electron density contours
10-77
CHAPTER 10
---
K= O- I
K= 3+-9
FEB. - APR.
---
.....
0- I
of
60-
in,
'40
257 MHz
60
20
24
12 LT
<'/\2
12
18
NATAL, BRAZIL
'
I
60
40
,-
20
lB
24
12 LT
NATAL, BRAZIL
12
18
24
12 LT
80 _
80
60
60
60
HUANCAYO,PERU
~40
_
20
2
12
18
24
20
12 LT
'A
HUANCAYO,PERU
~40
12
18
24
12 LT
Figure 10-70. Seasonal patterns of occurrence of scintillation activity > 6 dB (S4 = 0.3) at 250 MHz for very quiet (Kp = 0- 1+)
+
(Kp = 3 - 9) magnetic conditions for Nov-Apr.
10-78
12 LT
20 -
12
.
_I
wu
24
'
60 40
0o
ACCRA, GHANA
ACCRA, GHANA
20 -
18
257
'
"
40-
80
80
40 -
12
^A
NOV. - JAN.
1. The equatorial anomaly has considerably higher electron density values in high sunspot number years than in
years of low solar activity,
2. The occurrence of maximum electron density for
anomaly latitudes is near sunset in the years of high sunspot
number and in the afternoon in years of low solar activity.
Thus the post sunset irregularity patches attain high AN
levels in the years of high solar flux.
K =OIK 3t-9
AUG. - OCT.
MAY - JULY
60
60-
257 MHz
257MHz
40
40 -
ARAGHANA
ACCRA,
12
24
18
12 LT
12
18
24
12 LT
60
60
60
NATAL, BRAZIL
z 40
40
NATAL, BRAZIL
40
20 -
20-
ACCRA, GHANA
20 -
20
\
GHANA
60
40 12
40.
20
12
24
18
..
18
PERU
HUANCAYO,
12 LT 40
6
HUANCAYO,PERU
18
24
20 -
24
40 12
12 LT
12
18
Figure 10-71. Seasonal patterns ofoccurrence of scintillationactivity > 6dB (S4 = 0.3) for very quiet (Kp
magnetic conditions for May-Oct.
12 LT
24
12 LT
3+
9)
CHAPTER 10
6O0km
60mSUNS
AURORAL
OVAL
PLASMAPA-USE
400km
pOLAR
HIGH SOLAR FLUX
'
JAN
FES.
/'2
wSOL A
LOW '# -"-,"
*6b0
800
700
SURISE
9)
NOV
DEC
MAR APR
MAY
JUN
CORRECTED GEOMAGNETIC
LATITUDE
JUL
Figure 10-72. Depiction of high latitude irregularities -22-02 LT. Sheetlike irregularities are seen in the auroral oval, rod irregularities at higher and lower latitudes.
SEP _
AUG
r-
OCT -
21
15
CGM
09
MEAN
SCINTILLATION
INDEX
(dB)
NARSSARSSUAQ
1968-1974Kpz4-9
15
LT
1981b]
SUNS
JAN
Perhaps the most consistent studies of long term behavior of scintillations have been made in the auroral zone, at
Alaskan longitudes and along the 70oW meridian.
Both the diurnal pattern of scintillation activity and the
seasonal behavior as observed from one site can be noted
in Figures 10-73 and 10-74. The data used for this long
term study [Basu and Aarons, 1980] were taken over a
period of 6 years from Narssarssuaq by observing 137 MHz
scintillations of the ATS-3 beacon; the propagation path
traversed the ionosphere at -63o CGL.
The long term study used for Figure 10-75 incorporated
data from three observatories (Narssarssuaq, Greenland; Goose
SUNRISE
DEC_NOV
-Bay,
\
DEC
ing 1968-1974.
FEB -
MAR
APR
0MAY -I
JUN -
|D65' MAY
JUL
JUL
_
_
/g
/.
THRUdULY
....
.
<
S~~p~l(~~
I ""-
60<-
>/d
6001~//
NARSSARSSUAO
\V60'-
AOOSEBA
OCT
55,
-15
21-121
9
+GM
09
03
MEAN SC NTILLATION
15---L_6
INDEX (dB)
SAGAmORE
HILL
oIo
00
2
MEAN Sl(dB)
-8 2A
27 PTS
24 LT
Kp 4 9
10-80
~~0
~ ~ ~ ~ 1~
1~|~
I~Ii
_
|
l
POKERFLAT
S~
[1~~~~
I7
3 -
1977- 1978
-F
6
G
_ 1~~~~~
4-_
\I /
55
60
2/ J.'
\\
ture-at least as observed from Alaska. Figure 10-78 illustrates the daytime increase with increasing latitude.
65
/11
70
DIPLAT-
75
80
deg
number.
l,
i 2o
____197850
lo_
|=
5.
/\1977
\_r
--
K 3 -9
430
1020
7
20-
o7
iso
/iiJ
5J--
L___
"2400
)
40b.t-
/-,
'
APR-OCT
oI,r
975
200
10-81
CHAPTER
10
PERCENTOCCURRENCE
GREATERTHAN lOdB
PERENTOCURRNC
GREATERTHAN1
-NO
i\DEC W NTER
h -2--0----
-JAN
e-
EABFA
a:2
HA
9 /L \)~20
~JUNSUMMER
ti;)2A
TA
MAR APRlNG
APR
MAY p20
R
cL
Bnantly
SAMER
40N
'"--'" ~~
5 NoOL
LCG
03
THULE 1979
K =0O 3
AJUL
_eoarcrate
SC_
......
24.0. ,
..CG 200
9
15
UT
S
W
1~ ~~8
>1
//
8, ,,
3," a
855
f;5
//
,,,
\ 2\./ t
\5, c80/
10.7.11
/?
'
C,, ,3
,'
Behavior
.
-6
'vz
:ometer-size
SUN-FISE
SUNSET
LAT.
~~~~~~~75
o~
~The
DAWN-DUSK
ARC DRiFT
Figure 10-81. Schematicof smallscaleanti-sunwardirregularitydrift and
the patch motion.
10-82
CsL
E(Xm,XgT,D,R) + M(X,,,T)
4. H(Xm,Tm,Kp,R) (10.68)
+0.12FD)+0.09As(l
+ 176FD)]
Tm = geomagnetic time,
Kp = planetary geomagnetic activity index.
The three terms in Equation (10.68) respectively describe
the strength of equatorial, midlatitude, and high-latitude
irregularities. The first two have not been tested extensively
against Wideband data but the high latitude term H has been.
The high-latitude term is based on the observation that
there often is a more-or-less abrupt boundary [Aarons et al.,
1969] between the midlatitude region of relatively smooth
ionosphere and the highlatitude scintillation region. It is
located, typically, equatorward of discrete-arc auroras in
the general vicinity of the diffuse auroral boundary.
+ 5. 2)]}2 10o3Kp(l
+0.FD)+O8As(l+ 2FD)]
FD = cos (DA + 0.5) + 0.2 cos [2(DA - 99)]
(10.69b)
Sagamore Hill (53 CGL intersection)
SI(dB) = 0.33 + 0.02(1 + 0.2FD) {1 + 1.2(1 - 0.01FD)
x cos (HL - 0.4 - 0.15 Kp) - 0.3
x cos [2(HL - 0.8)] - 0.1 cos [3(HL
+ 6. )]}[038Kp(1+03FD)+3
IAs(l-02FD)]
2
FD = cos (DA + 56) + 0.7 cos [2(DA - 143)].
(10.69c)
10-83
CHAPTER 10
10.8 IONOSPHERIC TIME DELAY
EFFECTS ON EARTH-SPACE
PROPAGATION
10-5
o
1
10-7
0-
Z 10-9
i
,
\/o
io0MHz
200
IGHz
4 5
FREQUENCY
TEC (s),
(10.70)
10-84
At
k
x TEC, t2
= x TEC,
(10.71)
or:
(10.72)
The value (At) is obtained from the difference of the simultaneous measurements of the total range, including ionospheric time delay, at the two frequencies, f1 and f2, since
the geometric distance is, of course, the same at all frequencies. The quantity f2/ (f, - f2) is called the ionospheric
scaling factor. For ratios of f2/f1 near unity, the required
precision of the differential measurement may be unreasonably large. A plot of this quantity, normalized by fl, is given
in Figure 10-83. In this derivation the contribution of re-
6!
5-
I/
24t/
<3
---
03
~0.4
05
06
07
08
RATIOOF SECONDARY
FREOUENCY f2
0.9
f1.0
PRIMARY
(10.73)
(10.74)
x 1016 el/m 2
TEC (cycles),
CHAPTER 10
quenciesfor ionosphericdifferentialcarrierphase measure-
A2
ments. For the pair of frequencies used by GPS, approximately 1.2 and 1.6 GHz, the differential carrier phase shift,
referenced to the lower frequency, is
hA4 =
1.34
L- x
10 - 7 X
4=
(4(u-
(m2 - 1)
x TEC (cycles),
0)
- (c-
1.34
k4L)=
10 - 7 X
3 x
,, +
2 x k,
TEC (cycles)
(10.76)
thus
10-86
2.68 x
=2.
f(f2
TEC (cycles)
TEC (cycles).
0
f2 )
When
f2
>>
f2
(10.77)
2.68 x 10-7f2m
A2 4
TEC (cycles).
For a carrier frequency of 100 MHz a modulation frequency of 1.93 MHz would be required to give 2ii of second
differential phase for a TEC value of 1018el/m.2 A value
of A2 of 2ii for 1018el/m2 is a reasonable compromise
between the requirement for minimizing chances of an ambiguity in absolute TEC and accuracy in measuring TEC
relative changes. The second difference of carrier phase has
been used with the DNA-002 satellite to make estimates of
the absolute value of TEC [Freemouw et al., 1978].
TEC (Hz).
dt
(10.78)
dt
looo
a 100
z
l
/0 '
0
<
I.0
\-
10.8.4 Faraday Polarization RotationWhen a linearly polarized radio wave traverses the ionosphere the wave undergoes rotation of the plane of polar-00
ization. At frequencies of approximately100 MHz and higher
the amount of this polarization rotation can be described by:
- s
2.36 x 10
= 2
f B cos ONdl.
500
IGHz
(10.79a)
where the quantity inside the integral is the product of electron density times the longitudinal component of the earth's
magnetic field, integrated along the radio wave path. Many
ionospheric workers have used this effect, named for Michael Faraday who first observed polarization changes in an
optical experiment, to make measurements of the TEC of
the ionosphere. Since the longitudinal magnetic field intensity changes much slower with height than the electron
density of the ionosphere, the equation can be rewritten as
( = K Bin x TEC,
.ol
(10.79b)
10.8.5 AngularRefraction
The refractive index of the earth's ionosphere is responsible for the bending of radio waves from a straight
line geometric path between satellite and ground. This angular refraction or bending produces an apparent higher
elevation angle than the geometric elevation. Millman and
Reinsmith [1974] have derived expressions relating the refraction to the resultant angular bending. Perhaps the easiest
expressions to use, as given by Millman and Reinsmith
[1974] relate the ionospheric range error to angular refraction:
AE
x-
AR
R
(10.80)
10-87
CHAPTER 10
.2
20
_=18
o1
VERTICAL
15
15
NO
IN
C3
~f.
oo
Z
N
L
=
<
10
15
20
25
IE
U
<
cos E1
:0
10
25
20
15
For low
elevation angles and satellites well above most
of theionization, R>rosinEo, and the angular refraction can
be
expressed
as:
AE=
cos
2hi Eo AR.
(10.81)
for a TEC
Typical values of elevation refraction error
of 1018el/m2 column are shown in Figure 10-85 for several
frequencies. Notethat, atthe lowest frequency, 100 MHz,
near the horizon the refraction is well over 1.5 degrees! The
curves shown
in Figure 10.85 have been constructed using
the approximationderived by Millman and Reinsmith [1974]
for low elevation angles given in Equation (10.81).
Generally, the range error itself is the main ionospheric
navigation systems,and elevation
for
problemadvanced
angle errors are insignificant. Satellite detection radar systems, on the other hand, do have the requirement to know
accurate pointing elevation angles for their large aperture
arrays, though generally the accurate tracking is done by
using range rate information, and elevation angle is of secondary importance as long as the beamwidth of the antenna
is large enough to see the target.
Errors in the azimuth of radio waves transmitted through
the ionosphere can also occur; they depend upon azimuthal
gradients in TEC which are generally small and which can
usually be neglected in practical cases.
10-88
At =
Af x TEC,
(10.82)
1700
1600-[
1600
1500-
1500
1400--\
1400-
13001200
1300
\
1100\
1000
900
900
800-
800
700
700
600-
600
500-
500
400-
400
300--
300
200-
200
100 I
100 -
,o
O-I'
104
MARCH 1980
(Fio. 7 167)
1200
1100 1000
BOULDER,COLORADO
105
el /cm
106
3
.2
.4
0 I I
.6
el/cm
.8
1.0
1.2
1.4
1.6
units of 106
Figure 10-86. Typicalprofileof electrondensityvs height.In (a) log Ne is plotted;in (b) Ne is plottedon a linear scale.
ionosondes have been used since the 1930s to make continuous, routine measurements of the density at the peak of
the F2 region, measured by ionosondes as foF2, and equated
to Nmaxby
(10.83)
10-89
CHAPTER
10
experimental, data based, or empirical model, in representing the actual worldwide foF2 is due to the large amount
of data available from ionosondes in many regions of the
world. Other characteristics of this model are discussed by
Dandekar 11982]and in Section 10.3 of this chapter.
For the TEC parameter, data availability have been, and
will likely continue to be, much more sparse. First, TEC
measurements have generally been calculated from measurements of Faraday polarization rotation using VHF signals of opportunity transmitted from geostationary satellite
telemetry transmitters. A few lunar reflected Faraday rotation measurements in the late 1950s and early 1960s and
the TEC obtained from a few low orbit satellites did not
contribute significantlyto our knowledge of world-wideTEC
behavior, at least not for modeling average ionospheric conditions. Only since the early to mid-1960s have TEC values
been obtained on a more-or-less regular basis. Even today
fewer than one dozen stations regularly contribute TEC data,
which can be used in TEC modeling purposes, to a world
data center.
Fortunately, most of the contribution toTEC comes from
near the F2 region density peak where models of foF2 are
available. These foF2 models can be combined with some
limited knowledge of topside ionospheric thickness obtained
from topside sounders, along with topside in situ density
measurements, to produce a complete ionospheric height
profile model. The most well known of these models is the
one by Bent [Llewellyn and Bent, 1973] which uses ITS78 coefficients for foF2, and topside exponentials for com-
puting TEC. A representation of world-wide average behavior of TEC is illustrated in Figure 10-87 for 2000 hours
UT. To first order the TEC contours shown in Figure 1087 move westward along magnetic, rather than geographic,
latitude lines, at the earth's rotation rate. The Bent model
was constructed using solar maximum data from the 19681969 period and had to be adjusted upward somewhat to
account for the much higher 1979-1980 solar maximum than
that of 1968-1969. This adjustment was necessary to adequately represent the actual TEC values from stations making observations in March 1980, which was near the maximum of the second highest solar cycle ever recorded in the
more than 200 year history of solar cycle observations. The
Bent model, appropriately adjusted for high solar cycle values, does however, represent fairly well the average behavior of TEC for many locations tested. Other worldwide
ionospheric electron density profile models from which average TEC can be obtained include ones by Ching and Chiu
[1973], and Chiu [1975], Kohnlein [1978], the 4-D model
[VonFlotow, 1978], and the International Reference Ionosphere (IRI)[Rawer, 1981]. The characteristics of some of
these models are described in Dandekar [1982] and in Section 10.3.
Other empirical models of TEC have been developed
directly from TEC data alone, though these have necessarily
been limited in temporal and geographic extent by the available data base. These include models of TEC over Europe
and the Mediterranean [Klobuchar, 1973] for low and medium activity portions of the 11 year solar cycle, and a
406"~~~~~~~~
CONTOURS
OF TOTAL ELECTRON
CONTENTUNITS
MARCH 1980
167)OF 10- el/r
0LATITUDE-400
I___
___
30
010-080
60
9090
-
120
150
L 180
LONGITUDE
COLUMN
Il
Il
I09
120
100
90
60
2000 UT
Figure 10-87.
10-90
30
of curves of diurnal changes in TEC for a northern midlatitude station for twelve months during a solar maximum
period is shown in Figure 10-88. The standard deviation
from monthly mean diurnal behavior is approximately 20%25%, during the daytime hours when the absolute TEC
values are greatest. Figure 10-89 shows the standard deviation from monthly average TEC behavior for the midday hours for a number of stations during the solar maximum
period 1968-1969. Again 20%-25% is a good value for the
standard deviation from the monthly average behavior. The
standard deviation is somewhat higher during the nighttime
hours, but the absolute TEC values are much lower during
these periods.
If a satellite ranging system has error requirements such
that it must correct for monthly average ionospheric time
delay, but still can tolerate the approximate 20%-25% variability of TEC from monthly average conditions, approximately 70%-80% of the ionospheric effect on the system
can be eliminated by the use of an average TEC model such
as the one constructed by Bent [Llewellyn and Bent, 1973].
If the system only requires an approximate 50% rms correction of the ionospheric time delay, the algorithm developed by Klobuchar [ 1975 can be used. On the other hand,
if corrections for some portion of the remainder of the ionospheric time delay are required, after a state of the art TEC
20
JANUARY1979
MARCH
FEBRUARY
APRIL
60-
20
'MAY
JULY
AUGUST
E
1
JUNE
60-
20-
1 20
SEPTEMBER
8g
AIo
i58o
L- +2
OCTOBER
NOVEMBER
14
-t
'
18
1 ,J4 4
DECEMBER
:.
11
TOTALEQUIVALENTVERTICALELECTRONCONTENT
FROMHAMILTON(ATS-5)
(TIME UT)
Figure 10-88
Monthly overplots of TEC diurnal curves for Hamilton, Mass. for 1979.
10-91
CHAPTER 10
80
PERCENTAGE
STANDARD
12 - 16 LT
DEVIATIONS
TEC
,
EDMONTON 1969
RBERYSTWYTH1969
HRMILTON 1969
60 -
STRNFORD 1969
HONOLULU 1969
50 ,
L
HONG KONG
40 -
1968
30
20
10
JRN
FEB
MAR
RPR
MAY
JUN
JUL
RUG
SEP
OCT
NOV
DEC
Figure 10-89. Percentage standard deviations for daytime TEC from the stations indicated.
model, such as the Bent one, has been used to take out the
monthly mean TEC, then the short term (a few hours) temporal variability as well as the geographic variability, of
TEC must be considered.
)0-5]
(10.84)
10-92
Z8
.
LL
zM
.
00
-_
.8
LE
r 0.7
SE
W
.6
t
1000
2000
3000
4000
DISTANCE
IONOSPHERIC
LATITUDE
5000
(km)
SEPARATION
A
70N W LONGITUDE
<
ZH
a)20
r~~Z
0
Z
bJ
UI
00
~~~~1.0
U.
.._
:t
bJ3Z
cc
.4
wH0
n
U-
z0
>
0 'l~~~~~~0
.8
.0
0
Er
:
0
Zz
w
E
.2
1000
3000
2000
IONOSPHERIC DISTANCE (km)
4000
5000
10-93
CHAPTER10
MAY 1980
SEPTEMBER
JUNE
JULY
AUGUST
OCTOBER
NOVEMBER
DECEMBER
FEBRUARY
MARCH
40 i
1801
JANUARY191
.4
APRIL
TOTAL EQUIVALENT
VERTICALELECTRON
CONTENT
FROM ASCENSION
ISLAND(SIRIO)
TIME(UT)
Figure 10-91. Monthly overplots of TEC diurnal curves for Ascension Island, May 1980-April 1981.
SEPTEMBER
AUGUST
JULY
OCTOBER
xXXxx
XX
NO DATA
'""'-~ j
xXF
X
s ~ ~ ~~J-~
_ O2
JANUARY
DECEMBER
6-1 NOVEMBER
-3
,FEBRUARY
CL
uI
o
SAO XXX
PAUL,
0
0
12
18
24
I2
18
24
BRAZIL
I0
6
12
lB 24
xxxSAD PAULO,BRAZIL
12
18
24
LOCAL TIME
Figure 10-92. Monthly average plasmaspheric electron content vs. local time for Aberystwyth, Wales (dashed line) and for Hamilton, Mass. (solid line).
Also plotted are values from Kiruna, Sweden for October 1975 and from Sao Paulo, Brazil for May 1975.
specific large increases in TEC, likely due to auroral precipitation, are not individually predictable, but may be statistically characterized as a function of magnetic activity.
In the polar cap region a negligible amount of TEC data
exists. The absolute TEC values are probably lower in this
region than in the midlatitudes, and the variability of the
polar cap TEC is probably very high.
10.9.5 ProtonosphericElectronContent
Most of the available TEC data has been taken by measurements of Faraday rotation of single frequency radio waves
transmitted from geostationary satellites to ground observers. The electron content obtained from Faraday rotation
observations, while made from radio waves transmitted from
satellites at geostationary satellite height above the earth's
surface, only includes the contribution of electrons up to
heights of approximately 2000 km. This is because the integrated product of the longitudinal component of the earth's
magnetic field times the electron density, above approximately 2000 km, is negligible. The only measurements of
the additional contribution of electrons above the Faraday
maximum height have been made using signals from an
ionospheric beacon on the geostationary satellite ATS-6.
Davies [19801 has reviewed the overall results of the ATS6 experiment. A summary of typical protonospheric electron
content data is shown in Figure 10-92 taken from Klobuchar
et. al. [1978]. Note that the protonospheric values are fairly
low in absolute value.
During the nighttime hours when the ionospheric TEC
is low, the protonospheric contribution may become a fairly
large percentage of the total number of electrons between
a satellite at geostationary height and an observer on, or
near the earth's surface. Unfortunately, no protonospheric
electron content data are available during solar maximum
conditions.
CHAPTER 10
of a few percent of the background TEC. A 10% ionospheric
disturbance with respect to the background TEC is uncommon, while a 1% TEC perturbation is common. Titheridge
[1968] and Yeh [1972] have made studies of the statistics
of traveling ionospheric disturbances (TIDs), in TEC for
midlatitude regions.
A system that requires correction for the rate of change
of TEC cannot rely on models of TEC to provide reliable
information on short term rate of change of TEC information, and can use available TID information only in a statistical manner. The only recourse for a system significantly
affected by rate of change of TEC is to use a dual frequency
measurement technique to directly measure the ionospheric
contribution to range rate.
10.9.7 Conclusions
10.10.1 Chemical Releases
There are at least three categories of systems potentially
affected by ionospheric time delay. For the first category
of user the potential systems effects may be small, at least
under any naturally occurring worst case ionospheric conditions. In the second category, a user may require a nominal
correction for average ionospheric time delay, but is able
to tolerate the 20%-25% standard deviation from average
conditions. He should expect at least a 50% correction for
ionospheric time delay effects using a relatively simple time
delay algorithm, and up to 70%-80% for a state of the art,
fairly complex model. These model corrections can be improved by the use of actual ionospheric measurements within
a reasonable temporal and spatial frame. For the third category of user ionospheric model corrections, even updated
with near-real-time measurements, may not be sufficient to
correct for ionospheric time delay, and the system must then
make its own ionospheric correction. Fortunately, the ionosphere is a dispersive medium and the use of identical
modulation on two, widely-spaced frequencies will allow a
direct measurement to be made of ionospheric range delay.
Two coherently-derived carrier frequencies may be used to
obtain accurate time rate of change information for TEC.
Details of measuring ionospheric effects directly by a system's use of multiple frequencies are available in Burns and
Fremouw [1970].
ATS-3
E 20
..
SKYLAB
LAUNCH
16 -
o
,,12
s,
<
a
~-F~~~~~~~~~~~
/
W
~-
of/
aJ
a | 1. l,,,,,
O0
03
I,,,,,,,,
06
09
I,,
12
5I
14 MAY1973
18
Figure 10-93. TotalElectronContent(TEC)datausedto detectthe"SKYLAB effect" on 14 May 1973. The dashed curve gives the
anticipated diurnal TEC behavior based upon a monthly
K, = 2.4 x 10
O+ + H2
--
OH
+ H
cm3/s
(10.85)
(10.86)
O + + CO 2 -O2
+ CO
Once a molecular ion is formed, its dissociative recombination with an ambient electron occurs rapidly,
d,
H2 O' + +
e -OH
+ H a, = 3.0 x 10-7 cm3/s
(10.88)
,,
-- H + O
(10.89)
2+ 0 + +
e-,0
periments.
K3 = 1.2 x 10 - 9cm 3 /s
(10.87)
OH + + e
a,=
2.0 x
10-7
cm3 /s
(10.90)
and hence an "ionospheric hole" is formed. A review of
rocket induced ionospheric disturbances has been given by
Mendillo [1981].
Figure 10-94 contains a schematic showing the many
physical and chemical processes associated with artificiallyinduced depletions in the F region. It should be noted that
CHAPTER 10
EXPANUSION
SF
by
IVPQH|
CHEMISTRYRRECELC
CHEICAL
XCITAF
ELECTRONS
sl#
ION-MOLECUL
IEXCINION
Of
LL
, -.ItCP
_ NEUTRAL
RISR
IE
EC!.s,
RCOE.
RECOMBINATION
|
RAPID INITIAL
EXPANS
HANCED
|O
PHOTO LELCIRN
*SCAPL 10
CON)UIGAT
IONOFSPHElLRE
l
JTC
C R IRONA
PLASMA
FLO
INTODEPLTE
OFC
|TRIGCERING
i~P(EQUATORIAL)
PLASMA
~~~~ACOUSTIC
~~GRA~VlI~~TY
BEACONSATELHLIESS_
FLUX TURN
IPROTONOSPHRICI
-i
THEORNE-
-A.~
\7/
BEACONSATELLITES
S4nv
C
PLASIA
HOLES
AT\~
v////
POIINCOHERENT
IN- SITU
PROBES
~/'PLASMA
DIAGNOST
SDIAGNOSTICS
INCOHERENT
SCATTER
DIAGNOSTICS
GROUND-BASED -- AIRBORNE -Figure 10-94.
Schematic summary of possible rocket effluent effects upon the upper atmosphere
1985.
mission scheduled for Spring-Summer
absorption of the heater wave near the altitude of HF reflection. The following manifestations ofplasma instabilities
have been observed: (1) artificially created spread F; (2)
strongly enhanced radio wave absorption; (3) the creation
of field-aligned density irregularities which scatter (this phenomenon is sometimes called field-aligned scatter of FAS)
an incident HF, VHF, or UHF wave with virtually no frequency change and make certain types of scatter communication circuits possible; (4) scattering process in which
the frequency of the scattered wave differs from the frequency of the incident wave by roughly the frequency of
10-98
IN-SPACE
B2:56:05
C02:58:36
E.2:59:35
F-3:00:5
G.:03:42
panding6300A airglowcloudproduced
byexcitedoxygen
atomscreatedfromtherecombination
of freeelectrons
and
molecularions (021 , OHt, H2 0z ) produced
by exhaust
molecules
(CO2, H2, H20) andambientatomicions(0O).
Times are a.m., PST [Mendilloand Baumgardner,1982].
The long-wavelength (-1 km) field-aligned irregularities giving rise to artificial spread F [Utlaut et al., 1970;
Utlaut and Violette, 1972; Wright, 1973] could not, however, be explained in terms of the above instability process.
The causative mechanism for the generation of long wavelength irregularities remained obscure for quite a while and
exhibitscintillations[Rufenach,1973;PopeandFritz, 1974;
Bowhill, 1974]. Radio star scintillation measurements at 26
MHz during ionospheric modification indicated the presence
of either rapid and random or deep long-period (-5 mins)
fluctuations. In order to avoid some of the difficulties of
ACCELERATED ELECTRONS
(SOME ESCAPE TO CONJUGATE
HEMISPHERE)
at
HF REFLECTION ALTITUDE
ELECOeV S
TO
-J Wt
300 -
(-km's THICK)
<lOOm THICKSLABOF
... :'{'.::!
4'.'
.t-
*lr.:
r,4F.gAf/
I.
/~
AIRGLOW EXCITATION:
O :'-N2(,029
/
(H10's km THICK)
/.
'
5
fHF
Figure 10-96. Effects produced by ground-based transmitter of power aperture of the order of 104Mwm2 in the 4-12 MHz frequency range. Energy
deposited in the ionospheric plasma alters both the thermal and nonthermal properties of its charged particle population. Controlled
experiments have applications to aeronomy, chemical rates, atomic cross sections, communications, and a number of areas of plasma
10-99
CHAPTER 10
REFERENCES
Aarons, J., "A Descriptive Model of F-Layer High-Latitude
Irregularities as Shown by Scintillation Observations,"
J. Geophys. Res., 78: 7441-7450, 1973.
Aarons, J., "Equatorial Scintillations: A Review," IEEE
Trans. Antennas Propagat., AP-25, 729, 1977.
Aarons, J., "Global Morphology of Ionospheric Scintillations," Proc. IEEE, 70: 360-378, April 1982.
Aarons, J. and R.S. Allen, "Scintillation Boundaries During
Quiet and Disturbed Magnetic Conditions," J. Geophys.
Res., 76: 170, 1971.
Aarons, J. and E. Martin, "The Effects of the August 1972
Magnetic Storms on Ionospheric Scintillations," Radio
Sci., 10: 547-554, 1975.
Aarons, J., J.P. Mullen, and H.E. Whitney, "The Scintillation Boundary," J. Geophys. Res., 74: 884-889, 1969.
Aarons, J., E. MacKenzie, and K. Bhavnani, "High-Latitude Analytical Formulas for Scintillation Levels," Radio Sci., 15: 115-127, 1980a.
Aarons, J., J.P. Mullen, J.P. Koster, R.F. daSilva, J.R.
Medeiros, R.T. Medeiros, A. Bushby, J. Pantoja, J.
Lanat, and M.R. Paulson, "Seasonal and Geomagnetic Control of Equatorial Scintillations in Two Longitudinal Sectors," J. Atmos. Terr. Phys., 42: 861-866,
1980b.
Aarons, J., E. MacKenzie, and S. Basu, "Microwave Equatorial Scintillation Intensity During Solar Maximum,"
Radio Sci., 16: 939-945, 1981a.
Aarons, J., J.P. Mullen, H. Whitney, A. Johnson, and E.
Weber, "VHF Scintillation Activity Over Polar Latitudes," Geophys. Res., Lett., 8: 277-280, 1981b.
AGARD, "Modern HF Communications," AGARD Lecture
Series 127, (avail. through NTIS, Springfield, Va) 1983.
Akasofu, S.-I., "Midday Auroras and Polar Cap Auroras,"
Geofysiske Publikasjoner (Norske Videnskaps-Akademi
i Oslo), 29: 73-85, 1972.
Albrecht, H.J., (ed.), Artificial Modification of Propagation
Media, AGARD Conf. Proc. 192, Brussels, 26-29 April
1976.
Al'pert, Ia. L., "Ionospheric Propagation of Radio Waves,"
Acad. of Sci. of USSR, Moscow, Section on Long Waves
(avail) as Translation T5-60, National Bureau of Standards, Boulder, 1960.
Ames, J.W. and R.D. Egan, "Digital Recording and ShortTerm Prediction of Oblique Ionospheric Propagation,"
IEEE Trans. Antennas Propag. AP-15: 382, 1967.
Ames, J.W., R.D. Egan and G.F. McGintie, "Short-Term
Prediction of HF Communication Circuit Performance"
in Ionospheric Forecasting, AGARD Conf. Proc., 49,
AD700896, 1970.
Anderson, D.N., "Modelling the Ambient Low Latitude FRegion Ionosphere-A Review,"J. Atmos. Terr. Phys.,
43: 753, 1981.
Anderson, D.N. and P.A. Bernhardt, "Modeling the Effects
of an H Gas Release the Equatorial Ionosphere," J.
Geophys. Res., 83: 4777, 1978.
Appleton, E.V. and W.J.G. Beynon, "The Application of
Ionospheric Data to Radio-Communication Problems:
Part 1," Proc. Phys. Soc., 52: 518, 1940.
10-100
Baggaley, W.J., "BackscatterObservationsof F-RegionFieldAligned Irregularities During the I.Q.S.Y.," J. Geophys. Res., 75: 152-158, 1970.
Bain, W.C., "Models of the Mid Latitude D Region at
Noon" in Medium, Long and Very Long Wave Propagation, AGARD Conf. Proc., 305: 34-1-34-8
ADA 113969, 1982.
Balsley, B.B. and W.L. Ecklund, "VHF Power Spectra of
the Radar Aurora," J. Geophys. Res., 77: 4746-4760,
1972.
Balsley, B.B., W.L. Ecklund, and R.A. Greenwald, "VHF
Doppler Spectra of Radar Echoes Associated with a
Visual Auroral Form: Observations and Implications,"
J. Geophys. Res., 78: 1681-1687, 1973.
Bannister, P.R., "Variations in Extremely Low Frequency
Propagation Parameters," J. Atmos. Terr. Phys., 37:
1203-1210, 1975.
Bannister, P.R., "Extremely Low Frequency (ELF) Propagation," NUSC Scientific and Engineering Studies, Naval Underwater Systems Center, New London, Conn.,
February 1980.
Bannister, P.R., "Overview of ELF Propagation," in Medium, Long, and VeryLong Wave Propagation, AGARD
Conf. Proc., 305: 10-1-10-13 ADA113969, 1982.
Barghausen, A.L., J.W. Finney, L.L. Proctor, and L.D.
Schultz, "Predicting Long-term Operational Parameters
of High-Frequency Sky-Wave Telecommunication Systems," ESSA Tech. Rept. ERL 110-ITS 78, 1969.
Barr, R., "The Effect of Sporadic E on the Nocturnal Propagation of ELF Radio Waves," J. Atmos. Terr. Phys.,
39: 1379-1387, 1977.
Barry, G.H., "A Low Power Vertical Incidence lonosonde"
IEEE Trans. Geosci. Elect., GE-9(2): 86-89, 1971.
Basu, S. and Su. Basu, "Modeling of Equatorial Phase and
Amplitude Scintillations from OGO-6 and AE Irregularity Data, Low Latitude Aeronomical Processes,"
COSPAR Symposium Series, 8: 187-199, Pergamon Press
1980.
Basu, S. and H.E. Whitney, "The Temporal Structure of
Intensity Scintillation Near the Magnetic Equator," Radio Sci., 18: 263-271, 1983.
Basu, S., Su. Basu, A.L. Johnson, J.A. Klobuchar, and
C.M. Rush, "Preliminary Results of Scintillation Measurements Associated with Ionosphere Heating and Possible Implications for the Solar Power Satellite," Geophys. Res. Lett., 7: 609, 1980a.
Basu, S., Su. Basu, J.P. Mullen, and A. Bushby, "Longterm 1.5 GHz Amplitude Scintillation Measurements at
the MagneticEquator," Geophys. Res. Lett., 7: 259-262,
1980b.
Basu, S., J.P. McClure, Su. Basu, W.B. Hanson, and J.
Aarons, "Coordinated Study of Equatorial Scintillation
and in-situ and Radar Observations of Nighttime F Region Irregularities," J. Geophys. Res., 85: 5119-5130,
1980c.
Basu, Su. and J. Aarons, "The Morphology of High-Latitude VHF Scintillation Near 70W," Radio Sci., 15:
59-70, 1980.
CHAPTER 10
Simultaneous Magnetic Disturbances," J. Atmos. Terr.
Phys., 27: 1151-1158, 1965.
Brooks, D., "Observations of Auroral Echoes Following
Worldwide Storm Sudden Commencements," J. Atmos.
Terr. Phys., 28: 103-109, 1966.
Buchau, J., W.N. Hall, B.W. Reinisch, and S. Smith,
"Remote ionospheric Monitoring from Effects of the
Ionosphere on Space and Terrestrial Systems," Proc.
IES (Ionospheric Effects Symposium), edited by John
H. Goodman, U.S. Government Printing Office, Washington, D.C., 1978.
Buchau, J., B.W. Reinisch, E.J. Weber, and J.G. Moore,
"Structure and Dynamics of the Winter Polar Cap F
Region," Radio Sci., 19: 995-1010, 1983.
Budden, K.B., Radio Waves in the Ionosphere, Cambridge
University Press, Cambridge, UK, 1961.
Buneman, O., "Scattering of Radiation by the Fluctuations
in a Non-Equilibrium Plasma," J. Geophys. Res., 67:
2050, 1962.
Burke, W. (ed.), Proceedings of International Symposium
on Active Experiments in Space, 23-29 May 1983, Alpbach, Austria, European Space Agency Report ESASP-195, Paris, 1983.
Burns, A.A. and E.J. Fremouw, "A Real-Time Correction
Technique for Transionospheric Ranging Error," IEEE
Trans. Antennas Propag., AP-18:(6), Nov. 1970.
Burrows, M.L., ELF Communications Antennas, Peter Peregrinus Ltd., Stevenage, England, 1978.
Carlson, H.C. and L.M. Duncan, "HF Excited Instabilities
in Space Plasmas," Radio Sci., 12: 1001, 1977.
Carlson, H.C., W.E. Gordon, and R.L. Showen, "High Frequency Induced Enhancements of the Incoherent Scatter
Spectrum at Arecibo,"J. Geophys.Res., 77: 1242, 1972.
CCIR (Comite Cosultatif International des Radiocommunications = International Radio Consultative Committee), "World Distribution and Characteristics of Atmospheric Radio Noise," Report 332, Doc. Xth Plenary
Ass., International Telecommunications Union (ITU),
Geneva, 1963.
CCIR, "CCIR Atlas of Ionospheric Characteristics," Report
340 and 340-2, Doc. XIth Plenary Ass., ITU, Geneva,
1966.
CCIR, "Interim Method for Investigating Sky-Wave Field
Strength and Transmission Loss at Frequencies Between
the Approximate Limits of 2 and 30 MHz," Report 252-2,
Doc. XIlth Plenary Ass., ITU, Geneva, 1970.
CCIR, "Second CCIR Computer-Based Interim Method for
Estimating Sky-Wave Field Strength and Transmission
Loss Between 2 and 30 MHz," Suppl. Rept. 252-2, Doc.
XIVth Plenary Ass., ITU, Geneva, 1978.
CCIR, "Prediction of Sky-Wave Field Strength Between
150 and 1500 kHz," CCIR Recommendation453-3, Doc.
XIV Plenary Ass., ITU, Geneva, 1978.
CCIR, "Special Properties of the High Latitude Ionosphere
Affecting Radio-communications,"Study Program 25B/6,
Doc. 6/1005-E, 14 Oct 81, XVth Plenary Ass., ITU,
Geneva, 1981a.
CCIR,"Short-Term Forecasting of Operational Maximum
Usable Frequencies and Total Electron Content," Study
Program 27C/6, Doc. 6/1017-E, 14 Oct 81, XVth Plenary Ass., ITU, Geneva, 1981b.
10-102
10-103
CHAPTER 10
Forbes, J.M., "Upper Atmosphere Modifications Due to
Chronic Discharges of Water Vapor from Space Vehicle
Exhausts," in Space Systems and their Interactions with
Earth's Space Environment edited by H.B. Garrett and
C.P. Pike, American Institute of Aeronautics and Astronautics, New York, 1980.
Franklin, C.A. and M.A. MacLean, "The Design of Swept
Frequency Topside Sounders," Proc. IEEE, 57(6):
897-929, 1969.
Fremouw, E.J., "Geometrical Control of the Ratio of Intensity and Phase Scintillation Indices," J. Atmos. Terr.
Phys., 42: 775-782, 1980.
Fremouw, E.J. and J.F. Bates, "Worldwide Behavior of
Average VHF-UHF Scintillation," Radio Sci., 6:
863-869, 1971.
Fremouw, E.J. and C.L. Rino, "Modelling of Transionospheric Radio Propagation," RADC-TR-76-35, Rome
Air Development Center, Griffiss AFB, NY, 1976.
Fremouw, E.J., C.L. Rino, A.R. Hessing, and V.E. Hatfield, "A Transionospheric Communication Channel
Model," Quart. Tech. Rep. 7, SRI International (Stanford Research Institute), Menlo Park, Calif., 1977.
Fremouw, E.J. and C.L. Rino, "A Signal-Statistical and
MorphologicalModel of Ionospheric Scintillation," Proc.
AGARD Conf. Operational Modelling of the Aerospace Propagation Environment, Ottawa, Canada,
1978.
Fremouw, E.J., R.L. Leadabrand, R.C. Livingston, M.D.
Cousins, C.L. Rino, B.C. Fair, and R.A. Long, "Early
Results from the DNA Wideband Satellite ExperimentComplex Signal Scintillation," Radio Sci., 13: 167-187,
1978.
Fremouw, E.J., R.C. Livingston, and D.A. Miller, "On the
Statisticsof ScintillatingSignals," J. Atmos. Terr. Phys.,
43: 717-731, 1980.
Fremouw, E.J. and J.M. Lansinger, "Recent High Latitude
Improvement in a Computer Based Scintillation Model,"
Proc. Ionospheric Effects Symposium, 3rd, 1981, Alexandria, Va., edited by J.M. Goodman, Naval Research Lab., Washington, D.C., 1982.
Frisius, J., A. Ehmert, and D. Stratmann, "Effects of High
Altitude Nuclear Tests on VLF-Propagation," J. Atmos.
Terr. Phys., 26: 251-262, 1964.
Galejs, J., Terrestrial Propagation of Long Electromagnetic
Waves, Pergamon Press, Oxford, UK, 1972.
Gallagher, P.B. and R.A. Barnes, "Radio-Frequency Backscatter of Artificial Electron Clouds," J. Geophys. Res.,
68: 2987-3010, 1963.
Gautier, T.N. and D.H. Zacharisen, "Use of Space and
Time Correlations in Short-Term Ionospheric Predictions," Conference Record, 1st Annual IEEE Communications Convention, 671-676, 1965.
Goldman, M.V., "Field Aligned Instability Due to Stimulated Scattering of Intense Radio Waves from Diffusion
Quasimodes," Radio Sci., 9: 1077, 1974.
Goodman, J.M., "A Survey of Ionospheric Models. A Preliminary Report on the Development of an Ionospheric
Model Thesaurus and Users Guide," NRL Memorandum
Report 4830, Naval Research Lab., Washington, D.C.,
1982.
Gordon, W.E., "Incoherent Scattering of Radio Waves by
10-104
10-105
CHAPTER 10
AGARD Conf. Proc., 305: 28-1-28-17, ADA 113969
1982.
Kofsky, I.,"Further Evaluationsof Infrared SimulationData,"
Tech. Rept. DNA Contract 001-81-C-003, Photo Metrics, Inc., Woburn, Mass., 1981.
Kohnlein, W., "Electron Density Models of the Ionosphere," Rev. Geophys. Space Phys., 16: 341, 1978.
Kossey, P.A. and E.A. Lewis, "Detection and Classification
of Solar X-Ray Flares Using VLF Phase and Amplitude
Information," AFCRL-TR-74-0468, ADA003383 1974.
Kossey, P.A., E.A. Lewis, and E.C. Field, "Relative Characteristicsof TE/TM Waves Excited by Airborne VLF/LF
Transmitters," in Medium, Long, and Very Long Wave
Propagation, AGARD Conf. Proc., 305: 19-1-19-10
ADA 113969 1982.
Kossey, P.A., J.E. Rasmussen, and E.A. Lewis, "VLF
Pulse Ionosounder Measurements of the Reflection
Properties of the Lower Ionosphere," COSPAR, Lower
Ionosphere Structure, 133-138, Akademie Verlag, Berlin, 1974.
Kossey, P.A., J.P. Turtle, R.P. Pagliarulo, W.I. Klemetti,
and J.E. Rasmussen, "VLF Reflection Properties of the
Normal and Disturbed Polar Ionosphere in Northern
Greenland," Radio Sci., 18: 907-916, 1983.
Koster, J.R., "Some Measurements of the Irregularities Giving Rise to Radio-Star Scintillations at the Equator," J.
Geophys. Res., 68: 2579-2590, 1963.
Kouba, J. "A Review of Geodetic and Geodynamic Satellite
Doppler Positioning, Rev. Geophys. and Space Phys.
21:27-40, 1983.
Kouris, S.S. and L.M. Muggleton, "Diurnal Variation in
the E Layer Ionization," .. Atmos. Terr. Phys., 35:
133-139, 1973a.
Kouris, S.S. and L.M. Muggleton, "A Proposed Prediction
Method for Monthly Median foE," CCIR Report 2522, Contrib. No. 6/3/07 to Interim Working Party 6/3,
1973b.
Krause, G.E., R.J. D'Accardi, and E.L. Roswell III, "Field
Test of a Near-Real-TimeIonosphericForecasting Scheme
(60 km)," U.S. Army ECOM Rept. 3345, 1970.
Krause, G.E., R.J. D'Accardi, and E.L. Roswell 111,"Field
Test of a Near-Real-TimeIonosphericForecasting Scheme
(200 km)," U.S. Army ECOM Rept. 4144, 1973a.
Krause, G.E., R.J. D'Accardi, and E.L. Roswell III, "Field
Test of a Near-Real-TimeIonosphericForecasting Scheme
(500 km)," U.S. Army ECOM Rept. 4145, 1973b.
Lakshmi, D.R., S. Aggarwal, P.K. Pasricha, and B.M.
Reddy, "HF Communication Problems at Low Latitudes
due to Steep Spatial and Temporal Gradients," STPP,
4: D2-58 (see R.F. Donnelly, ed.), 1980.
Larson, A.C. and J.C. Hodges, "VHF/UHF Auroral Radar
Measurements," Stanford Research Institute International, Calif., June 1967.
Leadabrand, R.L., A.G. Larson, and J.C. Hodges, "Preliminary Results on the Wavelength Dependence and
Aspect Sensitivity of Radar Auroral Echoes Between 50
and 3000 MHz," .1. Geophys. Res., 72: 3877-3887,
1967.
Leadabrand, R.L., J.C. Schlobohm, and M.J. Baron, "Simultaneous Very High Frequency and Ultra High Fre-
10-106
quency Observations of the Aurora at Fraserburgh, Scotland," J. Geophys. Res., 70: 4235-4284, 1965.
Lee, M.C., A. DasGupta, J.A. Klobuchar, S. Basu, and
S. Basu, "Depolarization of VHF Geostationary Satellite
Signals near the Equatorial Anomaly Crest," Radio Sci.,
17: 399, 1982.
Leftin, M., "Numerical Representation of Monthly Median
Critical Frequencies of the Regular E Region (FOE),"
OT Rept. 76-88k, U.S. Office of Telecommunications,
Boulder, Colo., 1976.
Leftin, M., S.M. Ostrow, and C. Preston, "Numerical Maps
of foEs for Solar Cycle Minimum and Maximum,"
ESSA Tech Rept. ERL 73-ITS-63, Boulder, Colo.,
1968.
Leiphart, J.P., R.W. Zeek, L.S. Bearce, and E. Toth, "Penetration of the Ionosphere by Very-Low-Frequency Radio Signals-Interim Results of the LOFTI I Experiment,"
Proc. IRE, 50: 6-17, 1962.
Lewis, E.A. and P.A. Kossey, "POWERFLUX I: A Method
of Estimating Wave Intensities at Large Distances from
Ground-Based Low Frequency Transmitters," AFCRLTR-75-0338, ADA014787, 1975.
Lewis, E.A., J.E. Rasmussen, and P.A. Kossey, "Measurements of Ionospheric Reflectivity from 6 to 35 kHz,"
J. Geophys. Res., 78: 3903-3911, 1973.
Lewis, R.L., "The Wave-Hop fields for an Inclined Dipole
Over a Spherical Earth with an Anisotropic Ionosphere,"
Tech Rept. OT/ITS RR 5, Institute for Telecommunication Sciences, Boulder, Colo., 1970.
Lied, F., (ed.), H.F. Radio Communications, Technivision,
Maidenhead, UK, 1967.
Lincoln, J.V. and R.O. Conkright (eds.), "International
Reference Ionosphere-IRI 79," Report UAG-82, World
Data Center A for Solar-Terrestrial Physics, Boulder,
Colo., 1981.
Llewellyn, S.K. and R.B. Bent, "Documentation and Description of the Bent Ionospheric Model," AFCRL-TR73-0657, AD 772733, 1973.
Lloyd, J.S., G.W. Haydon, D.L. Lucas, and L.R. Teters,
Estimating the Performance of Telecommunication Systems Using the Ionospheric Transmission Channel, National Telecommunications and Information Administration, Boulder, Colo., 1978.
Lyakhova, L.N., "Short-Range Forecasting of the State of
the Ionosphere," Trudy IZMIRAN, 17: 240, 1960.
Lyakhova, L.N. and L.I. Kostina, "On Quantitative Ionospheric Forecasting," Geomag. Aeron., 13: 50, 1973.
Maeda, K.I., "Study on Electron Density in the Lower
Ionosphere,"J. Geomag. Geoelect., 23: 133-159, 1971.
Malik, C. and J. Aarons, "A Study of Auroral Echoes at
19.4 Megacycles per Second," J. Geophys. Res., 69:
2731-2736, 1964.
Marmo, F.F., L.M. Ashcenbrand, and J. Pressman, "Artificial Electron Clouds-1, Planet," Space Sci., 1: 277,
1959.
Martin, E. and J. Aarons, "F-Layer Scintillation and the
Aurora," J. Geophys. Res., 82: 2717, 1977.
McNamara, L.F., "Ionospheric Predictions on Transequatorial Circuits," Proc. IREE (Aust.), 35: 17, 1974a.
McNamara, L.F., "Unusual Nighttime Propagation Modes
Cape Circuits,"
Iono-
Millington, G., "Ground-Wave Propagation Over an Inhomogeneous Smooth Earth," Proc. IEE, 96, Part III:
53-64, 1949.
Millman, G.H. Modern Radar, ed. R.S. Berkowitz, Wiley,
New York, 1965.
Millman, G.H., "An Evaluation of HF Ionospheric Backscatter Echoes," G.E. Tech. Information Series Rept.
No. R75EMH19, Nov 1975.
Millman, G.H. and C.D. Bell, "Ionospheric Dispersion of
an FM Electromagnetic
Naval Research
10-107
CHAPTER 10
VLF Fields Excited by Dipoles of Arbitrary Orientation
and Elevation," Radio Sci., 35: 1445-1452,
1970.
Radio Communica-
Sounder Circuit,"
K., International
Reference
Ionosphere-IRI-79,
CHAPTER 10
at AFGWC," STPP, 1: 355 (see R.F. Donnelly, ed.),
1979.
Thomson, J.J., Conduction of Electricity Through Gases,
p. 321, Cambridge University Press, London, 1906.
Thrane, E.V. (Chairman), "D-Region Predictions," STPP,
2: 573 (see R.F. Donnelly, ed.), 1979.
Tinsley, B.A., A.B. Christensen, J.A. Bittencourt, P.D.
Angreji, and H. Takahashi, "Excitation of Oxygen Permitted Emissions in the Tropical Nightglow," J. Geophys.Res., 78: 1174, 1973.
Titheridge, J.E., "Periodic Disturbances in the Ionosphere,"
J. Geophys. Res., 73: 243-252, 1968.
Trizna, D.B. and J.M. Headrick, "Ionospheric Effects on
HF Over-the Horizon Radar," in Effects of the Ionosphere on Radiowave Systems, Proc. Ionospheric Effects
Symposium 1981,Alexandria, Va. edited by J.M. Goodman, Naval Research Lab., Washington, D.C., 1982.
Tsunoda, R.T., "On the Spatial Relationship of 1-mEquatorial Spread F Irregularities and Plasma Bubbles," J.
Geophys. Res. 85: 185, 1980.
Tsunoda, R.T., R.C. Livingston, J.P. McClure, and W.B.
Hanson, "Equatorial Plasma Bubbles" Vertically Elongated Wedges from the Bottomside F-Layer," J. Geophys. Res., 387: 9171, 1982.
Turner, J.F., "The Development of the Ionospheric Index
T," Ionospheric Prediction Science Series R Reports,
IPS-R12, Sydney, Australia, 1968.
Turner, J.F. and P.J. Wilkinson, "A Weekly Ionospheric
Index," STPP, 1: 279 (edited by R.F. Donnelly), 1979.
UAG Report 10, Atlas of Ionograms, (avail. World Data
Center A, U.S. Dept. of Commerce, NOAA/EDIS,
Boulder, Colo.), 1970.
Probing in Geophysics,
UAG Report 23a, URSI Handook of Ionogram Interpretation and Reduction, Revision of Chapters 1-4, (avail.
and Reduction,
Waves
per Second to 500 Kilocycles," NBS Circular 574, National Bureau of Standards, Washington, D.C., U.S.
Government Printing Office, 1956.
Winkler, C., "Radio Wave Guidance at VHF Through Equatorial Plasma Bubbles," J. Atmos. Terr. Phys., 43: 307,
1981.
Wright, J.W., "Kinesonde Observations of Ionosphere Modification by Intense Electromagnetic Fields from Platteville, Colorado," J. Geophys. Res., 78: 5622, 1973.
Wright, J.W. and A.K. Paul, "Toward Global Monitoring
of the Ionosphere in Real Time by a Modern Ionosonde
Network: The Geophysical Requirements and Technological Opportunity," NOAA Special Report, Space Environment Lab., July 1981.
Wright, J.W. and M.L.V. Pitteway, "Real-Time Data Acquisition and Interpretation Capabilities of the Dynasonde, 1. Data Acquisition and Real Time Display, 2.
Determination of Magnetosonic Mode and Echolocation
using a Small Spaced Receiving Array," Radio Sci., 14,
815-835, 1979.
Wright, J.W. and M.L.V. Pitteway, "Data Processing for
the Dynasonde: The Dopplionogram,"J. Geophys. Res.,
87: 1589, 1982.
Yeh, K.C., "Travelling Ionospheric Disturbances as a Diagnostic Tool for Thermospheric Dynamics," J. Geophys. Res., 77: 709-719, 1972.
Yeh, K.C. and C.H. Liu, "Ionospheric Effects on Radio
Communication and Ranging Pulses," IEEE Trans. Antennas Propag., AP-27(6), Nov 1979.
Yeh, K.C. and C.H. Liu, "Radio Wave Scintillations in the
Ionosphere," Proc. IEEE, 70, 1982.
Yeh, K.C., H. Soicher, and C.H. Liu, "Observations of
Equatorial Ionospheric Bubbles by the Radio Propagation Method," J. Geophys. Res., 84: 6589, 1979.
Zacharisen, D.M., "Space-Time Correlation Coefficients
for use in Short-Term Ionospheric Predictions," NBS
Report 8811, 1965.
Zevakina, R.A. E.V., Lavrova, and L.N. Lyakhova, Manual on Short-Term Predictions of Ionospheric Geomagnetic Storms and Radio-Propagation Forecasting, (text
10-111