Energy mapping of Jupiter’s auroral electrons from Juno/UVS data using a new H 2 UV emission model

Context. Juno, which studies the Jovian system, continues to expand our knowledge of Jupiter’s magnetosphere and its environment. Thanks to onboard instruments such as Jupiter Energetic Particle Detector Instrument (JEDI) and Jovian Auroral Distributions Experiment (JADE), in situ measurements have allowed us to derive a realistic representation of charged particle energy distributions precipitating in the auroral regions. Because of the distance between Juno’s measurement location and the position of impact of the charged particles, where auroral emissions are produced, these energetic distributions of magnetospheric particles are likely to be affected by various phenomena such as wave-particle interactions on their way from Juno to the atmosphere. These processes can accelerate or decelerate the particles, changing their average energies. Hence, the energy distributions of particles measured at Juno’s altitude are likely different from those at auroral altitudes. Aims. In this study we develop a UV emission model, combined with an electron transport model, that allows us to relate the auroral emission spectra of H 2 molecules with the energy distribution of impinging electrons. Methods. Thanks to observations of the Jovian aurora by the Ultraviolet Spectrograph (UVS) on board Juno, we determined the characteristic energies of electrons precipitating in auroral regions during perijove 32. We modeled the relationship between color ratio (CR) and the characteristic energy of precipitating electrons. Initially, we considered mono-energetic electron fluxes. In a second step, we considered fluxes governed by a kappa distribution. Results. We derived characteristic energy maps for electrons precipitating in Jupiter’s auroral regions. In comparison with similar previous studies based on Space Telescope Imaging Spectrograph on board Hubble Space Telescope (HST/STIS) observations, we find that modeling the CR with a mono-energetic distribution leads to a systematic underestimation of the average energy of electrons precipitating in the auroral regions by a factor of 3–5. Conclusions. In this study we show that it is possible to derive a more realistic estimate of electron energy flux distributions at auroral altitudes.


Introduction
Observations of Jupiter in the UV range have revealed the existence of extremely bright polar auroral phenomena.The first observational evidence of these phenomena was provided by observations of UV emissions from atomic hydrogen (H Lymanα emission) and molecular hydrogen (Lyman band emissions: B 1 Σ + u → X 1 Σ + g and Werner bands: C 1 Π u → X 1 Σ + g ) during the flyby of the Voyager 1 spacecraft (Broadfoot et al. 1979).In addition, Jovian aurorae have been extensively studied by the International Ultraviolet Explorer (IUE) for approximately two decades (Clarke et al. 1980;Livengood et al. 1992;Gladstone & Skinner 1989;Harris et al. 1996).These observations, in the mid-UV range (between 120 nm and 170 nm), have allowed the characterization of the power of the aurorae, and also the study of their structure, variability, and intensity.Thanks to these spectral measurements, the first models of UV auroral emissions were developed (Yung et al. 1982;Gladstone & Skinner 1989).This demonstrated that the direct excitation of molecular hydrogen by electrons and absorption of CH 4 below 140 nm could, overall, reproduce UV auroral spectra well in the range of 120-170 nm.
The Hubble Space Telescope (HST) has dramatically contributed to the study of Jovian aurorae, thanks to observations of auroral structures by the Faint Object Camera (FOC; Dols et al. 1992;Gérard et al. 1993Gérard et al. , 1994;;Prangé et al. 1998) and Wide Field Planetary Camera 2 (WFPC2; Clarke et al. 1996Clarke et al. , 1998;;Grodent et al. 1997), as well as UV spectral observations by the Space Telescope Imaging Spectrograph (STIS) on board HST (HST/STIS; Gustin et al. 2002).Despite these numerous studies, the morphology of the Jovian aurora remains very difficult to describe in an exhaustive way as their structure is complex and includes spatially and temporally variable substructures (see Grodent 2015).However, these various observations have made it possible to characterize, in a simple way, the morphology of Jupiter's aurorae: the main auroral emissions form a partially closed oval, with highly variable structures in the polar region inside the main emissions and equatorward emission structures outside, including the footprints of Io, Europa, and Ganymede.
The auroral emission on Jupiter occurs due to the interaction between its magnetosphere and atmosphere.This interaction leads to the precipitation of energetically charged particles along magnetic field lines.In these regions, electrons are the primary species that precipitate, as stated by Rego et al. (2001).In the UV domain [80 nm, 180 nm], Jupiter's auroral spectral emission is dominated by H Lyman-α emission and the de-excitation of H 2 molecules by electronic transitions: B1 Σ + u −→ X 1 Σ + g (Lyman bands) and the R 1 , P, and Q branches of the C 1 Π u −→ X 1 Σ + g transition (Werner bands).Moreover, this spectral emission is strongly influenced by the presence of hydrocarbons in Jupiter's atmosphere, including methane, which mainly absorbs UV photons below 140 nm.
The main studies resulting from the various spectral observations of Jupiter's auroral regions are those concerning the energy characterization of the electrons precipitating in these regions.Although numerous probes have studied Jupiter's magnetosphere at the spacecraft altitude, such as Voyager (Russell 1993), Ulysses (Zarka 1998), Galileo, Cassini-Huygens (Hansen et al. 2004), and the New Horizon flyby (Krupp 2007), measurements and observations of magnetospheric plasma have never characterized the energy distributions of the electrons that precipitate immediately above the polar aurora's altitudes.There are complex processes, such as acceleration by inertial Alfvén waves (Hess et al. 2010(Hess et al. , 2013;;Saur et al. 2018) and by whistler waves (Elliott et al. 2018), and ion and electron inverted-V structures (Mauk et al. 2017(Mauk et al. , 2018;;Clark et al. 2017), that alter the energy flux distribution of electrons between the altitude of measurement and the high-latitude ionospheric regions.Thus, to measure the shape of the energy flux distributions of electrons precipitating in auroral regions, only low altitude measurements can be effective.
For the present study, we used a combination of spectral observations and modeling to study these characteristics.The method we adopted is phenomenological; it was first proposed by Yung et al. (1982) and is based on the far-UV (FUV) color ratio (CR) of the auroral H 2 emission spectrum.Our aim is to take advantage of the wavelength-dependent absorption of auroral emission by hydrocarbons such as CH 4 , which is the third most abundant molecule in Jupiter's stratosphere (e.g., Moses et al. 2005;Hue et al. 2018).Thus, by considering the unabsorbed part of the H 2 emission spectrum, we can derive a CR that allows us to characterize the energy distributions of the electron fluxes that precipitate in these regions, as well as their characteristic energies.In the case of absorption by methane, this ratio is defined by CR = I(155 nm-162 nm) I(123 nm-130 nm) (Gustin et al. 2013), where I(λ min − λ max ) = λ max λ min I dλ and I is the spectrum flux intensity and where the ranges [123-130 nm] and [155-162 nm] represent, respectively, a range absorbed by CH 4 and the unabsorbed part of the spectrum.This allows us to infer, for a fixed emission angle, that over the range [123-130 nm] increasing spectral absorption means that electrons penetrate deeper into the atmosphere before being thermalized, and are therefore more energetic.This method is very advantageous since it does not require the use of absolute spectra of auroral emission as CR can be measured using only arbitrary units within the same spectra.
The relationship between the CR and the characteristic electron energy E 0 (CR(E 0 )) is monotonic, and is modeled using a combination of an electronic transport model and a H 2 UV emission model in the auroral regions.This supposes modeling the excitation of H 2 molecules by electron collisions, before calculating their de-excitation from rovibrational levels producing UV emissions.Each modeled spectrum allows the CR to be linked to the characteristic energy of the electrons precipitated in the atmospheric model.Finally, by varying the characteristic energy of the electrons in the transport model, we build the relationship CR(E 0 ), which is compared with the observed CR to estimate E 0 .
Several studies (e.g., Trafton et al. 1994Trafton et al. , 1998;;Grodent et al. 2001;Gustin et al. 2002Gustin et al. , 2016;;Ajello et al. 2005;Gérard et al. 2014) have used this technique to characterize the electron energy in auroral regions using HST observations.In most of these studies, the auroral electron transport was modeled by initial mono-energetic, Maxwellian, or kappa phenomenological flux distributions (see the example by Gustin et al. 2016).However, despite the high spectral resolution of the HST/STIS observations exploited in previous studies, the signal-to-noise ratio (S/N) was limited and the spatial coverage of the aurora was partial and highly dependent on the planet's tilt axis.Thus, it is only since Juno's arrival (Bolton et al. 2017) and the Ultraviolet Spectrograph (UVS) observations (Gladstone et al. 2017;Bonfond et al. 2017) that we have had full access to Jupiter's northern and southern local time.In addition, UVS observations are highly spatially resolved near Jupiter's closest approach, with a spectral resolution of around 0.6-2.4nm.This allows correct spectral sampling with a better S/N than the HST/STIS observations.
In this study, we map for the first time the characteristic energy of electrons precipitating in Jupiter's auroral regions using Juno/UVS observations.We have developed a new UV emission model of H 2 , inspired by the models of Dols et al. (2000), Gustin et al. (2002) and Menager (2011), in a more optimized version that takes into account nine H 2 electronic states including cascade excitation and auto-absorption in the Lyman and Werner bands.This model is now available for the community, and can be used in every electron transport model.The excited states of H 2 are calculated through the outputs of our TransPlanet electronic transport model (Stamnes & Rees 1983b;Gronoff 2009;Menager 2011;Benmahi 2022).Additionally, for the CR modeling, we modeled the electronic transport using mono-energetic initial electron flux distributions and a kappa distribution (Coumans et al. 2002;Scherer et al. 2018) derived from Jupiter Energetic Particle Detector Instrument (JEDI) measurements obtained during the first 20 perijoves (PJs; Salveter et al. 2022) of the Juno mission.
The outline of our study is as follows.We first describe the electronic transport model and the UV emission model.In the second step, we describe the Juno/UVS observations and explain the mapping of the characteristic energy method.Finally, we present our results and discussions before concluding.

Electron transport model: TransPlanet
To simulate electron precipitation in Jupiter's atmosphere, we used the TransPlanet model developed in collaboration with the A26, page 2 of 20 Benmahi, B., et al.: A&A, 685, A26 (2024) Institut de Planétologie et d'Astrophysique de Grenoble (IPAG).This transport code was first created by Lilensten et al. (1989) and was modified and improved by Blelly et al. (1996) for application to terrestrial cases (Simon et al. 2007).The model was diversified and adapted to several planets over the years.The core of this algorithm was used in Trans-Mars (Witasse et al. 2002(Witasse et al. , 2003;;Simon et al. 2009;Nicholson et al. 2009), Trans-Venus (Gronoff et al. 2007(Gronoff et al. , 2008)), Trans-Titan (Lilensten et al. 2005a,b;Gronoff et al. 2009a,b), Trans-Uranus, Trans-Jupiter (Menager et al. 2010), Aeroplanets (Gronoff et al. 2012a(Gronoff et al. ,b, 2014)), and recently Trans-Planet (Benmahi 2022).The Trans* code core is divided into two parts: a kinetic part that calculates the interaction of precipitating electrons with atmospheric particles and a fluid part that is implemented only in the TRANSCAR and TRANS4 versions (Lathuillere et al. 1997;Simon et al. 2007), which uses a 13-moment fluid closure description, calculating, among other parameters, the number density, velocity, heat flux, and plasma temperature (of electrons and ions).Thus, in the Trans* version we used for this study, there is no fluid part.Compared with the various existing electron transport codes, Trans* allows multi-stream modeling of electron transport with electron scattering over a wide range of electron energies and pitch angles.The electronic transport model we use in this study is detailed in Appendix A. The uncertainties in that class of models has been studied in Gronoff et al. (2012a,b) and highlighted that the cross sections are one of the major sources of uncertainties.Efforts have been made to improve cross section datasets through the Atomic and Molecular Cross section for Ionization and Aurora Database (ATMOCIAD; Gronoff et al. 2021) and also through comparison with the experimental data (Wedlund et al. 2011).
For this study, we modeled electron transport taking into account only magnetospheric electron precipitation.Secondary electrons resulting from ionization by solar UV radiation are neglected as their penetration capacity in Jupiter's atmosphere is low and considerably above the homopause of the hydrocarbons considered, such as CH 4 , C 2 H 2 , or C 2 H 6 .
The atmospheric model of the auroral region we used to model electronic transport is described in Grodent et al. (2001) and presented in Fig. 1.This model is 1D, and takes into account the majority of neutral species (H, H 2 , He, and CH 4 ) that predominate in Jupiter's atmosphere.It extends from the tropopause at a pressure of ∼100 mbar (altitude ∼100 km above the cloud level) to the upper thermosphere (altitude ∼2300 km above the cloud level), corresponding to a pressure of ∼10 −9 mbar.The initial density of thermalized electrons considered in the model is that obtained by Hinson et al. (1998) from radio occultations during the Voyager 2 flyby.Because of the limited data available, the initial electron temperature is thought to be similar to the temperature of the neutral atmosphere.In addition, as the atmospheric model used is 1D, we do not take into account any spatial or temporal variability in the abundance of neutral species in auroral regions, particularly the variability of CH 4 .Hence, since methane is the only tracer used in this study to model CR, any variability in its abundance can influence the CR.Thus, in this study we consider a homogeneous and steady chemical composition in Jupiter's whole aurora, which probably represents a significant approximation.
The electron-matter interactions considered in the physics of electron transport are the elastic interactions in which total kinetic energy is conserved, and the inelastic interactions of electrons with atmospheric particles illustrated in Table 1.The elastic cross sections e + H (Kingston & Walters 1980), e + He (Porter et al. 1987), e+H 2 (Muse et al. 2008), and e+CH 4 Fig. 1.Atmospheric model described by Grodent et al. (2001), which considers only the neutral compounds (H, H 2 , He, CH 4 , and C 2 H 2 ) that predominate in Jupiter's atmosphere.For electronic transport modeling, only H, H 2 , He, and CH 4 compounds were considered.For the UV emission model, spectral absorption by CH 4 and C 2 H 2 was taken into account.(Davies et al. 1989) considered are measured in different energy ranges, which may not entirely cover the energy grid 2 ranges needed for the electron transport modeling.For cross sections that do not cover the entire energy grid considered, we used power-law extrapolations to fill the gap.Above 400 eV the cross section σ(E) is considered to be proportional to E −0.65 (Wedde & Strand 1974), and above 2.2 keV the cross section is proportional to E −1 (Rees 1989).For inelastic cross sections we used the same approach as for elastic collisions.The power-law evolution of these cross sections at high energies makes it possible to use a decreasing logarithmic extrapolation to cover the entire range of the modeling energy grid.
As described in Appendix A, by solving the Boltzmann equation we calculate the electron flux F(z, E) as a function of altitude and electron energy.This flux results from the interaction between magnetospheric electrons precipitating in auroral regions and the neutral atmospheric particles considered in our model.To model auroral emission by H 2 , we used F(z, E) as an initial condition in our UV emission model.

H 2 UV emission model
The R and P branches of the Lyman band of H 2 (B 1 Σ + u −→ X 1 Σ + g ) correspond to the group of rovibrational transitions 2 For kappa distributions, we simulate electron transport in the energy range [1 eV, 1 MeV].For mono-energetic distributions, the electron transport modeling energy ranges are defined by [1 eV, E 0 ], where E 0 is the characteristic or average energy of the mono-energetic distribution considered.
Reactions Products that produce spectral lines with wavelengths in the range [80-190 nm].For the Werner band, in addition to the R and P branches, there is a third branch, the Q branch, corresponding to the rovibrational transitions that produce spectral lines in the spectral range [80-160 nm].There are also other transitions in the UV spectrum of H 2 , at shorter wavelengths, arising from the excited levels According to the atmospheric model used (molecular abundances and thermal profile), we begin by calculating the number density n(z, X, v i , J i ) at altitude z of the H 2 ground state levels.Thus, assuming that neutral species are thermalized in the atmosphere, the population of the ground state of H 2 follows the Boltzmann distribution given by where n H 2 (z) (cm −3 ) is the density of H 2 at altitude z; X, v i , and J i are respectively the ground electronic level n i , and the vibrational and rotational quantum numbers; g I (i) and E i are respectively the degree of degeneracy of the i state and its energy; k B is the Boltzmann constant; and T (z) the temperature at altitude z.In the denominator, the sum is made over all the rovibrational ground state levels of H 2 .
An H 2 molecule can be excited into a n j , v j , and J j state by various processes.It can be excited directly by absorbing a photon or by collision with an electron or other atmospheric particles.It can also be excited in this state by cascade de-excitation from higher states.Unlike models of H 2 UV auroral emission that use the Born 3 approximation to calculate the excitation rates of the different excited states of H 2 (e.g., Waite et al. 1983), in our model we calculate the excitation rates of the considered electronic levels through electronic transport by modeling e − + H 2 → H * 2 + e − collisions.Transitions from the EF, GK, and H H states to the X 1 Σ + g ground state are forbidden due to the selection rule on g −→ g transitions, and thus a non-negligible part of the B and C states are populated by these transitions (Liu et al. 2002).In this model we also take into account the excitation of H 2 to the EF, GK, and H H states, as well as the cascade populating of the B and C states.Excitation by other collisional processes with neutral particles is neglected because the atmospheric temperature is not high enough to produce UV emission from collisions of H 2 molecules with neutral particles (e.g., H 2 + H 2 = H * 2 + H 2 ).Thus, the volume excitation rate (cm −3 s −1 ) of a rovibrational state j is a linear combination of the direct excitation rate g direct and the cascade excitation rate g cascad and is given by (2)

Direct excitation rate
The direct excitation rate of H 2 by electron collisions is described by the following formula: Here the indices j and i are used to identify the upper quantum state and the ground state, respectively; is the electron flux at altitude z with energy between E and E + dE and is modeled by electron transport; σ i j (E) (cm 2 ) is the excitation cross section of the j level from the i level by collision with an electron of energy E (see Liu et al. 1998), and is given by where a 0 is the Bohr radius, 2J j +1 (dimensionless) is the oscillator strength of the transition between the i and j levels, and where A ji (s −1 ) is the Einstein factor of the j → i transition; R y = m e e 2 8h 2 ϵ 2 0 is the Rydberg constant; E i j is the energy of the transition from the i state to the j state; E is the energy of the incident electron; and x = E E i j .The coefficients C k and α were obtained experimentally by Liu et al. (1998Liu et al. ( , 2003) ) by fitting the excitation functions of the transitions X 2).This parameter was measured for the ungerade (odd) and gerade (even) levels of the excited states B and C (Liu et al. 1998(Liu et al. , 2003)).However, these factors were not measured for the B ′ , B ′′ , D, D ′ excited levels, and in this study, following Menager (2011), we consider that these coefficients are also valid for all ungerade states.

Cascade excitation rate
Electron collisions populate the g states including the EF, GK, and H H states from the X 1 Σ + g ground state.However, since g −→ g are forbidden dipolar transitions, the process of populating the g levels differs from that of the u levels.This populating is described by the same relationship (see formula (3)) as that for the u states, but with a different cross section given by A26, page 4 of 20 The function F(x) describes the excitation of the EF state from the ground state.This function is given by the formula where C m and α are described in Table 2 for transitions to g states.F(x) was measured only for the excited EF state.For the GK and H H states, Liu et al. (2002) suggested using the same excitation function, but multiplying it by a scaling factor to take into account the excitation efficiency of the different g states.This factor is 0.8 for GK and 0.35 for H H.
FC(v i , J i , n k , v k , J k ) are the Franck-Condon factors that describe the overlap of the wave functions of the i and k states and depend on the quantum numbers v i and J i for the ground state and on n k , v k , and J k for the upper g state.The Franck-Condon factors follow the selection rules ∆J = 0, ±2 and were calculated by Hervé Abgrall and Evelyne Roueff for Liu et al. (2002).
Finally, the function S r (J i , J k ) represents the rotational terms that were calculated by Abgrall et al. (1999) and adopted by Liu et al. (2003), and are given by the following formula: Here δ is the Kronecker parameter and β is an anisotropy parameter for which Liu et al. (2003) recommended a value of 0.6.Cascade excitation of the B and C states mainly increases emission from the low vibrational levels of B and accentuates emission from the C state by a smaller proportion.The cascade excitations of the B ′ , B ′′ , D, and D ′ states are not taken into account due to the lack of appropriate data.

Volume emission rate
The discrete volume emission rate η (cm −3 s −1 ) of a transition from a state j to a state i at altitude z is given by where A tot j (s −1 ) is the total Einstein factor of the upper level ( j) and is given by with A cont j the probability that level j transmits into the continuum.This results in the dissociation of H 2 into two fragments with kinetic energy E c , whose expression is given by where is the differential probability of dissociation of the j state into two fragments of kinetic energy E c and quantum number J i .
The Einstein factors for the B → X, C → X, B ′ → X, and D → X transitions were calculated by Abgrall et al. (1994) and are available in the MOLAT4 database.Those for transitions B ′′ → X and D ′ → X were obtained by personal communication from Abgrall and Roueff (Menager et al. 2010).The differential probability of dissociation A E c j→X,E c ,J i (E c ) are also available in the MOLAT database and were calculated by Abgrall et al. (1997) only for transitions from the B, C, B ′ , and D states.
The differential volume emission rate in the continuum η cont λ (cm −3 s −1 nm −1 ) comes from excited states above the dissociation threshold of H 2 .Its intensity at wavelength λ and altitude z after the dissociation of molecules from excited state j to dissociated state i is given by where A λ (n j , v j , J j → v i , J i ) is the differential probability of dissociation of state j, as a function of wavelength λ, and is obtained directly from A E c j→X,E c ,J i (E c ).Quantum transitions n j , v j , J j → n i , v i , J i can be represented by the wavelength λ which corresponds to the energy of each transition.Thus, by substitution, for discrete transitions we can write η(z, n j , v j , J j → X, v i , J i ) = η discr (z, λ), and for continuum transitions we can write

Auto-absorption
When an emitted photon from an excited state of H 2 is energetic enough to excite an H 2 molecule initially in the ground state, it can be reabsorbed by H 2 to emit another photon of lower energy.This self-absorption therefore tends to attenuate the UV emission spectrum toward short wavelengths (below 120 nm) and amplify it toward low-energy wavelengths.To take account of this phenomenon in the model we describe here, we used the results of Jonin et al. (2000), who experimentally studied the UV spectrum of H 2 in the wavelength range [90 nm; 120 nm].Thus, considering the volume emission rate η j→i (z 0 ) at altitude z 0 by the transition between quantum levels i and j, the volume emission rate reaching altitude z is given by where κ j→i (z) is a dimensionless attenuation factor depending on the extinction coefficient ϵ ji , the column density ζ i (z, z 0 ) of the quantum state i, and is given by 1.44  (13) with where λ i j is the wavelength of the j → i transition given in centimeters.
Thanks to this approach, self-absorbed photons are redistributed to lower-energy transitions using branching ratios calculated with the appropriate Einstein factors.In addition, as most of the continuum emission of the excited states of H 2 that we consider in this study occurs above 120 nm, continuum self-absorption is not taken into account in this model.

Synthetic spectrum
The flux I(λ) (cm −2 s −1 nm −1 sr −1 ) of UV emission from the atmosphere, in the θ direction, at infinite spectral resolution and without taking into account absorption by hydrocarbons (e.g.; CH 4 , C 2 H 2 , C 2 H 6 , ...) is given by where η tot λ (z, λ) (cm −3 s −1 nm −1 ) is the total differential volume emission rate at altitude z, which is the linear combination of the differential volume emission rate in the continuum (η cont λ (z, λ) and the differential volume emission rate of discrete transitions (η discr λ (z, λ)) in (cm −3 s −1 nm −1 )).
The thermal agitation of H 2 molecules having a mass m H 2 leads to the discrete emission lines with a spectral broadening defined by In our atmospheric model, at altitudes of around 400 km above the cloud level (average altitude of the auroral emission peak according to Bonfond et al. 2015) where the average temperature is around 600 K (Grodent et al. 2001), the mean spectral broadening of H 2 emission at 140 nm is around ∆λ ∼ 0.002 nm.This Doppler broadening is well below the broadening of the instrumental spectral resolution.In addition, since optical depth is very low above 400 km altitude, all spectral lines are optically thin.Thus, to calculate the differential volume emission rate of discrete transitions η discr λ (z, λ), we used this mean Doppler broadening for all the spectral lines considered in our model at all altitudes.Thus, considering the Doppler profile f ∆λ (λ) (nm −1 ) with a full width at half maximum ∆λ representing the mean Doppler broadening in the atmospheric model used, and using the volume emission rate of discrete transitions η discr (z, λ), we have Finally, the synthetic spectrum I synth (λ) is calculated, taking into account the instrumental resolution ∆λ ′ with where g ∆λ ′ (λ) is a Gaussian function with full width at half maximum ∆λ ′ and λ ′ is a dummy variable.

Hydrocarbon absorption
When magnetospheric electrons penetrate deeply enough below the homopause, the emission produced by the de-excitation of H 2 is attenuated at certain wavelengths by absorption from hydrocarbons.The attenuated spectral emission may be given by where .. is the total optical depth of hydrocarbons considered in the atmospheric model used; z 0 is the minimum altitude; and is the optical depth of species m, whose number density is given by n m (z), whose absorption cross section is given by σ m (λ), and where z ′ is a dummy variable.In this study we used the cross sections of CH 4 (Au et al. 1993;Kameta et al. 2002;Lee et al. 2001) and C 2 H 2 (Cooper et al. 1995;Nakayama & Watanabe 2004;Wu et al. 2001) measured experimentally in the UV range.

Comparison with laboratory spectra and validation
To validate the H 2 UV emission model, we compared our simulations with experimental results from the study of Liu et al. (1995).In this study, the authors used a 100 eV mono-energetic electron beam for bombarding H 2 molecules and measured the UV emission.The spectral resolution was of about 0.0125 nm in the spectral range (114-170 nm).In addition, H 2 molecules with a density of n H 2 = 4.55 × 10 12 cm −3 were placed in a cell at a pressure of around 0.4 µbar and a temperature of 300 K.
To validate our calculations, we compared relative spectra and neglected absolute intensities between the modeled and observed spectra in order to avoid the quantity of H 2 used in the Liu et al. (1995) experiment.To do this, we considered a thin atmosphere with a temperature of 300 K, in which we precipitated 100 eV mono-energetic electrons (best fit A) in order to match the experimental conditions of Liu et al. (1995) as faithfully as possible.We then followed the same approach by precipitating a Maxwellian electron (best fit B) beam with an average energy of 100 eV.In Fig. 2, we represent this comparison graphically in four spectral ranges (panels i to iv), each spanning 1 nm in order to evaluate the differences at very high spectral resolution.
For the rest of the spectrum, our model is in very good agreement with the experimental results.We note that the amplitudes of some spectral lines are not perfectly reproduced by the model.These differences are minimal and are generally caused by the inhomogeneity of the electron energy spectrum exciting H 2 molecules.In reality, the incident electron beam is not perfectly mono-energetic, and not perfectly Maxwellian either.As shown in the study by Dols et al. (2000), the variation in the energy of the electrons exciting the H 2 molecules may have a strong effect on some particular spectral lines.This influences their widths and amplitudes.
We also note a spectral shift of around 0.005 nm at some wavelengths between model and observation (e.g., between 123.8 and 124 nm in panel i, between 121.7 and 122 nm in panel ii, and between 120.4 and 121 nm in panel iii).According to Dols et al. (2000), this shift is caused by thermal expansion of the structure of the spectrometer's sensor during the measurements by Liu et al. (1995).

Juno/UVS observations
The aim of the Juno mission, launched in August 2011, is to study the planet Jupiter and its environment (Bolton et al. 2017).Its insertion into a highly elliptical polar orbit was achieved on July 5, 2016, and its first PJ was carried out on August 27, 2016.Since then the spacecraft has made several dozen PJs, with a periodicity of around 53.5 days during the nominal mission (until PJ37), leading to close flybys of the polar regions allowing us to study the interaction of the magnetosphere with the Jovian atmosphere.Juno hosts several scientific instruments, including the UltraViolet Spectrograph (UVS; Gladstone et al. 2017).The UVS is specifically designed to study Jupiter's atmosphere and auroral emissions in the extreme UV (EUV) and FUV domains.The wavelengths ranging from 68 nm to 210 nm are dispersed over a 256 spatial channel × 2048 spectral channel sensor (Davis et al. 2011;Greathouse et al. 2013;Gladstone et al. 2017).The spectrometer's slit has a dog-bone shape, and is oriented parallel to the axis of the spacecraft's rotation.This slit has a field of view at the edges of 2.55 • × 0.2 • and a spectral resolution of around 1.9-3.0nm, and a field of view at the center of 2 • × 0.025 • with a spectral resolution of ∼1.3 nm (Greathouse et al. 2013).Juno is a spin-stabilized probe with a period of about 30 s.As a result, the UVS slit is scanned across the sky to measure the UV emission spectrum in its field of view including the emission spectrum from Jupiter's poles (Bonfond et al. 2017).Each detected photon is associated with ancillary information including (latitude, longitude), x and y coordinates on the UVS detector, wavelength, emission angle from the planet.Counts recorded by UVS are converted into physical flux units using the instrument effective area derived from thousands of stellar observations during regular calibration phases (Hue et al. 2019(Hue et al. , 2021)).This photon list is rearranged in latitude-longitude-wavelength data cubes for each hemisphere.Latitude and longitude is sampled every 1 • and we used a 0.1 nm spectral sampling that fulfills the Nyquist criterion based on UVS spectral PSF.In addition, to increase the S/N of the UV emission spectra, we only map photons measured by the two UVS wide slits, for which the spectral resolution is around 2.1 nm (Greathouse et al. 2013), and discarded the photons coming from the narrow slit.
For the present study, we used spectral data obtained during PJ32, from 2021-Feb.-2116:23:09 UTC to 2021-Feb.-2122:38:45 UTC, for a total acquisition time of 6.25 h.This dataset includes 1.17 h acquisition time for the northern polar region and 4.04 h for the southern polar region, with about 1 h between acquisitions at the two regions (Juno's passage over the equator).
A26, page 7 of 20  The acquisition time at the south pole is naturally longer, due to the inclination of the semimajor axis of Juno's orbit to Jupiter's equatorial plane.During each 30 s spin, the UVS field of view intercepts Jupiter.Thus, each point on Jupiter or on the sky is looked at with an exposure of ∼18 ms during one spacecraft spin.
Figure 3 shows the integrated non-absorbed (N.A.) UV emission over the northern and southern polar regions.In order to isolate the auroral photons from the solar emission backscattered by the Jovian atmosphere, we established a selection criterion for pixels within the aurora (see Fig. C.1).We only selected pixels corresponding to UV emission spectra with a signal-to-noise ratio (S/N) ≥ 3. To evaluate the S/N, we consider the average of the unabsorbed part of the UV emission spectrum of H 2 in the spectral range [155 nm; 162 nm] to define the signal.For the noise, we estimated it within the same spectral range by subtracting the average signal value and calculating the standard deviation.
The UVS observations are co-added over the acquisition time into a large spectral datacube (latitude vs. longitude vs. wavelength), from which the CR CR = I(155 nm-162 nm) I(125 nm-130 nm) is then calculated, characterizing the absorption of the UV emission spectrum by CH 4 .It is important to note that we cannot use the initially defined wavelength range at the denominator (i.e., 123-130 nm) because of the uncertain calibration due to the detector degradation, which is due to gain sag on and around Lyman-α in the wide slit region of the detector.Hence, we start at 125 nm instead.Polar maps of the CR in the northern and southern hemispheres are shown in Fig. 4. As a consequence of this different wavelength range, the minimum CR in these maps is about 1.8, which is higher than the minimum CR ∼ 1.1 observed by Gustin et al. (2013) and corresponds to an unabsorbed UV emission spectrum.Regarding the maximum CR value, we have an overall CR max ≈ 30 for both poles.

Method
To map the characteristic energy of primary electrons precipitating in auroral regions, we modeled the CR using the TransPlanet electronic transport model combined with the H 2 UV emission model.For these simulations, we considered the Grodent et al. (2001) atmospheric model (an atmosphere of H, H 2 , He, and CH 4 ) with an altitude range from 100 km (∼1 mbar) to 2300 km (∼5.3×10 −12 bar) above the cloud level (see Fig. 1).To simulate the electron transport, we used two types of initial electron flux distribution.First, we used a mono-energetic distribution Φ(E, z max ) characterized only by a characteristic energy E 0 and given by Φ(E, z max ) ∼ δ(E − E 0 ), where z max represents the altitude at which the initial electron flux is injected into the atmosphere.In a second step, we used a kappa-type distribution Φ(E, z max ) ∼ f κ (E, ⟨E⟩) (Coumans et al. 2002) characterized by an average energy ⟨E⟩ and a κ parameter governing the logarithmic gradient of the distribution toward high energies.The kappa distribution used in this study is given by where Q 0 is the total energy flux and ⟨E⟩ is given as a function of the characteristic energy E 0 by the expression ⟨E⟩ = 2E 0 κ κ−2 , with E 0 representing the energy of the maximum amplitude of the distribution.
In Fig. 5, we give some examples of kappa electron flux distributions.Unlike the Maxwellian distributions used in previous electron transport models (e.g., Gustin et al. 2016), this distribution extends to higher energies.The value of kappa (κ = 2.5) that controls the amplitude of the distribution toward high energies and that we use in this study was derived from observed electron fluxes by Juno/JEDI during the first 20 PJs (Salveter et al. 2022).Hence, this kappa value gives rise to a realistic distribution of electron energy flux precipitating in Jupiter's auroral regions.
After modeling the electron transport, the resulting Φ(E, z) electron flux is used to calculate the excitation rates of H 2 Fig. 6.Optical depth calculated over the atmospheric column for CH 4 and C 2 H 2 .The transparent green and cyan bands represent the absorption spectral ranges used for the CR calculations.For CH 4 the absorption spectral range is considered between 125 nm and 130 nm and for C 2 H 2 between 150 and 153 nm.The transparent magenta band is the non-absorbed (N.A.) spectral range over which hydrocarbon absorption was assumed to be negligible.before simulating the atmospheric spectral emission in a given θ direction.Depending on the characteristic energy of the initial precipitated electron flux distributions, the electrons penetrate to varying depths into the atmosphere.As a result, the UV spectral emission of H 2 is absorbed to varying degrees by hydrocarbons below the homopause.
In the UV emission model, we considered absorption by CH 4 , which absorbs mainly below 140 nm (see Fig. 6 for the optical depth of CH 4 and C 2 H 2 ).We also included absorption by C 2 H 2 , whose distribution profile is shown in Fig. 1.This second hydrocarbon absorbs mainly in the spectral range [150-153 nm] and in other small ranges below 140 nm (see Fig. 6).This means that longward of 145 nm, only C 2 H 2 can attenuate the UV emission spectrum.However, short of this wavelength, CH 4 is the major absorber and C 2 H 2 has a weak influence on the amplitude of the UV spectrum.
During PJ32, the close flyby of Jupiter's poles is such that each point in both polar zones (north or south) is observed with a different viewing angle.Accordingly, each spectrum is measured with a different median θ emission angle (see Fig.

C.2).
In the next step, for each characteristic energy, we modeled the spectral emission with a spectral resolution of 2.1 nm (comparable to UVS spectra), varying the emission angle between 0 • and 80 • .This allowed us to establish the relationship CR(E 0 , θ), which links the CR, the emission angle, and the characteristic energy of the initial electron flux distribution injected at the top of the modeled atmosphere.With the assumption that the CR(E 0 , θ) function is monotonic, we could then invert this relationship to calculate E 0 using the CR observed by UVS and the emission angle maps.
In this study we distinguished the CR(E 0 , θ) relationship in the northern and southern hemispheres according to the magnetic dip angle in these auroral regions (see Appendix D.1).Thus, for each type of electron flux distribution and for each polar region, the relationship CR(E 0 , θ) is modeled and fitted5 using the formula D.2 (see Appendix D.2).For the case of mono-energetic initial electron flux distribution, the CR(E 0 , θ) relationship we modeled is shown in gold in the left panel of  of the relationship CR(E 0 , θ) in green grid lines6 , and we evaluated the uncertainty of this fit by plotting the absolute differences between the fit and the model for the north and south poles.The mean absolute difference is around 0.6, which is comparable to the uncertainty of the CR(E 0 , θ) relationship that we modeled for the case of the mono-energetic distribution (see Table D.1).
For the case of the kappa distribution, the results are shown in the right panel of

Results and discussion
In comparison with the CR(E 0 ) relationship modeled by Gustin et al. (2016) by precipitating a Maxwellian electron flux, our results show that the CR increases about two times faster in the case of a mono-energetic distribution.In their model, Gérard et al. (2014) also modeled the CR(E 0 ) relationship when precipitating mono-energetic electron fluxes.For a fixed θ emission angle, our results show that the CR increases about three times faster than that of Gérard et al. (2014).This can be explained by the fact that the modeling of auroral spectral emission in these previous studies is different from our case.In our UV emission model we calculate the volume emission rate (VER) by considering all the excitation rates of all the rovibrational levels (v, J) of the H 2 electronic states that we take into account.However, in the studies by Gérard et al. (2014) and Gustin et al. (2016), the UV emission of H 2 molecules is modeled using the volume emission rate (VER) obtained directly from the cross sections of the inter- measured by Dalgarno et al. (1999).The unabsorbed spectrum is obtained by multiplying the VER by a synthetic spectrum of H 2 (see Eq. ( 6) in Gustin et al. 2016).The absorbed spectrum is calculated taking into account the abundance of hydrocarbons in the atmosphere in the same way as explained previously.Thus, the differences in calculated CRs are mainly due to very small differences in the VER peak altitude obtained from these different studies.In addition, the electronic transport model used in these previous studies is based on Monte Carlo simulations.Compared with our transport model, which uses a radiative transfer solver, this could introduce additional discrepancies in the results.In comparison with the studies that modeled the CR(E 0 ) relationship, we still have a monotonically increasing CR as a function of characteristic energy.

Mapping the characteristic energy using the mono-energetic initial distribution of the electron flux
Using the modeled CR(E 0 , θ) relationship (shown in Figs.D.2 and D.3), we inverted Eq. (D.2) to determine E 0 (CR, θ).We obtained characteristic energy maps E 0 from the observed PJ32 CR and emission angle maps shown in Figs. 4 and C.2, respectively.
Figure 7 shows the resulting characteristic energy maps for the mono-energetic case in the northern (left panel of Fig. 7) and southern (right panel of Fig. 7) auroral regions.The average estimated uncertainty on the characteristic energy is derived from the average uncertainty on the CR(E 0 , θ) modeling presented in Table D.1.For the mono-energetic distribution, it is around 6 keV for both auroral regions.At the north pole the estimated maximum characteristic energy is around (150 ± 6) keV in the polar emission region.There are also other peaks with characteristic energies of up to (120 ± 6) keV in the injection zone of the main emission region near the pole.For the south pole the estimated maximum characteristic energy is around (170 ± 6) keV in the main emission region, with secondary peaks ranging from (90 ± 6) keV to (150 ± 6) keV and spread throughout the auroral arc.
In the polar emission zone the characteristic energy is low, with a peak of around ( 60 polar emission in the south was not very bright.This may be due to exceptionally low-energy electron precipitation or it could be due to exceptionally low electron flux in general, which does not produce strong UV emission.For other PJs (e.g., PJ4, PJ5, PJ8, PJ13, PJ14), the polar emission region in the south pole was relatively bright, and its spectral emission was largely detectable by UVS (Greathouse et al. 2021).
We also observed that there is no similarity between the characteristic energies of the auroral regions connected by the magnetic field lines between the north and south poles.This seems, at first order, to indicate a different electronic precipitation between the two auroral regions.However, the overall range of energies along the auroral ovals appears to be similar.
In our modeling of the CR(E 0 , θ) relationship we did not take into account the width of the horizontal extension of the auroral emission resulting from main emissions and Io's magnetic footprint (Bonfond et al. 2009).Thus, the energies inferred in the very narrow regions bordering the auroral ovals are not usable.In consequence, the results presented here are not valid for the Io footprint region.In addition, we did not take into account the SIII longitudinal motion of Io's footprint during the acquisition times in the northern and southern auroral regions.This could significantly decrease the emission and the CR of Io's footprint.Furthermore, according to Bonfond (2010), the average energies of electrons precipitating in this region are around 1-2 keV, which, according to our model, cannot be detected.For electrons with energies below 5 keV, the CR ratio produced is minimal because they cannot penetrate deeply enough below the CH 4 homopause.

Mapping the characteristic energy using kappa initial distribution of the electron flux
In the same way as for the case of the mono-energetic distribution, we also derived characteristic energy maps for kappa distributions of electrons precipitating in the auroral regions.Figure 8 displays the characteristic energy maps for the northern and southern polar regions (left and right panels of Fig. 8, respectively).The mean energy uncertainty of these maps is estimated around 4 keV for both auroral regions.At the north pole the maximum characteristic energy is around (40 ± 4) keV in the polar emission region.In the main emission regions we obtained a maximum characteristic energy of around (10 ± 4) keV.At the south pole we obtained peaks between 30 and 50 keV in the main emission oval.
The characteristic energy maps in Figs.7 and 8 cannot be directly compared.In the case of a mono-energetic distribution, the characteristic energy is identical to the average energy of the distribution, whereas for the kappa distribution E 0 represents the energy of the distribution peak.In the energy maps shown in Fig. 8, each pixel represents a kappa distribution described by a characteristic energy E 0 and a parameter κ = 2.5.Therefore, only the average energy of the kappa distribution can be used for comparison with a mono-energetic distribution.For a Kappa distribution, the conversion from characteristic energy to mean energy can be obtained from the relation ⟨E⟩ = 2E 0 κ κ−2 .Figure 9 represents the corresponding mean energy maps of electrons precipitating in the auroral regions for the case of a kappa distribution in the north (left panel of Fig. 8) and in the south (right panel of Fig. 8).The mean energy ⟨E⟩ being proportional to 2 κ κ−2 , for a parameter κ = 2.5 its value is therefore ten times larger than the characteristic energy.The uncertainty on the energies is also ten times larger.Comparison with the maps in Fig. 7 shows that the spatially averaged mean energy of electrons precipitating in auroral regions is significantly underestimated when the CR(E 0 , θ) relationship is modeled by mono-energetic distributions.However, for a fixed emission angle (θ = 0 • ) we compared the CR(⟨E⟩) relationship, at low energies, between the two cases of electron flux distribution precipitating in auroral regions (see panel c in Fig. D.4).We found that the CR(⟨E⟩) relationship obtained for the case of a mono-energetic distribution overestimates the mean electron energy below ⟨E⟩ = 90 keV.Above 90 keV, the energy of precipitating electrons is underestimated by the mono-energetic distribution assumption.This value of ⟨E⟩ = 90 keV, which represents the intersection between the two CR(⟨E⟩) relationships (obtained for the case of a mono-energetic electron flux distribution and for the case of a kappa distribution), seems to be A26, page 11 of 20 Benmahi, B., et al.: A&A, 685, A26 (2024) Fig. 9. Maps of the average energy ⟨E⟩ defined in Eq. ( 19) and given by the formula ⟨E⟩ = 2E 0 κ κ−2 .These maps are calculated for the north (left panel) and south (right panel) poles directly from the maps in Fig. 8.The iso-energy lines are defined as for Figs.Previous studies that modeled the relationship between the mean energy of precipitating electrons and CR in Jupiter's auroral regions obtained mean energy maps that differ from our results.In the case of Gérard et al. (2014), the mean-energy maps were based on HST observations from January 2014 and considered mono-energetic distributions of electrons precipitating in the atmospheric model described by Grodent et al. (2001).Gérard et al. (2014) obtained mean energy peaks of up to 500 keV in the northern polar emission region, which is comparable to the energy peaks we obtained in the case of a kappa distribution (see Fig. 9).However, this agreement should be qualified by the fact that these observations were obtained almost 10 yr apart and may be very different.Moreover, auroral fluctuations in UV-auroral emission brightness can be significant even over relatively short periods of time.Similarly, the Gustin et al. (2016) energy measurements cannot be directly compared with our results since they used observations obtained several years after PJ32.
Furthermore, we compared the average energies obtained in our study with the in situ measurements conducted by Mauk et al. (2020) using the JEDI instrument during PJs 4, 6, 7, and 10 of the Juno mission.However, the JEDI measurements by Mauk et al. (2020) for each PJ cover only a minimal fraction of the polar regions associated with the intersection of Juno's footprint and the auroral emission oval in the northern and southern auroral regions.We evaluated that the average energies of electrons precipitating in these auroral regions, measured by JEDI at Juno's altitudes, range between 150 and 300 keV.In comparison with the results of our study, where we consider kappa distributions for precipitating electrons, we observe energies distributed between 100 and 200 keV in approximately 70% of the auroral region (either north or south), with peaks reaching up to 600 keV in the remainder of these zones.Thus, the mean energies derived at auroral altitudes are generally of the same order as the results obtained by Mauk et al. (2020).However, at this stage, it is premature to draw a definitive conclusion regarding the comparison between our results and the measurements of electron energy distributions conducted by Mauk et al. (2020).The auroral emission maps used to derive the mean energies of precipitating electrons are integrated over several hours, whereas JEDI measurements are almost instantaneous.Moreover, given the rapid dynamics of polar regions, including phenomena like short-lived bright flares and local temporal variations, it is challenging to go beyond a comparison of orders of magnitude.This comparison thus requires a more in-depth analysis, which will be the subject of further study.
In the present study, and in previous similar studies (e.g., Trafton et al. 1994Trafton et al. , 1998;;Gustin et al. 2002Gustin et al. , 2016;;Ajello et al. 2005;Gérard et al. 2014), the CR(⟨E⟩) relationship is modeled by considering a 1D atmosphere model assuming a constant homopause throughout the auroral regions, which is a very approximate hypothesis.Hydrocarbon abundances in auroral regions are expected to be influenced by the precipitation of magnetospheric charged particles.Recent observations of these abundances (e.g., Sinclair et al. 2018) demonstrate that the spatial distribution of the main hydrocarbons in the auroral region is inhomogeneous.The maximum molar fraction variability factor in the auroral region is around 1.2 for C 2 H 2 , 1.1 for C 2 H 6 , and up to 1.3 for C 2 H 4 at 0.01 mbar pressure only in the north pole (see Sinclair et al. 2018).For the southern auroral region the variability is also significant, and differs from that of the north.As CH 4 photolysis is one of the sources of production of these different hydrocarbons, this suggests that the altitude of its homopause is also variable in the auroral regions.This may impact the estimates of the average energy of electrons precipitating in these regions.To evaluate this impact, we modeled the CR(⟨E⟩) relationship using two different CH 4 abundance profiles based on the A and C eddy diffusion models of Moses et al. (2005) and Hue et al. (2018).These profiles are not representative of CH 4 auroral abundance, but have a higher homopause compared with the atmospheric model of Grodent et al. (2001, Benmahi, B., et al.: A&A, 685, A26 (2024) for the average energy map determination may also be significant.However, this issue is beyond the scope of this article, and is left to a future investigation.

Conclusions
The present study is a further step forward in the investigation of electron energies precipitating in Jupiter's auroral regions.Inspired by Dols et al. (2000), Gustin et al. (2002) and Menager (2011), we developed a new model of UV emission from H 2 in Jupiter's auroral regions and adapted it to Juno/UVS observations.This model is more complete than previous studies.It takes into account nine electronic states of the H 2 molecule, and we combined it with the TransPlanet electronic transport model (Stamnes & Rees 1983b;Simon et al. 2007;Gronoff 2009;Menager 2011;Benmahi 2022).The H 2 auroral emission model was validated by simulating UV emission spectra at very high spectral resolution and comparing it with the results of Liu et al. (1995).The modeled UV spectra between 125 and 170 nm are compared with some of the emission spectra observed by UVS.The fit results are in good agreement with the observations (see Figs. B.1 and B.2) except in the wavelength range between 140 nm and 150 nm where the modeled spectra are more intense because we only take into account absorption by CH 4 and C 2 H 2 .Between 150 and 153 nm, absorption by C 2 H 2 is not strong enough for a good fit of the UVS spectra shown in Figs.B.1 and B.2.This suggests that the C 2 H 2 abundance profile we used in this study is underestimated in Jupiter's auroral regions.We did not include absorption by C 2 H 6 because we do not have an auroral abundance profile for this chemical species.However, as demonstrated in Gustin et al. (2016), C 2 H 6 absorbs in the interval [140 nm-150 nm] and influences the amplitude of the spectrum when electrons reach the homopause, typically for average energies above 20 keV.
Thanks to Juno/UVS observations during PJ32, we mapped the CR at the north and south polar regions and used them to map the energy of precipitating primary electrons.The relationship CR(E 0 , θ) was modeled taking into account the emission angle at each observed point in both auroral regions.We used the atmosphere model of Grodent et al. (2001) and considered a constant homopause throughout the auroral regions.
In the northern and southern auroral regions, using the JRM33 magnetic field model (Connerney et al. 2022), we found that the magnetic dip angle varies between ψ ∼ 60 • and ψ ∼ 80 • .As the penetration depth of electrons precipitating in these regions is influenced by the magnetic dip angle, this has an impact on the CR.We modeled the CR(E 0 ) relationship for small variations in ψ and found a small influence.This allowed us to consider only the median value of ψ in each auroral region for our modeling, in order to distinguish between the modeling of the CR(E 0 ) relationship in the north and south.
Modeling of the CR(E 0 , θ) relationship was carried out by considering mono-energetic initial electron flux distribution and kappa distribution (κ = 2.5).This allowed us to compare results from previous works to those obtained with a more realistic broadband population.We found that when considering a mono-energetic distribution, the average energy of electrons precipitating in auroral regions is globally underestimated by a factor of 3-4.We also found that at low energies (below ∼100 keV), the CR(⟨E⟩) relationships intersect (see panel c in Fig. D.4).Thus, below ⟨E⟩ = 90 keV, the average electron energy is overestimated if the mono-energetic hypothesis is used to infer the mean energy of a broadband population.Above 90 keV, the energy of electrons precipitating into auroral regions is underestimated.
Our results clearly demonstrate the importance of considering broadband distributions representative of the actual particle observations (e.g., Salveter et al. 2022) and modeled here as a kappa distribution when modeling the CR relationship as a function of the average energy of electrons precipitating in auroral regions.The average energies inferred by this method under the above-mentioned hypotheses for the atmospheric composition profile lie in the 300-500 keV range in the polar emission region of the north.In the main emission zones we found average energies up to 550 keV, with peaks along the auroral oval.In the outer emission regions the average energies lie between 5 and 50 keV.In the south the polar emissions are much fainter (see also Greathouse et al. 2021) with a mean energy peak of about 100 keV.In the main emission regions we found several average energy peaks from 150 keV to 600 keV.
Finally, thanks to this study, our work can be readily applied to mapping the average energy of auroral electrons for all Juno mission PJs.This will allow us to establish a temporal map of the electron energies precipitating in these regions.We also aim to improve this study by taking into account the meridional and latitudinal variabilities of the CH 4 homopause.For the case of mono-energetic electron flux distribution, these results demonstrate the importance of taking into account the geometry of the magnetic field lines at high energies because, typically above 120 keV (see Fig. D.1), the relationship CR(E 0 ) becomes dependent on the magnetic dip angle.For the kappa distribution case, the evolution of the CR(E 0 ) relationship is less influenced by the magnetic dip angle.D.2.Analytic formula and fit of the CR(E 0 , θ) relationship Following the study of Gustin et al. (2016), we considered that CR(E 0 ) follows a hyperbolic law at low energy and increases as a logarithmic law at high energy.Using our result from Fig.
D.1, we derive a phenomenological relation CR(E 0 ) in order to fit the CR modeling as a function of characteristic energy.This step allows us to obtain an analytical form of CR as a function of energy.In addition, thanks to a MCMC fit, it makes it possible to estimate the uncertainty on the modeled CR.Our analytical formulation of CR(E 0 ) for a fixed θ emission angle is given by where A is the minimum amplitude of the modeled CR; E c is a threshold energy; and B, C, D, α, and β are fit parameters that constrain the shape of the curve throughout the energy range.
To adjust these fitting parameters, we used the Python emcee package developed by Foreman-Mackey et al. (2013), which implements the MCMC method using the Metropolis-Hastings algorithm.The fitting configuration is characterized by 250 Markov chains and 2500 iterations.These two parameters were determined after several runs of the burn-in size determination.We found that the Markov chains converge after 500-1000 iterations on average for all fit parameters for the mono-energetic and kappa distributions.The choice of 2500 iterations ensures convergence of the Markov chains in all cases.
In Fig. D.1, each CR-energy relationship that we modeled (red dots for the north and blue dots for the south) was fitted using the formula D.1.The magenta and cyan colored envelopes represent the 1σ confidence band for each fit.This allowed us to estimate the uncertainty of the CR modeling as a function of energy.We found a mean CR uncertainty of around 0.5 for the case of a mono-energetic distribution and 0.2 for the case of a kappa distribution.
For a variable emission angle θ and a fixed characteristic energy, we found that the CR(θ) relationship follows a sinusoidal law.By taking into account the variability of the emission angle and the characteristic energy simultaneously, the CR(E 0 , θ) relationship is two-dimensional and is given by where δ and γ are additional fit parameters.As in the one-dimensional case (i.e., CR(E 0 )), we modeled the CR(E 0 , θ) relationship in a 2D (E 0 , θ) grid.This grid is defined by 14 characteristic energy points ranging from 1 keV to 220 keV for the case of a mono-energetic distribution, and 10 characteristic energy points ranging from 1 keV to 85 keV for the case of a kappa distribution.For each characteristic energy E 0 , we also considered an emission angle grid θ of ten points ranging from 0 • to 80 • .This corresponds to modeling 140 UV emission spectra, for each auroral region, to obtain a map of the CR(E 0 , θ) relationship for the case of mono-energetic distribution.For the case of kappa distribution, this corresponds to modeling 100 UV emission spectra, for each auroral region, to obtain the maps of the CR(E 0 , θ) relationship.
The electron transport modeling accuracy at each of the (E 0 , θ) grid points (illustrated above) is evaluated by a total energy conservation rate of the electrons precipitated in the model.Thus, for each electron transport simulation, the average rate of conservation of the total energy precipitated is about 99% for the case of a mono-energetic distribution and 99.5% for the case of a kappa distribution, which represents an energy loss, respectively, of 1% and 0.5%.All parameters were fitted in the same way as for the 1D case, using the same burn-in size for the MCMC method.In Figures D.2 and D.3, we represent in yellow the results of the CR(E 0 , θ) modeled relationship.To invert this phenomenological relationship, we used the formula D.2 and fit its parameters using the MCMC method explained in the text.In Table D.1, we presented the fit values of the parameters of the formula D.2 for each of the cases considered in our study.These parameters can be used directly in Eq.D.2 to map the characteristic energy of electrons precipitating in Jupiter's auroral regions for any PJ.∆ CR is the fitted noise of the modeled CR and corresponds to the uncertainty of CR(E 0 , θ).The average absolute difference between the fit and the model result is comparable to the uncertainty on the CR obtained by the MCMC method.show us the differences in the CR(⟨E⟩) relationships obtained for the case of a kappa distribution and the case of a mono-energetic distribution of precipitating electrons.Thus, using the Grodent et al. (2001) atmospheric model, we found that below 90 keV the CR(⟨E⟩) relationship obtained for the case of a mono-energetic distribution overestimates the average energy of precipitating electrons.Above 90 keV, the same relationship underestimates the average energy of electrons precipitating in auroral regions.
By analyzing the ratio of the maps of the mean energies derived from two different electron distributions, a kappa distribution and a mono-energetic one, we obtain the results presented in  where the ratio of mean energies approaches ∼1, the mean energies derived from the kappa distribution are around ∼90 keV.In areas where the ratio is less than 1, the mean energies associated with the kappa distribution are strictly less than 90 keV.

Fig. 2 .
Fig. 2. Examples of comparisons of modeled synthetic spectra with experimental spectra (red line) obtained by Liu et al. (1995) with a spectral resolution of 0.0125 nm.The blue line and the dashed green line represent the best fits obtained by precipitating a 100 eV mono-energetic (best fit A) and Maxwellian (best fit B) electron flux distributions.The four spectral windows shown have a width of 1 nm.In the ii window around 121.567 nm, the Lyman-α line was filtered out.

Fig. 3 .
Fig. 3. Integrated non-absorbed (N.A.) UV emission from Jupiter's auroral regions observed by Juno/UVS during PJ32 in the SIII jovicentric reference frame.The acquisition time over the northern polar region is 4282 s.At the southern polar region the acquisition time is 14 561 s.The plus sign (+) in red represents the average solar longitude during the selected acquisition times for the northern and southern polar regions individually.Thus, the Sun's longitudinal path between its mean position during acquisition at the northern hemisphere and the southern hemisphere is approximately 130 • westward.

Fig. 4 .
Fig. 4. Color ratio of Jupiter's auroral regions (left: northern hemisphere; right: southern hemisphere) observed during PJ32.The panels show the CR calculated for each pixel of the UV emission map as defined by Gustin et al. (2016) with CR = I(155 nm-162 nm) I(125 nm-130 nm) .
Fig. D.2 for the north pole, and separately in the left panel of Fig. D.3 for the south pole.In these same figures, we plot the fit

Fig. 7 .
Fig. 7. Characteristic energy maps E 0 obtained from the CR(E 0 , θ) relationships (see Figs. D.2 and D.3) modeled for the case of an initial monoenergetic electron flux distribution, and from the CR observed during PJ32, at the north (left panel) and south (right panel) poles.Iso-energy lines are defined for values of 1, 5, 10 keV, in steps of 20 keV between 10 and 300 keV, and then in steps of 100 keV between 300 and 900 keV.
Fig. D.2 for the north pole and in the right panel of Fig. D.3 for the south pole, using the same conventions.The average absolute difference between fits and models is around 0.5, which is comparable to the uncertainty of the modeled CR(E 0 , θ) relationship in this case.Whatever the initial electron flux distribution, our results shows that neglecting the emission angle leads to an underestimation of the characteristic energy.
Fig. 8. Characteristic energy maps E 0 obtained from the CR(E 0 , θ) relations (see Figs. D.2 and D.3) modeled for the case of an initial kappa electron flux distribution, and from the CR observed during PJ32, at the north (left panel) and south (right panel) poles.Iso-energy lines are defined in the same way as in Fig. 7.
7 and 8 and with the same color bar normalization.linked to the atmospheric model used and particularly to the CH 4 homopause considered (see panel c in Fig. D.4).
see panel a in Fig. D.4).Using this approach, the results obtained (Fig. D.4 shown in Appendix D.3) suggest that, depending on the CH 4 abundance profile used in our atmospheric model, the CR(⟨E⟩) relationship increases more rapidly or less rapidly as a function of ⟨E⟩, as shown in Fig. D.4.The consequences of this variability A26, page 12 of 20

Fig
Fig. D.2. 3D representation of the CR(E 0 , θ) modeled relationship for the case of a mono-energetic distribution (left panel) and for the case of a kappa distribution (right panel).Both relationships are calculated for the north pole (i.e., for ψ = 65.7 • ).The green grid lines represent the fit of the modeled CR(E 0 , θ) relationship.The colored surface (color bar from blue to red) represents the absolute difference between the modeled surface and the fit surface, with an average CR uncertainty of around 0.6 for the case of the mono-energetic distribution and an uncertainty of around 0.5 for the case of the kappa distribution.
Fig. D.3.3D representation of the CR(E 0 , θ) modeled relationship for the case of mono-energetic distribution (left panel) and kappa distribution, south pole (i.e., for ψ = 65.7 • ) and with the same conventions as Fig. D.2.The modeled CR uncertainty is estimated respectively at around 0.6 and 0.5 for both cases.

D. 3 .
Evolution of the CR(E 0 , θ) relationship using different CH 4 abundance profiles To evaluate the impact of methane distribution on the CR, we modeled the CR(⟨E⟩) relationship for an emission angle θ = 0 • and using two different CH 4 abundance profiles from the A and C eddy diffusion models ofMoses et al. (2005) andHue et al. (2018).These profiles are not representative of CH 4 auroral abundance, but have different homopauses compared to that of theGrodent et al. (2001) atmospheric model (see panel a in Fig. D.4).The CR(⟨E⟩) relationship was modeled for the cases of a kappa distribution and a mono-energetic distribution for electrons precipitating in auroral regions.The results obtained are shown in panels b and c in Fig. D.4.This shows us that the CR depends on the CH 4 abundance profile used in our atmospheric model.Thus, depending on the CH 4 homopause, the CR(⟨E⟩) relationship increases as a function of ⟨E⟩ more rapidly or less rapidly, as shown in Fig. D.4.The impact of this variability on the determination of the average energy map is also significant.The results shown in panel c in Fig. D.4 Fig. D.4.Model of the CR(⟨E⟩, θ) relationship for θ = 0 • using different CH 4 abundance profiles.In panel a, the red, blue, and green solid lines represent the CH 4 abundance profile from, respectively, the A and C eddy diffusion models of Moses et al. (2005) and Hue et al. (2018), and the Grodent et al. (2001) atmospheric model.In panel b is represented the modeled CR(⟨E⟩) relationship.The solid lines represent the CR modeled using kappa electron flux distribution and the dotted lines by using a mono-energetic flux distribution.The color conventions are the same as described for panel a.In panel c a zoomed-in version of panel b is shown to distinguish between the various CR(⟨E⟩) results at low energy.

Fig
Fig. D.5.Ratio of average energy map for the kappa-distribution case (Fig.9) to that of the monoenergetic case (Fig.7).Left: Northern auroral region.Right: Southern auroral region.

Table 2 .
Electronic excitation function parameters C k and α. always refers to the ground state and k to the upper g state (EF, GK or H H), x = E E ik , E is the energy of the incident electron, and E ik is the excitation threshold of the i −→ k transition.

Table D .
1. Fit parameters of Eq.D.2 for each of the cases considered in this study.