Planetary Ionospheres Research Articles

An Investigation on the Distribution of Martian Ionospheric Particles, Based on the Mars Atmosphere and Volatile Evolution (MAVEN)

In this paper, we use the key parameters data set of the Neutral Gas and Ion Mass Spectrometer from the Mars Atmosphere and Volatile Evolution (MAVEN) mission. The particle density profiles of electrons, CO2+/N2+, CO+, O2+, O+, NO+, O2 and O from 90 to 500 km have been deduced by adopting the Chapman modeling methodology. The correlation of the peak density/altitude with the solar zenith angle, the changes in the profile of the Martian ionosphere during solar flares, and the effects of Martian dust storms are analyzed. The results exhibit a positive/negative correlation between the peak density/altitude of the M2 layer and the solar zenith angle. Within the MAVEN observational record available, only three C-Class flares occurred on 26 August 2016, 29 November 2020, and 26 August 2021. The analysis reveals during these solar flare events, the electron density of the M2 layer above 200 km increases obviously. The peak density of M1 increases by 33.4%, 13.2% and 7.4%, while the peak height decreases by 0.1%, 10.2% and 4.4%, respectively. The Martian dust storm causes the peak height of the M2 layer to increase by 19.5 km, and the peak density to decrease by 4.2 × 109 m−3. Our study shows that the Martian ionosphere is similar to the Earth’s, which is of great significance for understanding the planetary ionosphere.

The Infrared Auroral Footprint Tracks of Io, Europa and Ganymede at Jupiter Observed by Juno‐JIRAM

AbstractThe electromagnetic coupling between the Galilean satellites at Jupiter and the planetary ionosphere generates an auroral footprint, which is detected with high spatial resolution in the infrared L band by the Jovian InfraRed Auroral Mapper (JIRAM) onboard the Juno spacecraft. We report the JIRAM data acquired since 27 August 2016 until 23 May 2022, which are used to compute the average position of the footprint tracks of Io, Europa and Ganymede. The result of the present analysis help to test the reliability of magnetic field models, to calibrate ground‐based observations and to highlight the variability in the footprint positions, which can be used to probe the plasma environment at the orbit of the satellites. The determination of the plasma properties around the moons is particularly relevant to complement the Juno flybys of the moons during its extended mission, and to support the future Juice and Europa Clipper missions. Lastly, we report no clear evidence of the auroral footprint of Callisto, which is likely due to a combination of its low expected brightness and its position very close to the main Jovian aurora.

Characteristic Timescales for the Dayside Martian Ionosphere: Chemistry, Diffusion, and Magnetization

Different boundaries could be defined in a planetary ionosphere where the dominant process in function switches from one to the other. Identifying these boundaries and understanding their variations are hence crucial for disentangling the complexity of the ionosphere. Focusing on Mars, we perform a data-driven analysis of various boundaries and the associated time constants based on the multi-instrument measurements made by the Mars Atmosphere and Volatile Evolution mission during six campaigns that sample broadly different internal and external conditions. The boundaries we investigate include the photochemical equilibrium (PCE) boundary, the magnetic frozen boundary, the ion collision boundary, and the ion gyration boundary. Our analysis reveals systematic solar cycle and diurnal variations in that all boundaries tend to be elevated at enhanced solar activity and on the dayside and duskside of Mars. The variations with the magnetic environment are not observed for all boundaries except for the PCE boundary that exhibits an obvious elevation in strongly magnetized regions. Finally, our analysis suggests interesting species-dependent variations of different boundaries. In particularly, the PCE boundary shows the largest variability among all, with reduced boundary locations for all terminal species (NO+, HCO+, , and H3O+) and one extra nonterminal species (CO2 +) owing to different chemical properties rendered by different ions.

Unified fluid theory of the collisional thermal Farley–Buneman instability including magnetized multi-species ions

This paper develops a unified linear theory of cross field plasma instabilities, including the Farley–Buneman, electron thermal, and ion thermal instabilities, in spatially uniform collisional plasmas with partially unmagnetized multi-species ions. Collisional plasma instabilities in weakly ionized, highly dissipative, weakly magnetized plasmas play an important role in the lower Earth's ionosphere and may be of importance in other planetary ionospheres, stellar atmospheres, cometary tails, molecular clouds, accretion disks, etc. In the Earth's ionosphere, these collisional plasma instabilities cause intense electron heating. In the solar chromosphere, they can do the same—an effect originally suggested from spectroscopic observations and modeling. Based on a simplified 5-moment multi-fluid model, the theoretical analysis presented in this paper produces the linear dispersion relation for the combined Thermal Farley–Buneman Instability with an important long-wavelength limit analyzed in detail. This limit provides an easy interpretation of different instability drivers and wave dissipation. This analysis of instability, combined with simulations, will enable us to better understand plasma waves and turbulence in these commonly occurring collisional space plasmas.

Experimental Study of the Response of Sheared E × B Flow to Varying Ion–Neutral Collisions

Partially ionized plasma is a common occurrence in astrophysical and space environments. The emergence and development of plasma instabilities are significantly impacted by the inelastic collisions between the ions and neutrals in the partially ionized plasma, such as the charge exchange. In this study, the effect of the ion–neutral collisions on the sheared E × B flow was experimentally investigated. In the weak collision range, the shear-driven plasma instability, such as Kelvin–Helmholtz instability, was excited by the velocity-sheared flow. However, increasing ion–neutral collisions resulted in a decrease in the magnitude of the sheared E × B flow due to charge exchange–induced drag forces. Consequently, the Kelvin–Helmholtz instability is suppressed, and the Rayleigh–Taylor instability is triggered. The underlying mechanism was elucidated through experimental findings and numerical analysis. The result of this study proposes that a transition between the two modes occurred with increasing ion–neutral collision strength. It could be applied to the study of the solar chromosphere and prominence and planetary ionospheres, where plasma is partially ionized and the sheared E × B flow is often encountered.

Simulating Secondary Electron and Ion Emission from the Cassini Spacecraft in Saturn’s Ionosphere

The Cassini spacecraft’s Grand Finale flybys through Saturn’s ionosphere provided unprecedented insight into the composition and dynamics of the gas giant’s upper atmosphere and a novel and complex spacecraft–plasma interaction. In this article, we further study Cassini’s interaction with Saturn’s ionosphere using three-dimensional particle-in-cell simulations. We focus on how electrons and ions, emitted from spacecraft surfaces due to the high-velocity impact of atmospheric water molecules, could have affected the spacecraft potential and low-energy plasma measurements. The simulations show emitted electrons extend upstream along the magnetic field, and for sufficiently high emission rates, charge the spacecraft to positive potentials. The lack of accurate emission rates and characteristics, however, makes differentiation between the prominence of secondary electron emission and ionospheric charged dust populations, which induce similar charging effects, difficult for Cassini. These results provide further context for Cassini’s final measurements and highlight the need for future laboratory studies to support high-velocity flyby missions through planetary and cometary ionospheres.

What Is the Altitude of Thermal Equilibrium?

AbstractThermal equilibrium in planetary atmospheres occurs at altitudes where the ion, electron, and neutral temperatures are equal. Thermal equilibrium is postulated to occur in the collision‐dominated ionosphere. This postulated altitude is above the lower boundary of all empirical models of planetary ionospheres. Physics‐based model predictions of the altitude cannot be validated due to a lack of adequate simultaneous observations of temperature profiles. This study presents temperature profiles from simultaneous observations on Atmosphere Explorer–C below 140 km and quiet‐time neutral observations from Thermosphere Ionosphere Mesosphere Energy and Dynamics/Global UltraViolet Imager over Millstone Hill. These are compared with profiles from physics‐based models with a discussion of their respective limitations. We conclude that there does not yet exist a quantitative understanding of the ion, electron, and neutral thermalization processes in low‐altitude planetary ionospheres. Progress on this topic requires an adequate database of simultaneous ion, electron, and neutral temperature profiles in the 110–140 km altitude range.

Impact of Heavy Ions on the Structure of Current Sheets in the Gravity Field of Exoplanets and Stars

Boundary current sheets (CSs) can be formed in collisionless space plasmas in the environment of exoplanets and cold stars. Usually they represent curved surfaces carrying the electric current analogous to the well-known planetary ionospheres, magnetopauses, or stellar coronas surrounding celestial bodies. At smaller local scales, some of them can be imagined as planar current layers of a finite scale located parallel to the surface of a celestial object and, correspondingly, perpendicular to the direction of the gravitational force. In some cases, this force crossed with magnetic field can influence the dynamics of charged particles in CSs and substantially change the structure of both the current layer and the magnetic field. We have generalized our prior model, taking into account a multi-ion plasma composition and a magnetic field configuration with a shear. It is shown that, due to the drift motion of plasma particles in the crossed gravitational and magnetic fields, the structure of CSs becomes more complex, accruing asymmetric and shifted profiles of the current and plasma densities that depend on dominating current carriers and the characteristics of the magnetic shear. We discuss possible applications of the results to the interpretation of observations of boundary layers in different space plasmas.

The three-fluid generalized Ohm's law: A theoretical study

The two-fluid generalized Ohm's law (GOL) is based on the assumption that plasma is composed of only protons and electrons. The three-fluid GOL is obtained theoretically for the three-fluid plasma consisting of heavy ions, light ions, and electrons, which prevails in planetary ionospheres and magnetospheres. Three inertial lengths corresponding to the three-scale diffusion region in the three-fluid magnetic reconnection are derived. The ion inertial lengths and reconnection rate as well as the Hall magnetic and electric fields are modified due to the two-step decoupling process of ions. Our results provide a framework to extend the reconnection theory for even more ion species.

A Significant Sudden Ionospheric Disturbance Associated with Gamma-Ray Burst GRB 221009A

We report the detection of a significant sudden ionospheric disturbance in the D-region of Earth’s ionosphere (∼60–100 km), which was associated with the massive γ-ray burst GRB 221009A that occurred on 2022 October 9. We identified the disturbance over northern Europe—a result of the increased ionization by X- and γ-ray emission from the GRB-using very low frequency radio waves as a probe of the D-region. These observations demonstrate that an extra-galactic GRB (z ∼ 0.151) can have a significant impact on the terrestrial atmosphere and illustrates that the Earth’s ionosphere can be used as a giant X- and γ-ray detector. Indeed, these observations may provide an insight into the impacts of GRBs on the ionospheres of planets in our solar system and beyond.

In‐Situ Measurements of Ion Density in the Martian Ionosphere: Underlying Structure and Variability Observed by the MAVEN‐STATIC Instrument

AbstractMeasurement of the dense cold thermal plasma in planetary ionospheres via orbiting spacecraft is challenging because ion energies are small (0–4 eV), densities can vary by four orders of magnitude, composition varies with altitude, spacecraft charging varies in time and must be measured very accurately, and instrumental effects (e.g., detector dead‐time and background) can be significant. The SupraThermal And Thermal Ion Composition instrument team has recently released a new set of data products that contain density moments of the primary ion species at Mars, including those derived at periapsis, subject to the full suite of calibration factors required. This article discusses the challenges associated with deriving these densities and provides examples of the key caveats that users of the data should be aware of. A preliminary statistical study of this new data set focuses on the structure and variability of Mars' ionosphere, demonstrating that solar zenith angle effects, the crustal magnetic fields, and electron precipitation on the nightside, drive the strongest structural features, consistent with photochemical theory and previous studies. Dayside ionospheric density profiles are highly repeatable below altitudes of 200 km, marking the region where photochemistry and collisions dominate. In the upper dayside ionosphere (altitudes >300–400 km) changes in the solar wind dynamic pressure on timescales of Mars Atmosphere and Volatile EvolutioN's orbit (hr) drive the largest (factors of 1–3) variability in ionospheric density. In contrast variability in ionospheric density peaks between 150 and 250 km altitude on the nightside (factors of 1–2), consistent with electron precipitation driving ionization in this region.

Formation of Nonlinear Stationary Structures in Ionospheric Plasma

Solar radiation, along with cosmic radiation, ionizes the Earth's atmosphere and creates a dense layer known as the ionosphere. By considering weakly relativistic degenerate plasma in the planetary ionosphere, we have studied the formation and nature of the solitary structure, electrostatic double layers (DLs), and so on. As we have considered weak relativistic degeneracy, only electrons get accelerated and are the key to stationary structures. We have implemented Sagdeev's pseudopotential method, standard Gardner equation, and accordingly identified regimes where the solitary formation and DLs may be observed. We have studied the parametric influences on solitons and DLs. Furthermore, we extended our investigation to the oscillatory Rossby solitons in the ionospheric plasma. The results may help interpret many high-energy atmospheric observations in the ionospheric plasma.

Species-dependent solar rotation effects on the Martian ionosphere

ABSTRACT Atmospheric escape is a central issue in understanding how Mars lost the habitability and it is critically controlled by the link between the atmosphere and the interplanetary space, i.e. the ionosphere. The predominant driver of the Martian ionospheric variability is the solar extreme ultraviolet radiation. To date, how the ionospheric electrons respond to this radiation has been well investigated, but the responses of various ion species are to be understood. Based on a multi-instrument data set from the Mars Atmosphere and Volatile Evolution, we report direct observations of quasi-periodical variations of photoelectrons and ions, with a common period of ≈27.2 d and almost in exact phase with the solar rotation. These diverse variations present remarkably different amplitudes. The ion density variation has a much larger amplitude than the solar flux variation and the electron density variation has a smaller amplitude. For comparison, the amplitude of the photoelectron intensity variation is nearly identical to that of the solar flux variation. The species-dependent solar rotation effects provide a good diagnostic of the upper atmospheric and ionospheric chemistry, urging reconsiderations of the solar-driven composition and variability of any planetary ionosphere.

Effect of Meteoric Ions on Ionospheric Conductance at Jupiter

AbstractIonospheric Pedersen and Hall conductances play significant roles in electromagnetic coupling between the planetary ionosphere and magnetosphere. Several observations and models have suggested the existence of meteoric ions with interplanetary origins in the lower part of Jupiter’s ionosphere; however, no models have considered the contributions of meteoric ions to ionospheric conductance. This study is designed to evaluate the contribution of meteoric ions to ionospheric conductance by developing an ionospheric model combining a meteoroid ablation model and a photochemical model. We find that the largest contribution to Pedersen and Hall conductivities occurs in the meteoric ion layer at altitudes of 350–600 km due to the large concentration of meteoric ions resulting from their long lifetimes of more than 100 Jovian days. Pedersen and Hall conductances are enhanced by factors of 3 and 10, respectively, in the middle‐ and low‐latitude and auroral regions when meteoric ions are included. The distribution of Pedersen and Hall conductances becomes axisymmetric in the middle‐ and low‐latitude regions. Enhanced axisymmetric ionospheric conductance should impact magnetospheric plasma convection. The contribution of meteoric ions to the ionospheric conductance is expected to be important only on Jupiter in our solar system because of Jupiter’s intense magnetic and gravitational fields.

Conductivities of Titan's Dusty Ionosphere

AbstractTitan's ionosphere hosts a globally distributed non‐trivial dusty ion‐ion plasma, providing an environment for studies of dusty ionospheres that is in many aspects unique in our solar system. Thanks to the Cassini mission, Titan's ionosphere also features one of the largest dusty plasma data sets from 126 flybys of the moon over 13 years, from 2004 to 2017. Recent studies have shown that negatively charged dust dramatically alters the electric properties of plasmas, in particular planetary ionospheres. Utilizing the full plasma content of the moon's ionosphere (electrons, positive ions, and negative ions/dust grains), we derive the electric conductivities and define the conductive dynamo region. Our results show that using the full plasma content increases the Pedersen conductivities at ∼1,100–1,200 km altitude by up to 35% compared to the estimates using only electron densities. The Hall conductivities are in general not affected but several cases indicate a reverse Hall effect at ∼900 km altitude (closest approach) and below. The dayside conductivities are shown to be factor ∼7–9 larger than on the nightside, owing to higher dayside plasma densities.

The Jovian Ionospheric Alfvén Resonator and Auroral Particle Acceleration

AbstractThe ionospheric Alfvén resonator (IAR) is a structure formed by the rapid decrease in the plasma density above a planetary ionosphere. This results in a corresponding increase in the Alfvén speed that can provide partial reflection of Alfvén waves. At Earth, the IAR on auroral field lines is associated with the broadband acceleration of auroral particles, sometimes termed the Alfvénic aurora. This arises since phase mixing in the IAR reduces the perpendicular wavelength of the Alfvén waves, which enhances the parallel electric field due to electron inertia. This parallel electric field fluctuates at frequencies of 0.1–20.0 Hz, comparable to the electron transit time through the acceleration region, leading to the broadband acceleration. The prevalence of such broadband acceleration at Jupiter suggests that a similar process can occur in the Jovian IAR. A numerical model of Alfvén wave propagation in the Jovian IAR has been developed to investigate these interactions, indicating that the IAR resonant frequencies are in the same range as those at Earth. This model describes the evolution of the electric and magnetic fields in the low‐altitude region close to Jupiter that is sampled during Juno's perijove passes. In particular, the model relates measurement of magnetic fields below the ion cyclotron frequency from the MAG and Waves instruments on Juno and electric fields from Waves to the associated parallel electric fields that can accelerate auroral particles.

Incorporating Inner Magnetosphere Current-driven Electron Acceleration in Numerical Simulations of Exoplanet Radio Emission

Abstract We present calculations of auroral radio emission for an Earth-like planet produced by field-aligned current (FAC) driven electron acceleration using a coupled global magnetohydrodynamic (MHD) and inner magnetosphere model, extending the capabilities of previous works which focus solely on the direct transmission of magnetic energy between the stellar wind and ionosphere. Magnetized exoplanets are expected to produce radio emission via interaction between the host star’s stellar wind and planetary magnetosphere-ionosphere system. The empirically derived Radiometric Bode’s Law (RBL) is a linear relation between the magnetic solar wind power and total emitted radio power from magnetized Solar System planets, and is often extrapolated to extreme exoplanet systems. It has been shown that the magnitudes of the FACs coupling the stellar wind to planetary ionospheres are likely to be significantly limited (often referred to as ionospheric saturation), resulting in an estimated radio power up to several orders of magnitude less than that predicted by RBL. In this paper, we demonstrate the significance of intense, sporadic FACs, driven by nightside magnetic reconnection and inner magnetosphere plasma flow, to the total radio power produced by wind–ionosphere interaction in terrestrial planets. During periods of strong stellar wind variability, the contribution from these secondary currents can be over an order of magnitude greater than the primary current systems that previous models describe. The results highlight the role of the variability of the stellar wind on the magnitude and location of the resulting emission, subsequently affecting the conditions for detectability.

Perspective on Energetic and Thermal Atmospheric Photoelectrons

Atmospheric photoelectrons are central to the production of planetary ionospheres. They are created by photoionization of the neutral planetary atmosphere by solar EUV and soft X-ray irradiance. They provide the energy to heat the thermosphere. Thermalized photoelectrons permeate magnetospheres creating polarization electric fields and plasma waves as they interact with ions to maintain charge neutrality. Energetic photoelectrons (>1 eV) have a distinctive energy spectral shape as first revealed in data from the Atmosphere Explorer satellites. Energetic photoelectrons escaping the ionosphere follow local magnetic fields illuminating the planet's magnetic topology. Current models using state-of-the-art EUV observations accurately capture their production and transport. However, in spite of 60 years of space research the electron thermalization processes occurring below 1 eV at low altitudes in planetary thermospheres are not understood quantitatively. Results from event analysis of data from the Mars Atmosphere and Volatile Evolution (MAVEN) mission are not consistent with current models of photoelectron thermalization. The lack of quantitative understanding reflects the complexity of the physics and the lack of a large data base of simultaneous neutral, ion, and electron densities and temperatures in lower planetary thermospheres.

The global potential energy surface of extrapolated to the complete basis set limit

It has been established that the proton transfer dynamics of play an important role in determining the chain reaction of phosphorus-containing compounds in the planetary ionosphere. This work presents an accurate global potential energy surface (PES) of for the first time by fitting extensive ab initio energies from the aug-cc-pV(T, Q)Z level of theory using the multi-reference configuration interaction method including Davidson correction, and extrapolating the points to the complete basis set limit. The spectral parameters of PH+(A 2Δ) and are shown to be in agreement with the data available in previous literature. Meanwhile, a detailed study of the topographical features of the global PES could be used as a reliable photolytic kinetic theory for the reaction. Furthermore, to demonstrate the validity of the new PES, we have explicitly taken into account the reaction, and assessed its feasibility in terms of reaction dynamics by calculating the integral cross-section via the time-dependent wave packet and quasi-classical trajectory approaches. The consequent results indicate that the new PES is suitable for thermochemical reactions.

Cross‐Shock Electrostatic Potentials at Mars Inferred From MAVEN Measurements

AbstractA bow shock is generated by the interaction of the solar wind with the planetary global dipole field (e.g., Earth), or with (mainly) the planetary ionosphere (e.g., Mars). The cross‐shock potential has been well studied at Earth but not yet for Mars. We infer and approximate the peak in the frame invariant de Hoffmann‐Teller shock potential profile at Mars (ϕ) with data from the Mars Atmospheric and Volatile EvolutioN (MAVEN) mission. We find that ϕ and its ratio to the solar wind ram ion energy (Eram) vary similarly against solar zenith angle (SZA, a proxy for the angle between the solar wind flow and the shock normal) and magnetic latitudes. Our results also reveal no significant dependence of the shock potential on parameters such as the angle between the upstream interplanetary magnetic field (IMF) and the shock normal and plasma beta of upstream solar wind. There is a somewhat positive correlation with the magnetosonic Mach number and the magnetic amplification ratio across the shock. We also find a solar cycle effect on the shock location, closer to the planet near the solar minimum, as expected. Lastly, similarities and differences of cross‐shock potentials at Mars and Earth are discussed. Characterizing electron energization and high‐altitude ion loss at Mars is influenced by the bow shock and thereby the work here.