MHD Wave Theory in the Magnetosphere A Tutorial for Ground Magnetometer Analysts

1. MHD Basics — Frozen-In Flux

In a highly conducting plasma (like the magnetosphere), magnetic field lines are "frozen in" to the plasma. Where the plasma goes, the field goes. This means disturbances propagate as waves on the magnetic field, not through the field. The magnetic field is the medium.

This frozen-in condition holds whenever the magnetic Reynolds number R_m = μ₀σVL ≫ 1, which is satisfied everywhere in the magnetosphere except at reconnection sites (Section 7). It means that plasma elements that start on the same field line stay on that field line as the system evolves.

Three Key Speeds

The dynamics of a magnetized plasma are controlled by three characteristic wave speeds:

The key insight: V_A varies enormously across the magnetosphere because B decreases as r⁻³ (dipole) while ρ increases with L-shell (more plasma at larger distances). This variation in V_A is what creates the rich wave structure you observe.

2. The Three MHD Wave Modes

Three fundamental wave types exist in a magnetized plasma. Understanding these is essential because each mode creates different ground signatures with different spatial coherence properties:

1. Alfvén Wave (Shear / Transverse)

Field lines bend transversely but do not compress. Like shaking a rope sideways. Propagates only along B (the background magnetic field). Speed = V_A. The wave carries a field-aligned current (FAC) that closes through the ionosphere — this is what your magnetometers detect as a magnetic perturbation in the N-S direction (the ΔB_X component from the Hall current).

2. Fast Magnetosonic Wave (Compressive)

Field lines AND plasma compress together. Propagates in all directions (isotropic in the limit V_A ≫ C_s). Speed = V_f ≈ V_A. This is the "sound wave" of MHD. Because it propagates isotropically, it can carry energy across field lines — crucial for coupling different L-shells and creating global signatures.

3. Slow Magnetosonic Wave (Compressive)

A compressive wave guided along B at speed below V_A. Less important for ground signatures in the Pc3–Pc5 bands because it is strongly damped in the magnetosphere. Included for completeness.

3. Alfvén Waves — The Plucked String

An Alfvén wave is a transverse oscillation of a magnetic field line, exactly like a vibrating guitar string. The physics is a direct analogy:

In the magnetosphere, field lines connect the northern and southern ionospheres. The ionospheres act as partial reflectors: the conductive ionosphere shorts out the wave's electric field, reflecting the Alfvén wave back along the field line. Standing waves form — just like standing waves on a guitar string with fixed ends.

The fundamental frequency of a field line is:

where L_path is the total path length from ionosphere to ionosphere along the field line. Higher harmonics: f_n = n · f₁.

4. Field Line Resonances (FLR) — Why Pc5 is Localized

Each magnetic field line has its own eigenfrequency, determined by its length and the Alfvén speed along it. Because both B and ρ vary along the field line, the eigenfrequency must be computed as an integral:

But the key scaling is simple: longer field lines (larger L) have lower eigenfrequencies.

Field Line Resonance (FLR) is what happens when a broadband fast-mode wave (Section 2) encounters a region where the local eigenfrequency matches the wave frequency. Energy transfers from the fast mode (compressive, global) to the Alfvén mode (transverse, localized). The field line resonates.

But here is the critical point: each field line resonates at its own frequency. At a given frequency f, only the L-shell where f_eigen(L) = f will resonate. Neighboring L-shells have different eigenfrequencies and do not resonate at that frequency. The result is a perturbation that is highly localized in latitude — extending only ΔL ≈ 0.5–1 R_E around the resonant shell.

5. Cavity/Waveguide Modes — Why Pi2 Coherence Recovers

The magnetosphere also acts as a waveguide or cavity for fast-mode (compressive) waves. The boundaries are regions where V_A changes sharply, creating partial reflections:

• Outer boundary: The magnetopause, where the solar wind's low V_A meets the magnetosphere's high V_A.
• Inner boundary: The plasmapause (L ≈ 4), where a sharp density increase causes V_A to drop abruptly. Also, the upper ionosphere serves as a lower boundary.

Fast-mode waves bounce between these boundaries, setting up standing patterns — cavity modes. These are the MHD analog of acoustic resonances in a room. Their periods are determined by the cavity dimensions and the fast-mode speed:

The essential difference from FLR: cavity modes are global. The entire magnetosphere between the boundaries participates in the same oscillation. All ground stations within the cavity see the same phase of compression and rarefaction. This is why your coherence recovers at Pi2 frequencies.

Your Pc5 coherence dip is therefore the transition zone: from global convection (high coherence at low-f) through localized FLR (low coherence at Pc5) back to global cavity modes (recovered coherence at Pi2).

6. Magnetosphere Breathing — Solar Wind Pressure Pulses

The most basic MHD response of the magnetosphere to external forcing: when solar wind dynamic pressure increases, the magnetopause is compressed inward; when it decreases, the magnetopause expands outward. This is a fast-mode compression/rarefaction, and it creates the largest-scale perturbations observable on the ground.

Sudden commencements (SC): A sharp increase in solar wind pressure at storm onset drives rapid magnetopause compression. A fast-mode pulse propagates inward through the magnetosphere at V_f ≈ 1000 km/s. Since the IMAGE network spans ~1000 km and V_f is ~1000 km/s, the transit time across the network is ~1 second. This means SC signals arrive at all stations nearly simultaneously → extremely high coherence.

Continuous pressure fluctuations in the solar wind drive continuous fast-mode compressions at the magnetopause. These are the primary driver of the high coherence at frequencies below the FLR band (< 1.7 mHz in your data).

7. Field Line Decoupling During Reconnection

Reconnection is where the frozen-in condition (Section 1) breaks down. At an X-point (or X-line), the resistivity is locally enhanced, allowing field lines to change their topology. This is not merely a theoretical concept — it is a wave physics event that changes the entire character of MHD wave propagation.

Dayside Reconnection

When the IMF (interplanetary magnetic field) is southward, it is antiparallel to the Earth's northward field at the subsolar magnetopause. At the X-line:

• A closed field line (both ends on Earth) splits into two open field lines (one end on Earth, one end in the solar wind)
• These open lines get swept poleward and tailward by the solar wind flow
• Open flux accumulates in the tail lobes

Nightside (Tail) Reconnection

In the magnetotail, the reverse process:

• Two open field lines approach the tail X-line
• They reconnect → one closed line (snaps earthward as a dipolarization front) + one disconnected IMF line (ejected tailward as a plasmoid)

Wave Physics Consequences

An open field line has no southern reflection point → no standing wave → no FLR eigenfrequency. When reconnection opens or closes field lines, it changes the Alfvén travel time discontinuously. The wave physics of the magnetosphere changes fundamentally. This is why FLR behavior changes dramatically between quiet and disturbed conditions.

8. Plasma Sheet Dynamics — Substorm Field Evolution

Growth Phase

During the substorm growth phase (typically 30–60 minutes):

• Open flux accumulates in the tail lobes (from dayside reconnection, Section 7)
• The lobes expand, squeezing the plasma sheet thinner
• The cross-tail current intensifies to maintain pressure balance: J × B = ∇p
• The magnetic field becomes increasingly stretched — tail-like rather than dipolar

Expansion Onset

When the tail becomes too thin and stretched, instability triggers:

• A Near-Earth Neutral Line (NENL) forms at ~20 R_E
• Reconnection begins, and a dipolarization front surges earthward at ~200–400 km/s
• Stretched field lines snap back to dipolar configuration (think of a rubber band released)
• The cross-tail current is disrupted, and current is redirected as FACs into the ionosphere — the substorm current wedge (SCW)
• The SCW produces the large ΔB_X perturbations your magnetometers record during substorms

Pi2 Onset Signature

The formation of the SCW and the dipolarization front launch waves in the Pi2 band (40–150 s). These are the Pi2 onset pulsations used universally to time substorm onset. They are a mix of cavity modes (Section 5) and transient Alfvén oscillations of the newly dipolarized field lines.

9. The 12-Hour View — How a Ground Station Sees the Magnetosphere

A ground magnetometer at fixed geographic coordinates rotates through magnetic local time (MLT) over 24 hours. At different MLT, it samples completely different parts of the magnetospheric current system:

The station sees a completely different magnetosphere every 6 hours. The magnetic perturbation switches character: from compression-dominated (noon) to current-sheet-dominated (midnight) to shear-dominated (dusk). This is why diurnal variation in spatial statistics is not merely a nuisance — it reflects genuine physical changes in which current system is overhead.

10. Your Four Regimes Explained — The Complete Physical Picture

Now we tie everything together. Your cross-spectral analysis found four coherence regimes. Here is the MHD wave physics behind each one:

Regime 1: >10 min (<1.7 mHz) — HIGH Coherence (0.85–0.95)

These are convection-driven perturbations from the Dungey cycle (Section 7). The electrojet is a single large-scale current sheet spanning the entire auroral zone. The current flows E-W for thousands of kilometers. All stations within the electrojet see essentially the same current → high coherence, long L (~2000–5000 km). The spatial pattern is strongly E-W extended, explaining your low-frequency E-W anisotropy.

Regime 2: 3–10 min (1.7–5.6 mHz) — DIP (0.55–0.75)

This is the Pc5 / FLR regime (Section 4). Each field line resonates at its own frequency, determined by its L-shell. Neighboring L-shells resonate at slightly different frequencies → perturbations are localized in latitude (ΔL ~ 0.5–1 R_E) → low coherence, short L (~800–1500 km). The L(f) minimum is here. FLR drives FACs that are N-S aligned (localized filaments in the ionosphere), beginning the transition from E-W to N-S anisotropy.

Regime 3: 40s–3 min (5.6–25 mHz) — RECOVERY (0.65–0.76)

Pi2 cavity/waveguide modes (Section 5). Fast-mode waves bouncing between magnetopause and plasmapause create global standing patterns. All stations see the same cavity mode → coherence recovers, L increases. The Z component does NOT recover — these are horizontal (X, Y) oscillations. The cavity modes are the MHD "drumhead" of the magnetosphere.

Regime 4: <40s (>25 mHz) — HIGH (0.72–0.84)

Pc3 upstream waves. Ion cyclotron waves generated in the foreshock region, where reflected solar wind ions gyrate and interact with the incoming solar wind. The wave frequency is set by the ion cyclotron frequency in the solar wind: f ≈ 0.015 × B_IMF(nT) Hz. These waves propagate through the magnetosheath and magnetopause into the dayside magnetosphere. They have a single, spatially extended source region → high spatial coherence across the ground network. The Pc3 signal is strongest at noon MLT (closest to the foreshock).