11.3 Critical Frequency, Maximum Usable Frequency, and HF Propagation Paths

11.3 Critical Frequency, Maximum Usable Frequency, and HF Propagation Paths

Radio waves travelling through the ionosphere do not "bounce" off a sharp metallic ceiling; rather, they undergo continuous refraction—a gradual, progressive bending within an electrically charged plasma. By understanding the mathematical laws of plasma refraction, the geometry of oblique incidence, and the solar and geomagnetic disturbances that perturb the upper atmosphere, radio amateurs can accurately predict band openings, select optimal operating frequencies, and exploit complex propagation paths across the globe.

---

Physics of Ionospheric Refraction and Snell's Law

The ionosphere behaves as a dispersive, non-magnetic dielectric plasma. When an electromagnetic wave of angular frequency $\omega = 2\pi f$ traverses this plasma, its alternating electric field ($E(t) = E_0 e^{j\omega t}$) exerts an oscillating Lorentz force on the free electrons. Because an electron has a microscopic mass ($m_e \approx 9.109 \times 10^{-31}\text{ kg}$) compared to a heavy positive ion, it accelerates rapidly in response to the electric field:

$$m_e \frac{d^2 x}{dt^2} = -e E_0 e^{j\omega t} \implies x(t) = \frac{e E(t)}{m_e \omega^2}$$

This oscillating displacement of charge produces an internal macroscopic polarization density ($\vec{P} = -N_e e \vec{x}$), which opposes the applied electric displacement field:

$$\vec{D} = \varepsilon_0 \vec{E} + \vec{P} = \varepsilon_0 \left( 1 - \frac{N_e e^2}{\varepsilon_0 m_e \omega^2} \right) \vec{E}$$

The Plasma Frequency ($f_p$)

The natural resonant frequency of the electron plasma is defined as the plasma frequency ($f_p$, or Langmuir frequency):

$$f_p = \frac{1}{2\pi} \sqrt{\frac{N_e e^2}{\varepsilon_0 m_e}}$$

Substituting physical constants—electron charge $e = 1.602 \times 10^{-19}\text{ C}$, electron mass $m_e = 9.109 \times 10^{-31}\text{ kg}$, and permittivity $\varepsilon_0 = 8.854 \times 10^{-12}\text{ F/m}$:

$$f_p \approx \sqrt{80.6 \times N_e} \approx 8.98 \sqrt{N_e} \approx 9 \sqrt{N_e} \quad (\text{Hz})$$

where $N_e$ is the free electron density in $\text{electrons per cubic meter}$ ($\text{m}^{-3}$).

The Plasma Refractive Index ($n$)

The relative dielectric permittivity ($\varepsilon_r$) and refractive index ($n = \sqrt{\varepsilon_r}$) of a collisionless, unmagnetised plasma are:

$$n = \sqrt{1 - \frac{f_p^2}{f^2}} = \sqrt{1 - \frac{80.6 N_e}{f^2}}$$

                 Local angle definitions: flat-layer sketch
                              normal
                                |  / ray direction
                                | /
                                |/ theta_i (from normal)
Base of layer ------------------o-------------------------
                               /
                              /
                             /
Ground ------------------TX o---------------------------->
                            takeoff angle Delta above horizon

theta_i is measured at layer entry, not from the ground horizon.
As electron density rises, the ray bends away from the local normal.

Applying Snell's Law of Refraction at the boundary of successive plasma layers of increasing electron density:

$$n_1 \sin \theta_1 = n_2 \sin \theta_2$$

As the radio wave enters higher altitudes where electron density ($N_e$) increases, the refractive index ($n$) drops below $1.0$. Because $n_2 < n_1$, the refracted angle must increase ($\sin \theta_2 > \sin \theta_1$). The wave bends away from the normal (closer to the horizontal). If the electron density is sufficiently high or the incidence angle is sufficiently oblique, the wave reaches the horizontal ($\theta = 90^\circ$, $\sin \theta = 1.0$) and bends continuously downward back toward the Earth—appearing to undergo total reflection.

---

Critical Frequency ($f_c$) and the Secant Law (MUF)

Critical Frequency ($f_c$ or $f_oF_2$)

The Critical Frequency ($f_c$, designated $f_oF_2$ for the ordinary wave of the F2 layer) is defined as the highest frequency that will be returned to Earth when a radio wave is transmitted at vertical incidence (launched straight up at a takeoff angle of $\Delta = 90^\circ$, corresponding to an angle of incidence $\theta_i = 0^\circ$ relative to the normal).

At vertical incidence, the wave will penetrate the layer unless the refractive index drops to zero ($n = 0$) at the apex of its trajectory. Setting $n = 0$ in the refractive index formula yields:

$$\sqrt{1 - \frac{f_c^2}{f_c^2}} = 0 \implies f_c = f_p(\text{max}) \approx 9 \sqrt{N_{\text{max}}}$$

• If $f \le f_c$, the vertically directed wave is reflected back to Earth.
• If $f > f_c$, the refractive index never reaches zero; the wave penetrates the layer and escapes into outer space.

Critical frequencies are measured continuously worldwide by automated pulse transmitters called ionosondes, which sweep from $1\text{ to }20\text{ MHz}$ to produce real-time ionograms.

Oblique Incidence and the Secant Law

When a radio wave strikes the ionosphere at an oblique angle of incidence ($\theta_i > 0^\circ$, launched at a low antenna takeoff angle $\Delta$), it requires far less refractive bending to be returned to Earth than a vertically directed wave.

At the apex of an oblique ray's path, the wave travels horizontally: $\theta_r = 90^\circ$ and $\sin \theta_r = 1.0$. Applying Snell's Law between free space ($n = 1$) at the base of the layer and the reflection apex ($n$):

$$1 \times \sin \theta_i = n \sin(90^\circ) = n$$

Squaring both sides and substituting the plasma refractive index ($n^2 = 1 - f_c^2 / f^2$):

$$\sin^2 \theta_i = 1 - \frac{f_c^2}{f^2}$$

Using the trigonometric identity $\cos^2 \theta_i = 1 - \sin^2 \theta_i$:

$$\cos^2 \theta_i = \frac{f_c^2}{f^2} \implies f = \frac{f_c}{\cos \theta_i} = f_c \sec \theta_i$$

This fundamental formulation is known as the Secant Law for the Maximum Usable Frequency (MUF):

$$\text{MUF} = \frac{f_c}{\cos \theta_i} = f_c \sec \theta_i$$

Because the secant of an oblique angle is always greater than unity ($\sec \theta_i > 1.0$), the Maximum Usable Frequency for an oblique path is always significantly higher than the vertical critical frequency ($f_c$).

The Geometric MUF Factor ($M$-Factor)

Because the Earth and the ionosphere are concentric spheres, the maximum angle of incidence ($\theta_i$) achievable for an F2 ray launched tangent to the Earth's horizon (takeoff angle $\Delta = 0^\circ$) is geometrically limited to approximately $71^\circ\text{ to }74^\circ$.

Evaluating the secant at this grazing angle:

$$\sec(70.5^\circ) = \frac{1}{\cos(70.5^\circ)} \approx 3.0$$

The ratio of the MUF over a standard $3000\text{ km}$ path to the vertical critical frequency is termed the oblique factor or $M$-factor ($M(3000)F_2$):

$$\text{MUF}(3000\text{ km}) = f_oF_2 \times M(3000)F_2 \approx 3.0 \times f_oF_2$$

Step-by-Step Worked Calculation: MUF and Oblique Incidence

Problem: A local ionosonde station records an F2-layer critical frequency of $f_oF_2 = 7.20\text{ MHz}$. An amateur operator wishes to contact a station at a distance corresponding to an ionospheric angle of incidence of $\theta_i = 62.0^\circ$.

1. Calculate the Maximum Usable Frequency (MUF) for this circuit.
2. Assess which listed bands meet this model’s upper refraction limit, then identify what else is needed to establish usable communication.

Step 1: Compute the MUF via the Secant Law:

$$\cos(62.0^\circ) \approx 0.4695 \implies \sec(62.0^\circ) = \frac{1}{0.4695} \approx 2.130$$ $$\text{MUF} = f_c \sec \theta_i = 7.20\text{ MHz} \times 2.130 \approx 15.34\text{ MHz}$$

Step 2: Assess the Upper Refraction Limit:

• $7.0\text{ MHz}$: Below $15.34\text{ MHz}$, so refraction is possible in the stated model. Absorption, noise and the required radiation angle must still permit reception.
• $14.0\text{ MHz}$: Also below the calculated MUF and potentially useful, provided it is above the circuit's LUF and the antenna/path support that mode.
• $18.1\text{ MHz}$ and $21.0\text{ MHz}$: Above the calculated MUF, so this modelled F2 path does not return them.

Thus MUF alone cannot establish that either lower band is open. A usable frequency must also satisfy the receiving system's SNR requirement.

---

Optimum Working Frequency (OWF) and Lowest Usable Frequency (LUF)

Optimum Working Frequency (OWF / FOT)

Electron density and the upper supported frequency vary with time, location and disturbances. Operation very near the upper limit is vulnerable to penetration fades when that limit falls.

For reliable circuit planning, a predicted Optimum Working Frequency (OWF/FOT) is commonly the lower-decile MUF: the upper refraction limit is expected to support that frequency on about $90\%$ of days at the specified time. A familiar amateur-exam approximation is:

$$\text{OWF} \approx 0.85 \times \text{MUF}$$

This rule gives a $15\%$ frequency margin, rather than a guaranteed maximum signal or immunity to fading. A frequency near the upper part of the usable window often has lower D-region absorption, but it must remain above the LUF and below the supported upper limit. Use path forecasts and current conditions when choosing the operating frequency.

Lowest Usable Frequency (LUF)

While the upper limit of the HF communication window is governed by refraction mechanics (MUF), the lower limit is governed by absorption and noise. This lower boundary is the Lowest Usable Frequency (LUF).

Non-deviative absorption in the D-layer scales inversely with frequency squared ($A \propto 1/f^2$). As an operator shifts to lower frequencies, signal attenuation increases sharply. The LUF is defined as the frequency below which the received signal power drops below the minimum signal-to-noise ratio (SNR) required for intelligible demodulation.

Unlike the MUF (which is dictated solely by ionospheric physics), the LUF is system-dependent:

• Increasing transmitter output power (e.g. from $100\text{ W}$ to $1500\text{ W}$) lowers the LUF.
• Utilizing high-gain directional beam antennas (such as Yagis) lowers the LUF.
• A suitable narrow-band or coded mode can lower the SNR needed for a contact. For example, FT8 occupies about $50\text{ Hz}$, but its quoted sensitivity near $-21\text{ dB}$ uses a $2500\text{ Hz}$ noise reference. Compare mode sensitivities using the same reference bandwidth; occupied bandwidth and reference noise bandwidth are different quantities.
• High ambient electromagnetic interference (EMI) or atmospheric thunderstorm static raises the LUF.

                      The HF Usable Propagation Window

   Signal Strength / Path Loss
        ^
        |           ========================
        |          /                        \
        |         /   USABLE HF PASSBAND     \   Rapid Penetration
        |        /   (LUF <= f <= MUF)        \  into Outer Space
        |       /                              \
        |      /          * OWF (0.85 MUF)      \
        +-----+----------------------------------+---------------------> Frequency
             LUF                                MUF
        (D-layer Absorption                    (Refraction Limit,
         & Noise Floor Limit)                   f_c * sec θ_i)

The Operational Window and Circuit Blackout

For any skywave path, reliable communication is strictly confined to the frequency passband between LUF and MUF:

$$\text{LUF} \le f_{\text{operating}} \le \text{MUF}$$

During intense space weather events, D-layer absorption can skyrocket while geomagnetic turbulence depresses the F2 critical frequency. Under these conditions, the LUF rises while the MUF falls. If $\text{LUF} > \text{MUF}$, the usable window closes entirely, resulting in a complete circuit blackout.

---

Skywave Path Geometry: Skip Distance, Skip Zone, and Multihop

               Coverage along a path at one operating frequency

TX o========= ground-wave coverage =========a.........b==== skywave
                                            a: ground-wave limit; b: first skywave return
   |<------------------------------------------------>|
   Distance TX to b = skip distance
                                            |<------->|
                                            a to b = skip zone

The skip zone exists when ground-wave and ordinary skywave coverage
do not meet at the selected frequency. Sketch is not to scale.

Core Geometric Definitions

1. Virtual Height ($h'$): When pulse sounders measure reflection altitude via time-of-flight ($\Delta t$), they assume the wave travels at the constant speed of light ($c$) along straight lines: $h' = c \Delta t / 2$. However, because the wave actually slows down ($v_g < c$) as it refracts through the dense plasma, the virtual height is always greater than the true physical reflection height ($h' > h_{\text{true}}$). A typical daytime F2 virtual height is $300\text{ to }350\text{ km}$.
2. Takeoff Angle (Elevation Angle $\Delta$): The angle above the local horizon at which an antenna launches its RF energy. Lower takeoff angles ($5^\circ\text{ to }15^\circ$) strike the ionosphere at more oblique angles (larger $\theta_i$), yielding higher MUFs and greater hop distances.
3. Skip Distance ($d_{\text{skip}}$): The shortest surface distance from the transmitter to an ordinary returned skywave at the selected frequency. If the frequency exceeds the layer's vertical critical frequency, rays above a critical elevation penetrate instead of returning. At a suitable lower frequency, near-vertical rays can return and there need be no skip zone.
4. Skip Zone: The region between the useful ground-wave limit and the nearest skywave return for that frequency. The wanted station may be absent by those ordinary modes there; lowering frequency can bring skywave coverage closer to the transmitter.

Single-Hop Geometric Limits

Due to the spherical curvature of the Earth and the finite altitude of the ionospheric layers, single-hop skywave propagation has strict theoretical distance limits:

• E-Layer ($h \approx 100\text{ km}$): The maximum single-hop distance for a wave launched tangent to the horizon ($\Delta = 0^\circ$) is approximately $2000\text{ km}$.
• F2-Layer ($h \approx 320\text{ km}$): The maximum single-hop distance for a tangent ray is approximately $4000\text{ km}$.

To bridge intercontinental distances exceeding $4000\text{ km}$—such as an Australia-to-Europe circuit ($16,000\text{ km}$)—the wave must propagate via multihop paths (e.g. 4 or 5 successive F2 hops, bouncing alternately between the ionosphere and the Earth's surface) or chordal-hop propagation (successive ionospheric refractions without an intervening ground reflection, reducing ground-reflection loss).

Near Vertical Incidence Skywave (NVIS)

While long-distance DX demands low takeoff angles, tactical regional communication ($0\text{ to }400\text{ km}$) requires an entirely different approach: Near Vertical Incidence Skywave (NVIS).

• Choose a frequency above the circuit's LUF and below the vertical critical frequency for the relevant region. The $3.5\text{ MHz}$ or $7\text{ MHz}$ bands are common choices when conditions support them. Use an antenna with strong high-angle radiation, typically around $70^\circ$–$90^\circ$ elevation.
• Returning high-angle rays can fill the ordinary skip zone and cover nearby and regional stations, often out to a few hundred kilometres. This can bypass terrain that obstructs a direct path, but usable coverage still requires adequate signal strength and a supporting ionosphere.

---

Fading Mechanisms and Space Weather Disturbances

Fading (QSB) Mechanisms

Signal fading—termed QSB in amateur radio Q-code—arises from three primary physical mechanisms:

1. Multipath Fading: The transmitted signal reaches the receiver over multiple concurrent paths of varying length (e.g. a 1-hop path interfering with a 2-hop path, or low-angle and high-angle Pedersen rays). As ionospheric winds shift path lengths by fractions of a wavelength, the arriving signals alternate between constructive and destructive phase interference, causing deep signal swings.
2. Selective Fading: Different frequency components within a single modulated signal experience different phase shifts across the multipath channels. In an AM or SSB signal, the carrier or individual audio sidebands may cancel independently, causing the audio to sound hollow, garbled, or severely distorted.
3. Polarisation Fading: Due to continuous variations in Faraday rotation, the arriving wave's polarisation vector rotates across the receiving antenna's aperture, inducing periodic signal nulls as it passes through cross-polarisation.

Space Weather Phenomena: SIDs and Geomagnetic Storms

                       Solar disturbance timing

Solar flare -> X-rays reach Earth in about 8.3 minutes
            -> increased daylight D-region absorption
            -> severity ranges from degradation to HF blackout

Earth-directed CME -> plasma and magnetic field arrive later
                   -> roughly 15 hours for the fastest, days for slower events
                   -> possible geomagnetic storm
                   -> path-dependent MUF changes and auroral effects

Photon travel time is not an 8-minute warning after a flare is seen
at Earth: its light and X-rays arrive at approximately the same time.

1. Sudden Ionospheric Disturbances (SID / Mögel-Dellinger Effect)

• Cause: A solar flare emits enhanced X-rays; the $0.1$–$0.8\text{ nm}$ band is used in operational flare classification.
• Arrival Time: Travelling at the speed of light, the X-ray blast strikes Earth in $8.3\text{ minutes}$.
• Mechanism: Hard X-rays penetrate deeply down to the D-region ($60\text{ to }90\text{ km}$), producing catastrophic ionization and driving electron densities up by orders of magnitude.
• Impact: Increased D-region absorption degrades HF skywave on sunlit paths. Effects range from reduced signal strength or limited outages to a widespread blackout in strong events. Lower frequencies are usually affected more and recover later. Duration and geographic extent depend on the flare, path and frequency; not every SID blacks out the entire dayside.

2. Coronal Mass Ejections (CMEs) and Geomagnetic Storms

• Cause: A massive bubble of billions of tons of magnetized solar plasma (protons and electrons) is blasted outward from the Sun.
• Arrival Time: Only an Earth-directed CME can directly affect Earth's environment. The fastest can arrive in roughly $15$–$18\text{ hours}$; slower events take several days. A fixed $20$–$48\text{-hour}$ interval does not cover the range.
• Mechanism: The arriving plasma and embedded magnetic field can disturb the magnetosphere and upper atmosphere. Changes in composition, circulation and electric fields alter F-region electron density and therefore HF propagation.
• Impact on Propagation:
  • MUF Changes: A storm can initially enhance the supported frequency on some paths and later depress it; other paths show a decrease directly. A depressed MUF can close a previously usable high-HF circuit, but no fixed percentage drop or worldwide all-band outcome follows from the event label alone.
  • Auroral Flutter: Signals passing across polar or high-latitude paths scatter off turbulent auroral plasma columns, introducing rapid multi-frequency Doppler shifts that convert clear CW signals into raspy, hissing noise.

3. Space Weather Forecasting Indices

Amateur operators evaluate real-time space weather data using two primary geomagnetic indices:

• $K$-Index and Planetary $K_p$-Index: A quasi-logarithmic index from $0\text{ to }9$ measuring geomagnetic field fluctuations over a $3\text{-hour}$ synoptic interval:
  • $K_p=0$–$3$: Low to unsettled geomagnetic activity. Other factors still determine which HF paths are usable.
  • $K_p=4$: Active conditions; increased variability may affect propagation.
  • $K_p=5$: Minor geomagnetic storm (G1), with possible effects especially at higher latitudes.
  • $K_p=6$: Moderate storm (G2); high-latitude HF can fade.
  • $K_p=7$, $8$ and $9$: Strong, severe and extreme activity (G3–G5). HF may become intermittent or unusable in affected areas; these indices do not guarantee a global blackout.
• $A$-Index and Planetary $A_p$-Index: A linear scale from $0\text{ to }400$ representing the 24-hour daily average of geomagnetic activity derived from the eight 3-hour $K$-indices. $A_p < 7$ denotes quiet conditions, whereas $A_p > 30$ warns of significant geomagnetic storm activity.