Detection of lightning in Saturns Northern Hemisphere

During Cassini flyby of Saturn at a radial distance 6.18R_s (Saturn Radius), a signal was detected from about 200 to 430 Hz that had the proper dispersion characteristics to be a whistler. The frequency-time dispersion of the whistler was found to be 81 Hz1/2s. Based on this dispersion constant, we determined, from a travel time computation, that the whistler must have originated from lightning in the northern hemisphere of Saturn. Using a simple centrifugal potential model consisting of water group ions, and hydrogen ions we also determine the fractional concentration and scale height that gave the best fit to the observed dispersion.

💡 Research Summary

The paper reports the detection and analysis of a whistler‑type electromagnetic wave recorded by the Cassini spacecraft during its October 28 2004 flyby of Saturn at a radial distance of 6.18 Rₛ and a spacecraft latitude of 12.31°. The signal spanned roughly 200–430 Hz and lasted about three seconds. Because whistler waves travel faster at higher frequencies, the observed frequency‑time structure displayed the classic dispersion where the high‑frequency components arrived first. By plotting arrival time versus inverse square root of frequency (the Eckersley law) the authors obtained a straight line whose slope corresponds to a dispersion constant D ≈ 81 Hz·s¹ᐟ².

To interpret this dispersion, the authors employed the theoretical framework originally developed by Helliwell (1965). In the limit where the electron cyclotron frequency ω_c greatly exceeds the electron plasma frequency ω_p, the group refractive index simplifies to n_g ≈ ½ ω_p/ω_c, leading to a travel‑time expression D = ½ c (ω_p/ω_c) L, where L is the integral of the magnetic‑field‑line path from the lightning source to the spacecraft. The magnetic field was modeled as a tilted dipole (Parks 1991), giving B(θ) = B₀ (1 + 3 sin²θ)¹ᐟ²⁻¹. Using the measured field of 89.513 nT at the spacecraft latitude, the equatorial field B₀ was inferred to be 73 nT, which yields the local cyclotron frequency.

The plasma density model assumes that beyond ~2.3 Rₛ the dominant force shaping the plasma distribution is the centrifugal force of the co‑rotating magnetosphere. The authors derived a Boltzmann‑type density profile for each ion species, n_i = n₀ exp