pith. sign in

arxiv: 2604.12729 · v1 · submitted 2026-04-14 · 🌌 astro-ph.EP · physics.space-ph

Constraining Electron-Impact Ionization of O₂ Through UV Aurora Observations at Ganymede

Stefan Duling , Joachim Saur , Darrell Strobel , Philippa Molyneux , Jamey R. Szalay , Thomas K. Greathouse This is my paper

Pith reviewed 2026-05-10 14:25 UTC · model grok-4.3

classification 🌌 astro-ph.EP physics.space-ph
keywords Ganymedeelectron-impact ionizationOI 1356 ÅUV auroraO2 atmosphereJuno UVSionospheric outflowsurface ice erosion
0
0 comments X

The pith

UV aurora brightness at Ganymede directly converts to electron-impact ionization rates ten times higher than photoionization.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a way to turn measurements of faint ultraviolet light from oxygen atoms into estimates of how many ions electrons create when they strike Ganymede's thin oxygen atmosphere. Cross-section data show that for every 10 to 60 excitations that produce the 1356-angstrom line, one ionization occurs, no matter the electron energy, which caps the conversion uncertainty at less than a factor of six. When this relation is applied to Juno spacecraft images of the moon's auroral ovals and background regions, the resulting global ionization map exceeds what sunlight alone can produce by at least an order of magnitude. The finding matters because it supplies the first observation-based total for how many ions escape the moon and how quickly its surface ice is lost.

Core claim

We present a novel approach to quantify electron-impact ionization rates directly through OI 1356 Å emission brightness observations. The analysis of measured cross sections reveals that the ionization-to-excitation ratio is limited to 10-60 over all electron energies, reducing the uncertainty of estimating ionization rates to a factor less than 6. We apply this method to Juno UVS observations of Ganymede's aurora. We find that the OI 1356 Å brightness of the auroral ovals is well described by 3-5° latitude wide Gaussian distributions centered on the open-closed field line boundary, with an average peak of 120 R. The average brightness outside the ovals in the polar and equatorial background

What carries the argument

The ionization-to-excitation ratio for the OI 1356 Å line, bounded between 10 and 60 for all electron energies, which converts observed brightness into ionization rate.

If this is right

  • Electron-impact ionization produces at least ten times more ions than photoionization across Ganymede's entire atmosphere.
  • Total global ionization reaches 1.3-7.6×10^{26} s^{-1}, with column rates of roughly 5×10^9 cm^{-2}s^{-1} inside the auroral ovals and 3×10^8 cm^{-2}s^{-1} outside.
  • Transport dominates ion loss, so radio-occultation electron densities match model predictions only when outflow is included.
  • O2+ outflow totals 0.1-2×10^{26} s^{-1} or 0.5-11 kg s^{-1}, eroding surface ice at 0.03-0.5 cm per million years.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same brightness-to-ionization conversion could map ionization on other icy moons that show oxygen aurora without needing separate electron-density models.
  • Repeated UV observations during different Jupiter plasma conditions would reveal how variable the ionization rate is on short timescales.
  • Higher ionization strengthens the expected link between magnetospheric electrons and surface sputtering, suggesting faster redistribution of ice components than current models assume.

Load-bearing premise

The ionization-to-excitation ratio for OI 1356 Å stays between 10 and 60 for every electron energy present at Ganymede and that other processes add little to the observed emission.

What would settle it

A spacecraft measurement of both OI 1356 Å brightness and local electron density or ionization rate during a Ganymede flyby that falls outside the predicted range by more than a factor of six.

Figures

Figures reproduced from arXiv: 2604.12729 by Darrell Strobel, Jamey R. Szalay, Joachim Saur, Philippa Molyneux, Stefan Duling, Thomas K. Greathouse.

Figure 1
Figure 1. Figure 1: a) Cross sections for electron-impact ionization and OI 1356 ˚A dissociative excitation of O2 as function of electron energy from Lindsay & Mangan (2003) (blue), Schram et al. (1965) (dashed blue) and Kanik et al. (2003) (green). The measurements are extrapolated using Bethe-Oppenheimer relations (dotted). Excitation cross sections for atomic oxygen (black) are from Julienne & Davis (1976). b) Ratio of O2 … view at source ↗
Figure 2
Figure 2. Figure 2: a) OI 1356 ˚A emission brightness as observed by Juno UVS during the PJ34 flyby. Reflected sunlight has been subtracted. The color bar is saturated, 56 pixels have values above 300 R, 15 pixels above 500 R. The orange line shows the terminator. The right color bar shows electron-impact ionization rates of O2 for an ionization-to-excitation ratio of βion/1356 = 40, characterizing measured electron distribut… view at source ↗
Figure 3
Figure 3. Figure 3: OI 1356 ˚A brightness as function of latitudinal distance ∆θOCFB,fit from the fitted OCFB (Figure 2b). The left side shows the northern hemisphere, the right side the southern hemisphere. a) Number of pixels in each 1° wide bin of latitude. b): Brightness of each pixel sorted by 1° wide bins of latitude. c): OI 1356 ˚A emission brightness averaged over all longitudes as function of ∆θOCFB,fit (red, blue). … view at source ↗
Figure 4
Figure 4. Figure 4: O + 2 density at Ganymede’s surface, modeled from the chemical equilibrium between ionization and recombination of O2, using the electron-impact ionization rates from Figure 2c. Photoionization is included in the sunlight region, which is shown by the terminator in orange. White numbers show modeled density values in cm−3 . For comparison, the red numbers show electron densities detected by Galileo (G8) an… view at source ↗
Figure 5
Figure 5. Figure 5: Percentages of ion loss by transport. The figures show how much of the O+ 2 produced by electron-impact and photoionization is transported out of an ionospheric column between Ganymede’s surface and 200 km altitude as function of surface densities n0 and scale heights He− . The top row uses a realistic ionization-to-excitation ratio of βion/1356 = 40. Figure a considers ionization rates consistent with 8 R… view at source ↗
Figure 6
Figure 6. Figure 6: Lower limits of electron density and downward energy fluxes required to excite O2 OI 1356 ˚A emissions of certain brightnesses. The maximum possible excitation rate coefficient of a monoenergetic distribution at 185 eV was used for calculation (k1356,max = 4 × 10−9 cm3 s −1 ). Associated column electron-impact ionization rates were calculated with βion/1356 = 46, which corresponds to 185 eV electrons. The … view at source ↗
read the original abstract

While photoionization rates of Ganymede's O$_2$ dominated atmosphere are well constrained, the contribution of electron-impact ionization is rather uncertain. Previous quantitative estimates have relied on assumptions about densities and energy distributions of precipitating electrons, or on rare spacecraft measurements that cannot be unambiguously mapped to the regions of ionization. In this study, we present a novel approach to quantify electron-impact ionization rates directly through OI 1356 \r{A} emission brightness observations. The analysis of measured cross sections reveals that the ionization-to-excitation ratio is limited to 10-60 over all electron energies, reducing the uncertainty of estimating ionization rates to a factor less than 6. We apply this method to Juno UVS observations of Ganymede's aurora. We find that the OI 1356 \r{A} brightness of the auroral ovals is well described by 3-5{\deg} latitude wide Gaussian distributions centered on the open-closed field line boundary, with an average peak of 120 R. The average brightness outside the ovals in the polar and equatorial background regions is ~8 R. From these observations, we derive a global map of electron-impact ionization rates, which are at least an order of magnitude higher than photoionization rates. The estimated total global ionization rate is 1.3-7.6$\times$10$^{26}$ s$^{-1}$, with average column rates of ~5$\times$10$^{9}$ cm$^{-2}$s$^{-1}$ in the ovals and ~3$\times$10$^{8}$ cm$^{-2}$s$^{-1}$ in the background regions. Comparison of radio occultation measurements with predicted electron densities indicates that transport processes are the dominant loss mechanism in Ganymede's ionosphere. The rate of ionospheric outflow of O$_2^+$ is 0.1-2$\times$10$^{26}$ s$^{-1}$ or 0.5-11 kg s$^{-1}$, indicating 0.03-0.5 cm Myr$^{-1}$ erosion of Ganymede's surface ice.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper introduces a method to constrain electron-impact ionization rates of O2 in Ganymede's atmosphere by bounding the ionization-to-excitation ratio for OI 1356 Å (10-60 across all electron energies) from measured cross sections, then scaling Juno UVS auroral brightness observations. This produces a global ionization rate map with total rates 1.3-7.6×10^{26} s^{-1} (ovals ~5×10^9 cm^{-2}s^{-1}, background ~3×10^8 cm^{-2}s^{-1}), exceeding photoionization by an order of magnitude, and infers transport-dominated ion loss with O2+ outflow 0.1-2×10^{26} s^{-1} and surface ice erosion 0.03-0.5 cm Myr^{-1}.

Significance. If the ratio bound holds for Ganymede's electron spectra, the work supplies an observationally anchored ionization estimate that avoids assumed electron distributions or sparse in-situ data, directly linking UV brightness to ionization and loss processes. The factor-of-<6 uncertainty reduction and explicit comparison to photoionization rates represent a concrete advance for modeling Ganymede's ionosphere and atmospheric erosion.

major comments (2)
  1. [cross-section analysis and abstract] The ionization-to-excitation ratio bound of 10-60 (abstract and cross-section analysis section) is load-bearing for every quantitative result. The manuscript must demonstrate that this ratio remains within the stated limits for the specific electron energy distributions precipitating at Ganymede (typically 10-1000 eV), including an explicit plot or table of the ratio versus energy, and must rule out significant 1356 Å contributions from dissociative excitation of O2, radiative cascades, or optical-depth effects that would violate the direct brightness-to-rate conversion.
  2. [results and discussion] § on global rate derivation: the total ionization rate (1.3-7.6×10^{26} s^{-1}) and column rates are obtained by direct scaling of observed 1356 Å brightness by the ratio range. Any excursion of the ratio outside 10-60 for the actual precipitating spectrum scales all reported rates and the order-of-magnitude dominance claim by more than the stated factor of 6; the paper should therefore provide a sensitivity test using plausible Ganymede electron spectra rather than the 'all energies' bound alone.
minor comments (2)
  1. [observations and analysis] The abstract states that ovals are 'well described by 3-5° latitude wide Gaussian distributions' with peak 120 R; the full text should report the fitting procedure, reduced-χ² values, and uncertainties on width and amplitude.
  2. [observations and analysis] Background brightness is given as '~8 R'; clarify whether this is an average over specific latitude bands and whether it includes any residual airglow or instrumental background after subtraction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have helped us improve the clarity and robustness of our analysis. We address each major comment point by point below.

read point-by-point responses

  1. Referee: [cross-section analysis and abstract] The ionization-to-excitation ratio bound of 10-60 (abstract and cross-section analysis section) is load-bearing for every quantitative result. The manuscript must demonstrate that this ratio remains within the stated limits for the specific electron energy distributions precipitating at Ganymede (typically 10-1000 eV), including an explicit plot or table of the ratio versus energy, and must rule out significant 1356 Å contributions from dissociative excitation of O2, radiative cascades, or optical-depth effects that would violate the direct brightness-to-rate conversion.

    Authors: We agree that the applicability of the 10-60 bound to Ganymede's electron energies requires explicit demonstration. In the revised manuscript we add a new figure plotting the ionization-to-excitation ratio versus electron energy (1 eV to 5 keV) derived directly from the laboratory cross-section data. For the 10-1000 eV range relevant to Ganymede, the ratio remains between 12 and 55. We have also expanded the cross-section analysis section with a dedicated paragraph addressing dissociative excitation of O2, radiative cascades, and optical-depth effects. Under Ganymede's low column densities the 1356 Å optical depth is <0.1, and literature on Ganymede aurora indicates that direct excitation dominates; these points are now quantified to support the direct brightness-to-rate conversion. revision: yes

  2. Referee: [results and discussion] § on global rate derivation: the total ionization rate (1.3-7.6×10^{26} s^{-1}) and column rates are obtained by direct scaling of observed 1356 Å brightness by the ratio range. Any excursion of the ratio outside 10-60 for the actual precipitating spectrum scales all reported rates and the order-of-magnitude dominance claim by more than the stated factor of 6; the paper should therefore provide a sensitivity test using plausible Ganymede electron spectra rather than the 'all energies' bound alone.

    Authors: We thank the referee for highlighting the value of a spectrum-specific test. We have added a sensitivity analysis in the results section using representative Ganymede precipitating electron spectra (power-law distributions between 10-1000 eV constrained by Juno JEDI and UVS observations). For these spectra the effective ratio lies between 18 and 42, narrowing the global ionization rate to 2.0-5.5×10^{26} s^{-1}. The order-of-magnitude dominance over photoionization is preserved. The revised text, abstract, and discussion now incorporate this test and the updated uncertainty range. revision: yes

Circularity Check

0 steps flagged

No significant circularity; rates derived from independent cross sections and observations

full rationale

The central derivation scales observed OI 1356 Å brightness by an ionization-to-excitation ratio (10-60) obtained from analysis of external measured cross sections. This ratio is not fitted to the Juno UVS data or defined from the paper's own observations, and no equation reduces the output ionization rates (global map, total rate 1.3-7.6×10^{26} s^{-1}, column rates) to a self-referential input. No load-bearing self-citation or ansatz is invoked for the ratio or the scaling step. The chain remains self-contained against external lab data.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that OI 1356 Å emission is produced solely by electron-impact excitation of O2 with a known bounded ratio to ionization, and that observed brightness maps directly to local ionization without transport or other source contamination.

axioms (2)
  • domain assumption OI 1356 Å emission brightness is produced primarily by electron-impact excitation of O2
    Invoked to link brightness observations to ionization rates via the cross-section ratio
  • domain assumption The ionization-to-excitation ratio remains between 10 and 60 across relevant electron energies
    Derived from measured cross sections; used to bound the conversion factor

pith-pipeline@v0.9.0 · 5722 in / 1499 out tokens · 28734 ms · 2026-05-10T14:25:35.239353+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    Alday, J., Roth, L., Ivchenko, N., Retherford, K.D., Becker, T.M., Molyneux, P., & Saur, J. (2017). New constraints on Ganymede’s hydrogen corona: Analysis of Lyman-αemissions observed by HST/STIS between 1998 and 2014,Planetary and Space Science,vol. 148, pp. 35–44, ISSN 0032- 0633, 10.1016/j.pss.2017.10.006 Allegrini, F., Bagenal, F., Ebert, R.W., Louar...

    work page doi:10.1016/j.pss.2017.10.006 2017

  2. [2]

    108(E11), ISSN 0148-0227, 10.1029/2000je001423 Kanik, I., Trajmar, S., & Nickel, J.C

    Absolute emission cross sections of the OI(130.4 nm) and OI(135.6 nm) lines,Journal of Geophysical Research: Planets,vol. 108(E11), ISSN 0148-0227, 10.1029/2000je001423 Kanik, I., Trajmar, S., & Nickel, J.C. (1993). Total electron scattering and electronic state excitations cross sections for O2, CO, and CH4,Journal of Geophysical Research: Planets,vol. 9...

    work page doi:10.1029/2000je001423 1993

  3. [3]

    108(E11), ISSN 0148-0227, 10.1029/2000je001422 Marconi, M

    Kinetic energy distributions of fast oxygen atoms,Journal of Geophysical Research: Planets,vol. 108(E11), ISSN 0148-0227, 10.1029/2000je001422 Marconi, M. (2007). A kinetic model of Ganymede's atmosphere,Icarus,vol. 190(1), pp. 155–174, ISSN 0019-1035, 10.1016/j.icarus.2007.02.016 Marzok, A., Schlegel, S., Saur, J., Roth, L., Grodent, D., Strobel, D.F., &...

    work page doi:10.1029/2000je001422 2007