Ionospheric Physics of Radio Wave Propagation

Edwin C. Jones, MD, PhD (AE4TM)
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996-1200

Introduction

A basic understanding of the ionosphere properties is paramount to understanding HF radio propagation over distance. This page is dedicated to providing a simplified overview for some of the important ionosphere mechanics used in propagation research as well as discussing limitations due to the complex dynamics of the ionosphere. Because the physics of ionospheric radio wave propagation is an extremely vast and complex topic, to describe the physics of ionospheric radio wave refraction in a short and concise manner without overwhealming the reader requires that certain assumptions be made. Although these assumptions may appear to oversimplify the problem at first, the predictions made by the final solution of the Equation of Motion do adequately describe the observations made by radio wave propagation experiments; therefore, the following physics has formed the basis of understanding of this problem in research groups around the world. In my summary, I have included supporting data for these conclusions that were obtained by recent HF radio testing using the PTC-II HF controller to derive HF radio skip distance over a range of HF frequencies. Finally, in the physics below, Gaussian units will be utilized and all vectors will be denoted in bold.

In order to properly view the mathematical equations, you will need to keep the fonts setting in your preference menu set to the document-specified fonts position. Furthermore, you must have the Symbols and Map Symbols fonts installed to your FONT directory, i.e. WINDOWS/FONT directory for Windows users. If you are unable to properly view the equations simply save the following two links to your Operating System FONT directory:

Symbols & Map Symbols.

Simplified Model of Propagation in the Ionosphere

The simplest approach to describing radio wave propagation is to solve for the index of refraction h = (m e)1/2, where m = magnetic permeability (1.25664 x 10-6 H m-1) and e = dielectric constant. The index of refraction, in turn, describes the relationship between the angles of incidence and refraction through Snell's Law

Sin i / Sin r = h' / h.

This is shown graphically in the figure below which shows an incident wave k striking a plane interface between different media, giving rise to a reflected wave k" and a refracted wave k'.

Refracted and Reflected Waves

Since the magnetic permeability is constant in the ionosphere, our goal is to now solve for the dielectric constant from the Equation of Motion, where the dielectric constant is defined as the ratio of the strength of an electric field in a vacuum to that in the ionosphere.

We now consider a simple problem of a tenuous electron plasma of uniform density trapped in a strong, static, and uniform magnetic induction Bo. If we assume that the transverse radio waves propagate parallel to Bo, the Equation of Motion for electrons trapped in this ionospheric plasma is given by

m d2x/dt2 - e/c Bo X dx/dt = -e E e-iwt ,

where the influence of the B field of the transverse wave has been neglected compared to the static induction Bo and the electron charge is given by -e. It is now customary to describe the electric field component of the radio waves as circularly polarized which implies

E = (e1 +_ ie2) E

and a similar expression for x. Since Bo is orthogonal to both e 1 and e2, the cross product in our Equation of Motion has components only in the directions e1 and e2; therefore, the transverse components decouple. This leads to a steady-state solution given by

x = e E / m w (w _+ wB) ,

where wB is the frequency of precession of a charged particle in a magnetic field which is given by

wB = e Bo / mc ~ 6 x 106 sec-1 (in the earth's magnetic field of Bo=0.3 gauss).

The frequency dependence of our steady-state solution can be determined by transforming our Equation of Motion to a coordinate system precessing with frequency wB about the direction of Bo. If the static magnetic field is neglected, the force on the electrons has an effective frequency (w +_ wB), depending on the sign of the circular polarization.

The steady-state solution implies a dipole moment for each electron and yields, for a bulk sample, the dielectric constant of the ionosphere

e_+ = 1 - {wp2 / w(w _+ wB)},

where the upper sign corresponds to a positive helicity wave (left-handed circular polarization in optics notation), while the lower sign corresponds to negative helicity. Furthermore, w is the frequency of our radio wave of interest and wp is the plasma frequency of the ionosphere and it is given by

wp2 = 4p NZ e2 / m ,

where NZ is the density of electrons per unit volume. For propagation antiparallel to the magnetic field Bo, the signs are reversed. Furthermore, for propagation in directions other than (anti)parallel to the static field Bo, it is straight forward to show that, if terms of order wB2 are neglected compared to w2 and wwB, the dielectric constant is still given by e_+ above.

In this simplified problem of ionospheric radio wave propagation, we see the essential characteristic that waves of right-handed and left-handed circular polarizations propagate differently. In other words, the ionosphere has both birefringent and anisotropic scattering properties.

For the earth's ionosphere, the density of free electrons ranges between 104 and 106 electrons/cm3, corresponding to a plasma frequency on the order of wp ~ 6x106 - 6x107 sec-1. This along with the precession frequency wB given above implies a wide interval in frequency within the HF spectrum where one state of circular polarization cannot propagate within the medium at all; rather, this wave is totally reflected back towards the earth. The other state of polarization is partially transmitted. Thus, when a linearly polarized wave is incident on a plasma, the reflected wave will be elliptically polarized, with its major axis generally rotated away from the direction of polarization of the incident wave.

The propagation of HF radio waves off the ionosphere is explicable in terms of these ideas, but the presence of several layers of plasma with densities and relative positions varying with height and time makes this problem considerably more complicated than this simplified model implies. Therefore, most ionospheric research groups seek to understand the variation of electron density versus vertical height above the earth's surface. The electron densities at various heights can be inferred by studying the reflection of radio pulses transmitted vertically upwards (known as ionograms). Such studies have shown that the number n0 of free electrons per unit volume slowly increases with height within a given layer of the ionosphere where it reaches a maximum before falling off abruptly with further increases in height as shown in the figure below.

Vertical Radio Wave Reflections

A pulse of given frequency w1 enters a layer without reflection because of the slow change in n0. When the density n0 is large enough, wp(h1) ~ w1. At this point, the dielectric constant vanishes and the pulse is reflected vertically back to the earth. The electron density NZ required to induce this vertical reflection for a pulse of frequency w is found by combining the equations above to give

NZ = (m / 4pe2) (w2 _+ wwB).

Defining the minimum frequency at which a radio wave just penetrates a specific layer of ionization as the critical frequency fc for that layer, the above equation can be rearranged into a more useful form. Using MKS units and frequencies given in Hertz, we have two modes contributing to vertical reflections off a layer of ionization of electron density NZ. The critical frequencies of these two modes are given by the following:

fc = 8.98*(NZ)1/2 for the ordinary mode and

fc = 8.98*(NZ)1/2 + 0.5*Bo(e/m) for the extraordinary mode.

By measuring the time interval between the initial transmission and reception of the reflected signal, the height h1 corresponding to this electron density can be found. Therefore, by varying the frequency f and studying the change in time intervals (ionograms), the electron density as a function of height can be determined. If the frequency f is too high, the index of refraction remains finite and no vertical reflection occurs. Therefore, the minimum frequency at which the vertical reflections just disappear determines the maximum electron density in a given ionospheric layer. Shown below is an example of an ionogram typically recorded at various ionospheric research stations around the world.


An Ideal Ionogram
Click on image for additional information on interpreting ionograms

The maximum frequency that vertically reflects off the F2-layer of the ionosphere at a given time and place is known as foF2. The extraordinary mode reflecting off the F2-layer is known as fxF2. Similarly, the maximum frequencies that vertically reflect off the E-layer and F1-layer are known as foE and foF1, respectively. Real-time foF2 data are readily available from several government agencies in the US and Australia. The two plots below are real-time foF2 plots of the USA and Australia that are provided by the Australian Government.


Real Time foF2 Soundings in North America

Source: http://www.ips.gov.au/

Real Time foF2 Soundings in Australia

Source: http://www.ips.gov.au/

Unfortunately, vertical reflections typically disappear at frequencies higher than the ~30m range so that the more popular HF DX bands (10m-20m) cannot be adequately characterized by these vertical ionograms. For these higher frequencies, oblique ionograms are generally utilized to yield empiric data on the refractive properties of the ionosphere and the virtual altitudes of these reflections. However, such data are generally limited since the goal of most ionospheric research studies are aimed at determining the variation of electron density with height. The figure below illustrates the minimum ground distance needed (blind zone) to establish a digital communications link via reflection off the ionosphere for these higher frequencies. These data were recently obtained by connecting to my Nashville, TN pactor base station from my mobile pactor station during two recent trips to Key West, Florida. The data clearly show that as the frequency increases, the blind zone increases with no vertical reflections occurring at frequencies higher than ~7.5 MHz (SFI=179) or ~14MHz (SFI=261). Finally, the data indicate that the blind zone decreases with higher solar flux (SFI) and increases with frequency.

Skip Distance vs Frequency

Whistlers

The frequency dependence of the dielectric constant e-(w) at low frequencies is responsible for a peculiar magnetospheric propagation phenomenon known as "whistlers". As w->0, e-(w) tends to positive infinity as e-~wp2/wwB. Propagation occurs, but with a wave number

k_~(wp/c){w/wB}1/2 .

This corresponds to a highly dispersive medium. Recalling that energy transport is governed by the group velocity vg=[dw/dk]0 leads to a solution for the group velocity in the magnetosphere at the MF range as

vg(w)~2vp(w)~2c{(wwB)1/2/wp} .

This equation indicates that pulses of radiation at different frequencies travel at different speeds, in particular, the lower the frequency, the slower the speed. Whistlers occur when thunderstorms in one hemisphere generates a wide spectrum of frequencies, some of which propagate along the dipole lines of the earth's magnetic field described above with the higher frequency components reaching the antipodal point first and the lower ones later. These whistlers generally occur at frequencies below 100kHz and when received by a radio receiver sound like a whistle dropping in frequency. By assuming a distance of 104 km to a lightening discharge, the time scale for these whistlers occurs on the order of seconds. An example of a whistler can be heard at the following link: Whistler.wav file.


RF Skin Effect

The skin effect penetration depth is defined as the distance to which a radio wave can penetrate into a conductive medium (metal, salt water, ionosphere, etc.) leaving only 37% of its initial intensity. As predicted by Maxwell's Equations below, rf energy decays exponentially when it encounters a conductive medium.

From Maxwell's Equations,

, x H = [4pJ + dD/dt] / c

In an electron gas, E = J / s; therefore, this becomes

, x H = [4psE + dD/dt] / c,

where H is the magnetic field in matter, B is the magnetic induction, c is the velocity of light, J is the volume current, D is the electric field displacement, s is the electrical conductivity, and E is the electric field.

Faraday's law of induction states

, x E = - (dB/dt) / c

Combining the above two equations gives

, x (, x E) = ,(, . E) - ,2E = - ,2E = - [d(, x B)/dt] / c

,2Eo + (4pms / c2) dE / dt = 0 .

Recognizing that this is the diffusion equation, write

E(t) = Eo e-iwt
,2Eo + t2Eo = 0, then
t2 = (4psmw / c2) i
t = (1 + i) [2psmw]1/2 / c = (1 + i) / d .

So the skin depth d is given by

d = c / [2psmw]1/2 ,

where s is the electrical conductivity, m is the permeability, and w is the angular frequency.

In MKSA units, the skin depth d is given by

d = {2 / mws}1/2 .

Recognizing that the relative RF intensity at the skin depth is 37% (1/e) of the incident intensity and that this is equivalent to 8.68dB, the attenuation of radio waves in dB/ft can be solved for various conductivities s. In the figure below, the attenuation of sea water (s=4 Siemens) is shown in blue as a function of frequency. In addition, typical ionospheric conductivities(*) lie between 1x10-7Siemens<s<1x10-4Siemens; therefore, the attenuation per foot is shown as red and green for the upper (daytime) and lower (nighttime) conductivity limits for the ionosphere, respectively.


Attenuation of RF Passing Through Conductive Media

In sum, the figure above indicates that as the frequency increases the RF attenuation per unit length increases. Furthermore, the RF attenuation also increases with electrical conductivity s. The implications of these findings include the existence of discrete ranges of frequencies that can reflect off the ionosphere as well as the ability to communicate through sea water using extremely low frequencies (ELF).

_____________________________________________________________________________________

(*), From statistical mechanics, the electrical conductivity of a free electron gas is given by

s = ne2t / m ,

where e is the electron charge, n is the density of electrons, m is the electron mass, and t is the average time between electron collisions. Furthermore, the average time between electron collisions t can be calculated from

t = lmfp (m / 3kT)1/2,

where lmfp is the mean free path between electron collisions, k is the Boltzmann constant, and T is the absolute temperature in Kelvins.

Furthermore, the conductivity of the ionosphere depends of the angle between the incidence of rf with the earth's magnetic field. However, because this effect amounts to less than one order of magnitude difference between the perpendicular and parallel directions, this effect was ignored in the simplified model above.




For more information on radio wave propagation see:

The Basics of Radio Wave Propagation.


AE4TM Propagation Research Page


AE4TM Ed's HomePage


Government Ionospheric Research Sites


NGDC Ionospheric Physics Group

NGDC Real Time Ionograms

Earth's Conjugate Auroras

NOAA Space Environment Center

IPS Radio and Space Services

Propagation Forecast via hfradio.org

Background to Ionospheric Sounding

The Source of Energy for Auroras

Antarctica: Whistler Sounds




7/1999-6/2004 Edwin C. Jones





Counter