Io's auroral emissions via global hybrid plasma simulations

Štverák, Štěpán; Trávníček, Pavel; Šebek, Ondřej; Herčík, David

doi:https://doi.org/10.5194/angeo-2021-10

Preprints

https://doi.org/10.5194/angeo-2021-10

Preprints

12 Feb 2021

Review status: this preprint has been withdrawn by the authors.

Io's auroral emissions via global hybrid plasma simulations

Štěpán Štverák^1,2, Pavel Trávníček³, Ondřej Šebek¹, and David Herčík¹ Štěpán Štverák et al.

Show author details

¹Institute of Atmospheric Physics, CAS, Prague, Czech Republic
²Astronomical Institute, CAS, Ond ˇ rejov, Czech Republic
³Space Sciences Laboratory, UCB, Berkeley, CA, USA

Received: 28 Jan 2021 – Accepted for review: 09 Feb 2021 – Discussion started: 12 Feb 2021

Abstract. We tackle the Io's aurora source and topology by carrying out a set of global hybrid simulations of Io's interaction with the plasma torus under different model geometry and background conditions. Based on the simulated results, we compute the photon emission rates above the Io's surface and present the resulting images from a virtual telescope and topological maps showing the distribution of the emission sources across the moon's surface. This allows us to compare the structure of the aurora with the real observations and conclude on the different assumptions. We found a reasonable agreement with the real observations in the case of non-collisional background electron populations. From the comparison of the local magnetic field topology with model aurora structures, we also infer that an induced dipole feature is more probable to play a role in the interaction of Io with the Jovian magnetosphere. In addition we also examine the potential contribution of energetic electron beams, being observed in the Io's wake region, to the overall auroral emissions.

This preprint has been withdrawn.

Download & links

Preprint (PDF, 2325 KB)

Withdrawal notice
This preprint has been withdrawn.
Preprint (2325 KB)

Download & links

This preprint has been withdrawn.

Preprint (2325 KB)
BibTeX
EndNote

Štěpán Štverák et al.

Interactive discussion

Status: closed

RC1:
'Comment on angeo-2021-10', Anonymous Referee #1, 23 Feb 2021

see attached pdf file.

Citation: https://doi.org/10.5194/angeo-2021-10-RC1
- AC1: 'Reply to Editor on RCs', Stepan Stverak, 19 May 2021
  
  Dear Editor,
  Two very detailed referee's reports were provided to the manuscript "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. with a common recommendation for major revisions. Although the received comments pose a rather strong criticism of the manuscript, we are very thankful for such feedback showing us the weaknesses of the performed analysis and the way forward how the paper can be still significantly improved. Since the performed analysis was partly limited on the already existing simulation data set (Šebek et al., 2019), we see a necessity to perform as a first step a new set of dedicated simulation runs to correctly address and analyse the role of the individual plasma processes on the Io's interaction with the magnetospheric environment and thus on the generation of the aurora emissions. This is, however, a highly time demanding task (order of months) and we will therefore rather ask to withdraw the recent submission and later resubmit the fully reworked manuscript once available. We believe that some of the (minor) comments could be either argued or implemented into the manuscript already now, anyway, we prefer to do so after a careful revision of the whole manuscript and solving the major revisions first.
  Kind regards,
  Štěpán Štverák (on behalf of all co-authors)
  
  Citation: https://doi.org/10.5194/angeo-2021-10-AC1
RC2:
'Comment on angeo-2021-10', Anonymous Referee #2, 14 Apr 2021

The paper "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. used hybrid

simulations to investigate the aurora morphology of Jupiter moon Io. Unfortunately, there are numerous shortcomings and misconceptions as well as a lack of quantitative results in my view. Furthermore, I doubt the use of the approach of hybrid simulations (with an isothermal electron fluid and treating ions as particles) to study the issue of Io's aurora. The aurora is excited by collisions between electrons and neutrals in atmospheric regions and hardly/not affected by ion kinetic effects.

I will discuss some of the shortcoming below but these are not all. (I do not itend to be overly critical and try to keep it to the facts.)

The treatment of the electrons with two populations of old (upstream torus) and new (from ionization of atmosphere) electrons is oversimplifying the complex issue of the exchange of energy between the electrons. The study distinguishes between two cases where these two populations exchange energy instantaneously or not at all. In reality, the electron heat conductivity along the magnetic field in the plasma (like the electric conductivity) is almost infinity except where the neutral density becomes comparable (Banks, P.M., Kockarts, G., 1973. Aeronomy, volume A. Academic Press, San Diego CA). This means that the electron energy distribution around Io will be complex and one also needs to considers that electrons far away can exchange energy with the close electrons as long as they are on one magnetic field line.

The simplified treatment leads to wrong distribution of the electron energy around Io and thus wrong aurora morphology.

It is unclear how the scaling of the simulated system in the hybrid model is set up. In their previous studies (Sebek et al. 2016, 2019), Io was downscaled by a factor of 10 in order to make the hybrid approach possible. How is his treated here and how does it work to compare aurora images on real scales?

(And what did you assume for the atmospheric scale height (L119)?)

There are significant general mismatches between the simulated aurora presented here and the observations. The authors yet cherry-pick specific features and make hand-waiving arguments (with strange confusion of longitude and latitude!) about similarities and even conclude about presence of induction in an ocean.

A few examples for mismatches:

The simulated results maximize in the region upstream of Io. (Which the authors also state, e.g., in line 2017: "that most of the emissions are produced in the upstream hemisphere"). In the observations the emissions on the upstream central hemisphere are, however, very faint. There are a few FUV observations with a useful geometry for illustrate this. For example, in Figure 3 of Roth et al. 2014, bottom panel first row, the OI13565 A image marked with "38" in the lower left corner. This view shows the sub-jovian hemisphere and upstream is to the right (similar but vertically flipped view as VIEW-B in this paper). Here it is obvious that there are hardly any detectable emissions upstream of the moon. In the simulation here, in contrast, (VIEW-B in Figure 2 here) the upstream emissions is much much brighter than everything else. (Only the beam-simulated wake feature in the simulation reaches somewhat similar brightness).

The shape of the flank emissions on the anti-Jovian side is very extended in the simulation, but is just as confined as on the sub-jovian flank in the observations. The authors mention this, but they do not discuss this.

The overall mismatch definitely prevents conclusions on details like the role of induction in the simulation.

The paper completely ignores our work published in Roth et al. 2014 on the aurora observations. The Roth et al. 2014 study used all existing FUV HST observations of the aurora for the first time and presented simple mathematical equations that describe the aurora distribution in all these images with one model.

The Roth et al. 2014 study thus exactly provides high-level data in a format to be used in studies numerical simulations of Io's interaction like the one presented!

One issue related to the negligence of the Roth et al. (2014) results is the assumption that Io's FUV aurora contains a feature in the wake, denoted here and before "wake emission". In the Roth et al. 2014 study, it is explained why none of the emissions in the data can unambiguously assigned to be located in the wake and that, in contrary, features that appear to be in the wake in the 2d images originate from the flanks of Io.

The assumption of aurora by an electron beam is obsolete. There is no FUV wake feature.

The literature review in the introduction and entire study is incomplete and outdated, missing some key findings and papers. This becomes already evident by the years of the cited publication, all of which date more than 6 years back although there were relevant papers in more recent years.

One example in L36: The discussion of the volcanic vs sublimation supply totally ignores the latest literature and findings: Tsang et al. 2016 J. Geophys. Res. Planets 121, and de Kleer et al. 2019 Icarus 317, to name only two important new results.

Also, there is no reference at all to results of the Juno mission which is now almost half a decade at Jupiter. For example, they authors cite Bonfond et al. 2008 on the Io footprint morphology but there have been high-resolution images by Juno which revealed numerous new features in the footprint (Mura et al. 2020). In addition, there are several follow-up papers by Bonfond in Io's footprint morphology, which should be mentioned.

The study of Roth et al. 2017 is described in a misleading and actually wrong way, both regarding the results and the method.

It is a data analysis study where new observations of the aurora are presented and then analyzed with various methods considering different interior effects. One of the methods include MHD simulations to generally investigate the behaviors of tangential magnetic field lines. This does not make it an "MHD study", as claimed in several places.

The result is indeed that a magma ocean is inconsistent with the data, but also that induction in Io's core is well consistent. So the sentence "They showed that magnetic field without contribution of induced dipole field corresponds better to assumed locations of auroral equatorial spots." is just wrong.

Most importantly, the Roth et al. 2017 study uses the time-variability of the aurora as one of the two main diagnostics for magnetic induction. The lack of a *temporal phase shift* is the main argument against the presence of induction near the surface. This is completely ignored here.

Some line-by-line comments:

L173: "we note that local rate of any typical emission observed from Io ... would have similar structure and would scale according to magnitude of efficiency of given excitation reaction."

That is not correct. Some of the emissions related to allowed electronic transitions in the atoms and can be optically thick and proper radiative transfer modeling would be needed. Therefore, many studies use the forbidden transitions, e.g., OI1356 A.

L262: "...are to certain extent in fair agreement with the real observations "

I disagree and think it requires any kind of quantiative comparison (not only visual, although I think even the visual comparison shows mostly obvious diasgreements) to claim any agreement.

L284-L286: The spot latitudes derived in Retherford et al. 2000 are off for some images due to observational limitations and unfortunate conditions/geometries. The Roth et al. 2017 study used targeted new observations in 2013 and 2014 intentionally set up to optimize the geometries in order to properly address the issue of the spot locations. It is really a sneaky (and wrong) argument to use these old results (from basically same authors actually as in Roth et al. 2017) to spread doubts about the Roth et al. 2017 results.

And, btw, Retherford et al. 2000 is here cited with "e.g.". What are the other studies *quantifying* scattered spots?

L297: "are here located at rather low latitudes while inclusion of the induced dipole moves the spots to φ ≈ −30â¦ "

Here the authors confuse latitude and longitude. The angle phi derived in Retherford et al. 2000 is a LONGITUDE.

Citation: https://doi.org/10.5194/angeo-2021-10-RC2
- AC2: 'Reply to Editor on RCs', Stepan Stverak, 19 May 2021
  
  Dear Editor,
  Two very detailed referee's reports were provided to the manuscript "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. with a common recommendation for major revisions. Although the received comments pose a rather strong criticism of the manuscript, we are very thankful for such feedback showing us the weaknesses of the performed analysis and the way forward how the paper can be still significantly improved. Since the performed analysis was partly limited on the already existing simulation data set (Šebek et al., 2019), we see a necessity to perform as a first step a new set of dedicated simulation runs to correctly address and analyse the role of the individual plasma processes on the Io's interaction with the magnetospheric environment and thus on the generation of the aurora emissions. This is, however, a highly time demanding task (order of months) and we will therefore rather ask to withdraw the recent submission and later resubmit the fully reworked manuscript once available. We believe that some of the (minor) comments could be either argued or implemented into the manuscript already now, anyway, we prefer to do so after a careful revision of the whole manuscript and solving the major revisions first.
  Kind regards,
  Štěpán Štverák (on behalf of all co-authors)
  
  Citation: https://doi.org/10.5194/angeo-2021-10-AC2

Interactive discussion

Status: closed

RC1:
'Comment on angeo-2021-10', Anonymous Referee #1, 23 Feb 2021

see attached pdf file.

Citation: https://doi.org/10.5194/angeo-2021-10-RC1
- AC1: 'Reply to Editor on RCs', Stepan Stverak, 19 May 2021
  
  Dear Editor,
  Two very detailed referee's reports were provided to the manuscript "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. with a common recommendation for major revisions. Although the received comments pose a rather strong criticism of the manuscript, we are very thankful for such feedback showing us the weaknesses of the performed analysis and the way forward how the paper can be still significantly improved. Since the performed analysis was partly limited on the already existing simulation data set (Šebek et al., 2019), we see a necessity to perform as a first step a new set of dedicated simulation runs to correctly address and analyse the role of the individual plasma processes on the Io's interaction with the magnetospheric environment and thus on the generation of the aurora emissions. This is, however, a highly time demanding task (order of months) and we will therefore rather ask to withdraw the recent submission and later resubmit the fully reworked manuscript once available. We believe that some of the (minor) comments could be either argued or implemented into the manuscript already now, anyway, we prefer to do so after a careful revision of the whole manuscript and solving the major revisions first.
  Kind regards,
  Štěpán Štverák (on behalf of all co-authors)
  
  Citation: https://doi.org/10.5194/angeo-2021-10-AC1
RC2:
'Comment on angeo-2021-10', Anonymous Referee #2, 14 Apr 2021

The paper "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. used hybrid

simulations to investigate the aurora morphology of Jupiter moon Io. Unfortunately, there are numerous shortcomings and misconceptions as well as a lack of quantitative results in my view. Furthermore, I doubt the use of the approach of hybrid simulations (with an isothermal electron fluid and treating ions as particles) to study the issue of Io's aurora. The aurora is excited by collisions between electrons and neutrals in atmospheric regions and hardly/not affected by ion kinetic effects.

I will discuss some of the shortcoming below but these are not all. (I do not itend to be overly critical and try to keep it to the facts.)

The treatment of the electrons with two populations of old (upstream torus) and new (from ionization of atmosphere) electrons is oversimplifying the complex issue of the exchange of energy between the electrons. The study distinguishes between two cases where these two populations exchange energy instantaneously or not at all. In reality, the electron heat conductivity along the magnetic field in the plasma (like the electric conductivity) is almost infinity except where the neutral density becomes comparable (Banks, P.M., Kockarts, G., 1973. Aeronomy, volume A. Academic Press, San Diego CA). This means that the electron energy distribution around Io will be complex and one also needs to considers that electrons far away can exchange energy with the close electrons as long as they are on one magnetic field line.

The simplified treatment leads to wrong distribution of the electron energy around Io and thus wrong aurora morphology.

It is unclear how the scaling of the simulated system in the hybrid model is set up. In their previous studies (Sebek et al. 2016, 2019), Io was downscaled by a factor of 10 in order to make the hybrid approach possible. How is his treated here and how does it work to compare aurora images on real scales?

(And what did you assume for the atmospheric scale height (L119)?)

There are significant general mismatches between the simulated aurora presented here and the observations. The authors yet cherry-pick specific features and make hand-waiving arguments (with strange confusion of longitude and latitude!) about similarities and even conclude about presence of induction in an ocean.

A few examples for mismatches:

The simulated results maximize in the region upstream of Io. (Which the authors also state, e.g., in line 2017: "that most of the emissions are produced in the upstream hemisphere"). In the observations the emissions on the upstream central hemisphere are, however, very faint. There are a few FUV observations with a useful geometry for illustrate this. For example, in Figure 3 of Roth et al. 2014, bottom panel first row, the OI13565 A image marked with "38" in the lower left corner. This view shows the sub-jovian hemisphere and upstream is to the right (similar but vertically flipped view as VIEW-B in this paper). Here it is obvious that there are hardly any detectable emissions upstream of the moon. In the simulation here, in contrast, (VIEW-B in Figure 2 here) the upstream emissions is much much brighter than everything else. (Only the beam-simulated wake feature in the simulation reaches somewhat similar brightness).

The shape of the flank emissions on the anti-Jovian side is very extended in the simulation, but is just as confined as on the sub-jovian flank in the observations. The authors mention this, but they do not discuss this.

The overall mismatch definitely prevents conclusions on details like the role of induction in the simulation.

The paper completely ignores our work published in Roth et al. 2014 on the aurora observations. The Roth et al. 2014 study used all existing FUV HST observations of the aurora for the first time and presented simple mathematical equations that describe the aurora distribution in all these images with one model.

The Roth et al. 2014 study thus exactly provides high-level data in a format to be used in studies numerical simulations of Io's interaction like the one presented!

One issue related to the negligence of the Roth et al. (2014) results is the assumption that Io's FUV aurora contains a feature in the wake, denoted here and before "wake emission". In the Roth et al. 2014 study, it is explained why none of the emissions in the data can unambiguously assigned to be located in the wake and that, in contrary, features that appear to be in the wake in the 2d images originate from the flanks of Io.

The assumption of aurora by an electron beam is obsolete. There is no FUV wake feature.

The literature review in the introduction and entire study is incomplete and outdated, missing some key findings and papers. This becomes already evident by the years of the cited publication, all of which date more than 6 years back although there were relevant papers in more recent years.

One example in L36: The discussion of the volcanic vs sublimation supply totally ignores the latest literature and findings: Tsang et al. 2016 J. Geophys. Res. Planets 121, and de Kleer et al. 2019 Icarus 317, to name only two important new results.

Also, there is no reference at all to results of the Juno mission which is now almost half a decade at Jupiter. For example, they authors cite Bonfond et al. 2008 on the Io footprint morphology but there have been high-resolution images by Juno which revealed numerous new features in the footprint (Mura et al. 2020). In addition, there are several follow-up papers by Bonfond in Io's footprint morphology, which should be mentioned.

The study of Roth et al. 2017 is described in a misleading and actually wrong way, both regarding the results and the method.

It is a data analysis study where new observations of the aurora are presented and then analyzed with various methods considering different interior effects. One of the methods include MHD simulations to generally investigate the behaviors of tangential magnetic field lines. This does not make it an "MHD study", as claimed in several places.

The result is indeed that a magma ocean is inconsistent with the data, but also that induction in Io's core is well consistent. So the sentence "They showed that magnetic field without contribution of induced dipole field corresponds better to assumed locations of auroral equatorial spots." is just wrong.

Most importantly, the Roth et al. 2017 study uses the time-variability of the aurora as one of the two main diagnostics for magnetic induction. The lack of a *temporal phase shift* is the main argument against the presence of induction near the surface. This is completely ignored here.

Some line-by-line comments:

L173: "we note that local rate of any typical emission observed from Io ... would have similar structure and would scale according to magnitude of efficiency of given excitation reaction."

That is not correct. Some of the emissions related to allowed electronic transitions in the atoms and can be optically thick and proper radiative transfer modeling would be needed. Therefore, many studies use the forbidden transitions, e.g., OI1356 A.

L262: "...are to certain extent in fair agreement with the real observations "

I disagree and think it requires any kind of quantiative comparison (not only visual, although I think even the visual comparison shows mostly obvious diasgreements) to claim any agreement.

L284-L286: The spot latitudes derived in Retherford et al. 2000 are off for some images due to observational limitations and unfortunate conditions/geometries. The Roth et al. 2017 study used targeted new observations in 2013 and 2014 intentionally set up to optimize the geometries in order to properly address the issue of the spot locations. It is really a sneaky (and wrong) argument to use these old results (from basically same authors actually as in Roth et al. 2017) to spread doubts about the Roth et al. 2017 results.

And, btw, Retherford et al. 2000 is here cited with "e.g.". What are the other studies *quantifying* scattered spots?

L297: "are here located at rather low latitudes while inclusion of the induced dipole moves the spots to φ ≈ −30â¦ "

Here the authors confuse latitude and longitude. The angle phi derived in Retherford et al. 2000 is a LONGITUDE.

Citation: https://doi.org/10.5194/angeo-2021-10-RC2
- AC2: 'Reply to Editor on RCs', Stepan Stverak, 19 May 2021
  
  Dear Editor,
  Two very detailed referee's reports were provided to the manuscript "Io’s auroral emissions via global hybrid plasma simulations" by Štverák et al. with a common recommendation for major revisions. Although the received comments pose a rather strong criticism of the manuscript, we are very thankful for such feedback showing us the weaknesses of the performed analysis and the way forward how the paper can be still significantly improved. Since the performed analysis was partly limited on the already existing simulation data set (Šebek et al., 2019), we see a necessity to perform as a first step a new set of dedicated simulation runs to correctly address and analyse the role of the individual plasma processes on the Io's interaction with the magnetospheric environment and thus on the generation of the aurora emissions. This is, however, a highly time demanding task (order of months) and we will therefore rather ask to withdraw the recent submission and later resubmit the fully reworked manuscript once available. We believe that some of the (minor) comments could be either argued or implemented into the manuscript already now, anyway, we prefer to do so after a careful revision of the whole manuscript and solving the major revisions first.
  Kind regards,
  Štěpán Štverák (on behalf of all co-authors)
  
  Citation: https://doi.org/10.5194/angeo-2021-10-AC2

Štěpán Štverák et al.

Data sets

Io hybrid simulations Štěpán Štverák http://doi.org/10.5281/zenodo.4472602

Štěpán Štverák et al.

Viewed

Total article views: 741 (including HTML, PDF, and XML)

HTML	PDF	XML	Total	BibTeX	EndNote
511	216	14	741	7	5

HTML: 511
PDF: 216
XML: 14
Total: 741
BibTeX: 7
EndNote: 5

Views and downloads (calculated since 12 Feb 2021)

Month	HTML	PDF	XML	Total
Feb 2021	194	29	5	228
Mar 2021	180	7	0	187
Apr 2021	33	9	3	45
May 2021	30	11	5	46
Jun 2021	11	3	0	14
Jul 2021	12	7	0	19
Aug 2021	7	2	0	9
Sep 2021	9	15	0	24
Oct 2021	10	74	0	84
Nov 2021	11	46	0	57
Dec 2021	8	6	1	15
Jan 2022	6	7	0	13

Cumulative views and downloads (calculated since 12 Feb 2021)

Month	HTML	PDF	XML	Total
Feb 2021	194	29	5	228
Mar 2021	180	7	0	187
Apr 2021	33	9	3	45
May 2021	30	11	5	46
Jun 2021	11	3	0	14
Jul 2021	12	7	0	19
Aug 2021	7	2	0	9
Sep 2021	9	15	0	24
Oct 2021	10	74	0	84
Nov 2021	11	46	0	57
Dec 2021	8	6	1	15
Jan 2022	6	7	0	13

Viewed (geographical distribution)

Total article views: 662 (including HTML, PDF, and XML) Thereof 662 with geography defined and 0 with unknown origin.

Country	#	Views	%

Latest update: 22 Jan 2022

Short summary

We present a study of Io's aurora patters based on global numerical simulations of the moon's interaction with the Jovian magnetospheric plasma torus. We demonstrate how different configurations of plasma interactions with the neutral atmosphere affect the resulting aurora's morphology. Comparison to real remote observations indicate existence of Io's induced magnetic dipole and low collisionality of thermal electrons. The conducted research was supported by Czech Science Foundation.


Total:	0
HTML:	0
PDF:	0
XML:	0