# Photochemistry-emission coupled model for Europa and Ganymede

Gaël Cessateur, Mathieu Barthelemy, Isabel Peinke

► To cite this version:

Gaël Cessateur, Mathieu Barthelemy, Isabel Peinke. Photochemistry-emission coupled model for Europa and Ganymede. Journal of Space Weather and Space Climate, 2016, 6, pp.A17. 10.1051/swsc/2016009. insu-02938731

HAL Id: insu-02938731

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

Submitted on 15 Sep 2020

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
# Photochemistry-emission coupled model for Europa and Ganymede

Gaël Cessateur<fnref n="1" />,<fnref n="2" />,*, Mathieu Barthelemy<fnref n="3" />, and Isabel Peinke<fnref n="3" />

Received 26 May 2015 / Accepted 26 January 2016

## ABSTRACT

In the frame of the JUICE mission, preliminary studies of the Jupiter's icy moons, such as Ganymede and Europa, are mandatory. The present paper aims at characterizing the impact of the solar UV flux and its variability on their atmospheres. The solar UV radiation is responsible for the photoionization, photodissociation, and photoexcitation processes within planetary atmospheres. A 1-D photoabsorption model has been developed for different observational geometries, on the basis of a neutral atmospheric model. Considering various production and loss mechanisms but also the transport of oxygen atoms, we estimate the red and green line emissions from photo impact-induced excitation only. These dayglow emissions can represent few percent of the global airglow emission, mainly dominated by electron-induced excitation in auroral regions. For limb viewing conditions, red line emission is bright enough to be detected from actual spectrometers, from 338 R to 408 R according to the solar activity. This is also the case for the green line with 8 R at limb viewing. Considering a different neutral atmosphere model, with an O₂ column density 50% more important, leads to a 14% increase in the red line emissions for limb viewing close to the surface. This difference could be important enough to infer which neutral model is the most likely. However, uncertainties on the solar UV flux might also prevent to constrain the O₂ column density when using ground-based observations in the visible only. The impact of solar flares on the red line emissions for Europa has also been investigated within a planetary space weather context.

**Key words.** Jovian moons – Aeronomy – Airglow – UV flux – Solar activity

## 1. Introduction and motivations

Planetary space weather is currently an area of high scientific interest (Lilensten et al. 2014), especially in the frame of the exploration of planetary environments such as the Jovian icy moons. In the frame of the JUICE (JUpiter ICy moons Explorer) mission, studies of such objects are of particular importance especially to plan instrument performances (Grasset et al. 2013). The scientific incentive for the JUICE mission is to study not only Jupiter, but also Europa and Ganymede, which are the arch type of icy moons. Both planetary objects are of particular interest for their potential habitability. This concept refers to whether environmental conditions are available that could support life and does not relate to whatever life presently exists or has existed on those worlds. JUICE will address this question not only by constraining the volume of water in the Jovian system, but also by realizing an inventory of elements on the surfaces as well as in the atmospheres of the icy moons. The need for accurate space weather studies of the exospheric variability regarding these two moons is therefore reinforced. An important challenge is thus to find the relevant observable quantities to characterize the space weather variations. In this paper, special focus will be given to the various emission and quenching processes within Ganymede's and Europa's atmospheres, regarding different levels of solar activity. As outcomes, several key physical parameters, namely atmospheric emissions from the visible spectral range, will be estimated.

The tenuous atmospheres of both icy moons are examples of surface-bounded atmospheres. Observations, combined with

numerous models of the neutral atmospheres for Europa and Ganymede, give valuable information about the complexity of such surface-bounded atmospheres. While Ganymede and Europa share similar characteristics, each of them has its own specific features. Observations from the *Hubble Space Telescope* (HST) have confirmed the presence of an O₂ atmosphere for both icy moons (Hall et al. 1998). The ionosphere has been detected by in-situ observations by Galileo (Gurnett et al. 1996; Kliore et al. 1997). The surface of both moons is widely covered by water ice (Calvin et al. 1995). The main processes responsible for the formation of such exospheres are direct sputtering by heavy ions (e.g. S⁺) located within Jupiter's strong magnetic field as measured by Galileo (Paranicas et al. 1999). This leads to the release of water molecules which can subsequently get radiolytically dissociated by sputtering resulting in the formation of new molecules (e.g. O₂, H₂), finally ejected from the surface into the exosphere. At Ganymede, sublimation of water molecules has an equivalent contribution to the exosphere generation for equatorial regions where the temperature of the surface is of about 150 K. The solar UV flux also contributes to the exosphere generation through the Photon Stimulated Desorption mechanism, but with an efficiency lower than that due to sputtering about 1.5 orders of magnitude (Plainaki et al. 2010). Based on these physical processes, numerical models show that the exospheres are mainly dominated by O₂ and H₂ for Europa (Shematovich et al. 2005; Smyth & Marconi 2006; Plainaki et al. 2013) with an O₂/H₂O ratio of about 15 at 0.04 *R*<sub>E</sub> and for polar region for Ganymede (Marconi 2007; Turc et al. 2014) with an O₂/H₂O ratio of about 130 at 0.04 *R*<sub>G</sub>. For equatorial regions at

<footnote n="1">¹ Royal Belgian Institute for Space Aeronomy, Ringlaan-3-Avenue Circulaire, 1180 Brussels, Belgium
*Corresponding author: gael.cessateur@aeronomie.be</footnote>

<footnote n="2">² Physikalisch-Meteorologishes Observatorium Davos, World Radiation Center (PMOD/WRC), Dorfstrasse 33, 7260 Davos Dorf, Switzerland</footnote>

<footnote n="3">³ UJF-Grenoble 1/CNRS-INSU, Institut de Planétologie et d'Astrophysique de Grenoble (IPAG), UMR 5274, 38041 Grenoble, France</footnote>
Ganymede, Marconi (2007) proposed an H₂O- and H₂-dominated atmosphere with a very low density of O₂ (with an O₂/H₂O ratio of about 10⁻⁴ at 0.04 R<sub>G</sub>). Recent numerical models by Plainaki et al. (2015) suggest however an O₂ atmosphere more homogeneous for both polar and equatorial regions within the first 100 km from the surface: the O₂ density is up to two orders of magnitude higher than the ones for polar regions derived from Marconi (2007) and Turc et al. (2014) at 0.05 R<sub>G</sub>. The water is still the dominant atmospheric species for equatorial regions at Ganymede. This discrepancy is explained by taking the multiple bouncing of the released O₂ molecules into account (see e.g. Cassidy et al. 2007). At higher altitudes, however, an inhomogeneous atmosphere is still expected following the inhomogeneity in the surface release, the sublimation being an efficient process at the equator. More detailed observations of Europa's and Ganymede's environment are thus needed in order to provide more accurate atmospheric models.

Ganymede is also unique for its intrinsic magnetic field with a dipole configuration (Kivelson et al. 1996). Interaction with the Jupiter magnetic field leads to open field lines in the polar regions, and closed ones for latitudes below 45°. Precipitating electrons along the open field lines have an impact on the atmospheric O₂, which in turn leads to auroral emissions in the UV as observed by the Hubble Space Telescope (HST; Hall et al. 1998; Feldman et al. 2000; McGrath et al. 2013). The observations of the auroral ovals as observed in the UV might also be used to characterize the subsurface ocean (Saur et al. 2015). Numerous modeling studies have showed that electron impact is the dominant dissociation mechanism of O₂ (Saur et al. 1998; Eviatar et al. 2001; Payan et al. 2015), leading to the generation of atomic oxygen in the O(⁵S) and O(³S) states which produce subsequently emission in the UV at 135.6 and 130.4 nm, respectively.

Observations of emissions in the visible are also good parameters for determining the different energetic inputs coming in upper planetary atmospheres, i.e. solar UV photons and precipitating particles, with the green (557.7 nm) and red-doublet (630, 636.4 nm) lines coming from the deactivation of oxygen atoms in the O(¹S) and O(¹D) states, respectively. These two visible lines have been indeed largely observed in the case of the Earth (see e.g. Witasse et al. 1999; Thirupathaiah & Singh 2014). More generally, these emissions are direct signatures of the presence of oxygen-bearing atmospheric species, such as O₂ or H₂O. Regarding Europa and Ganymede, observations of visible emissions are rather scarce. For Europa, observations in the visible have been reported by Porco et al. (2003): using the ISS narrow angle camera onboard Cassini with a broadband filter (200–1050 nm), emissions with a peak brightness of about 10 kR, with noise of the order of 5 kR, have been observed. Cassidy et al. (2008) have attributed these emissions majorly to the electron impact-induced excitation of sodium leading to emissions at 589 nm (Na(²S)-Na(²P) transition). The contribution of oxygen atoms to these visible emissions, including the red and green lines, is deduced to be limited with brightness values of about 100 R, and mostly at 777.4 nm and 844.6 nm from O(⁵P)-O(⁵S) and O(³P)-O(³S) transitions, respectively. Moreover, no emissions have been reported for Ganymede during the same Cassini observation campaigns. The relative importance of photochemical processes compared to electronic impact processes on the visible dayglow emissions has never been properly estimated for Europa and Ganymede. Electron impact is the dominant dissociation mechanism, but photodissociation is an important loss process and an important supply of oxygen to the neutral torus,

potentially including oxygen atoms in the O(¹S) and O(¹D) states, as discussed by Shematovich et al. (2005).

Observations of these visible emissions can help to constrain neutral atmospheric profiles. However, limb in-situ observations from an orbiter are necessary for measuring the green and red-doublet emissions. Earth-based observations provide only limited ability (e.g. with HST observations; Roth et al. 2014) to look at the vertical profile of the different volatile species of Europa and Ganymede or at the variation with latitude. Indeed, Earth-based observations suffer from too much scattered light. The first motivation of this paper is to provide some emission predictions to support Europa and Ganymede observations in the frame of the JUICE mission, by considering the photo impact-induced excitation only.

In a more general space weather frame, the second motivation is following Barthelemy & Cessateur (2014) studies on the solar UV flux accuracy to demonstrate on a case study, i.e. Europa and Ganymede, how the knowledge of the solar UV flux is of primary importance regarding the visible red and green dayglow emission estimations especially in the case of solar eruptive events. Since those estimations are directly linked to both the neutral atmospheric model of the planetary object and the solar UV flux, a flawed knowledge of the last one leads to an important uncertainty in the deduced atmospheric compositions. Adding to this, estimations of photochemical processes are relevant in the framework of the JUICE mission to properly disentangle the origin of potential detection of red and green line emissions.

The paper is organized as follows. In Section 2, the photolysis model within Europa's atmosphere is discussed, using two different solar activity levels as well as different observational geometries. In Section 3, dayglow emission for the red and green lines will be estimated based on several production and loss mechanisms of the atomic O(¹D) and O(¹S) states and by taking the transport of neutral oxygen atoms into account. The term dayglow is rather ambiguous since it can refer to the atmospheric emissions due to both solar UV photons and particle precipitations. In this paper, we consider the term dayglow referring to atmospheric emissions from photo impact-induced only. In Section 4, estimations of the red and green line dayglow for Ganymede will also be presented. Dayglow emission for Ganymede will be compared to auroral emissions in order to estimate the contribution of photon-induced to electron-induced emissions in Section 5. The importance of the chosen solar UV flux for dayglow modeling is discussed in Section 6, and a summary is given in Section 7.

## 2. Photolysis in Europa's atmosphere

### 2.1. Neutral atmospheric models

Two neutral atmospheric models will be used throughout this paper, the model of Smyth & Marconi (2006) and the one of Shematovich et al. (2005) (hereafter named Models 1 and 2, respectively):

- Model 1 from Smyth & Marconi (2006) is a two-dimensional kinetic model providing the density profiles for the water group species (H₂O, OH, O₂, O, H₂, and H). Europa's atmosphere is globally dominated by O₂ at lower altitudes, and by H₂ at higher altitudes, as illustrated in Figure 1.
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Radial density profiles for the water group species of Europa's atmosphere of Smyth & Marconi (2006) (thick lines) and of Shematovich et al. (2005) (dotted lines).</figcaption></figure>

– Model 2 is taken from Shematovich et al. (2005), who proposed a series of 1-D collisional kinetic models which differ in the initial $O_2$ and $H_2$ source rates and loss processes considered. Among these different models, we choose Model F, in which the $O_2$ surface source rate depends only on sputtering processes, while the $H_2O$ surface source rate depends both on sputtering and sublimation processes (with a surface temperature of about 100 K).

The key difference between these two neutral atmospheric models lies in the $O_2$ and O vertical profiles. In terms of column density, Model 2 has greater values compared to Model 1 with an increase of about 50% and 290%, respectively. This $O_2$ enrichment takes place for altitudes below 50 km. Model 2 does not provide the density profiles for $H_2$ and H. Comparison between atmospheric models will thus be performed when focusing on the processes involving $O_2$ and O, namely dayglow emission of atomic oxygen (see Sect. 3.3). The column densities estimated on the basis of Models 1 and 2 are also shown in Table 1. Let us also note that these column densities are higher compared to the ones presented by Smyth & Marconi (2006) and Shematovich et al. (2005). For numerical reasons, both atmospheric models have been interpolated till the ground, which leads to an increase in the column density values. It is important to mention that the following results regarding both photoionization and photoexcitation processes within Europa's atmosphere are only related to the column density values provided in Table 1, and not to the values provided by Smyth & Marconi (2006) and Shematovich et al. (2005).

## 2.2. Atmospheric response to solar UV flux

Four different cases with the geometries shown in Figure 2 are considered, where a spherical symmetry regarding the atmosphere is assumed. Case 1 represents the spacecraft above the equator looking to the nadir with a zero solar zenith angle. Case 2 corresponds to a nadir view at a latitude of 82°. Cases 3 and 4 simulate the observation of an orbiter which would look 10 km and 100 km above the surface of the planet, respectively. The line of sight crosses the whole atmosphere from beginning to end, with the closest altitude regarding Europa's surface being 10 km and 100 km. A line of sight abscissa is defined as the distance from the tangent point along the line of sight. The top of the atmosphere, closer to the observing

<table><caption>**Table 1.** Column density (in $cm^{-2}$) for the considered species based on both neutral atmospheric models.</caption><thead><tr><th>Species</th><th>Model 1 (Smyth & Marconi 2006)</th><th>Model 2 (Shematovich et al. 2005)</th></tr></thead><tbody><tr><td>O<sub>2</sub></td><td>6.2 × 10<sup>14</sup></td><td>9.48 × 10<sup>14</sup></td></tr><tr><td>H<sub>2</sub></td><td>7.4 × 10<sup>13</sup></td><td>—</td></tr><tr><td>H<sub>2</sub>O</td><td>5.5 × 10<sup>12</sup></td><td>5.48 × 10<sup>12</sup></td></tr><tr><td>OH</td><td>2.03 × 10<sup>11</sup></td><td>8.96 × 10<sup>10</sup></td></tr><tr><td>O</td><td>1.14 × 10<sup>12</sup></td><td>4.48 × 10<sup>12</sup></td></tr><tr><td>H</td><td>6.99 × 10<sup>10</sup></td><td>—</td></tr></tbody></table>

orbiter, is 2022 km (1920 km in case 4) and –2022 km (resp. –1920 km) for the farther point (corresponding to an altitude of 1000 km relative to Europa's surface). The atmosphere has been projected along the line of sight. The altitude regarding Europa's surface is calculated for each point of the line of sight, and the corresponding atmosphere is attributed.

Photoionization processes mainly involve the extreme UV (EUV, 10–121 nm), while photoexcitation and photodissociation processes involve the far UV (FUV, 122–200 nm) and middle UV (MUV, 200–300 nm). The solar UV flux as measured by recent space missions will be first used in order to estimate dayglow emission. Two different solar activities are considered, a moderate one (F10.7 = 83) and a strong one (F10.7 = 195). The data in the 27–115 nm range is measured by the EUV Grating Spectrograph (EGS), which is part of the Solar Extreme Ultraviolet Experiment (SEE) onboard TIMED (Woods et al. 2005). The 115–250 nm range is covered by the SOLar Stellar Irradiance Comparison Experiment (SOLSTICE) instrument onboard the SOLar Radiation and Climate Experiment (SORCE; Rottman 2005). Data in the 1–27 nm range come from the XUV Photometer System (XPS) instrument onboard SORCE, measuring the irradiance with several large passbands. The measurements are then processed with an algorithm using CHIANTI spectral models (Woods et al. 2008). A composite spectrum, with 1-nm bin resolution, is used, based on the observations of those three instruments. As those missions are Earth-located, a simple dilution law ($\lambda_{\text{dil}} = 1/2$) is used to estimate the flux at the Jovian system. In Section 6, different values for the solar UV fluxes will be used for inferring the sensibility of the respective upper atmospheric emissions.
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Modeling geometries. Cases 1 and 2 correspond to a nadir view for both equatorial and polar regions at noon. Cases 3 and 4 correspond to a limb view, with the closest altitude to Europa's surface being, respectively, 10 and 100 km. The red lines illustrate the line of sight within the atmosphere.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Photochemical production rates for all ions species and the electron production rate for nadir viewing at the equator (case 1) for high solar activity.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Photochemical production rates for all ions species and the photo-electron production rate for limb viewing (case 3) for high solar activity. The left y-axis corresponds to the line of sight abscissa, and the right y-axis indicates the corresponding altitude.</figcaption></figure>

For modeling the atmospheric response to the solar UV flux, a Beer-Lambert law is considered. The solar UV flux at each altitude z, $F(\lambda, z)$, is computed within the Europa's atmosphere, such as:

$$F(\lambda, z) = F_{\infty}(\lambda)e^{-\tau_{tot}}, \quad (1)$$

where $F_{\infty}$ represents the flux at the top of the atmosphere, i.e. the incoming solar UV flux, and $\tau_{tot}$ represents the optical depth, computed from:

$$\tau_{tot} = \sum_{k} \sigma_{tot}^{k}(\lambda) \int_{z}^{\infty} n_{k}(z') dz', \quad (2)$$

where $\sigma_{tot}^k$ are the wavelength-dependent total photoabsorption cross-sections for the six species, indexed k. $n_k$ is the neutral density at altitude z for each species k. These values are taken from the neutral atmospheric model.

The partial photoionization and photoexcitation cross-sections ($\sigma_a^k$) for several states (indexed a) are considered to calculate ion production rates (for $H_2O^+$, $OH^+$, $O_2^+$, $O^+$, $H_2^+$, and $H^+$), the electron production rate and production rates

for different atomic oxygen states (such as $O(^3P)$, $O(^1D)$, and $O(^1S)$). The production rate, $P_k^a$, is computed as follows:

$$P_k^a(z) = \int_{\lambda} n_k(z) \sigma_a^k(\lambda) F(\lambda, z) d\lambda. \quad (3)$$

The total and partial photoabsorption cross-sections used throughout this study have been summarized by Cessateur et al. (2012a).

## 2.3. Photoionization

The results of the photoionization in case 1 for a high solar activity (F10.7 = 195) are displayed in Figure 3. Very similar profiles have been obtained for case 2, with nadir observations at 82°, despite a different geometry. In both cases 1 and 2, the maximum of ionization occurs at the surface, with an electron production rate of about $8.6 \text{ cm}^{-3} \text{ s}^{-1}$. The predominant photochemically produced ion in this region is $O_2^+$, whereas at higher altitudes (above roughly 250 km), it becomes $H_2^+$.

Figure 4 displays the results of the photoionization in case 3, also for a high solar activity. The electron and ion productions are maximal where the atmosphere is the densest,
<table><caption>**Table 2.** Thermal electron density (column density in cm⁻²) integrated along the line of sight using the atmospheric model 1 for the four different cases and for two different solar activities.</caption><thead><tr><th></th><th>Solar activity<br>(F10.7)</th><th>Total thermal electron density,<br>or TEC (cm<sup>-2</sup>)</th></tr></thead><tbody><tr><td rowspan="2">Case 1</td><td>83</td><td>5.29 × 10<sup>10</sup></td></tr><tr><td>195</td><td>6.49 × 10<sup>10</sup></td></tr><tr><td rowspan="2">Case 2</td><td>83</td><td>5.24 × 10<sup>10</sup></td></tr><tr><td>195</td><td>6.44 × 10<sup>10</sup></td></tr><tr><td rowspan="2">Case 3</td><td>83</td><td>4.89 × 10<sup>11</sup></td></tr><tr><td>195</td><td>6.04 × 10<sup>11</sup></td></tr><tr><td rowspan="2">Case 4</td><td>83</td><td>1.63 × 10<sup>11</sup></td></tr><tr><td>195</td><td>1.66 × 10<sup>11</sup></td></tr></tbody></table>

i.e. at 10 km above the surface, with an electron production rate of about 6.1 cm⁻³ s⁻¹. The profiles look symmetric because the extinction is very small since Europa's atmosphere is very thin. However, the electron production is different at the entrance (at 2022 km) compared to the exit (−2022 km) by about 18%. Similar profiles have been obtained for case 4, but with significant lower photochemical production rates. At 100 km above the surface, the electron production rate reaches 0.17 cm⁻³ s⁻¹.

From the electron production rate, the thermal electron population within Europa's atmosphere can be estimated. In order to retrieve the electron density from the production rate, a steady state has been assumed, where production mechanisms are equivalent to loss mechanisms. According to Figures 3 and 4, there are two dominant ions with O₂⁺ and H₂⁺, regardless of the altitude. The following loss reactions H₂⁺ + e⁻ → 2H, and O₂⁺ + e⁻ → 2O are taken into account, with coefficient rates of 2.3 × 10⁻⁷ (300/Tₑ)⁰.⁴ cm³ s⁻¹ and 2 × 10⁻⁷ (300/Tₑ)⁰.⁷ cm³ s⁻¹, respectively (Florescu-Mitchell & Mitchell 2006). Tₑ represents the temperature of thermalized electrons. These reactions are known as the principal means of ion recombination within ionospheres (Banks & Kockarts 1973) and are more likely to occur than the simple recombination reaction (resulting in H₂ and O₂). The electron temperature is necessarily larger than or equal to the neutral temperature. As an approximation, the same temperature, as for H₂ provided by Smyth & Marconi (2006), of about 1000 K is considered. Table 2 summarizes the integrated electron densities (or column densities) along the line of sight, obtained with this rather simple chemical model. However, these values represent some lower limits, since neither electron precipitation, ion-neutral charge exchange, nor secondary photoionization within Europa's atmosphere has been considered. Some other loss processes should also be taken into account such as ion escape. However, Cessateur et al. (2012a) have applied the same model to estimate the electron density within Ganymede's atmosphere. The fact that a good agreement with in-situ measurements performed by Galileo (Gurnett et al. 1996) was found in that paper provides evidence for the validity of the applied methodological approach. The solar activity over the 11-year cycle leads to a total thermal electron density variability of 23% for cases 1–3, while it is only 2% for case 4.

## 2.4. Photodissociation and photoexcitation

Excitation production rates for several states of the atomic oxygen are now computed, namely the O(³P), O(¹D), and O(¹S), in

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

<figure><img src="image_7.png" /><figcaption>**Fig. 5.** Photochemical production rates of O(¹D) and O(¹S) from neutral species for equatorial viewing in case 1 (a) and limb viewing in case 3 (b), for high solar activity.</figcaption></figure>

the same way as for the photoionization. Those results are displayed for cases 1 and 3 in Figure 5, for high solar activity. Profiles for case 2 are very similar to case 1. In case 4, similar profiles to case 3 have been calculated but with significant lower production rates. While the O(³P) production rate maximum is 32 cm⁻³ s⁻¹ for case 3, it is only 0.8 cm⁻³ s⁻¹ for case 4. As for photoionization processes, atomic oxygen production rates maximize where the atmosphere is densest, i.e. close to the ground. For case 1, our model suggests a ratio of about 0.20 and 1.2 × 10⁻³ for O(¹D)/O(³P) and O(¹S)/O(³P), respectively. Similar values have been obtained for all other cases.

## 3. Dayglow emission modeling

From these previous results, estimation of the dayglow emission for several spectral lines can be made. Model 1 (i.e. Smyth & Marconi 2006) is first used for estimating dayglow emission. Comparisons between atmospheric models will be performed in Section 3.3. Two line emissions are of interest here: the radiative transition from O(¹D) to O(³P), which is responsible for the formation of the red line emission at 630 nm and the deactivation of the O(¹S) state leading to the green line emission at 557.7 nm. Both emission measurements are of primary importance since the green to red lines ratio is commonly used for determining the abundances of CO and CO₂ relative to the H₂O abundance in planetary and cometary atmospheres. Cessateur et al. (2012a) have shown that emission from the deactivation of OH(A²Σ⁺) to OH(X²Π) is very small in the
<table><caption>**Table 3.** Production and loss reactions of O(¹D) and O(¹S).</caption><thead><tr><th>Reaction</th><th>Rate<sup>a</sup></th><th>Reference</th></tr></thead><tbody><tr><td>(H<sub>2</sub>O, OH, O<sub>2</sub>, O) + hν → (O(¹D), O(¹S))</td><td>Model</td><td>This work</td></tr><tr><td>O<sub>2</sub><sup>+</sup> + e<sup>-</sup><sub>th</sub> → (O(¹D), O(¹S)) + O(³P)</td><td>2 × 10<sup>-7</sup> (300/T<sub>e</sub>)<sup>0.7</sup></td><td>Florescu-Mitchell & Mitchell (2006)</td></tr><tr><td>(1) O(¹D) + O<sub>2</sub> → O(³P) + O<sub>2</sub></td><td>k<sub>1</sub> = 3.2 × 10<sup>-11</sup></td><td>Atkinson et al. (2004)</td></tr><tr><td>(2) O(¹D) + O → O(³P) + O(³P)</td><td>k<sub>2</sub> = 2.5 × 10<sup>-11</sup></td><td>Kalogerakis et al. (2009)</td></tr><tr><td>(3) O(¹D) + e<sup>-</sup><sub>th</sub> → O(³P) + e<sup>-</sup><sub>th</sub></td><td>k<sub>3</sub> = 8.3 × 10<sup>-10</sup></td><td>Pavlov et al. (1999)</td></tr><tr><td>(4) O(¹D) + H<sub>2</sub>O → OH + OH</td><td>k<sub>4</sub> = 2.2 × 10<sup>-10</sup></td><td>Schmidt et al. (1988)</td></tr><tr><td>(5) O(¹D) → O(³P) + hν (630 nm)</td><td>A<sub>1</sub> = 6.478 × 10<sup>-3</sup></td><td>Fischer & Tachiev (2004)</td></tr><tr><td>(6) O(¹D) → O(³P) + hν (634.4 nm)</td><td>A<sub>2</sub> = 2.097 × 10<sup>-3</sup></td><td>Fischer & Tachiev (2004)</td></tr><tr><td>(7) O(¹S) → O(¹D) + hν (557.7 nm)</td><td>A<sub>3</sub> = 1.26</td><td>Wiese et al. (1996)</td></tr><tr><td>(8) O(¹S) → O(³P) + hν (297.7 nm)</td><td>A<sub>4</sub> = 0.134</td><td>Slanger et al. (2006)</td></tr><tr><td>(9) O(¹S) + O<sub>2</sub> → O(³P) + O<sub>2</sub></td><td>k<sub>5</sub> = 3.2 × 10<sup>-13</sup></td><td>Stuhl & Welge (1969)</td></tr><tr><td>(10) O(¹S) + H<sub>2</sub>O → OH + OH</td><td>k<sub>6</sub> = 3 × 10<sup>-10</sup></td><td>Zipf (1969)</td></tr></tbody></table>

case of Ganymede. Since the H<sub>2</sub>O and O<sub>2</sub> vertical density profiles are of the same order of magnitude within Europa's atmosphere, these emissions will not be computed here.

### 3.1. O(¹D) and O(¹S) emission lines

In order to compute aeronomic emissions, Cessateur et al. (2012a) assumed that radiative transition was the most effective way to deactivate excited species. Therefore different other loss mechanisms have been neglected. Several possible production and loss mechanisms for O(¹D) and O(¹S) have been identified by many studies (see e.g. Thirupathaiah & Singh 2014). In this current paper, significant chemical reactions are considered in order to provide accurate dayglow emission estimates. For the case of Europa's atmosphere, the different loss mechanisms as well as the production mechanisms are summarized in Table 3. The values of the rate constants and transition probabilities are also listed.

#### 3.1.1. Production mechanisms of O(¹D) and O(¹S)

The photoexcitation of the water group species ($P_{PE}$) is an important production pathway of the atomic O(¹D) and O(¹S) states. Figure 5 provides the radial profile of the total production rates for O(¹D) and O(¹S). The production rates here only include the impact of the solar UV flux on the water group species excluding other known processes. The dissociative recombination, ($P_{DR}$), of O<sub>2</sub><sup>+</sup> is also considered, i.e.,

$$O_2^+ + e_{th}^- \rightarrow O^+ + O(3P), \quad (4)$$

where O<sup>*</sup> = O(¹D), O(¹S), and O(⁵S). This reaction is generally an important source of atomic oxygen emissions (Banks & Kockarts 1973), and in particular for Europa (Nagy et al. 1998). The electrons involved in this reaction are mostly thermalized. As an approximation, a steady state is considered, where production and loss terms for the O<sub>2</sub><sup>+</sup> ion are supposed to be equal. Production rates of O<sub>2</sub><sup>+</sup> for several orbital configurations have been previously computed (see Figs. 3 and 4). The branching ratio for inferring the proportion of O(¹D), O(¹S), and O(³P) resulting from the dissociative recombination is provided by Kella et al. (1997). This allows us to derive the atmospheric emissions. Dissociative recombination of H<sub>2</sub>O<sup>+</sup> is neglected in the case of Europa, since the global water concentration considered here with Models 1 and 2 is particularly low. However, the existence of plumes in Europa as recently observed with HST could locally

increase the H<sub>2</sub>O and H<sub>2</sub>O<sup>+</sup> concentrations (Roth et al. 2014). The deactivation of O(¹S) by radiative decay at 557.7 nm leads also to the production of the O(¹D) state (see reaction (7)).

#### 3.1.2. Loss mechanisms of O(¹D) and O(¹S)

Several significant chemical reactions and radiative transitions are taken into account as loss mechanisms of O(¹D) and O(¹S), listed as reactions (1)–(10) in Table 3. For reactions (1), (2), and (4), the concentration profiles for O<sub>2</sub>, O, and H<sub>2</sub>O are taken from the Europa's atmospheric model (see Fig. 1). The thermal electronic density for reaction (3) is deduced from the thermal electronic production rate computed in the previous Section 2.3. Figure 6 displays the different loss rates for the O(¹D) state, for cases 3 and 4. Starting with case 3, the collisional quenching of the excited oxygen atoms by O<sub>2</sub> is a very effective loss process for altitudes lower than 80 km. The loss rate reaches 5.6 × 10<sup>-3</sup> s<sup>-1</sup> close to the ground. The collisional quenching by H<sub>2</sub>O shall also be considered but plays a minor role here with a loss rate of about 1.6 × 10<sup>-4</sup> s<sup>-1</sup> close to the ground. Regardless of the neutral and electronic temperatures, reactions (2) and (3) are not efficient for quenching the produced O(¹D) with maximal loss rate values of about 1.52 × 10<sup>-6</sup> s<sup>-1</sup> and 5.2 × 10<sup>-6</sup> s<sup>-1</sup>, respectively. Above 80 km, loss processes are mainly dominated by radiative processes leading to the red line emission at 630 nm. For case 4, with limb viewing 100 km above the surface, the O<sub>2</sub> concentration is not high enough in order to allow collisional quenching to be an effective loss process. Along the line of sight for case 4, deactivation of the O(¹D) goes through the radiative process with red line emission.

### 3.2. Red and green line dayglow emission estimations

In order to infer the red and green line emissions, the number densities of O(¹D) and O(¹S) within Europa's atmosphere have first to be computed. To do so, we solve the continuity equation for steady-state equations and spherical symmetry, such as

$$\frac{1}{r^2} \frac{d}{dr} (r^2 N_i v_i) = P_i - L_i, \quad (5)$$

where N<sub>i</sub> are the number densities for the excited atomic oxygen states, i.e. O(¹D) and O(¹S). P<sub>i</sub> and L<sub>i</sub> are the production and loss rates, respectively, as described in the previous sections. The radial velocity profiles, v<sub>i</sub>, for neutral oxygen

<footnote><sup>a</sup> Units are cm<sup>3</sup> s<sup>-1</sup> for two-body reactions, s<sup>-1</sup> for photolysis reactions and A<sub>i</sub> are in s<sup>-1</sup>.</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig. 6.** Total loss rate for O(¹D) for limb viewing at 10 km (a) and 100 km (b) above Europa's surface in black. Partial loss rates for radiative decay (in green thin line for reaction (5), in green dotted line for reaction (6)), for collisions with O₂ (in blue), and for collisions with H₂O (in red) are also indicated.</figcaption></figure>

atoms have been taken from Shematovich et al. (2005): below roughly 23 km, velocities are negative, i.e. the oxygen atoms are precipitating toward the surface, with velocities of around 100 m s⁻¹ at 10 km above the surface. For altitudes above 23 km, velocities are positive, i.e. the oxygen atoms are escaping. The transport will not impact the O(¹S) density because its lifetime is only 0.8 s. Oxygen atoms in the O(¹D) state with a lifetime of about 130 s, which represents a distance of 10–15 km at the considered velocity, will however be impacted. The computed densities are displayed in Figure 7 for case 2 (nadir viewing at 82°). For the sake of comparison, the number density of O(¹D) without transport is also represented.

The difference arising from the inclusion of the transport is significant through the whole line of sight. The O(¹D) number density is divided by 2 for altitudes lower than 15 km. Indeed, because of a high velocity, a large proportion of oxygen atoms in the O(¹D) state will impact the surface before radiating at 630 nm, which causes this depletion. As a boundary value, the O(¹D) number density is set to 0 at the surface. For altitudes above 100 km, where neutral oxygen atoms velocities are positive, the number density is shifted toward higher altitudes. Around 23 km, the transition zone occurs where the oxygen atoms have a mean velocity close to 0 which explains the local increase of the density with values close to photochemical equilibrium conditions (i.e. without transport). Generally, taking the transport into account leads to an

<figure><img src="image_9.png" /><figcaption>**Fig. 7.** O(¹D) number densities from photo impact-induced excitation, with (red line) and without transport (black line) for case 2.</figcaption></figure>

<table><caption>**Table 4.** Red and green line emissions in Rayleigh from photo impact-induced excitation, using the atmospheric model 1 for the four different cases and for two different solar activities.</caption><thead><tr><th></th><th>Solar<br>(F10.7)</th><th>Activity O(¹D⁻³P)<br>(redline)</th><th>O(¹S⁻¹D)<br>(greenline)</th></tr></thead><tbody><tr><td rowspan="2">Case 1</td><td>83</td><td>28.5</td><td>0.41</td></tr><tr><td>195</td><td>34.6</td><td>0.62</td></tr><tr><td rowspan="2">Case 2</td><td>83</td><td>28.3</td><td>0.41</td></tr><tr><td>195</td><td>34.2</td><td>0.61</td></tr><tr><td rowspan="2">Case 3</td><td>83</td><td>338.3</td><td>5.45</td></tr><tr><td>195</td><td>408.5</td><td>8</td></tr><tr><td rowspan="2">Case 4</td><td>83</td><td>46</td><td>0.3</td></tr><tr><td>195</td><td>55.6</td><td>0.32</td></tr></tbody></table>

O(¹D) depletion for altitudes below 50 km, and an O(¹D) enrichment for higher altitudes.

The red and green line emissions are directly deduced from the number densities by multiplying $N_i$ by the corresponding Einstein transition probability, i.e. $A_1$ and $A_3$ for the red line at 630 nm and the green line at 557 nm, respectively. The results for Europa's atmospheric emissions are summarized for the different orbital configurations and for different solar activities in Table 4. Similar values have been obtained for both nadir cases from 28 R to 34 R according to the solar activity. Let us precise that taking the transport into account leads to an emission decrease of about 23% for nadir cases because of the local depletion at lower altitudes and the loss by impact with the surface. For the limb viewing case 3, values are considerably higher compared to cases 1 and 2 only because of the longer line of sight within Europa's atmosphere, from 338.3 R to 408.5 R according to the solar activity. For case 3, taking the transport into account leads to an emission decrease of about 30.8%. For case 4, at 100 km above the surface, the emissions vary from 46 R to 55.6 R according to the solar activity. Red line emissions in case 4 are increased by a factor of about 125% compared to calculations without transport, which is explained by the O(¹D) enrichment at high altitudes.

The Jovian system is undoubtedly affected by the solar activity at any time scale (see Jackman & Arridge 2011).
<figure><img src="image_10.png" /><figcaption>**Fig. 8.** Total loss rate for O<sup>1</sup>D) for limb viewing 10 km above Europa's surface in black, using the atmospheric model 2. Partial loss rates for radiative decay (in green thin line for reaction (5), in green dotted line for reaction (6)), for collisions with O<sub>2</sub> (in blue) and for collisions with H<sub>2</sub>O (in red) are also depicted.</figcaption></figure>

Dayglow emissions are also impacted by the solar activity, but the effect is not linear. Both quiet and active solar activity conditions have been tested with F10.7 equal to 83 and 195, respectively. Dayglow emissions are increased by about 21% for every case considered here in this study. Actual spectrometers are nowadays sensitive enough to be able to detect such emissions for limb viewing (i.e. cases 3 and 4). Nadir cases are however not suitable for observations because of the reflected solar illumination on the surface. Therefore relative density measurements can be skewed by up to 21% if solar activity is neglected.

The green line emission at 577.7 nm is estimated in the same way. The radiative decay is very high ($A_3 = 1.26 s^{-1}$) compared to the loss rate resulting from the collisional reactions ($7.4 \times 10^{-5} s^{-1}$ and $2.72 \times 10^{-7} s^{-1}$ for reactions (9) and (10), respectively, at the ground). In the case of Europa's atmosphere, radiative decay is the most effective way to deactivate the O<sup>1</sup>S state of the atomic oxygen, with a quenching ratio of 0.9. The emissions by photo impact-induced excitation only, summarized in Table 4, are small, with a maximum of 8 R in case of the limb viewing (case 3) for a high solar activity. For case 4, the green line is pretty weak, close to 0.3 R, which might be difficult to detect. The green line emissions from photo impact-induced excitation are weak and will be discarded in the following.

### 3.3. Comparison with other Europa atmospheric models

Observations of dayglow emissions represent valuable information in order to constrain atmospheric species densities for planetary atmospheres. Dayglow emissions for the red and green lines have been estimated using the neutral atmospheric Model 1 (i.e. from Smyth & Marconi 2006). We consider from now on Model 2. The concentration profiles of O<sub>2</sub>, O, H<sub>2</sub>O, and OH have been displayed in Figure 1. As a reminder, there is a ratio of about 1.5 and 4 for column densities of O<sub>2</sub> and O, respectively, between the two neutral atmospheric models. Since H<sub>2</sub> and H are not provided, the previous density profiles from Model 1 are used. This will not affect our conclusions, since these two species, especially H<sub>2</sub>, are more important for the photoionization processes.

The loss rate using the atmospheric Model 2 is displayed in Figure 8. The loss rate is dominated by collisional reaction

<table><caption>**Table 5.** Red line emissions in Rayleigh using the model of Shematovich et al. (2005). Differences with the model of Smyth & Marconi (2006) in brackets.</caption><thead><tr><th></th><th>Solar (F10.7)</th><th>Activity O<sup>1</sup>D<sup>-3</sup>P (redline)</th></tr></thead><tbody><tr><td rowspan="2">Case 1</td><td>83</td><td>34.3 (+20.3%)</td></tr><tr><td>195</td><td>41.5 (+19.9%)</td></tr><tr><td rowspan="2">Case 3</td><td>83</td><td>385.9 (+13.9%)</td></tr><tr><td>195</td><td>465 (+14%)</td></tr><tr><td rowspan="2">Case 4</td><td>83</td><td>35.9 (-21.9%)</td></tr><tr><td>195</td><td>43.5 (-23.2%)</td></tr></tbody></table>

with O<sub>2</sub> for altitudes lower than 30 km, and the maximal loss rate is $1.7 \times 10^{-2} s^{-1}$ at 10 km above Europa's surface. Above 30 km, radiative decay is the dominant loss reaction. Table 5 summarizes the red line emission obtained with the model of Shematovich et al. (2005). Despite a higher loss rate, a 14% increase of the red line emission is estimated for case 3. For observations at the equator, estimations have been increased of 20%. Red line emissions are not linearly dependent on the O<sub>2</sub> concentration, because of the loss processes involving O<sub>2</sub>. However, an estimated 22–23% decrease of the red line emission is expected for case 4, for the limb viewing 100 km above Europa's surface. This result has to be related to the distribution along the altitude in O<sub>2</sub> for both atmospheric models (Fig. 1).

The difference in terms of aerodynamic emission (+14% and -22% for cases 3 and 4, respectively) might be important enough to infer an accurate estimation of the O<sub>2</sub> and O concentrations within Europa's atmosphere, providing a good knowledge of the solar UV flux, of the cross-section data, and of the relative brightness and morphologically distinguishable measurements for disentangling electron-impact excited emissions. Limb observations are clearly the best option for doing so, since dayglow emission estimations are high enough to be detected by actual onboard spectrometers. The signal will also not be polluted by the solar reflection from the surface. However, as will be explained in Section 6, these results are strongly dependent on the accuracy of the solar UV flux.

### 4. Ganymede's dayglow emission

Since the *Moons and Jupiter Imaging Spectrometer* (MAJIS) onboard JUICE will be devoted to the observations of both Europa and Ganymede, proper comparisons should be performed (Grasset et al. 2013). As previously mentioned, Cessateur et al. (2012a) already provided dayglow emission for the red and green lines in the case of Ganymede, but only for moderate solar activity (F10.7 = 82). However, neither the collisional loss nor the transport of neutral oxygen has been taken into account. Indeed, radiative decay was the only loss mechanism considered. In the light of the presented model for Europa, O<sup>1</sup>D) loss through collision and transport cannot be neglected, and we provide here more realistic emission estimates for red line emissions from photo impact-induced excitation.

For the sake of comparison, the same orbital configurations as for Europa are considered. Ganymede has, however, a different neutral atmosphere as modeled by numerous studies (Marconi 2007; Turc et al. 2014; Plainaki et al. 2015) with two distinguished regions: a polar one with an O<sub>2</sub>- and H<sub>2</sub>-dominated atmosphere, and an equatorial one where H<sub>2</sub>O
<figure><img src="image_11.png" /><figcaption>**Fig. 9.** Total loss rate for O(¹D) for the equatorial region of Ganymede. Partial loss rates for radiative decay (in blue) and for collision with H₂O (in red) are also illustrated.</figcaption></figure>

and H₂ are the most predominant species. The same reactions as listed in Table 3 for the production and loss mechanisms for the polar region are considered. The same reactions can also be applied to the equatorial region, but the dissociative recombination of H₂O⁺ (H₂O⁺ + e⁻ₜh → 2H + O) must be added. A branching ratio of 0.35 for producing the O(¹D) state from the dissociative recombination has been measured by Rosen et al. (2000). To take the transport into account, the mean radial velocity for oxygen atoms as computed by Marconi (2007) is considered. For polar regions at Ganymede, the radial velocity profile is very similar to the one at Europa. For equatorial regions, however, oxygen atoms always have a negative velocity, meaning that they are precipitating toward the surface. While the velocity is rather small (about 10 m s⁻¹) at low altitudes, it reaches 900 m s⁻¹ above 400 km up to 1000 km.

The loss rates for Ganymede's polar region are very similar to the ones obtained for Europa (see Fig. 6). However, they are different for the equatorial region since H₂O is the most predominant species there. Collisional quenching of the O(¹D) state by H₂O is the predominant loss reaction, especially for altitudes lower than 150 km, with a loss rate of about 0.2 s⁻¹ at the ground. For higher altitudes, the radiative decay of O(¹D) is the most dominant loss process, as displayed in Figure 9.

Red line emission can be estimated for several orbital configurations as summarized in Table 6. A direct comparison with previous estimations from Cessateur et al. (2012a) is not realistic since collisional quenching and transport were not taken into account. For equatorial observations, red line emissions are estimated from 5.3 R to 7.1 R when taking the transport into account for the two solar activities which represents an emission increase of 40%, between moderate and high solar activities. Photoexcitation of the H₂O is dominated by the solar Lyman α line, with roughly 70% of the produced O(¹D) state. The solar UV flux at this particular wavelength should be known with great accuracy in order to constrain H₂O density profiles within Ganymede's atmosphere at the equator when derived from visible dayglow data. Values for the nadir configuration at high latitudes are higher, from 20.6 R to 25 R depending on the solar activity which induced a variability of about 20%. For a spacecraft looking at the limb, at 10 km above Ganymede's surface, red line emissions from 346.8 R

<table><caption>**Table 6.** Red line emission for Ganymede for cases 1–3 (in Rayleigh).</caption><thead><tr><th></th><th>Solar activity (F10.7)</th><th>This work</th></tr></thead><tbody><tr><td rowspan="2">Case 1</td><td>83</td><td>5.3</td></tr><tr><td>195</td><td>7.1</td></tr><tr><td rowspan="2">Case 2</td><td>83</td><td>20.6</td></tr><tr><td>195</td><td>25</td></tr><tr><td rowspan="2">Case 3</td><td>83</td><td>346.8</td></tr><tr><td>195</td><td>418.7</td></tr></tbody></table>

to 418.7 R are estimated for solar UV photon impact-induced excitation only, according to the solar activity.

## 5. Limitation of our approach

Estimation of the red line emissions for Europa and Ganymede constitutes in fact a degenerate case study. Indeed, additionally to the solar UV flux impact, electron precipitation is expected to play an important role in the production of the O(¹D) oxygen state. All the previous results are obtained without considering this additional source and could be considered as accurate in regions with no electron precipitations. As discussed previously, the relative importance of photochemical processes compared to electronic impact processes in the visible dayglow emissions has never been estimated in the cases of Europa and Ganymede. Electron impact, however, is considered to be the dominant dissociation mechanism and must be taken into account for assessing the global red line emission estimate for those moons. In order to assess the contribution of the photo-induced emissions, we consider here the auroral regions of Ganymede, which are then the most sensitive to electron impact.

As for the impact of the solar UV flux, we consider an electron distribution (in energy) convolved with the appropriate electron impact cross-section, as provided by the ATMOCIAD database (Gronoff et al. 2012). A Beer-Lambert law has been used (see Eqs. (1)–(3) in Sect. 2.2). This gives directly an estimate of the reaction rate associated to the electron impact reactions on O₂, for a considered altitude in the atmosphere of the moon. The resulting red line emission is integrated along the line of sight. We consider here only case 2 for providing some estimates in the case of auroral regions. Based on Galileo data, we considered the electron distribution numerically described from 100 eV to few keV (see Eq. (2) from Paranicas et al. 1999) in the case of Ganymede. The cross-section for the O(¹D) oxygen state production requires to extend the energy range down to 2 eV. This is not necessary for the the O(³S) and O(⁵S) states, respectively, leading to the UV emission at 130.4 and 135.6 nm since their electronic impact cross-section cut-off energies are close to 100 eV. As reported by Paranicas et al. (1999), there are some concerns about the electron distribution for lower energies.

Since electron densities depend on local plasma properties, we use the UV emission at 135.6 nm as a proxy for evaluating visible emissions. For auroral regions in the case of Ganymede, UV emissions can reach from 50 to 300 R (see Fig. 3 from Feldman et al. 2000). These latter values are thus used to calibrate our electron distribution. In the case of 300 R for the UV emission, we estimate airlgow emission at 630 nm to be close to 1 kR for the orbital configuration 2, i.e. a nadir observation at high latitudes. In the previous section, we have estimated a
<table><caption>**Table 7.** Emission estimates at 135.6 nm and 630 nm, in Rayleigh, induced by electron impact for auroral regions at Ganymede (case 2) for an invariable neutral atmospheric model. The contribution (in %) of the photo-induced emission, estimated at 25 R (see Table 6) compared to the total emission at 630 nm, is also indicated.</caption><thead><tr><th>O(⁵S-³P) (135.6 nm)</th><th>O(¹D-³P) (630 nm)</th><th>Contribution (%)</th></tr></thead><tbody><tr><td>50</td><td>200</td><td>11</td></tr><tr><td>300</td><td>1000</td><td>2.4</td></tr></tbody></table>

dayglow of 25 R for high solar conditions (see Table 6), which represents 2.4% of the emission induced by electron impact only. At high latitude, it is possible to find regions with emission around 50 R. It could mean less precipitating electrons or a different atmospheric composition. As stated by Feldman et al. (2000), various combinations of O₂ column density and electron temperature are indeed possible. If keeping the same atmospheric model, the 630 nm emission induced by electron impact represents 200 R. In this case, visible emissions from photo-impact represent thus 11% in the case of Ganymede (see Table 7). These estimations are valid only for auroral regions in the case of Ganymede. In the case of Europa, we can expect a similar photon/electron-induced excitation ratio as a first approximation.

# 6. Discussions about the solar UV flux

The importance of an accurate knowledge of the solar UV flux and its variability for space weather purposes has been previously debated (Lilensten et al. 2008; Barthelemy & Cessateur 2014). On the basis of several examples, we will show how important an accurate knowledge of solar UV fluxes in terms of dayglow emission modeling is, especially during solar eruptive events. For the particular case of icy satellites, these emissions are not only powerful tools for constraining planetary atmospheres, but are also an important way to infer information regarding the interaction between the satellite and its environment. The impact of solar flares on Europa's atmosphere will first be discussed, for pointing out that Earth-based observations are not always a good solution in order to estimate the solar UV flux throughout the heliosphere. We will discuss some solar UV fluxes, defined as reference spectra, which lead to different dayglow emission estimations.

## 6.1. Solar flare events

In the frame of planetary space weather, the impact of solar flares on dayglow emissions is considered. During such events, the solar flux can increase up to 1000% for the XUV part, up to 100% for the EUV spectral range between 20 and 120 nm, and up to 10–30% for the FUV part between 120 and 200 nm. Both photoionization and photoexcitation processes will thus be directly affected. 4 X-class flares are taken into consideration here, listed in Table 8, for which the solar UV flux has been actually measured by TIMED. Figure 10 displays the solar UV flux for the X17 flare occurring on October, 28th 2003. This was one of the most powerful flares ever measured. The solar UV flux for high solar activity (F10.7 = 195) is also displayed for comparison. Earth-based observations are transposed at the Jovian system in order to take center-to-limb position effects into account of solar flares. It is quite well known that flare locations on the solar surface play an important role

<table><caption>**Table 8.** Red line emission estimations for case 3, for several solar flares events described by their dates and locations on the solar surface as well as their power densities.</caption><thead><tr><th>Events dates</th><th>Class</th><th>Location on solar surface</th><th>Red line emission (in Rayleigh)</th></tr></thead><tbody><tr><td>January, 15th 2005</td><td>X-1.2</td><td>Center</td><td>409.8</td></tr><tr><td>January, 20th 2005</td><td>X-7</td><td>Limb</td><td>411.3</td></tr><tr><td>October, 28th 2003</td><td>X-17</td><td>Center</td><td>644.6</td></tr><tr><td>November, 4th 2003</td><td>X-17.4</td><td>Limb</td><td>519</td></tr><tr><td>F10.7 = 195</td><td>—</td><td>—</td><td>408.5</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Solar UV fluxes at the Jovian system, for the X17 flare occurring on October, 28th 2003 (black curve), and for high solar activity conditions (F10.7 = 195, red curve).</figcaption></figure>

while looking at the impact on planetary atmospheres, especially in upper atmospheres (Qian et al. 2010). Considering that Jupiter and Earth are very far from opposition at some periods, such an extrapolation can generate some errors in the solar flux estimation.

Since aeronomic emissions are the most important and not polluted by solar reflection from the surface, only case 3 is considered here. Model 1 as neutral atmosphere is used. Table 8 presents our results regarding the different X-class flares for red line emission estimations. For the sake of comparison, estimation for a high level solar activity is also indicated. Not surprisingly, red line emissions increase when considering more powerful flare events. However, there is no simple relationship between the power density of flares and exospheric processes. Two X-17 flares lead to different red line emission and TEC values. Both flares are identical in terms of power density but occur at different locations, one near the disk center and the other near the limb as viewed from the Earth by TIMED. Most emissions in the XUV part of the spectrum ($\lambda \le 20$ nm) are optically thin and not absorbed by the solar atmosphere, while the EUV and FUV parts of the spectrum are optically thick, which leads to a smaller radiation enhancement. Photo-production of the O(¹D) state from O₂ is mainly induced by solar UV photons from 115 nm to 177 nm which are effectively absorbed by the solar atmosphere.

This location effect and the associated radiative transfer effect within the solar atmosphere are even more interesting when considering two flares with different power densities and locations. The X-1.2 flare located close to the disk center
leads to red line emissions of 409.8 R. An X-7 flare (nearly seven times more powerful in the x-ray domain) located near the limb induces a red line emission of about 411.3 R. Flare location is clearly an important factor in determining dayglow emissions for planetary space weather purposes. Using Earth-based measurements for flare events would not be a solution, especially when the Earth-Sun-Planet angle is too important. In the case of planetary space weather, this present study strongly outlines the need for distributed STEREO A & B like platforms or in-situ observations of the solar UV flux during flare events. For the X-17 solar flare, red line emissions for case 3 are predicted to be of 644.6 R, which represents an increased emission compared to quiet high solar conditions by about 58%. The impact of the flare location is again nicely illustrated with the X-17.4 solar flare with a red line estimation of about 519 R, 20% less than the X-17 solar flare located at the center of the Sun.

## 6.2. Modeling using reference spectra

Some reference spectra in the UV have been defined for aeronomical purposes. Two spectra, for moderate and high solar activities, are used for inferring rate constants regarding photoabsorption processes for several species (Huebner et al. 1992). Two more recent reference spectra have also been proposed for aeronomical purposes, based on the result of the Atlas mission (Thuillier et al. 2004) for two solar activities (F10.7 = 78 and F10.7 = 192). These two spectra will thus be used for assessing dayglow emission, which will be compared to the previous estimations obtained with recent measurements of the solar UV flux. This comparison makes sense since solar activity levels have been chosen to match the Atlas reference spectra ones.

Using both neutral atmospheric models of Smyth & Marconi (2006) (Model 1) and Shematovich et al. (2005) (Model 2), Table 9 summarizes the estimated red line emissions for observations at the limb. Similar red line emission estimations shall be expected when taking two solar fluxes representing the same level of solar activity into account. There is surprisingly a difference of about 6%. For a high solar activity level, red line emissions are quite variable, from 408.5 R to 433.1 R, in the case of the first neutral atmospheric model. Such a discrepancy is obviously coming from differences within the solar UV flux: the solar Lyman α line as well as the continuum between 130 and 180 nm are more intense within the Atlas spectrum compared to SORCE/SOLSTICE and TIMED/SEE observations (with a ratio of about 1.12 and 1.02, respectively). Solar indices, such as the F10.7, used here to compare the different solar UV fluxes, actually do not reflect the level of variability in the UV particularly well (Dudok de Wit et al. 2009). Some notable differences between two solar fluxes could thus be expected, even for very similar F10.7 values. Finally, EUV and FUV fluxes are measured for two distinct solar cycles (cycle 22 for Atlas, cycle 23 for TIMED and SORCE observations). Both solar cycles are quite different in terms of irradiance amplitude which could lead to differences within the solar UV flux (Ermolli et al. 2013).

The difference induced by the uncertainties in the solar UV flux might have some implications in planetary atmosphere modeling. Indeed, using two different solar UV fluxes for the same solar conditions as characterized by F10.7 leads to a difference of 6–7% on the red line emission estimations. The difference arising from using two neutral density profile models leads to an uncertainty of about 14% for limb observations

<table><caption>**Table 9.** Red line emission estimations (in Rayleigh) obtained with the Atlas spectrum, compared to those obtained with SORCE and TIMED observations for case 3.</caption><thead><tr><th></th><th>Model 1</th><th>Model 2</th></tr></thead><tbody><tr><td>Atlas 1</td><td>433.1</td><td>493.6</td></tr><tr><td>SORCE/TIMED (F10.7 = 195)</td><td>408.5</td><td>465</td></tr><tr><td>Atlas 3</td><td>361.5</td><td>412</td></tr><tr><td>SORCE/TIMED (F10.7 = 83)</td><td>338.3</td><td>385.9</td></tr></tbody></table>

(see Tables 5 and 9). When considering the solar UV flux impact only, its uncertainties then lead to a difference which is half as important compared to the difference induced by the chosen atmospheric model. However, as a reminder, there is a ratio of about 1.5 for the $O_2$ column density between Model 1 and Model 2 ($N_{O_2}^2/N_{O_2}^1 = 1.5$). Let us consider now a hypothetical model 3, i.e. an average profile $O_2$ density between Model 1 and Model 2, with a ratio of about 1.2 for the $O_2$ column density compared to model 2 ($N_{O_2}^2/N_{O_2}^3 = 1.2$). Red line emissions for limb observations at 10 km above the surface are estimated at 389.6 R using the Atlas 3 spectrum for Model 3. A difference of 1% then arises between the red line estimation from Model 2 using the solar UV flux from SORCE/TIMED (385.9 R) and Model 3 with the Atlas 3 spectrum (389.6 R). It would then be difficult to distinguish which $O_2$ column density is the most likely between Model 2 and Model 3 (difference of 20% for the $O_2$ column density) based on our current knowledge of the solar UV flux used as reference spectra.

In regard to our current knowledge of the atmospheric density of the Jovian moons, the uncertainties induced either by the solar UV flux or by the considered atmospheric models are minor (Eviatar et al. 2001; Plainaki et al. 2013, 2015). However, regarding the future JUICE measurements it might be virtually impossible to disentangle atmospheric model errors, even for the regions with no electron precipitation, up to roughly 20% for the $O_2$ column density because of the solar UV flux uncertainties. This will also necessitate an accurate estimate of the emissions induced by particle precipitation. A proper evaluation of the electron impact-induced excitation leading to the red line emission will be evaluated in a forthcoming paper, using a transport model within Europa and Ganymede atmospheres.

Previous results suggest that we need accurate solar data, and this for the whole EUV/FUV part of the spectrum. For doing so, special efforts have to be made to assess the absolute values of the solar UV flux along its intrinsic variability for space weather purposes. Using reference spectra such as Atlas is not a recommended solution as previously illustrated. One other alternative would be to rely on solar UV models. Major efforts have been done toward developing robust models, based on solar proxies, such as the radio flux at 10.7 cm (F10.7) (Tapping & Detracey 1990) or the MgII core-to-wing index (Heath & Schlesinger 1986), which should assess the solar UV flux variability. However, it has been largely argued that no combination of existing indices can properly reconstruct the solar EUV/FUV irradiance on all time scales (see e.g. Dudok de Wit et al. 2009; Cessateur et al. 2011). This is quite prejudicial for properly characterizing dayglow emissions as previously seen, and more generally for retrieving upper atmosphere parameters (see e.g. Lilensten et al. 2007; Cessateur et al. 2012a). There is then no better alternative than direct measurements of the local solar UV flux in near-real time.
We then advocate that any planetary missions should embark on some instruments devoted to the characterization of the solar UV flux. This would certainly remove one part of the uncertainties within modeling processes. Cessateur et al. (2012b) have proposed an instrumental concept, particularly suited for planetary missions. With a low weight and using radiation-hard materials, a radiometer composed of five spectral bands can provide the solar UV flux with great accuracy, but for a non-flaring Sun. Further developments are clearly needed for assessing the solar UV flux during flare events with simple instruments such as radiometers. This instrument is suited not only in the case of a Jovian mission for assessing Europa's and Ganymede's atmospheres, but clearly for all planetary missions, where dayglow emissions are measured.

## 7. Conclusions

We have explored in this paper the effects of the solar UV flux and its variability on the dayglow of Europa and Ganymede. The thermal electron population, as well as the photochemical production rates along the atmosphere of different oxygen states, namely O(³P), O(¹D), and O(¹S), are the outcomes of the modeling. The continuity equation within the atmosphere of both moons has been solved in order to calculate the number densities of the excited states. Considering several orbital configurations, dayglow emission induced by solar UV photons only has been subsequently estimated for the red and green lines at 630 nm and 577.7 nm, respectively. Taking the velocity as modeled by Shematovich et al. (2005) and Marconi (2007) into account, the transport is a critical process with a significant impact on the radial distribution of O(¹D). As orbital configurations, we considered limb observations 10 km above the surface, and red line emissions have been estimated around 400 R for high solar conditions, using the neutral atmospheric model of Smyth & Marconi (2006). Those emissions, induced by photo-impact only, are however a secondary process compared to electron-impact excitation and can contribute few percent of the total visible emission in auroral regions of Ganymede.

Using planetary atmospheric emissions for assessing both atmospheric compositions and energetic inputs is a highly degenerate case study, especially within regions where both electron precipitation and the solar UV flux are responsible for the dissociation and excitation of the neutral species. A way to reduce the degeneracy of the system is then to separate electron and the solar UV flux problems. The importance of an accurate knowledge of the solar UV flux and its variability for space weather purposes has also been discussed. The impact of the solar flares on dayglow emission from Europa has been explored for the first time. Location effects of solar flares have been enlightened, which supports the idea that in-situ observations of the solar UV flux are mandatory when performing planetary space weather studies. Using reference spectra might also lead to different dayglow emission estimations, compared to the ones obtained with recent solar UV measurements for the same level of solar activity. This difference of about 6% might affect planetary atmosphere modeling. For regions without electronic precipitation, this could also prevent the discrimination of two distinctive neutral atmospheric models if their O₂ column densities differ by roughly 20%, when only using red line emissions observations.

This present study makes sense when considering the future JUICE mission outcomes. Indeed, in-situ measurement will provide better neutral and electron atmospheric densities

and secondary processes as photo impact-induced ionization and excitation could be of importance to correctly interpret the data. In order to do so, a necessary step will be to disentangle the different contributions, electronic from photo processes. Moreover, the variability of the solar UV flux, from quiet conditions to transient events such as solar flares, can lead to important red line emission variability at limb viewing. The difference, induced by the solar UV flux, could impact the resulting planetary atmosphere models, which outlines the need for direct measurements of the solar UV flux.

*Acknowledgements.* GC was supported by the Fonds de la Recherche Scientifique Grant PDR T.1073.14 Comparative study of atmospheric erosion and received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under Grant Agreement No. 284461 (eHEROES project, http://soteria-space.eu/eheroes/html/) and No. 313188 (SOLID project, https://projects.pmodwrc.ch/solid/). GC also acknowledges support from Swiss COST office (Grant No. C11.0135). The authors and the editor thank two anonymous referees for their assistance in evaluating this paper.

## References

Atkinson, R., D.L. Baulch, R.A. Cox, J.N. Crowley, R.F. Hampson, R.G. Hynes, M.E. Jenkin, M.J. Rossi, and J. Troe. Evaluated kinetic and photochemical data for atmospheric chemistry: volume I – gas phase reactions of Oₓ, HOₓ, NOₓ and SOₓ species. *Atmos. Chem. Phys.*, **4** (6), 1461–1738, 2004, DOI: 10.5194/acp-4-1461-2004.

Banks, P., and G. Kockarts. *Aeronomy, Part A*, Academic Press, New York, 1973, DOI: 10.1016/B978-0-12-077802-7.50002-6.

Barthelemy, M., and G. Cessateur. Sensitivity of upper atmospheric emissions calculations to solar/stellar UV flux. *J. Space Weather Space Clim.*, **4** (27), A35, 2014, DOI: 10.1051/swsc/2014033.

Calvin, W.M., R.N. Clark, R.H. Brown, and J.R. Spencer. Spectra of the icy Galilean satellites from 0.2 to 5 µm: a compilation, new observations, and a recent summary. *J. Geophys. Res.*, **100**, 19041–19048, 1995, DOI: 10.1029/94JE03349.

Cassidy, T.A., R.E. Johnson, M.A. McGrath, M.C. Wong, and J.F. Cooper. The spatial morphology of Europa's near-surface O₂ atmosphere. *Icarus*, **191**, 755–764, 2007, DOI: 10.1016/j.icarus.2007.04.033.

Cassidy, T.A., R.E. Johnson, P.E. Geissler, and F. Leblanc. Simulation of Na *D* emission near Europa during eclipse. *J. Geophys. Res. [Planets]*, **113**, E02005, 2008, DOI: 10.1029/2007JE002955.

Cessateur, G., T. Dudok de Wit, M. Kretzschmar, J. Lilensten, J. Hochedez, and M. Snow. Monitoring the solar UV irradiance spectrum from the observation of a few passbands. *A&A*, **528**, A68, 2011, DOI: 10.1051/0004-6361/201015903.

Cessateur, G., J. Lilensten, M. Barthélémy, T. Dudok de Wit, C. Simon Wedlund, G. Gronoff, H. Ménager, and M. Kretzschmar. Photoabsorption in Ganymede's atmosphere. *Icarus*, **218**, 308–319, 2012a, DOI: 10.1016/j.icarus.2011.11.025.

Cessateur, G., J. Lilensten, T. Dudok de Wit, A. BenMoussa, and M. Kretzschmar. New observation strategies for the solar UV spectral irradiance. *J. Space Weather Space Clim.*, **2**, A16, 2012b, DOI: 10.1051/swsc/2012016.

Dudok de Wit, T., M. Kretzschmar, J. Lilensten, and T. Woods. Finding the best proxies for the solar UV irradiance. *Geophys. Res. Lett.*, **36**, L10107, 2009, DOI: 10.1029/2009GL037825.

Ermolli, I., K. Matthes, T. Dudok de Wit, N.A. Krivova, K. Tourpali, et al. Recent variability of the solar spectral irradiance and its impact on climate modelling. *Atmos. Chem. Phys.*, **13**, 3945–3977, 2013, DOI: 10.5194/acp-13-3945-2013.

Eviatar, A., D.F. Strobel, B.C. Wolven, P.D. Feldman, M.A. McGrath, and D.J. Williams. Excitation of the Ganymede ultraviolet aurora. *Astrophys. J.*, **555**, 1013–1019, 2001, DOI: 10.1086/321510.
Feldman, P.D., M.A. McGrath, D.F. Strobel, H.W. Moos, K.D. Retherford, and B.C. Wolven. HST/STIS ultraviolet imaging of polar aurora on Ganymede. *Astrophys. J.*, **535**, 1085–1090, 2000, DOI: 10.1086/308889.

Fischer, C.F., and G. Tachiev. Breit-Pauli energy levels, lifetimes, and transition probabilities for the beryllium-like to neon-like sequences. *At. Data Nucl. Data Tables*, **87** (1), 1–184, 2004, DOI: 10.1016/j.adt.2004.02.001.

Florescu-Mitchell, A.I., and J.B.A. Mitchell. Dissociative recombination. *Phys. Rep.*, **430**, 277-374, 2006, DOI: 10.1016/j.physrep.2006.04.002.

Grasset, O., M.K. Dougherty, A. Coustenis, E.J. Bunce, C. Erd, et al. Jupiter ICy moons Explorer (JUICE): an ESA mission to orbit Ganymede and to characterise the Jupiter system. *Planet. Space Sci.*, **78**, 1–21, 2013, DOI: 10.1016/j.pss.2012.12.002.

Gronoff, G., C. Simon Wedlund, C.J. Mertens, and R.J. Lillis. Computing uncertainties in ionosphere-airglow models: I. Electron flux and species production uncertainties for Mars. *J. Geophys. Res. [Space Phys.]*., **117**, A04306, 2012, DOI: 10.1029/2011JA016930.

Gurnett, D.A., W.S. Kurth, A. Roux, S.J. Bolton, and C.F. Kennel. Evidence for a magnetosphere at Ganymede from plasma-wave observations by the Galileo spacecraft. *Nature*, **384**, 535–537, 1996, DOI: 10.1038/384535a0.

Hall, D.T., P.D. Feldman, M.A. McGrath, and D.F. Strobel. The far-ultraviolet oxygen airglow of Europa and Ganymede. *Astrophys. J.*, **499**, 475–481, 1998, DOI: 10.1086/305604.

Heath, D.F., and B.M. Schlesinger. The Mg 280-nm doublet as a monitor of changes in solar ultraviolet irradiance. *J. Geophys. Res.*, **91**, 8672–8682, 1986, DOI: 10.1029/JD091iD08p08672.

Huebner, W.F., J.J. Keady, and S.P. Lyon. Solar photo rates for planetary atmospheres and atmospheric pollutants. *Astrophys. Space Sci.*, **195**, 1–289, 1992, DOI: 10.1007/BF00644558.

Jackman, C.M., and C.S. Arridge. Solar cycle effects on the dynamics of Jupiter's and Saturn's magnetospheres. *Sol. Phys.*, **274**, 481–502, 2011, DOI: 10.1007/s11207-011-9748-z.

Kalogerakis, K.S., T.G. Slanger, E.A. Kendall, T.R. Pedersen, M.J. Kosch, B. Gustavsson, and M.T. Rietveld. Remote Oxygen Sensing by Ionospheric Excitation (ROSIE). *Ann. Geophys.*, **27** (5), 2183–2189, 2009, DOI: 10.5194/angeo-27-2183-2009.

Kella, D., L. Vejby-Christensen, P. Johnson, H. Pedersen, and L. Andersen. The source of green light emission determined from a heavy-ion storage ring experiment. *Science*, **276**, 1530–1533, 1997.

Kivelson, M.G., K.K. Khurana, C.T. Russell, R.J. Walker, J. Warnecke, F.V. Coroniti, C. Polanskey, D.J. Southwood, and G. Schubert. Discovery of Ganymede's magnetic field by the Galileo spacecraft. *Nature*, **384**, 537–541, 1996, DOI: 10.1038/384537a0.

Kliore, A.J., D.P. Hinson, F.M. Flasar, A.F. Nagy, and T.E. Cravens. The ionosphere of Europa from Galileo radio occultations. *Science*, **277**, 355–358, 1997, DOI: 10.1126/science.277.5324.355.

Lilensten, J., T. Dudok de Wit, P.-O. Amblard, J. Aboudarham, F. Auchre, and M. Kretzschmar. Recommendation for a set of solar EUV lines to be monitored for aeronomy applications. *Ann. Geophys.*, **25** (6), 1299–1310, 2007, DOI: 10.5194/angeo-25-1299-2007.

Lilensten, J., T. Dudok de Wit, M. Kretzschmar, P. Amblard, S. Moussaoui, J. Aboudarham, and F. Auchère. Review on the solar spectral variability in the EUV for space weather purposes. *Ann. Geophys.*, **26**, 269–279, 2008, DOI: 10.5194/angeo-26-269-2008.

Lilensten, J., A.J. Coates, V. Dehant, T. Dudok de Wit, R.B. Horne, F. Leblanc, J. Luhmann, E. Woodfield, and M. Barthélemy. What characterizes planetary space weather? *Astron. Astrophys. Rev.*, **22**, 79, 2014, DOI: 10.1007/s00159-014-0079-6.

Marconi, M.L. A kinetic model of Ganymede's atmosphere, *Icarus*, **190**, 155–174, 2007, DOI: 10.1016/j.icarus.2007.02.016.

McGrath, M.A., X. Jia, K. Retherford, P.D. Feldman, D.F. Strobel, and J. Saur. Aurora on Ganymede. *J. Geophys. Res. [Space Phys.]*., **118**, 2043–2054, 2013, DOI: 10.1002/jgra.50122.

Nagy, A.F., J. Kim, T.E. Cravens, and A.J. Kliore. Hot oxygen corona at Europa. *Geophys. Res. Lett.*, **25**, 4153–4155, 1998, DOI: 10.1029/1998GL000139.

Paranicas, C., W.R. Paterson, A.F. Cheng, B.H. Mauk, R.W. McEntire, L.A. Frank, and D.J. Williams. Energetic particle observations near Ganymede. *J. Geophys. Res.*, **104**, 17459–17470, 1999, DOI: 10.1029/1999JA900199.

Pavlov, A., N. Pavlova, and A. Drozdov. Production rate of O(1D), O(1S) and N(2D) in the subauroral red arc region. *Geomag. Aeron.*, **39**, 201205, 1999.

Payan, A.P., C.S. Paty, and K.D. Retherford. Uncovering local magnetospheric processes governing the morphology and variability of Ganymede's aurora using three-dimensional multifluid simulations of Ganymede's magnetosphere. *J. Geophys. Res. [Space Phys.]*., **120**, 401–413, 2015, DOI: 10.1002/2014JA020301.

Plainaki, C., A. Milillo, A. Mura, S. Orsini, and T. Cassidy. Neutral particle release from Europa's surface. *Icarus*, **210**, 385–395, 2010, DOI: 10.1016/j.icarus.2010.06.041.

Plainaki, C., A. Milillo, A. Mura, J. Saur, S. Orsini, and S. Massetti. Exospheric O₂ densities at Europa during different orbital phases. *Planet. Space Sci.*, **88**, 42–52, 2013, DOI: 10.1016/j.pss.2013.08.011.

Plainaki, C., A. Milillo, S. Massetti, A. Mura, X. Jia, S. Orsini, V. Mangano, E. De Angelis, and R. Rispoli. The H₂O and O₂ exospheres of Ganymede: the result of a complex interaction between the Jovian magnetospheric ions and the icy moon. *Icarus*, **245**, 306–319, 2015, DOI: 10.1016/j.icarus.2014.09.018.

Porco, C.C., R.A. West, A. McEwen, A.D. Del Genio, A.P. Ingersoll, et al. Cassini imaging of Jupiter's atmosphere, satellites, and rings. *Science*, **299**, 1541–1547, 2003, DOI: 10.1126/science.1079462.

Qian, L., A.G. Burns, P.C. Chamberlin, and S.C. Solomon. Flare location on the solar disk: modeling the thermosphere and ionosphere response. *J. Geophys. Res. [Space Phys.]*., **115** (A9), A09311, 2010, DOI: 10.1029/2009JA015225.

Rosen, S., A. Derkatch, J. Semaniak, A. Neau, A. Al-Khalili, et al. Recombination of simple molecular ions studied in storage ring: dissociative recombination of H₂O⁺. *Faraday Discussion*, **115**, 295–302, 2000, DOI: 10.1039/A909314A.

Roth, L., J. Saur, K.D. Retherford, D.F. Strobel, P.D. Feldman, M.A. McGrath, and F. Nimmo. Transient water vapor at Europa's south pole. *Science*, **343**, 171–174, 2014, DOI: 10.1126/science.1247051.

Rottman, G. The SORCE mission. *Sol. Phys.*, **230**, 7–25, 2005, DOI: 10.1007/s11207-005-8112-6.

Saur, J., D.F. Strobel, and F.M. Neubauer. Interaction of the Jovian magnetosphere with Europa: constraints on the neutral atmosphere. *J. Geophys. Res.*, **103**, 19947–19962, 1998, DOI: 10.1029/97JE03556.

Saur, J., S. Duling, L. Roth, X. Jia, D.F. Strobel, et al. The search for a subsurface ocean in Ganymede with Hubble Space Telescope observations of its auroral ovals. *J. Geophys. Res. [Space Phys.]*., **120**, 1715–1737, 2015, DOI: 10.1002/2014JA020778.

Schmidt, H., R. Wegmann, W. Huebner, and D. Boice. Cometary gas and plasma flow with detailed chemistry. *Comput. Phys. Commun.*, **49** (1), 17–59, 1988, DOI: 10.1016/0010-4655(88)902014-7.

Shematovich, V.I., R.E. Johnson, J.F. Cooper, and M.C. Wong. Surface-bounded atmosphere of Europa, *Icarus*, **173**, 480–498, 2005, DOI: 10.1016/j.icarus.2004.08.013.

Slanger, T.G., P.C. Cosby, B.D. Sharpee, K.R. Minschwaner, and D.E. Siskind. O(1S → ¹D,³P) branching ratio as measured in the terrestrial nightglow. *J. Geophys. Res. [Space Phys.]*., **111**, A12318, 2006, DOI: 10.1029/2006JA011972.

Smyth, W.H., and M.L. Marconi. Europa's atmosphere, gas tori, and magnetospheric implications. *Icarus*, **181**, 510–526, 2006, DOI: 10.1016/j.icarus.2005.10.019.

Stuhl, F., and K.H. Welge. Deactivation of O(1S) and O₂(b¹Σg⁺). *Can. J. Chem.*, **47** (10), 1870–1877, 1969, DOI: 10.1139/v69-306.
Tapping, K.F., and B. Detracey. The origin of the 10.7 cm flux. *Sol. Phys.*, **127**, 321–332, 1990, DOI: **10.1007/BF00152171**.

Thirupathaiah, P., and V. Singh. An updated model of atomic oxygen redline dayglow emission. *Adv. Space Res.*, **54** (6), 939–945, 2014, DOI: **10.1016/j.asr.2014.05.022**.

Thuillier, G., L. Floyd, T.N. Woods, R. Cebula, E. Hilsenrath, M. Hersé, and D. Labs. Solar Irradiance Reference Spectra. In: J.M. Pap, et al., Editors. *Solar Variability and its Effects on Climate. Geophysical*, Vol. 141 of Washington DC American Geophysical Union Geophysical Monograph Series, 171, 2004, DOI: **10.1029/GM141**.

Turc, L., L. Leclercq, F. Leblanc, R. Modolo, and J.-Y. Chaufray. Modelling Ganymede’s neutral environment: a 3D test-particle simulation. *Icarus*, **229**, 157–169, 2014, DOI: **10.1016/j.icarus.2013.11.005**.

Wiese, W.L., J.R. Fuhr, and T.M. Deters. Atomic transition probabilities of carbon, nitrogen, and oxygen: a critical data compilation. *Monograph No. 7 of the Journal of Physical and Chemical Reference Data*, American Chemical Society for the

National Institute of Standards and Technology, Washington DC, ISBN: 1-56396-602-6, 1996.

Witasse, O., J. Lilensten, C. Lathuillère, and P.-L. Bleilly. Modeling the OI 630.0 and 557.7 nm thermospheric dayglow during EISCAT-WINDII coordinated measurements. *J. Geophys. Res.*, **104**, 24639–24656, 1999, DOI: **10.1029/1999JA900260**.

Woods, T.N., F.G. Eparvier, S.M. Bailey, P.C. Chamberlin, J. Lean, G.J. Rottman, S.C. Solomon, W.K. Tobiska, and D.L. Woodraska. Solar EUV Experiment (SEE): mission overview and first results. *J. Geophys. Res. [Space Phys.]*, **110**, A01312, 2005, DOI: **10.1029/2004JA010765**.

Woods, T.N., P.C. Chamberlin, W.K. Peterson, R.R. Meier, P.G. Richards, et al. XUV Photometer System (XPS): improved solar irradiance algorithm using CHIANTI spectral models. *Sol. Phys.*, **250**, 235–267, 2008.

Zipf, E.C. The collisional deactivation of metastable atoms and molecules in the upper atmosphere. *Can. J. Chem.*, **47** (10), 1863–1870, 1969, DOI: **10.1139/v69-305**.

**Cite this article as:** Cessateur G, Barthelemy M & Peinke I. Photochemistry-emission coupled model for Europa and Ganymede. *J. Space Weather Space Clim.*, **6**, A17, 2016, DOI: **10.1051/swsc/2016009**.