# Water-Group Pickup Ions From Europa-Genic Neutrals Orbiting Jupiter

J. R. Szalay, H. T. Smith, E. J. Zirnstein, D. J. Mccomas, L. J. Begley, F.

Bagenal, P. A. Delamere, R. J. Wilson, P. W. Valek, A. R. Poppe, et al.

► To cite this version:

J. R. Szalay, H. T. Smith, E. J. Zirnstein, D. J. Mccomas, L. J. Begley, et al. Water-Group Pickup Ions From Europa-Genic Neutrals Orbiting Jupiter. Geophysical Research Letters, 2022, 49, 10.1029/2022GL098111. insu-03867498

HAL Id: insu-03867498

https://insu.hal.science/insu-03867498

Submitted on 23 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Geophysical Research Letters®

## RESEARCH LETTER

10.1029/2022GL098111

**Key Points:**

* First identification of H₂⁺ in the Jovian magnetosphere from 13 to 18 R<sub>J</sub>, with a density ratio of H₂⁺/H⁺ = 8 ± 4%

* H₂⁺ in Jupiter's magnetosphere is predominantly produced by Europa's neutral toroidal cloud

* Europa's total H₂ neutral loss rate is 1.2 ± 0.7 kg s⁻¹

**Supporting Information:**

Supporting Information may be found in the online version of this article.

**Correspondence to:**

J. R. Szalay,
jszalay@princeton.edu

**Citation:**

Szalay, J. R., Smith, H. T., Zirnstein, E. J., McComas, D. J., Begley, L. J., Bagenal, F., et al. (2022). Water-group pickup ions from Europa-genic neutrals orbiting Jupiter. *Geophysical Research Letters*, *49*, e2022GL098111. https://doi.org/10.1029/2022GL098111

Received 28 JAN 2022

Accepted 8 APR 2022

© 2022. The Authors.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

<figure><img src="image_2.png" /></figure>

## Water-Group Pickup Ions From Europa-Genic Neutrals Orbiting Jupiter

J. R. Szalay¹ Ⓓ, H. T. Smith² Ⓓ, E. J. Zirnstein¹ Ⓓ, D. J. McComas¹ Ⓓ, L. J. Begley¹ Ⓓ, F. Bagenal³ Ⓓ, P. A. Delamere⁴ Ⓓ, R. J. Wilson³ Ⓓ, P. W. Valek⁵ Ⓓ, A. R. Poppe⁶ Ⓓ, Q. Nénon⁷ Ⓓ, F. Allegrini⁵,⁸ Ⓓ, R. W. Ebert⁵,⁸ Ⓓ, and S. J. Bolton⁵ Ⓓ

¹Department of Astrophysical Sciences, Princeton University, Princeton, NJ, USA, ²The Johns Hopkins University Applied Physics Laboratory, Laurel, MD, USA, ³Laboratory for Atmospheric and Space Physics, University of Colorado Boulder, Boulder, CO, USA, ⁴Geophysical Institute, University of Alaska Fairbanks, Fairbanks, AK, USA, ⁵Southwest Research Institute, San Antonio, TX, USA, ⁶Space Sciences Laboratory, University of California, Berkeley, CA, USA, ⁷Institut de Recherche en Astrophysique et Planétologie, CNRS-UPS-CNES, Toulouse, France, ⁸Department of Physics and Astronomy, University of Texas at San Antonio, San Antonio, TX, USA

**Abstract** Water-group gas continuously escapes from Jupiter's icy moons to form co-orbiting populations of particles or neutral toroidal clouds. These clouds provide insights into their source moons as they reveal loss processes and compositions of their parent bodies, alter local plasma composition, and act as sources and sinks for magnetospheric particles. We report the first observations of H₂⁺ pickup ions in Jupiter's magnetosphere from 13 to 18 Jovian radii and find a density ratio of H₂⁺/H⁺ = 8 ± 4%, confirming the presence of a neutral H₂ toroidal cloud. Pickup ion densities monotonically decrease radially beyond 13 R<sub>J</sub> consistent with an advecting Europa-genic toroidal cloud source. From these observations, we derive a total H₂ neutral loss rate from Europa of 1.2 ± 0.7 kg s⁻¹. This provides the most direct estimate of Europa's H₂ neutral loss rate to date and underscores the importance of both ion composition and neutral toroidal clouds in understanding satellite-magnetosphere interactions.

**Plain Language Summary** Jupiter's moons Europa, Ganymede, and Callisto all have icy surfaces which interact with their local environments. From this interaction, water-group atoms and molecules are released from the icy surfaces and orbit Jupiter as a collection of material in “neutral toroidal clouds.” The material in these toroidal clouds interact with the local charged particle environment, where neutrals in the toroidal clouds can become charged and incorporated into Jupiter's charged particle environment. Here, we highlight observations of Jupiter's charged particle environment and present the first detections of H₂⁺ in this environment. These H₂⁺ ions are shown to be originally produced from H₂ lost from Europa and the abundance of detected ions allows us to determine Europa is losing 1.2 ± 0.7 kg s⁻¹ of neutral H₂. This provides the most direct estimate of Europa's H₂ neutral loss rate to date and underscores the importance of both ion composition and neutral toroidal clouds in understanding how satellites interact with their local charged particle environments.

## 1. Introduction

Each of the outer three Galilean satellites Europa, Ganymede, and Callisto have surficial ice and atmospheric water group species (H, H₂, O, OH, H₂O, and O₂). These species are abundant to varying degrees due to differing surface compositions and atmospheric source/loss processes (e.g., Ip, 1996; Johnson, 1990, 2004; Marconi, 2007; McGrath et al., 2004; Roth, 2021; Saur et al., 1998; Smyth & Marconi, 2006), where atmospheric losses could sustain neutral toroidal clouds. Of these three satellites, only Europa's neutral toroidal cloud has been inferred from observations (Lagg et al., 2003; Mauk et al., 2003) and is assumed to be produced by sublimation and sputtering from its icy surface material (Smyth & Marconi, 2006). Depending on the lifetimes of orbiting neutrals, these clouds may or may not completely populate all longitudes throughout their parent body's orbits. Hence, we refer to them as “neutral toroidal clouds” in lieu of the historically common “neutral tori” terminology.

As observed by the Galileo spacecraft, energetic neutral atoms originating just outside Europa's orbit, born via charge exchange between a cold neutral toroidal cloud and energetic charged particle environment, suggest the densities in Europa's neutral gas toroidal cloud are comparable to Io's neutral gas toroidal cloud (Mauk
et al., 2003, 2020). Unlike near Io's orbit, however, where the dense Io plasma population quickly ionizes neutrals, the less dense plasma environment at Europa and beyond is more conducive to longer neutral lifetimes (Bagenal & Dols, 2020). Separately, observations of depletions in energetic protons that would primarily reside at low magnetic latitudes (pitch angles near 90°) were proposed to be due to charge exchange with the very same neutral Europa toroidal cloud (Lagg et al., 2003). Neither observation was able to directly determine the composition of the neutral cloud, however, ultraviolet observations ruled out a dominantly oxygen-rich cloud (Hansen et al., 2005). Europa's neutral toroidal cloud was subsequently modeled using multiple atmospheric constraints, and a cloud predominantly composed was proposed to explain both indirect cloud observations (Smith et al., 2019; Smyth & Marconi, 2006). H₂ is understood to have the highest escape rate because of its enhanced surface abundance from radiolysis, large-scale height, light mass, and non-condensability as a diatomic species from Europa (Smyth & Marconi, 2006), Ganymede (Marconi, 2007), and Callisto (Mogan et al., 2021).

While remote observations did not detect an oxygen-rich toroidal cloud and subsequent modeling efforts have suggested a dominantly H₂-rich Europa cloud, H₂ had yet to be definitively measured and direct constraints on its abundances or loss rates were not possible given the prior observations. Here, we present direct observations of H₂⁺ pickup ions sourced from a neutral H₂ toroidal cloud in the Europa-Ganymede region. We use these observations to constrain the total abundances and loss rate of neutral H₂ from Europa. In Section 2, we discuss how and where H₂⁺ pickup ions are observed from the Juno spacecraft. In Section 3, we use these observations to determine numerical densities of both H⁺ and H₂⁺ and compare these to a plasma advection model to estimate loss rates from Europa in Section 4. We conclude in Section 5 with a discussion on implications for loss rates from the other Galilean satellites and future work that could build on these results.

## 2. Observations of H₂⁺

The Juno mission (Bolton et al., 2017) carries the Jovian Auroral Distributions Experiment (JADE; McComas et al., 2017), a suite of plasma instruments that includes an ion time-of-flight (TOF) mass spectrometer. JADE's ion sensor measures the flux of ions with energy-per-charge of 10 eV/Q to 46 keV/Q and a mass-per-charge of ~1–64 AMU/Q. Juno's eccentric, nearly polar orbit precesses such that its furthest distance (apojove) is sequentially more southward after each close approach (perijove). Therefore, the location where Juno intersects Jupiter's equatorial plane moves inward each orbit, enabling ion composition observations by JADE's TOF mass-spectrometry of the equatorial plasma environments at varying radial distances (Kim et al., 2020a, 2020b).

By late-2021 up to and including its 37th perijove, Juno transited the equatorial plasma environments at the orbital distances of Europa, Ganymede, and Callisto. A JADE signal consistent with H₂⁺ is present near the equatorial region inside ~20 R<sub>j</sub> (1 R<sub>j</sub> ≡ 71,492 km). Isolating H₂⁺ in the JADE data set requires removing multiple background and foreground sources that overlap with the H₂⁺ signal in the combined TOF versus energy space. Notably, the H₂⁺ feature overlaps with H⁺ in this space, where H⁺ is more abundant, hence correctly subtracting the proton foregrounds is critical. We apply three background subtractions, and additionally remove the signature of H⁺ foreground to isolate the count rates due to H₂⁺. To remove the H⁺ signal, we use an empirical response from a unique event where Juno was connected to Io's Main Alfvén Wing (Szalay et al., 2020). During this period, large proton abundances were observed and energized throughout the entire observable range across multiple Juno instruments (Clark et al., 2020; Sulaiman et al., 2020) enabling a reliable determination of the in-flight response to protons alone. The subtraction schemes are described in detail in Supporting Information S1.

We manually identified regions where H₂⁺ could be reasonably separated, producing a conservative subset of JADE data with high-fidelity H₂⁺ data. In the region past 10 R<sub>j</sub> near the magnetic equator, an unambiguous H₂⁺ signature is observed from ~13 to 18 R<sub>j</sub>. Outside these locations, H₂⁺ signatures are still observed, but cannot be fully separated from foregrounds and backgrounds with the current techniques.

Figure 1 shows these H₂⁺ TOF observations, where each of the top TOF spectra are the superposed average count rate spectra within ±1 R<sub>j</sub> from the magnetic equator for each 1 R<sub>j</sub> radial bin shown in the bottom panel. Times
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Count rates of H₂⁺ from 13 to 18 Rⱼ (top) Ion energy per charge as a function of time-of-flight (TOF) for all species with nominal background subtraction (upper row) and with H₂⁺ isolated (lower row) for all observations within an averaging window of ±1 Rⱼ from the magnetic equator. H₂⁺ is unambiguously observed when foregrounds and backgrounds are removed, in the vicinity of the second gray TOF trace for a mass per charge of 2. The horizontal line indicates the expected detected energy for a freshly created pickup ion at the spacecraft's position (bottom) Locations of each observation in cylindrical coordinates with the +z axis are aligned with Jupiter's magnetic dipole. H₂⁺ is observed with the strongest magnitude in the inner-most detected locations.</figcaption></figure>

used in these combined TOF spectrograms are given in Table S1 in Supporting Information S1. Without background and foreground H⁺ subtraction, the counts along the *M/Q* = 2 track are dominated by H⁺ and to varying degrees, a long TOF tail from heavy ions that manifests as a horizontal bar around 2–4 keV notable in the second TOF spectrogram (14–15 Rⱼ). As indicated in Figure 1, which shows the expected energy and TOF tracks for different species, protons are detected along two separate tracks due to the internal workings of the instrument. See Text S1 in Supporting Information S1 for additional description on how this is accounted for in the H₂⁺ identification scheme.

The H₂⁺ count rates occur near and mostly below the expected local corotation energy (Text S2 in Supporting Information S1), consistent with a population of pickup ions corotating with Jupiter's magnetosphere sourced from a neutral gas moving slowly with respect to local corotation speed. Pickup ion fluxes are expected to peak in the detected energy range of 100s eV to a few keV when Juno is within 10–30 Rⱼ. Therefore, the instrument's energy range of 10 eV/Q–46 keV/Q is sufficient to capture any pickup ions generated in this region. The bottom panel of Figure 1 shows where these observations occurred along with magnetic field lines every 2 Rⱼ in equatorial distance using the JRM09 + CAN models (Connerney et al., 1981, 2018). As shown here, the largest H₂⁺ count rates are observed near the magnetic equator and closest to Jupiter over this distance range.
### 3. Numerical Densities

While the TOF data set does not have any information on directionality, JADE observed the full-sky over each ~30s spin and we can calculate a partial numerical density from the count rates as a function of energy over the JADE energy bandpass. After all foregrounds and backgrounds are subtracted, we sum count rates over all TOFs corresponding to AMU/Q between 1.5 and 2.5 to determine a total count rate $R_{obs}$ as a function of energy. JADE instantaneously observes an angular range of 270° extending from the anti-sunward spin axis, such that for each spacecraft rotation it records counts from a total angular extent of 6π sr, double-counting half the sky. We must reduce $R_{obs}$ by an appropriate factor to determine the “true” average count rate $R = \eta R_{obs}$ corresponding to the 4π sr full sky (Text S3 in Supporting Information S1). Due to the instrument mounting and orbit geometry, for each $H_2^+$ observation by JADE, plasma corotating with Jupiter is predominantly observed on the hemisphere where JADE double-counts incident populations, which also gives improved counting statistics. Table S1 in Supporting Information S1 gives values of $\eta$; for nearly all data analyzed here, $\eta \approx 0.5$.

We convert count rate $R$ into phase space density $f$ via $f(v) = R / (G_{eff}^v(v)v^4)$ where $G_{eff}^v = G_{eff}^E/2$ is the energy dependent geometric factor (Kim et al., 2020a), with a factor of two between the energy geometric factor (e.g., McComas et al., 2021) and velocity geometric factor, and $v$ is the measured energy per charge converted to speed for a molecule weighing two proton masses. In turn, the number density derived from a one-dimensional phase space density is $n_{num} = 4\pi \int_0^\infty f(v)v^2 dv$. For JADE data with count rates in discrete energy bins, the numerical partial number density is given by $n_{num} = \frac{4\pi}{3} \sum (v_{i,max}^3 - v_{i,min}^3) f(v_i)$, where $i$ is the energy bin and each energy bin spans from $v_{i,min}$ to $v_{i,max}$ in velocity space and $v_{i,max} = v_{i+1,min}$. We calculate numerical densities below 10 keV/Q to avoid additional backgrounds above this energy. This energy range is broad enough to well-capture the $H_2^+$ features analyzed in this study. Finally, we must scale the numerical densities by a correction factor $n = \epsilon n_{num}$ to account for a minority portion of $H_2^+$ counts in TOF × E space that are over-subtracted in the vicinity of the $H^+$ fork. See Text S1 in Supporting Information S1 and Kim et al. (2020a) for additional description on the proton “fork.” Text S4 in Supporting Information S1 describes this scaling factor, where we use a range of $\epsilon = 1.0$, corresponding to no over-subtraction, to $\epsilon = 1.6$, a conservative upper bound on the correction factor.

Figure 2 shows these numerical densities for $H^+$ and $H_2^+$ along with average values every 1 R$_J$ for $\epsilon = 1.2$. We include the $H^+$ numerical densities to compare the method to model expectations, as there are no published estimates for the densities of $H_2^+$ pickup ions in the Jovian magnetosphere. The dark orange dashed line in Figure 2 shows the average expected Jovian proton density from an empirical model based on a reanalysis of Voyager data (Bodisch et al., 2017). Numerical densities for $H^+$ are commensurate with previous JADE TOF-derived densities in this region (Kim et al., 2020), and are also consistent with empirical model results for protons, lending confidence that the same numerical method applied on $H_2^+$ yields accurate estimates of the partial densities in these locations.

Using these methods, we find the $H_2^+$ ion densities monotonically decrease from 13 to 18 RJ, where $H_2^+/H^+ = 8 \pm 4\%$ for $\epsilon = 1.0–1.6$. Such a low density ratio is consistent with the lack of detection by Voyager or Galileo plasma observations, which did not have the sensitivity or TOF mass-spectroscopy capabilities that JADE has. We do not see any appreciable enhancement in the vicinity of Ganymede's orbit, although we did not include data taken very near to Ganymede itself. In the next section, we compare these densities with expectations from pickup ions generated by Europa's neutral toroidal cloud.

### 4. Comparison With a Europa Toroidal Cloud Source

For an inclusive range of plasma and solar ultraviolet conditions near Europa's orbit, electron impact ionization of $H_2$ neutrals, $H_2 + e^- \rightarrow H_2^+ + 2e^-$, is expected to be the dominant loss process for neutral $H_2$ (Smith et al., 2019). $H_2^+$ pickup ions generated by this process will be radially transported outwards. The integrated transport time from 6 to 10 R$_J$ is 10–40 days, while the transport time from 10 R$_J$ to the outer magnetosphere (~50 R$_J$) is on the order of a few days with radial outflow speeds approaching 10s of km s$^{-1}$ (Bagenal & Delamere, 2011). Additionally, a reanalysis of Galileo plasma data from PLS showed the radial gradient of flux tube content beyond 10 R$_J$ is significantly reduced, precluding significant diffusion (Bagenal et al., 2016). Hence, we follow the assertion
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Partial density of H₂⁺ pickup ions (gray) and H⁺ (orange), with average values every 1 Rⱼ for ε = 1.2. The blue curves show expected H₂⁺ pickup ion densities for a neutral Europa loss H₂ rate of 6.3 kg s⁻¹ (Smyth & Marconi, 2006). The purple curves show the best fit solution to the observed H₂⁺ partial densities. These observations are consistent with a Europa neutral H₂ loss rate of 1.2 ± 0.7 kg s⁻¹. The orange dashed line shows an empirical model for proton densities based on Voyager observations (Bodisch et al., 2017).</figcaption></figure>

that while slow diffusive transport is valid in the inner magnetosphere ($<$10 Rⱼ), transport in the middle magnetosphere is better modeled with advection (Ng et al., 2018).

At these low magnetic latitudes, we assume H₂⁺ pickup ions are transported radially outward via cylindrically symmetric advection. From conservation of mass, we solve the 1D time-independent continuity equation in cylindrical coordinates $\frac{1}{\rho}\frac{\partial(n\rho v_{\rho})}{\partial\rho} = P$, where ρ is cylindrical radial distance, n is the number density of H₂⁺, v_ρ is the radial transport speed, and P is the production rate of pickup ions. This is solved using explicit finite differencing via $F_{i+1} = F_i + P_i \rho_i \Delta\rho$, where F = nρv_ρ, and a small value of Δρ = 0.025 Rⱼ is chosen such that the results are insensitive to the exact value of radial step size. We use a radial transport speed profile (Bagenal & Delamere, 2011) corresponding to a total Jovian plasma source of 0.7 ton/s, which represents the typical magnetospheric configuration (Delamere et al., 2004; Nerney & Bagenal, 2020). For the pickup ion input source, we approximate the modeled H₂⁺ pickup ion input given in Figure 10 of Smyth and Marconi (2006) as $P_E(R) = P_{0E}e^{-\frac{(R-R_E)^2}{2\sigma_E^2}}$ where $P_{0E}$ is the production rate of H₂⁺ and $R_E = 9.4$ Rⱼ, and $\sigma_E = 0.4$ Rⱼ. We assume transport occurs within a vertical extent of 2 Rⱼ.

To link the production rate of H₂⁺ pickup ions to neutral H₂ loss from Europa, we use existing estimates that 50%–77% of neutral H₂ in the toroidal cloud is lost to H₂⁺ pickup ions (Smith et al., 2019). The purple curves in Figure 2 show the best fit advection solutions to the binned H₂⁺ data. Blue curves show expectations for the initial estimates of H₂⁺ production, discussed in the following section. As shown here, the data are consistent with an advecting Europa-genic toroidal cloud source of H₂⁺ with a monotonically decreasing radial density profile
<caption>**Table 1**</caption>

Estimates for the Rate of Neutral H₂ Mass Loss (kg s⁻¹) and Number Loss (s⁻¹) From Europa

<table><thead><tr><th>H₂ Loss (kg s⁻¹)</th><th>H₂ Loss (s⁻¹)</th><th>Method</th><th>Reference</th></tr></thead><tbody><tr><td>6.3</td><td>1.9 × 10<sup>27</sup></td><td>Sputtering and atmosphere model</td><td>Smyth and Marconi (2006)</td></tr><tr><td>1.5</td><td>4.6 × 10<sup>26</sup></td><td>Sputtering model</td><td>Plainaki et al. (2012)</td></tr><tr><td>0.7</td><td>2 × 10<sup>26</sup></td><td>Sputtering model</td><td>Cassidy et al. (2013)</td></tr><tr><td>1.9</td><td>5.7 × 10<sup>26</sup></td><td>Plasma-atmosphere charge exchange model</td><td>Dols et al. (2016)</td></tr><tr><td>1.8</td><td>5.5 × 10<sup>26</sup></td><td>Sputtering and atmosphere model</td><td>Vorburger and Wurz (2018)</td></tr><tr><td>1.2 ± 0.7</td><td>3.6 ± 2 × 10<sup>26</sup></td><td>Pickup ion observations and advection model</td><td>This work</td></tr></tbody></table>

past ~13 R<sub>J</sub>. These fits are consistent with a total H₂⁺ pickup ion production rate of 0.7 ± 0.3 kg s⁻¹ for the full range of ε = 1.0–1.6. Hence, dividing this range by 0.5–0.77, a neutral H₂ production rate within a 1-sigma range of 1.2 ± 0.7 kg s⁻¹ from the Europa toroidal cloud, lost from Europa's surface, is derived from the H₂⁺ pickup densities observed. Table 1 gives these estimates along with those in the literature to date.

**5. Discussion and Conclusions**

Using plasma observations from the JADE instrument onboard Juno, we detected a persistent population of H₂⁺ pickup ions in Jupiter's magnetosphere near the magnetic equator from 13 to 18 R<sub>J</sub>. Using 1D numerical moments, we determine densities of H₂⁺ to be in the range of 5 × 10⁻³ to 1 × 10⁻¹ cm⁻³ throughout this region. The densities monotonically decrease as a function of radial distance and do not show any clear enhancement near Ganymede's orbit. The overall radial profile and total densities are consistent with pickup ions generated by a Europa-genic neutral toroidal cloud with a loss rate of H₂ from Europa of 1.2 ± 0.7 kg s⁻¹. Hence, these observations confirm the existence and composition of Europa's neutral H₂ toroidal cloud, which is expected to be the dominant component (Smith et al., 2019). Figure 3 shows a summary of the Jovian magnetosphere highlighting the existence of H₂ neutrals and H₂⁺ pickup ions associated with Europa.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Overview of Jupiter's magnetosphere in the vicinity of the Galilean satellites. H₂⁺ pickup ions (blue) originate from Europa's neutral toroidal cloud (brighter blue near Europa). Io and Europa contribute plasma pickup ions of different compositions to Jupiter's magnetosphere. Alfvén wings connected to the moons due to their interaction with corotating plasma are also shown in gray.</figcaption></figure>
The total loss rate of neutral H₂ from Europa estimated here is a factor of 3–13 lower than the original estimated value of 6.3 kg s⁻¹ (Smyth & Marconi, 2006), where the blue curves in Figure 2 show corresponding H₂⁺ densities for this estimate. This comparison also shows that were the source rates from Europa to be this large, H₂⁺/H⁺ ≥ 30% such that H₂⁺ pickup densities would nearly rival those of protons in the magnetosphere. This is certainly not the case as previous measurements would have identified such a population (Voyager, Galileo, Juno, and New Horizons), and we find the density ratio to be H₂⁺/H⁺ = 8 ± 4%.

While we derive our estimates from direct observations of the expected primary loss mechanism for Europa's neutral toroidal cloud (pickup ions), the values calculated in this work of 0.5–1.9 kg s⁻¹ are commensurate with more recent estimates of 0.7–1.9 kg s⁻¹ (Cassidy et al., 2013; Dols et al., 2016; Plainaki et al., 2012; Vorburger & Wurz, 2018) based on constraints of the loss rate from the surface due to sputtering inputs (Table 1). Both the sputtering estimates and our pickup ion calculations here require a synthesis of data and modeling to derive estimates; however, given that the two very different methods of estimating Europa's neutral H₂ loss rate yield similar results, we can be reasonably confident that the loss rate is on the order of ~1 kg s⁻¹.

We refrain from estimating neutral H₂ densities. Linking pickup ion densities to neutral cloud densities requires an additional set of model-dependent assumptions above those already used in this analysis. Namely, the absolute lifetime of the H₂ neutrals needs to be estimated to constrain neutral densities and doing so would incorporate additional model-dependent uncertainty. Constraining neutral toroidal cloud densities using accurate lifetime estimates requires more complex modeling and is suggested as a future line of study.

Protons are also produced as pickup ions from Europa's neutral toroidal cloud, where the production rate of H⁺ is expected to be ~15 times lower than that for H₂⁺ (Smyth & Marconi, 2006). Therefore, we anticipate an input rate of protons to the Jovian magnetosphere from Europa to be less than 0.1 kg s⁻¹. This value is at least an order of magnitude less input than the 2.5–13 kg s⁻¹ required to sustain magnetospheric proton abundances inside 30 R<sub>J</sub> (Bodisch et al., 2017), further pointing to Jupiter as the likely source for magnetospheric protons, where Juno observations have shown proton outflow contributes at least 1–5 kg s⁻¹ to the magnetosphere (Szalay et al., 2021).

The largest source of uncertainty in our analysis lies in the use of a scaling factor to correct for over-subtraction of H₂⁺ counts (Text S4 in Supporting Information S1). Our conservative range of scaling factors leads to a large error bar in H₂ mass loss from Europa that could be narrowed with improved methods to isolate H₂⁺ count rates in JADE measurements. We also do not estimate densities in the vicinity of Europa's orbit in this study. Inside 13 RJ, particularly near Europa's orbit, two issues make measuring H₂⁺ challenging. First, JADE experiences particularly large fluxes of penetrating radiation, adding to the background that must be removed. Second, and more importantly, locally produced pickup ions are expected to be observed around 0.5 keV and below, where the overlap between the H⁺ “fork” in TOF and H₂⁺ in the JADE TOF data is particularly pronounced and ε may be significantly larger near Europa using the current method. Due to these reasons, we do not include observations in this region; however, future efforts involving forward-modeling and/or fitting may help reveal the H₂⁺ content in this region.

Ganymede is also a source of neutrals (Huang & Siscoe, 1987) and like Europa, is expected to predominantly shed H₂ to generate a population of co-orbiting neutrals (Marconi, 2007). However, such a neutral source would be more dispersed in a larger total volume as the moon orbits at 15 R<sub>J</sub> compared to Europa's 9.4 R<sub>J</sub> distance. Additionally, the radial transport speed is approximately an order of magnitude higher at 15 R<sub>J</sub> compared to 9.4 R<sub>J</sub> (Bagenal & Delamere, 2011), leading to shorter pickup ion lifetimes and densities. Hence even with equal neutral loss source rates from the two bodies, lower densities of neutrals and significantly lower peak densities of subsequent pickup ions would be expected near each of their orbits. While the monotonically decreasing densities as a function of radial distance are entirely consistent with a single source from Europa, we can very roughly constrain the contribution from Ganymede. Assuming Ganymede's neutral toroidal cloud has the same pickup

ion profile $P_G(R) = P_{0G}e^{-(R-R_G)/(2\sigma_G)}$, where $R_G = 15.0$ R<sub>J</sub> and $\sigma_R = \sigma_E = 0.4$ R<sub>J</sub>, we find source rates above ~50% of Europa's would present a detectable departure from the approximately power-law-like radial profile observed for H₂⁺, and therefore do not expect Ganymede's toroidal cloud H₂⁺ pickup ion production rate to be larger than 50% of Europa's. If the conversion of H₂ neutrals to H₂⁺ pickup ions near Ganymede's orbit is similar to Europa's orbit, H₂ neutral production at Ganymede would also be lower than ~50% of Europa's. This effect is even more dramatic for Callisto at 26.3 R<sub>J</sub>, which also likely loses neutrals (Mogan et al., 2021) to an even more dispersed toroidal
cloud. In the vicinity of Callisto's orbit, the current foreground and background removal techniques were not able to unambiguously identify a clear H₂⁺ signal, however, future methods may help constrain the H₂⁺ densities in this region and in the region near Europa.

These plasma measurements could directly tie into observations by the high-energy particle detector onboard Juno, JEDI (Mauk et al., 2017), which can search for pitch-angle depletions (e.g., Kollman et al., 2016; Lagg et al., 2003; Nénon & André, 2019) and cross-compare with inferred neutral densities from these pickup ion observations. As the observed H₂⁺ pickup ions likely originate from neutrals produced via radiolysis due to energetic particle bombardment and subsequent thermal desorption (e.g., Johnson et al., 2004), JEDI measurements of this charged particle environment could also be cross-compared to these pickup observations and loss rate estimates. Constraints such as these on satellite water species loss processes and composition provide key information on the chemistry and interaction with their local charged particle environments, where present and future spacecraft will continue to shed light on these important processes over a broad range of energies (Futaana et al., 2015; Grasset et al., 2013; Paranicas et al., 2021; Phillips & Pappalardo, 2014).

Thus, the JADE observations presented here for the first time directly measure ions originating from a neutral H₂ toroidal cloud at Jupiter, prove the cloud provides an additional plasma source in Jupiter's magnetosphere, and provide the most direct constraints on Europa's loss of neutral H₂ via observations of the neutral toroidal cloud's primary loss process: pickup ions. Future analyses of JADE H₂⁺ pickup ion data may further constrain loss rates and satellite-magnetospheric interactions at Europa, Ganymede, and Callisto. These results underscore the importance of both ion composition and neutral toroidal clouds in understanding the evolution of planetary bodies, not just in the Jovian system, but at all outer planetary systems (e.g., Richardson et al., 1986; Smith & Richardson, 2021), and at exoplanetary systems.

## Data Availability Statement

The JNO-J/SW-JAD-3-CALIBRATED-V1.0 data presented in this manuscript and Supporting Information S1, doi:10.1007/s11214-013-9990-9, can be obtained from the Planetary Data System (PDS) at https://pds-ppi.igpp.ucla.edu/mission/JUNO/JNO/JAD.

## References

* Bagenal, F., & Delamere, P. A. (2011). Flow of mass and energy in the magnetospheres of Jupiter and Saturn. *Journal of Geophysical Research*, **116**(A), A05209. https://doi.org/10.1029/2010JA016294

* Bagenal, F., & Dols, V. (2020). The space environment of Io and Europa. *Journal of Geophysical Research: Space Physics*, **125**(5), 8241–8257. https://doi.org/10.1029/2019ja027485

* Bagenal, F., Wilson, R. J., Siler, S., Paterson, W. R., & Kurth, W. S. (2016). Survey of Galileo plasma observations in Jupiter's plasma sheet. *Journal of Geophysical Research: Planets*, **121**, 871–894. https://doi.org/10.1002/2016JE005009

* Bodisch, K. M., Dougherty, L. P., & Bagenal, F. (2017). Survey of Voyager plasma science ions at Jupiter: 3. Protons and minor ions. *Journal of Geophysical Research: Space Physics*, **122**(8), 8277–8294. https://doi.org/10.1002/2017ja024148

* Bolton, S. J., Lunine, J., Stevenson, D., Connerney, J. E. P., Levin, S., Owen, T. C., et al. (2017). The Juno Mission. *Space Science Reviews*, **213**(1–4), 1–37. https://doi.org/10.1007/s11214-017-0429-6

* Cassidy, T. A., Paranicas, C. P., Shirley, J. H., Dalton, J. B. I., Teolis, B. D., Johnson, R. E., et al. (2013). Magnetospheric ion sputtering and water ice grain size at Europa. *Planetary and Space Science*, **77**, 64–73. https://doi.org/10.1016/j.pss.2012.07.008

* Clark, G., Mauk, B. H., Kollmann, P., Szalay, J. R., Sulaiman, A. H., Gershman, D. J., et al. (2020). Energetic proton acceleration associated with Io's footprint tail. *Geophysical Research Letters*, **47**, e2020GL090839. https://doi.org/10.1029/2020GL090839

* Connerney, J. E. P., Ağuna, M. H., & Ness, N. F. (1981). Modeling the Jovian current sheet and inner magnetosphere. *Journal of Geophysical Research*, **86**(A10), 8370–8384. https://doi.org/10.1029/JA086iA10p08370

* Connerney, J. E. P., Kotsiaros, S., Oliversen, R. J., Espley, J. R., Joergensen, J. L., Joergensen, P. S., et al. (2018). A new model of Jupiter's magnetic field from Juno's first nine orbits. *Geophysical Research Letters*, **45**, 2590–2596. https://doi.org/10.1002/2018GL077312

* Delamere, P. A., Steffl, A., & Bagenal, F. (2004). Modeling temporal variability of plasma conditions in the Io torus during the Cassini era. *Journal of Geophysical Research: Space Physics*, **109**(A), A10216. https://doi.org/10.1029/2003JA010354

* Dols, V. J., Bagenal, F., Cassidy, T. A., Crary, F. J., & Delamere, P. A. (2016). Europa's atmospheric neutral escape: Importance of symmetrical O₂ charge exchange. *Icarus*, **264**, 387–397. https://doi.org/10.1016/j.icarus.2015.09.026

* Futaana, Y., Barabash, S., Wang, X.-D., Wieser, M., Wieser, G. S., Wurz, P., et al. (2015). Low-energy energetic neutral atom imaging of Io plasma and neutral tori. *Planetary and Space Science*, **108**(C), 41–53. https://doi.org/10.1016/j.pss.2014.12.022

* Grasset, O., Dougherty, M. K., Coustenis, A., Bunce, E. J., Erd, C., Titov, D., et al. (2013). JUpiter ICy moons Explorer (JUICE): An ESA mission to orbit Ganymede and to characterise the Jupiter system. *Planetary and Space Science*, **78**, 1–21. https://doi.org/10.1016/j.pss.2012.12.002

* Hansen, C. J., Shemansky, D. E., & Hendrix, A. R. (2005). Cassini UVIS observations of Europa's oxygen atmosphere and torus. *Icarus*, **176**(2), 305–315. https://doi.org/10.1016/j.icarus.2005.02.007

* Huang, T. S., & Siscoe, G. L. (1987). Types of planetary tori. *Icarus*, **70**, 366–378. (ISSN 0019-1035). https://doi.org/10.1016/0019-1035(87)90142-4

* Ip, W. H. (1996). Europa's oxygen exosphere and its magnetospheric interaction. *Icarus*, **120**(2), 317–325. https://doi.org/10.1006/icar.1996.0052

* Johnson, R. E. (1990). *Charged particle interactions with atmospheres and surfaces*. Springer-Verlag.

## Acknowledgments

The authors would like to thank the many JADE and Juno team members that made these observations possible. We thank M. Imai for producing the JRM09 magnetic field line integrations, T. Kim for discussions on JADE TOF analysis, A. Sulaiman for discussions of neutral cloud comparisons, P. Swaczyna for discussions on pickup ions, B. Smith for his help creating the Jovian overview figure, and B. Bonfond for helpful comments on the overview figure. We acknowledge Juno's NASA contract NNM06AA75C and the NASA Juno Participating Scientist Program. A.R.P. acknowledges NASA NFDAP grant #80NSSC21K0823.
Johnson, R. E., Carlson, R. W., Cooper, J. F., Paranicas, C., Moore, M. H., & Wong, M. C. (2004). Radiation effects on the surfaces of the Galilean satellites. *Jupiter: The planet, satellites and magnetosphere*, 485–512.

Kim, T. K., Ebert, R. W., Valek, P. W., Allegrini, F., McComas, D. J., Bagenal, F., et al. (2020a). Method to derive ion properties from Juno JADE including abundance estimates for O⁺ and and S²⁺. *Journal of Geophysical Research: Space Physics*, 125(2), e026169. https://doi.org/10.1029/2018ja026169

Kim, T. K., Ebert, R. W., Valek, P. W., Allegrini, F., McComas, D. J., Bagenal, F., et al. (2020b). Survey of ion properties in Jupiter's plasma sheet: Juno JADE-I observations. *Journal of Geophysical Research: Space Physics*, 125(4), e27696. https://doi.org/10.1029/2019JA027696

Kollmann, P., Paranicas, C. P., Clark, G., Roussos, E., Lagg, A., & Krupp, N. (2016). The vertical thickness of Jupiter's Europa gas torus from charged particle measurements. *Geophysical Research Letters*, 43(1), 9425–9433. https://doi.org/10.1002/2016GL070326

Lagg, A., Krupp, N., Woch, J., & Williams, D. J. (2003). In-situ observations of a neutral gas torus at Europa. *Geophysical Research Letters*, 30(1), 1556. https://doi.org/10.1029/2003GL017214

Marconi, M. L. (2007). A kinetic model of Ganymede's atmosphere. *Icarus*, 190(1), 155–174. https://doi.org/10.1016/j.icarus.2007.02.016

Mauk, B. H., Clark, G., Allegrini, F., Bagenal, F., Bolton, S. J., Connerney, J. E. P., et al. (2020). Juno energetic neutral atom (ENA) remote measurements of magnetospheric injection dynamics in Jupiter's Io torus regions. *Journal of Geophysical Research: Space Physics*, 125(5), e27964. https://doi.org/10.1029/2020JA027964

Mauk, B. H., Haggerty, D. K., Jaskulek, S. E., Schlemm, C. E., Brown, L. E., Cooper, S. A., et al. (2017a). The Jupiter Energetic Particle Detector Instrument (JEDI) investigation for the Juno mission. *Space Science Reviews*, 213(1–4), 289–346. https://doi.org/10.1007/s11214-013-0025-3

Mauk, B. H., Mitchell, D. G., Krimigis, S. M., Roelof, E. C., & Paranicas, C. P. (2003). Energetic neutral atoms from a trans-Europa gas torus at Jupiter. *Nature*, 421(6), 920–922. https://doi.org/10.1038/nature01431

McComas, D. J., Alexander, N., Allegrini, F., Bagenal, F., Beebe, C., Clark, G., et al. (2017). The Jovian Auroral Distributions Experiment (JADE) on the Juno mission to Jupiter. *Space Science Reviews*, 213(1), 547–643. https://doi.org/10.1007/s11214-013-9990-9

McComas, D. J., Swaczyna, P., Szalay, J. R., Zirnstein, E. J., Rankin, J. S., Elliott, H. A., et al. (2021). Interstellar pickup ion observations halfway to the termination shock. *The Astrophysical Journal: Supplement Series*, 254(1), 19.

McGrath, M. A., Lellouch, E., Strobel, D. F., Feldman, P. D., & Johnson, R. E. (2004). Satellite atmospheres. In *Jupiter: The planet, satellites and magnetosphere*, pp. 457–483.

Mogan, S. R. C., Tucker, O. J., Johnson, R. E., Vorburger, A., Galli, A., Marchand, B., et al. (2021). A tenuous, collisional atmosphere on Callisto. *Icarus*, 368, 114597. https://doi.org/10.1016/j.icarus.2021.114597

Nénon, Q., & André, N. (2019). Evidence of Europa neutral gas Torii from energetic sulfur ion measurements. *Geophysical Research Letters*, 46(7), 3599–3606. https://doi.org/10.1029/2019GL082200

Nerney, E. G., & Bagenal, F. (2020). Combining UV spectra and physical chemistry to constrain the hot electron fraction in the Io plasma torus. *Journal of Geophysical Research: Space Physics*, 125(4), e27458. https://doi.org/10.1029/2019JA027458

Ng, C. S., Delamere, P. A., Kaminker, V., & Damiano, P. A. (2018). Radial transport and plasma heating in Jupiter's magnetodisc. *Journal of Geophysical Research: Space Physics*, 123(8), 6611–6620. https://doi.org/10.1029/2018JA025345

Paranicas, C., Szalay, J. R., Mauk, B. H., Clark, G., Kollmann, P., Haggerty, D. K., et al. (2021). Energy spectra near Ganymede from Juno data. *Geophysical Research Letters*, 48, https://doi.org/10.1029/2021GL093021

Phillips, C. B., & Pappalardo, R. T. (2014). Europa clipper mission concept: Exploring Jupiter's ocean moon. *Eos, Transactions American Geophysical Union*, 95(20), 165–167.

Plainaki, C., Milillo, A., Mura, A., Orsini, S., Massetti, S., & Cassidy, T. (2012). The role of sputtering and radiolysis in the generation of Europa exosphere. *Icarus*, 218(2), 956–966. https://doi.org/10.1016/j.icarus.2012.01.023

Richardson, J. D., Eviatar, A., & Siscoe, G. L. (1986). Satellite tori at Saturn. *Journal of Geophysical Research*, 91, 8749–8755. (ISSN 0148-0227). https://doi.org/10.1029/JA091iA08p08749

Roth, L. (2021). A stable H₂O atmosphere on Europa's trailing hemisphere from HST Images. *Geophysical Research Letters*, 48(2), e94289. https://doi.org/10.1029/2021GL094289

Saur, J., Strobel, D. F., & Neubauer, F. M. (1998). Interaction of the Jovian magnetosphere with Europa: Constraints on the neutral atmosphere. *Journal of Geophysical Research*, 103(E), 19947–19962. https://doi.org/10.1029/97je03556

Smith, H. T., Mitchell, D. G., Johnson, R. E., Mauk, B. H., & Smith, J. E. (2019). Europa neutral torus confirmation and characterization based on observations and modeling. *The Astrophysical Journal*, 871(1). https://doi.org/10.3847/1538-4357/aaed38

Smith, H. T., & Richardson, J. D. (2021). The 3D structure of Saturn magnetospheric neutral tori produced by the Enceladus plumes. *Journal of Geophysical Research: Space Physics*, 126. https://doi.org/10.1029/2020JA028775

Smyth, W. H., & Marconi, M. L. (2006). Europa's atmosphere, gas tori, and magnetospheric implications. *Icarus*, 181(2), 510–526. https://doi.org/10.1016/j.icarus.2005.10.019

Sulaiman, A. H., Hospodarsky, G. B., Elliott, S. S., Kurth, W. S., Gurnett, D. A., Imai, M., et al. (2020). Wave-particle interactions associated with Io's Auroral footprint: Evidence of Alfvén, ion cyclotron, and whistler modes. *Geophysical Research Letters*, 47(22), 261–310. https://doi.org/10.1029/2020GL088432

Szalay, J. R., Allegrini, F., Bagenal, F., Bolton, S. J., Bonfond, B., Clark, G., et al. (2020). A new framework to explain changes in Io's footprint tail electron fluxes. *Geophysical Research Letters*, 47(1), e89267. https://doi.org/10.1029/2020GL089267

Szalay, J. R., Allegrini, F., Bagenal, F., Bolton, S. J., Clark, G., Connerney, J. E. P., et al. (2021). Proton outflow associated with Jupiter's Auroral processes. *Geophysical Research Letters*, 48(1), 1–10. https://doi.org/10.1029/2020GL091627

Vorburger, A., & Wurz, P. (2018). Europa's ice-related atmosphere: The sputter contribution. *Icarus*, 311, 135–145. https://doi.org/10.1016/j.icarus.2018.03.022