Proton Aurora & Charge-Exchange Dynamics: How Precipitating Magnetospheric Protons and Hydrogen Doppler Shifts Forge Diffuse Polar Glows

Opening Scene: The Ghostly Sheet Beyond the Curtains

Stand atop an exposed plateau in the Norwegian Arctic during the deepest hours of the polar night, and the sky will teach you that not all light is born of the same fire.

The air is bitter, biting through wool and membrane at twenty-five degrees below zero. Every exhalation crystallises into microscopic ice needles that catch the faint glimmer of the stars. High above, the familiar drama of the aurora borealis unfolds in theatrical splendour: razor-sharp, emerald-green curtains pleat and snap across the zenith, their lower borders etched against the blackness with knife-edge precision. These ribboned folds ripple at frantic speeds, dynamic and volatile, dancing in direct communion with invisible forces.

       MAGNETIC FIELD LINES (B)
             |       |       |
[Discrete]   |  e-   |  e-   |   --> Electrons stay tightly bound to B-field
Curtain      |   \   |   /   |       Result: Razor-sharp, dancing curtains
             |    \  |  /    |
========================================================================
             |       |       |
[Proton]     |  H+   |       |   --> Proton captures electron (H+ -> H*)
Diffuse      |   \   |       |   --> Unbound neutral H* flies across field
Arc          |    -----> H*  |   --> Collisional stripping restores H+
             |       |    \  |   Result: Broad, amorphous, Doppler glow

Yet, if you turn your gaze equatorward, away from the furious green ribbons, you may notice something subtly uncanny. Spanning hundreds of kilometres across the southern horizon rests a vast, motionless luminous veil. It possesses neither folds nor rays; it exhibits no furious undulating hems. Instead, it hovers as a faint, velvety expanse—an eerie, diffuse emission of muted slate-blue and deep magenta that appears completely detached from the manic choreography overhead.

No wind stirs this ghostly sheet. It does not flicker. It remains stubbornly diffuse, defying the magnetic architecture that corrals the rest of the ionosphere into tight filaments. To the untrained eye, it seems like a high-altitude haze reflecting distant city lights. But you are observing one of the most sophisticated atomic relay races in space physics: the descent of energetic protons escaping the planet’s magnetic trap through a sequence of quantum charge-exchange disguise.

What's Actually Happening — Plain English First

To understand why this quiet glow behaves so differently from the rippling curtains of a classic aurora, we must first examine the invisible cage that governs our upper atmosphere: the Earth's geomagnetic field.

Most auroral displays are sculpted by energetic electrons. Electrons are extraordinarily light, negatively charged particles. As they cascade downward from the magnetosphere into the upper atmosphere, they are held in a vice-like grip by the geomagnetic field lines. The fundamental law of electrodynamics dictates that a charged particle traversing a magnetic field experiences a deflecting force—the Lorentz force—which compels it to spiral tightly around a single magnetic field line. Because an electron’s mass is so minuscule, the radius of its spiral (its gyroradius) is negligible, often measuring mere metres. Consequently, millions of descending electrons remain strictly confined to their respective magnetic tracks, striking atmospheric gases in narrow, crisp sheets. They paint the sky with the precision of a fine-tipped calligraphy brush.

Protons, by contrast, are the heavyweights of the cosmic plasma. A proton is the nucleus of a hydrogen atom—nearly two thousand times more massive than an electron. When magnetospheric storms accelerate these ions toward the Earth, they enter the thin upper thermosphere at velocities exceeding several thousand kilometres per second.

If they remained purely charged particles, they too would be forced to spiral along magnetic field lines, albeit in wider loops. But as a fast proton plunges into the thermospheric gas between 100 and 300 kilometres above our heads, it encounters neutral oxygen and nitrogen molecules. At these collision energies, a proton does something extraordinary: it snatches an electron from an unsuspecting atmospheric molecule.

                   THE CHARGE-EXCHANGE ENGINE

      Incoming Ion          Atmospheric Gas         Neutralized Atom
        [ H+ ]      +            [ M ]       --->        [ H* ]
    (Charged / Trapped)       (Target Gas)         (Neutral / Ballistic)
                                                           |
                                                   Unbound by Lorentz Force
                                                   Crosses Magnetic Lines
                                                           |
                                                   Radiates Photon (Ha / Hb)
                                                           |
      Restored Ion          Stripped Electron        Stripping Collision
        [ H+ ]      +            [ e- ]      <---   [ H ] + [ M ]
   (Re-trapped on B')

The moment the positively charged proton captures this electron, it undergoes a phase of atomic neutralization, transforming into a high-speed neutral hydrogen atom ($H^*$). Because it now carries a net electrical charge of zero, the Earth's magnetic field completely loses its grip on the particle. The Lorentz force drops instantly to zero:

$$\mathbf{F} = q(\mathbf{E} + \mathbf{v} \times \mathbf{B}) = 0$$

Free from the magnetic cage, the neutral hydrogen atom hurtles forward in a straight ballistic trajectory, sailing unimpeded across adjacent magnetic field lines.

As it travels, the newly formed, excited hydrogen atom sheds excess energy by emitting a photon of light—radiating characteristic spectral signatures in the Balmer series (such as Hydrogen-alpha at $656.3\text{ nm}$ and Hydrogen-beta at $486.1\text{ nm}$). Soon after, it slams into another neutral atmospheric atom, which violently strips the captured electron away. Instantly, the particle becomes a positive ion once more, snapping back under the control of the geomagnetic field—only now, it is gyrating around an entirely different field line kilometres away from where it began.

This rapid, cyclic alternation between charged ion and neutral ghost repeats hundreds of times for every single particle. It acts as an atomic random walk, scattering energy laterally across immense swaths of the thermosphere. The sharp, knife-edge brush of the electron curtain is thereby replaced by a wide, soft-focus atmospheric wash: the diffuse proton aurora.

The Science: Wave-Particle Resonances, Charge-Exchange Kinetics, and Doppler Shifts

For those who wish to delve into the rigorous physical mechanics, the life cycle of a proton aurora can be deconstructed into three distinct stages: magnetospheric pitch-angle scattering, atmospheric charge-transfer equilibrium, and Doppler-shifted optical emission.

1. Pitch-Angle Scattering via EMIC Waves

Protons responsible for auroral precipitation originate predominantly within the ring current and the plasma sheet of the inner magnetosphere, carrying kinetic energies typically between $10\text{ keV}$ and $100\text{ keV}$. Under quiescent conditions, these ions remain trapped within the geomagnetic dipole bottle, bouncing back and forth between magnetic mirror points in opposite hemispheres.

To precipitate into the upper atmosphere, the protons must have their pitch angles—the angle $\alpha$ between their velocity vector $\mathbf{v}$ and the local magnetic field vector $\mathbf{B}$—scattered into the loss cone ($\alpha < \alpha_c$). This pitch-angle diffusion is mediated by resonant wave-particle interactions with Electromagnetic Ion Cyclotron (EMIC) waves.

               MAGNETOSPHERIC LOSS CONE SCATTERING

       Trapped Trajectory            EMIC Wave Interaction
       ==================            =====================
          \        /                     \     ~~~~~ (EMIC Wave)
           \  /\  /                       \   /
            \/  \/                         \ /  Resonant Pitch-Angle Shift
       (Mirrors above atmosphere)           |   Enters Loss Cone (a < ac)
                                            v
                                     [ Thermosphere ]

When an energetic proton traverses an EMIC wave field, it satisfies the Doppler-shifted cyclotron resonance condition:

$$\omega - k_\parallel v_\parallel = \frac{\Omega_p}{\gamma}$$

where $\omega$ is the wave frequency, $k_\parallel$ is the parallel wavevector, $v_\parallel$ is the proton's parallel velocity along the magnetic field line, $\Omega_p = \frac{e B}{m_p}$ is the local non-relativistic proton gyrofrequency, and $\gamma = (1 - v^2/c^2)^{-1/2}$ is the relativistic Lorentz factor. This resonant interaction breaks the first adiabatic invariant ($\mu = \frac{m v_\perp^2}{2B}$), efficiently scattering ions into the loss cone and dumping them into the high-latitude thermosphere.

2. The Charge-Exchange and Stripping Equilibrium Engine

Upon reaching altitudes between $110\text{ km}$ and $250\text{ km}$, where neutral atmospheric densities ($n_n$) become substantial, the incoming proton beam initiates the charge-exchange cycle through two primary cross-sectional interactions with ambient neutrals ($M \in {N_2, O_2, O}$):

1. Electron Capture (Neutralization): $$H^+ + M \xrightarrow{\sigma_{10}} H^* + M^+$$
2. Electron Stripping (Ionization): $$H + M \xrightarrow{\sigma_{01}} H^+ + M + e^-$$

Here, $\sigma_{10}(E)$ denotes the electron capture cross section as a function of kinetic energy $E$, and $\sigma_{01}(E)$ represents the electron stripping cross section.

The spatial evolution of the downward neutral hydrogen flux ($\Phi_H$) and the proton ion flux ($\Phi_{H^+}$) as a function of altitude $z$ is governed by the following coupled differential equations:

$$\frac{d\Phi_H(z)}{dz} = \sigma_{10}(E) n_n(z) \Phi_{H^+}(z) - \sigma_{01}(E) n_n(z) \Phi_H(z)$$

$$\frac{d\Phi_{H^+}(z)}{dz} = -\sigma_{10}(E) n_n(z) \Phi_{H^+}(z) + \sigma_{01}(E) n_n(z) \Phi_H(z)$$

Because collision frequencies in the lower thermosphere ($z \approx 120\text{ km}$) are exceptionally high relative to the transit timescale of a several-keV particle, the beam rapidly reaches a state of charge-transfer equilibrium ($\frac{d\Phi_H}{dz} \approx 0$). Setting the derivatives to zero yields the fundamental equilibrium ratio of neutral to ionized hydrogen within the descending beam:

$$\frac{\Phi_H}{\Phi_{H^+}} = \frac{\sigma_{10}(E)}{\sigma_{01}(E)}$$

Because the neutral atoms are unconstrained by magnetic coordinates, they undergo transverse ballistic displacements of several tens to hundreds of kilometres before being re-ionized, entirely smoothing out any sharp spatial gradients.

       ALTITUDE vs. CHARGE-TRANSFER EQUILIBRIUM

  Altitude (km)
      ^
  300 |  Pure Ion Beam: Phi_H+ dominates (Low ambient density n_n)
      |   |
  200 |   |  Collision onset: Neutralization begins
      |   \
  130 |    \---> Equilibrium Established: Phi_H / Phi_H+ = sigma_10 / sigma_01
      |          Peak Balmer-beta and Lyman-alpha emission layer
  100 |          Beam Thermalization / Complete Energy Loss
      +------------------------------------------------------------> Flux Ratio

3. Optical Signatures and the Asymmetric Doppler Blueshift

When neutral hydrogen atoms emerge from the electron capture collision in an electronically excited state ($H^*$, with principal quantum number $n \ge 3$), they spontaneously decay to lower energy levels, generating the characteristic hydrogen auroral emissions:

• Lyman-$\alpha$ ($2p \rightarrow 1s$): $\lambda_0 = 121.57\text{ nm}$ (Far Ultraviolet)
• Balmer-$\alpha$ ($H_\alpha$, $3d/3p/3s \rightarrow 2p/2s$): $\lambda_0 = 656.28\text{ nm}$ (Red)
• Balmer-$\beta$ ($H_\beta$, $4d/4p/4s \rightarrow 2p/2s$): $\lambda_0 = 486.13\text{ nm}$ (Cyan-Blue)

Crucially, because the emitting hydrogen atoms are moving with massive downward velocities ($v_\parallel \approx 1,000\text{--}3,000\text{ km/s}$) directly toward a ground-based observer looking up along the magnetic field line (the magnetic zenith), the observed spectral lines do not appear at their rest wavelengths $\lambda_0$. Instead, they experience a substantial Doppler blueshift.

The non-relativistic wavelength shift $\Delta \lambda$ observed along the line of sight is expressed as:

$$\Delta \lambda = -\lambda_0 \frac{v_\parallel}{c}$$

where $c$ is the speed of light in vacuum ($2.998 \times 10^8\text{ m/s}$), and $v_\parallel$ is the component of particle velocity oriented toward the observer.

Worked Numerical Example: Clocking an Atmospheric Proton Beam

Let us calculate the theoretical Doppler shift and line profile displacement observed in the Hydrogen-beta ($H_\beta$) emission line by a ground-based spectrometer situated in Tromsø, Norway, pointing directly into the magnetic zenith during a moderate geomagnetic storm.

                  GROUND-BASED SPECTROMETER GEOMETRY

                           Descending H* Beam
                               \     |     /  (Velocity v_parallel ~ 2,200 km/s)
                                \    |    /
                                 v   v   v
                             [ Photon Emission ]  (Wavelength = lambda_0 + delta_lambda)
                                     |
                                     |
                                     v
                           [ Ground Spectrometer ]

Step 1: Establish the Physical Parameters

• Rest wavelength of $H_\beta$: $\lambda_0 = 486.13\text{ nm} = 4861.3\text{ \AA}$
• Mean kinetic energy of precipitating proton: $E_k = 25\text{ keV}$
• Proton rest mass: $m_p = 1.673 \times 10^{-27}\text{ kg}$
• Elementary charge: $e = 1.602 \times 10^{-19}\text{ C}$

Step 2: Compute the Particle Velocity

Converting the kinetic energy to Joules:

$$E_k = 25\text{ keV} \times (1.602 \times 10^{-16}\text{ J/keV}) = 4.005 \times 10^{-15}\text{ J}$$

Using the classical kinetic energy relation ($v \ll c$):

$$v = \sqrt{\frac{2 E_k}{m_p}} = \sqrt{\frac{2 \times 4.005 \times 10^{-15}\text{ J}}{1.673 \times 10^{-27}\text{ kg}}} = \sqrt{4.788 \times 10^{12}} \approx 2.188 \times 10^6\text{ m/s} = 2,188\text{ km/s}$$

Step 3: Determine the Doppler Wavelength Shift

Assuming an average pitch angle giving a line-of-sight velocity toward the zenith of $v_\parallel \approx 2.00 \times 10^6\text{ m/s}$:

$$\Delta \lambda = -486.13\text{ nm} \times \left(\frac{2.00 \times 10^6\text{ m/s}}{2.998 \times 10^8\text{ m/s}}\right) = -486.13 \times 0.006671\text{ nm} \approx -3.24\text{ nm}$$

Step 4: Interpret the Observed Spectral Profile

The resulting wavelength recorded by the spectrometer is:

$$\lambda_{\text{obs}} = \lambda_0 + \Delta \lambda = 486.13\text{ nm} - 3.24\text{ nm} = 482.89\text{ nm}$$

                SPECTRAL LINE INTENSITY PROFILE

  Intensity
      ^
      |                 Observed Peak
      |               (Blueshifted Wing)          Rest Line (lambda_0)
      |                     /\                          |
      |                    /  \                         |
      |                   /    \                        |
      |                 /        \                      |
      |               /            \                    |
      |             /                \                  |
      +------------+------------------+-----------------+--------> Wavelength
                482.0 nm           482.9 nm          486.1 nm
                                (-3.2 nm shift)

Because the beam contains a continuum of particle energies and pitch angles, ground instruments do not measure a single shifted spike. Instead, they observe a broad, distinctly asymmetric profile with an extended blue wing stretching up to $3\text{--}4\text{ nm}$ below the unshifted rest line.

Crucially, because there is no corresponding redshifted wing when looking into the zenith (protons do not travel backward into space while emitting), this asymmetric blueshift constitutes definitive, indisputable spectroscopic proof of descending extraterrestrial ions.

Practical Outdoor Guidance for the Nighttime Observer

While proton auroras lack the high-contrast brightness of their electron-driven counterparts, an experienced skywatcher, equipped with the right space weather metrics and optical techniques, can readily identify and document these profound magnetospheric phenomena.

1. What to Look for in the Sky

• The Equatorward Diffuse Belt: Proton precipitation typically maps to the equatorward boundary of the main auroral oval. If you are watching an auroral display from sub-auroral latitudes (such as Southern Scandinavia, Scotland, the northern United States, or southern New Zealand), look for an expansive, structureless luminous band located south of the main rayed curtains.
• Absence of Fine Morphology: Unlike electron arcs, which feature sharp ribbons, curls, and fast flickering rays ($< 1\text{ km}$ thickness), proton arcs appear as homogeneous sheets spanning $200\text{ to }500\text{ km}$ in latitude. They exhibit virtually no rapid motion, maintaining a steady, smouldering glow over tens of minutes.
• Subtle Chromatic Balance: To the dark-adapted naked eye, proton arcs often appear faint white, silver, or sub-visual due to human scotopic vision limits. However, long-exposure digital photography reveals a distinctive muted violet-magenta cast, resulting from the coexistence of blue $H_\beta$ ($486.1\text{ nm}$), deep red $H_\alpha$ ($656.3\text{ nm}$), and molecular nitrogen band emissions ($N_2^+$ First Negative).

                      AURORAL MORPHOLOGY MAP

        POLEWARD                                         EQUATORWARD
           |                                                  |
           v                                                  v
   +-----------------------+                         +------------------+
   | DISCRETE ELECTRON ARC |                         | PROTON DIFFUSE   |
   | - Razor-thin (<1 km)  |                         |   AURORAL BELT   |
   | - Rapid ray motion    |                         | - Broad (200-500km)|
   | - Vivid Green (557nm) |                         | - Smooth & still |
   | - High dynamic flux   |                         | - Muted Magenta  |
   +-----------------------+                         +------------------+

2. Space Weather Instruments and Telemetry to Monitor

Before heading out into the field, track real-time telemetry provided by authoritative space weather monitoring networks:

• Solar Wind Dynamic Pressure: Proton auroras intensify dramatically during Storm Sudden Commencements (SSCs), when an interplanetary shock wave compresses the dayside magnetosphere. Monitor real-time solar wind density and velocity via the NOAA Space Weather Prediction Center. A sharp spike in dynamic pressure ($P_{\text{dyn}} = \frac{1}{2}\rho v^2 > 5\text{ nPa}$) indicates immediate ring current compression and strong proton precipitation.
• The Sym-H / Dst Index: A rapidly dropping Dst index (exceeding $-50\text{ to }-150\text{ nT}$) signals energetic ion injection into the ring current, amplifying EMIC wave generation and widening the proton auroral belt equatorward.
• Global Space Weather Alerts: Reference coordinated alerts and forecasts from the Met Office Space Weather Operations Centre and the World Meteorological Organization Space Weather Portal.
• Satellite Far-Ultraviolet Imagery: Spaceborne imagers, such as the historic Far Ultraviolet (FUV) instrument on the NASA IMAGE Mission, utilized the Doppler-shifted Lyman-$\alpha$ line ($121.6\text{ nm}$) to map global proton auroral precipitation continuously across both the daylit and nightside polar hemispheres, free from contamination by discrete electron emissions.

                   OBSERVER DECISION MATRIX

             Is there a sharp solar wind dynamic pressure jump?
                                    |
                    +---------------+---------------+
                    | YES                           | NO
                    v                               v
         Monitor Equatorward Horizon       Await Substorm Expansion
         Look for Diffuse Belt             Watch Zenith for Discrete
         (Proton Aurora Signature)         Electron Curtains
                    |
                    v
   Use Tilted Narrowband H-beta Filter (486.1 nm)
   to Detect Characteristic Doppler Blueshift Wing

3. Optical Filtering Techniques for the Advanced Observer

Standard astrophotography filters often struggle to distinguish proton auroras from background skyglow and electron diffuse glow. Advanced observers and citizen scientists employ the tilt-filter technique: * Mount a narrow-band interference filter centered on the $H_\beta$ line ($486.13\text{ nm}$ with a full-width at half-maximum, $\text{FWHM} \approx 1.5\text{ nm}$) in front of a monochromatic sensor. * When mounted perpendicularly to the optical axis, the filter transmits the rest wavelength. * By tilting the filter by several degrees relative to the incoming optical path, the transmitted passband shifts toward shorter wavelengths (blue-shifting the filter response). * If the observed glow brightens when the filter is tilted toward $482\text{--}484\text{ nm}$, you have successfully isolated and confirmed the Doppler-blueshifted hydrogen emission of an active proton precipitation beam.

Key Takeaway for the Night Sky

The next time you stand beneath the Arctic or subpolar sky, remember that the quiet, motionless veil on the horizon is not an absence of activity, but the signature of a vast atomic random walk—a testament to charged ions masquerading as neutral ghosts as they plummet through our atmosphere.