Upper atmosphere

The upper part of the terrestrial atmosphere includes regions where ionization makes the response to electromagnetic waves depend strongly on mobile Electrons. The ionosphere describes that electrical structure; its boundaries need not coincide with the boundaries defined by temperature in the mesosphere and thermosphere.

The ionosphere is a region of the upper atmosphere of Earth with substantial partial ionization. Its mobile Electrons and ions constitute a plasma mixed with neutral gas. For radio communication, the crucial question is whether a chosen radio wave returns toward the ground, is absorbed, or passes through. A skywave can connect stations beyond their direct view across the curved surface of Earth, while wave interference between a ground wave and a skywave can make reception fade even when the radio transmitter remains steady.

The ionosphere is useful both as a route for radio communication and as a remote laboratory: varying the probing frequency changes how deeply the radio wave penetrates the plasma. The Appleton-Barnett radio interference experiment turned an unwanted fluctuation in reception into a measurement of the overhead route. Cold-plasma radio dispersion explains why this is a tunable medium, while ionospheric virtual height explains why a measured delay is not automatically an altitude.

The E layer is a main region of the ionosphere below the F layer. The 1924 Appleton-Barnett radio interference experiment concerns this lower region.^[ref]

The F layer is a higher main region of the ionosphere, above the E layer. Its different electron density and persistence can permit skywave propagation when a probing radio wave passes through the lower E layer.

A ground wave is a radio wave propagating near and along the surface of Earth. In a short-baseline comparison with a skywave, it supplies a reference path. Approximating its distance by the station separation and its speed by the speed of light is a modeling choice; ground properties and frequency can also affect its attenuation and wave phase.

A skywave is a radio wave returned toward Earth through interaction with the ionosphere. In a smoothly varying plasma, the ray can bend continuously through refraction before returning: there need be no sharp material surface at the apparent point of reflection. Different frequencies, launch directions and wave polarizations can have different paths and losses.

On 11 December 1924, Edward Victor Appleton and Miles Aylmer Fulton Barnett perform an experiment: the Bournemouth broadcasting radio transmitter slowly changes wavelength; an Oxford radio receiver records alternating maxima and minima from ground wave-skywave wave interference. The inferred equivalent height is about $90\mathrm{km}$. A separate comparison of simultaneous loop and vertical antenna signals tests the downward arrival direction. The E layer is established first; the higher F layer is found in winter 1926–27. These dates describe different findings.^[ref]

Why change the wavelength deliberately? Ordinary fading alone could arise from changing emission, receiver drift or a moving propagation path. A controlled sweep predicts a specific progression of phase difference between two arrivals. The directional antenna comparison supplies an independent kind of evidence for an overhead contribution. Agreement supports an atmospheric return; fitting a distance alone would not establish a unique electron density profile. In the model below, comparable received amplitudes make the interference fringes easier to see. A practical repetition would also calibrate the frequency sweep, track radio transmitter output and receiver gain, and compare repeated sweeps with fixed-frequency intervals to detect time drift. These are controls inferred from the measurement model, rather than additional procedures attributed to the historical account.

The radio receiver does not need to resolve two pulses to detect two paths. If one path takes longer, the two contributions reaching the receiver during a slow frequency sweep carry slightly different emission frequencies and can produce a slow beat in the detected amplitude. At strictly fixed frequency with stationary paths, the arrivals instead have a fixed phase difference: there is no continuing beat merely because one path is longer. This distinction connects the historical observation to ionospheric interference fringe spacing.

Consider an original two-path model with a locally flat ground baseline $d$ and a symmetric equivalent return height $h$. Replace the skywave by two straight segments traveling at the speed of light $c$, and replace the ground wave by a straight segment of length $d$. Their lengths and positive excess distance areThis is an equivalent geometric construction, not a metallic layer overhead. It neglects the curvature of Earth, refraction within the return region, ground-induced wave phase shifts and optical dispersion.

For received complex amplitudes $A_g$ and $A_s$, absorb their constant relative wave phase into ${\varphi }_0$ and choose their magnitudes nonnegative. An ideal detector reading proportional to squared received field hasAssume $\Delta L$ and ${\varphi }_0$ are independent of frequency across the sweep, with slowly changing amplitudes. Consecutive maxima differ by $2\pi$ in phase difference, soThe often-used approximation $\Delta f\simeq c/\Delta L$ therefore carries a frequency-independent-path assumption. Counting $N$ complete fringe periods between $f_1$ and $f_2$ gives $\Delta L\simeq Nc/|f_2-f_1|$; equivalently, $N=\Delta L|1/{\lambda }_2-1/{\lambda }_1|$ when the vacuum wavelength is $\lambda =c/f$. Count whole periods, not maxima plus minima as though each were a whole period.

Choose $d=150\mathrm{km}$ and $h=100\mathrm{km}$, with $c\simeq 3.00\times 10^8\mathrm{m}{\mathrm{s}}^{-1}$. Then $L_s=250\mathrm{km}$, $\Delta L=100\mathrm{km}$ and $\Delta f=3.00\mathrm{kHz}$. A $12.0\mathrm{kHz}$ sweep covers four full fringe periods. Conversely, the inferred equivalent height isThese chosen values illustrate the inversion; they are not a transcription of the 1924 measurements. With $A_g=1$ and $A_s=0.7$, the normalized reading ranges from $(1-0.7{)}^2=0.09$ to $(1+0.7{)}^2=2.89$: incomplete cancellation does not remove the distance information.

For a sweep with linearly changing frequency $f(t)=f_0+\dot{f}t$ and two constant delays ${\tau }_s>{\tau }_g$, the received instantaneous frequency difference has magnitude $|\dot{f}|({\tau }_s-{\tau }_g)$. In this model ${\tau }_s-{\tau }_g=\Delta L/c$, giving a beat frequency $|\dot{f}|\Delta L/c$. Equivalently, traversing each $\Delta f$ of the sweep crosses one fringe. This calculation requires a sufficiently slow sweep that propagation, tuning and amplitude changes can be treated consistently; it is not a beat between two independent transmitters.

More generally, measure the total phase difference $\Phi (f)$, including propagation and any reflection or apparatus contribution. Define the differential group delayFor a locally monotone wave phase with slowly varying derivative, the local spacing is $\Delta f\simeq 1/|\Delta {\tau }_g|$. If $\Phi =2\pi fD(f)/c+{\varphi }_r(f)$, this becomesHere $D(f)$ is a phase-equivalent optical path difference, not necessarily the geometric excess distance $\Delta L$. A drifting ionosphere, several simultaneous skywaves, changing wave polarization, a frequency-dependent reflection wave phase or significant loss can invalidate the simple constant-spacing inversion.

Take a cold, unmagnetized, collisionless plasma with uniform equilibrium electron density $n_e$, fixed neutralizing ions, and a small electromagnetic wave that is a transverse wave. Cold means that thermal-motion corrections are neglected. For electron mass $m_e$, elementary charge magnitude $e$ and electric field $\mathbf{E}=\mathrm{Re}(\hat{\mathbf{E}}{e}^{-i\omega t})$, the linear Newton's second law equation is $m_e\dot{\mathbf{v}}=-e\mathbf{E}$. Thus $\hat{\mathbf{v}}=-ie\hat{\mathbf{E}}/(m_e\omega )$ and the current density is $\hat{\mathbf{J}}=i{n}_e{e}^2\hat{\mathbf{E}}/(m_e\omega )$. Writing the induced electric polarization as $\hat{\mathbf{J}}=-i\omega \hat{\mathbf{P}}$ gives a negative electric susceptibility. With vacuum permittivity ${ϵ}_0$,The plasma frequency ${\omega }_p$ controls the response. For $\omega >{\omega }_p$, the propagating branch has real refractive index $0<n<1$. For $\omega <{\omega }_p$, this homogeneous model has an evanescent wave, not an ordinary propagating wave with a negative real refractive index. In a slowly varying vertical profile, a normally incident ray approaches a turning region where ${\omega }_p(z)=\omega$; the geometrical-ray approximation fails right at that region. At oblique incidence the conserved horizontal wave vector changes the turning condition. A smooth density gradient can therefore return a skywave without a mirror or a layer of metal.

For a homogeneous propagating section, direct differentiation of the dispersion relation givesA phase velocity greater than $c$ does not send information faster than $c$. Here a narrow wave packet has group velocity $v_g$, while the full causal response must still respect special relativity. Even this stripped-down plasma is strongly dispersive near its plasma frequency. In the actual ionosphere, collisions produce absorption, and the magnetic field of Earth changes the response with direction and wave polarization; a single real scalar $n$ cannot capture all of that behavior.

For a vertical outgoing and returning radio wave, define the virtual height $h^{\prime }=c{\tau }_g/2$ from its round-trip group delay ${\tau }_g$, after removing instrument delays. It is the height a fictitious reflector would have if both legs traveled at $c$. It differs from both the real turning height $z_t$ and a height fitted to an oblique constant-path interference fringe model.

An original solvable profile makes the difference explicit. Approximate the region below $z_0$ by $n=1$, and prescribe a linear rise of squared plasma frequency above it:Here the symbol $u_t$ denotes the penetration depth, and $a$ has units ${\mathrm{s}}^{-2}{\mathrm{m}}^{-1}$. In the propagating part of the ramp, $n=\sqrt{1-(z-z_0)/u_t}$. In the leading geometrical-ray delay model, with a frequency-independent turning connection wave phase, the round-trip group delay isThe endpoint singularity is integrable: substitute $s=\sqrt{1-u/u_t}$ to obtain $2u_t$. This approximation does not resolve the finite-wavelength structure at the turning region. It illustrates the effect of slow group velocity before the return, rather than adding extra geometric altitude.

As an independent phase-based derivation, the propagation contribution to the round-trip wave phase isIts derivative $d\Psi /d\omega =2(z_0+2{\omega }^2/a)/c$ agrees with the delay integral. Using the phase integral itself as a distance at speed $c$ would give $z_0+2u_t/3$ instead of $h^{\prime }$: wave phase, group delay and real turning height are three different quantities.

Choose $z_0=70\mathrm{km}$ and stipulate $f_p(z)=3\mathrm{MHz}\sqrt{(z-z_0)/(30\mathrm{km})}$, where $f_p={\omega }_p/(2\pi )$. ThenAt $f=3\mathrm{MHz}$, the real turning height is $100\mathrm{km}$ but the virtual height is $130\mathrm{km}$, corresponding to ${\tau }_g\simeq 0.867\mathrm{ms}$ with $c\simeq 3.00\times 10^8\mathrm{m}{\mathrm{s}}^{-1}$. These stipulated parameters describe neither a measured profile nor a separate physical reflector at $130\mathrm{km}$. Because the ramp has no density maximum, it also has no finite upper penetration threshold; a real bounded E layer or F layer requires a richer profile.