Part 7, Chapter 2

Earth's Magnetosphere

The terrestrial magnetosphere is a plasma laboratory where MHD, kinetic theory, and reconnection physics can be tested against in-situ spacecraft measurements. This chapter derives the key structures — magnetopause, bow shock, tail, and radiation belts — from first principles.

2.1 The Geomagnetic Dipole

To first approximation, Earth's magnetic field is a magnetic dipole with moment M_E ≈ 8.0 × 10²² A·m². In spherical coordinates (r, θ) the dipole field components are:

$$B_r = -\frac{2\mu_0 M_E}{4\pi r^3}\cos\theta, \qquad B_\theta = -\frac{\mu_0 M_E}{4\pi r^3}\sin\theta$$

The total field magnitude is:

$$B = \frac{\mu_0 M_E}{4\pi r^3}\sqrt{1 + 3\cos^2\theta}$$

At the equator (θ = π/2) on Earth's surface (r = R_E), this gives B₀ ≈ 3.1 × 10⁻⁵ T = 0.31 Gauss. The field line equation follows from dr/B_r = r dθ/B_θ:

$$r = L R_E \sin^2\theta$$

where L is the McIlwain L-parameter — the equatorial crossing distance in units of R_E. A field line with L = 6.6 maps to the geosynchronous orbit.

2.2 The Magnetopause: Chapman–Ferraro Model

The magnetopause is the boundary where the solar wind dynamic pressure balances the magnetic pressure of the compressed dipole field. The Chapman–Ferraro standoff distance is derived from pressure balance at the subsolar point.

Derivation: Magnetopause Standoff Distance

Step 1. The solar wind arrives with density ρ_sw and velocity v_sw. At the subsolar stagnation point, the total ram pressure is:

$$P_{dyn} = \frac{1}{2}\rho_{sw} v_{sw}^2$$

Step 2. Inside the magnetopause, the field is approximately a compressed dipole. At the subsolar point, the field is enhanced by a factor of 2 due to the image currents on the magnetopause (Chapman–Ferraro current sheet):

$$B_{MP} = 2B_{eq}(R_{MP}) = 2\frac{B_0 R_E^3}{R_{MP}^3}$$

where B₀ is the equatorial surface field.

Step 3. Setting magnetic pressure equal to ram pressure:

$$\frac{B_{MP}^2}{2\mu_0} = \frac{1}{2}\rho_{sw} v_{sw}^2$$

Step 4. Substituting B_MP and solving for R_MP:

$$\frac{4 B_0^2 R_E^6}{2\mu_0 R_{MP}^6} = \frac{1}{2}\rho_{sw} v_{sw}^2$$

$$\boxed{R_{MP} = R_E\left(\frac{2 B_0^2}{\mu_0 \rho_{sw} v_{sw}^2}\right)^{1/6}}$$

For typical solar wind conditions (n ≈ 5 cm⁻³, v ≈ 400 km/s), this gives R_MP ≈ 10 R_E. During intense solar storms, R_MP can be pushed inside geosynchronous orbit (6.6 R_E).

The magnetopause is not a rigid boundary. It supports surface waves (Kelvin–Helmholtz instabilities) and can be penetrated by reconnection when the IMF has a southward component.

2.3 The Bow Shock

The solar wind is supersonic (M_s ~ 8) and super-Alfvénic (M_A ~ 6), so a standing bow shock forms upstream of the magnetopause. The shock is collisionless — dissipation occurs through wave-particle interactions rather than binary collisions.

Derivation: Rankine–Hugoniot Jump Conditions

Conservation of mass, momentum, and energy across the shock (in the shock frame) gives:

$$[\rho v_n] = 0 \qquad \text{(mass)}$$

$$\left[\rho v_n^2 + p + \frac{B_t^2}{2\mu_0}\right] = 0 \qquad \text{(normal momentum)}$$

$$\left[\rho v_n \mathbf{v}_t - \frac{B_n \mathbf{B}_t}{\mu_0}\right] = 0 \qquad \text{(tangential momentum)}$$

$$[B_n] = 0, \qquad [v_n B_t - v_t B_n] = 0 \qquad \text{(Maxwell)}$$

For a perpendicular shock (B_n = 0) with Mach number M, the density compression ratio is:

$$\boxed{\frac{\rho_2}{\rho_1} = \frac{(\gamma+1)M^2}{(\gamma-1)M^2 + 2}}$$

For a strong shock (M → ∞) with γ = 5/3, the maximum compression ratio is ρ₂/ρ₁ = 4. The bow shock standoff distance is approximately:

$$\Delta_{BS} \approx 1.1 R_{MP}\frac{(\gamma-1)M^2 + 2}{(\gamma+1)(M^2-1)} \approx 3\,R_E$$

Between the bow shock and magnetopause lies the magnetosheath — a region of shocked, turbulent solar wind plasma with β ~ 1–10. The plasma is subsonic but still flows around the magnetopause and reconnects at the flanks.

2.4 Magnetic Reconnection at the Magnetopause

When the interplanetary magnetic field (IMF) has a southward component (B_z < 0), it is antiparallel to the geomagnetic field at the subsolar magnetopause. Magnetic reconnection opens the magnetosphere to the solar wind, driving magnetospheric convection.

Derivation: Sweet–Parker Reconnection Rate

Step 1. Consider a resistive current sheet of length L and thickness δ. Plasma flows in at velocity v_in (normal to B) and is ejected at velocity v_outalong the sheet.

Step 2. Mass conservation requires:

$$\rho v_{in} L = \rho v_{out} \delta$$

Step 3. The outflow is accelerated by the magnetic tension force. Momentum balance gives the outflow velocity as the Alfvén speed:

$$v_{out} = v_A = \frac{B}{\sqrt{\mu_0 \rho}}$$

Step 4. The magnetic flux is carried into the sheet by the inflow and dissipated by resistive diffusion. Steady-state induction equation in the sheet:

$$v_{in} B \sim \frac{\eta B}{\delta} \quad \Rightarrow \quad \delta = \frac{\eta}{v_{in}}$$

Step 5. Combining mass conservation and the diffusion relation:

$$v_{in} = \frac{v_A \delta}{L} = \frac{v_A \eta}{v_{in} L}$$

$$v_{in}^2 = \frac{v_A \eta}{L} \quad \Rightarrow \quad \boxed{\frac{v_{in}}{v_A} = \frac{1}{\sqrt{R_m}} = \frac{1}{\sqrt{S}}}$$

where S = Lv_A/η is the Lundquist number. For the magnetopause, S ~ 10¹², giving v_in/v_A ~ 10⁻⁶ — far too slow to explain observed reconnection rates.

Petschek Reconnection

Petschek (1964) showed that if the diffusion region is localized, standing slow-mode shocks extend from it, converting magnetic energy over a much larger region. The Petschek rate is:

$$\boxed{\frac{v_{in}}{v_A} \sim \frac{\pi}{8\ln S}}$$

For S = 10¹², this gives v_in/v_A ~ 0.01–0.1, consistent with observations of fast reconnection at the magnetopause (~0.1 v_A).

2.5 The Dungey Cycle

Dungey (1961) proposed a cycle of open magnetospheric convection driven by reconnection:

Dayside reconnection: IMF connects to geomagnetic field lines at the subsolar magnetopause
Convection: Open field lines are dragged tailward by solar wind flow over the polar caps
Tail reconnection: Open field lines reconnect in the magnetotail (x-line at ~20–25 R_E)
Return flow: Closed field lines convect sunward through the inner magnetosphere

The convection electric field associated with this cycle is:

$$\mathbf{E}_{conv} = -\mathbf{v}_{sw} \times \mathbf{B}_{IMF}$$

The total cross-polar-cap potential drop (transpolar potential) is:

$$\Phi_{PC} = \int E_{conv}\, dl \approx v_{sw} B_z^{(south)} \cdot d_{eff}$$

where d_eff ~ 5–10 R_E is the effective reconnection width. Typical values are Φ_PC ≈ 30–100 kV during quiet to moderate activity, rising to > 200 kV during storms.

Derivation: Reconnection-Driven Energy Input

The electromagnetic energy flux (Poynting flux) entering the magnetosphere through the open field line region:

$$P_{input} = \frac{1}{\mu_0}(\mathbf{E} \times \mathbf{B}) \cdot \hat{n}\, A$$

Using E = v_inB and A ~ πR²_eff:

$$\boxed{P_{input} \approx \frac{v_{in} B^2}{\mu_0}\pi R_{eff}^2 \sim 10^{11}\text{–}10^{13}\;\text{W}}$$

This power drives auroral precipitation, ring current injection, and Joule heating of the ionosphere. The epsilon parameter (Akasofu coupling function) parameterizes this:

$$\epsilon = \frac{4\pi}{\mu_0} v_{sw} B^2 l_0^2 \sin^4\!\left(\frac{\theta_{clock}}{2}\right)$$

where l₀ ≈ 7 R_E and θ_clock is the IMF clock angle in the GSM y-z plane.

2.6 The Magnetotail

Open field lines dragged tailward by the Dungey cycle form two tail lobes separated by a thin current sheet (the cross-tail current). The Harris current sheet is the canonical equilibrium model.

Derivation: Harris Current Sheet

We seek a 1D equilibrium with B_x(z) reversing sign across z = 0. The Harris solution is:

$$B_x(z) = B_0 \tanh\left(\frac{z}{\delta}\right)$$

Pressure balance requires p + B²/2μ₀ = const. Using the field profile:

$$p(z) = \frac{B_0^2}{2\mu_0}\text{sech}^2\left(\frac{z}{\delta}\right) + p_\infty$$

The current density from Ampère's law:

$$J_y = \frac{1}{\mu_0}\frac{dB_x}{dz} = \frac{B_0}{\mu_0 \delta}\text{sech}^2\left(\frac{z}{\delta}\right)$$

The half-thickness δ of the magnetotail current sheet is typically 0.5–2 R_E during quiet times, thinning to ~100 km before substorm onset.

Tail lobe field: The magnetic flux in each lobe is conserved. Pressure balance between the lobe magnetic pressure and solar wind gives the tail radius:

$$R_{tail} \approx R_{MP}\left(\frac{B_{tail}}{B_{MP}}\right)^{1/2}$$

with B_tail ≈ 20–30 nT, giving R_tail ≈ 20–30 R_E.

2.7 Radiation Belts & Trapped Particles

Energetic charged particles trapped in the dipole field undergo three quasi-periodic motions: gyration, bounce, and drift. Each is associated with an adiabatic invariant.

The Three Adiabatic Invariants

First invariant (μ): The magnetic moment, conserved when the field changes slowly compared to the gyroperiod:

$$\mu = \frac{mv_\perp^2}{2B} = \text{const}$$

Second invariant (J): The bounce action, conserved when the field changes slowly compared to the bounce period:

$$J = \oint p_\parallel\, ds = 2\int_{s_1}^{s_2} mv_\parallel\, ds = \text{const}$$

where s₁, s₂ are the mirror points along the field line.

Third invariant (Φ): The magnetic flux enclosed by the drift shell:

$$\Phi = \oint \mathbf{A} \cdot d\mathbf{l} = \text{const}$$

Derivation: Mirror Point Condition

A particle with pitch angle α₀ at the equator (where B = B₀) mirrors when all kinetic energy is in the perpendicular component (v_∥ = 0). Conservation of μ gives:

$$\frac{v_\perp^2(0)}{2B_0} = \frac{v^2}{2B_m}$$

Since v_⊥(0) = v sin α₀:

$$\boxed{\sin^2\alpha_0 = \frac{B_0}{B_m}}$$

Particles with pitch angles smaller than the loss cone angle α_LC will mirror below the atmosphere and be lost:

$$\sin^2\alpha_{LC} = \frac{B_0}{B_{atm}} = \frac{1}{L^3\sqrt{1 + 3(1-1/L)}} \approx \frac{1}{L^3}$$

At L = 6, α_LC ≈ 3.8° — the loss cone is very narrow, so most particles are stably trapped.

Derivation: Gradient-Curvature Drift

In a dipole field, trapped particles drift azimuthally due to the gradient and curvature of B. The combined drift velocity is:

$$\mathbf{v}_{gc} = \frac{m}{qB^3}\left(\frac{v_\perp^2}{2} + v_\parallel^2\right)\mathbf{B}\times\nabla B$$

Ions drift westward and electrons drift eastward, producing a net westward ring current. The total ring current energy content is related to the Dst index via the Dessler–Parker–Sckopke relation:

$$\boxed{\Delta B_{surf} = -\frac{\mu_0}{4\pi}\frac{2E_{RC}}{B_0 R_E^3} = \text{Dst}}$$

where E_RC is the total energy of the ring current. A Dst of −100 nT corresponds to E_RC ≈ 4 × 10¹⁵ J.

Bounce and drift periods: For equatorially mirroring particles (α₀ = 90°):

$$\tau_{bounce} \approx \frac{L R_E}{v}(3.7 - 1.6\sin\alpha_0), \qquad \tau_{drift} = \frac{2\pi q B_0 R_E^2}{3 L W}$$

where W is the particle kinetic energy. For 1 MeV protons at L = 4: τ_bounce ≈ 20 s, τ_drift ≈ 1000 s.

2.8 Magnetospheric Substorms

Substorms are quasi-periodic (2–4 hour) episodes of energy storage and release in the magnetotail. The substorm cycle consists of three phases:

Growth phase (~1 hour): Southward IMF drives dayside reconnection faster than nightside. Open flux accumulates in the tail lobes. The cross-tail current intensifies and the current sheet thins.

Expansion phase (~30 min): Reconnection onset in the near-Earth tail (x ~ 20–25 R_E) releases stored magnetic energy. Dipolarization front propagates earthward. A plasmoid is ejected tailward. Energetic particles are injected into the inner magnetosphere.

Recovery phase (~1 hour): Reconnection rate decreases. IMF turns northward. Magnetotail returns to pre-substorm configuration.

The energy released during a typical substorm is:

$$E_{substorm} \sim \frac{B_{tail}^2}{2\mu_0} \cdot V_{lobe} \sim 10^{15}\;\text{J}$$

This energy is partitioned roughly as: ~30% to auroral precipitation, ~30% to Joule heating in the ionosphere, ~20% to ring current enhancement, and ~20% to plasmoid ejection.

← Ch 1: Solar Wind Ch 3: Aurora →

Plasma Physics Course