Ionospheric Outflow Research Articles

Variability of Earth’s ionospheric outflow in response to the dynamic terrestrial exosphere

The most abundant neutral constituent in the exospheric region (i.e., beyond ≈ 500 km altitude) is the atomic hydrogen (H); however, its density distributions predicted by physics-based models have been challenged by satellite-based observations of its far ultraviolet emissions. This discrepancy may impact magnetospheric ions’ densities and velocities since numerous chemistry and ion-neutral coupling interactions rely sensitively on the underlying neutral hydrogen population. The Polar Wind Outflow Model a first-principled model for relevant ion species in the high-latitude ionosphere, is employed to investigate the role of neutral H on the ionospheric outflow. Specifically, variability in the outflow of ionospheric H+, He+, N+, and O+ as a response to systematic enhancement and depletion of H number densities were simulated. The altitude-dependent ion density and energy partition profiles vary with neutral H density, solar activities, and ion species. These findings suggest that the exosphere plays a crucial role in controlling the production and loss of ions through ionospheric chemistry, as well as the energy contributions by altering ion-neutral-electron collisions and the ambipolar electric field to the high-latitude ionospheric outflow. As a result, the escape rates of the ionospheric outflow are directly associated with exospheric distributions. This work potentially helps understand the dominant mechanisms of atmospheric escape, particularly during a hydrogen-rich early Earth’s and exoplanet’s atmosphere, which is known to play a significant role in understanding the evolution of Earth’s atmosphere.

Atmospheric ion escape and solar wind deposition as a function of planetary radius

ABSTRACT We explore the ability of an unmagnetized planet to retain an atmosphere as a function of its radius. We use a particle-in-cell hybrid code to simulate the global plasma interaction of unmagnetized terrestrial planets at 1 au under average solar wind conditions. We vary the radius of the planet $(R_\mathrm{ p})$ from Mars-sized ($3390 \ \mathrm{km}$) to super-Earth-sized ($9390 \ \mathrm{km}$). We inject hydrogen and oxygen ion outflows from the ionosphere and quantify how the ion escape, recirculation, solar wind deposition, and net atmospheric mass flux vary as a function of planetary radius. We find that as the radius and the corresponding ionospheric outflow rate are varied, the fraction of outflowing $\mathrm{ H^+}$ that escapes remains at $15.5\pm 1.0{{\ \rm per\, cent}}$, while the rest recirculates back towards the planet. The fraction of produced $\mathrm{ O^+}$ that escapes from a Mars-sized planet is $27\pm 1{{\ \rm per\, cent}}$, and decreases to $7\pm 1{{\ \rm per\, cent}}$ for super-Earth, suggesting that smaller planets are less able to retain heavy ions. We find, however, that larger planets have lower solar wind deposition fractions because their bow shocks are at greater distances from the surface of the planet. The ionospheric outflow rate at which mass deposition is equal to mass escape is found to be proportional to $R_\mathrm{ p}^2$. Lastly, we propose that the bulk gyration of the solar wind at the induced magnetosphere can lead to differential escape trajectories of light and heavy ions.

Differences in Ionospheric O+ and H+ Outflow During Storms With and Without Sawtooth Oscillations

AbstractPrevious simulations have suggested that O+ outflow plays a role in driving the sawtooth oscillations. This study investigates the role of O+ by identifying the differences in ionospheric outflow between sawtooth and non‐sawtooth storms using 11 years of FAST/Time of flight Energy Angle Mass Spectrograph (TEAMS) ion composition data from 1996 through 2007 during storms driven by coronal mass ejections. We find that the storm's initial phase shows larger O+ outflow during non‐sawtooth storms, and the main and recovery phases revealed differences in the location of ionospheric outflow. On the pre‐midnight sector, a larger O+ outflow was observed during the main phase of sawtooth storms, while non‐sawtooth storms exhibited stronger O+ outflow during the recovery phase. On the dayside, the peak outflow shifts significantly toward dawn during sawtooth storms. This strong dawnside sector outflow during sawtooth storms warrants consideration.

Interhemispheric Asymmetry in the Seasonal Ionospheric Outflow

AbstractA comprehensive statistical study is conducted on O+ and H+ outflows obtained from the TEAMS/FAST data during the 23rd solar cycle (1996–2007). The study investigates interhemispheric asymmetry in ionospheric outflows during local summer, winter, and equinox seasons. Data are classified into two distinct periods: the pre‐storm and geomagnetic storm phases. Numerous statistical asymmetries are identified. The findings indicate that the dayside cusp consistently demonstrates more outflow rates of O+ and H+ in the northern hemisphere than southern hemisphere during geomagnetic storms in all seasons as well as during the pre‐storm period in the summer season with the exception of H+ during summer storms. Conversely, the nightside O+ and H+ outflow rates are higher in the southern hemisphere during pre‐storm and storm periods in the summer season. Additionally, the dawnside and duskside outflow rates of O+ and H+ are predominantly stronger in the southern hemisphere.

Single‐Hemisphere Oxygen Outflow From Earth's Subauroral Zone

AbstractBesides the cusp, polar cap, and auroral oval, the nightside subauroral zone has also recently been reported as a source region of the ionospheric oxygen outflows. However, the detailed mass and energy sources of these ions remain open questions. Here, we address this issue from the perspective of the response of conjugate hemispheres. Investigation of Van Allen Probes data demonstrates a notable preference of oxygen outflows from the nightside subauroral zone from the sunlit hemisphere. This characteristic eliminates the possibility of nightside auroral precipitation playing a significant role, as it is more prominent in darkness. Instead, it highlights sunlight‐induced ionization as the mass source and enhanced plasma waves from the magnetotail as the energy source. The results presented here further support the nightside subauroral zone as an independent source of magnetospheric oxygen ions.

Beam‐Driven Electron Cyclotron Harmonic and Electron Acoustic Waves as Seen in Particle‐In‐Cell Simulations

AbstractRecent study has demonstrated that electron cyclotron harmonic (ECH) waves can be excited by a low energy electron beam. Such waves propagate at moderately oblique wave normal angles (∼70°). The potential effects of beam‐driven ECH waves on electron dynamics in Earth's plasma sheet is not known. Using two‐dimensional Darwin particle‐in‐cell simulations with initial electron distributions that represent typical plasma conditions in the plasma sheet, we explore the excitation and saturation of such beam‐driven ECH waves. Both ECH and electron acoustic waves are excited in the simulation and propagate at oblique wave normal angles. Compared with the electron acoustic waves, ECH waves grow much faster and have more intense saturation amplitudes. Cold, stationary electrons are first accelerated by ECH waves through cyclotron resonance and then accelerated in the parallel direction by both the ECH and electron acoustic waves through Landau resonance. Beam electrons, on the other hand, are decelerated in the parallel direction and scattered to larger pitch angles. The relaxation of the electron beam and the continuous heating of the cold electrons contribute to ECH wave saturation and suppress the excitation of electron acoustic waves. When the ratio of plasma to electron cyclotron frequency ωpe/ωce increases, the ECH wave amplitude increases while the electron acoustic wave amplitude decreases. Our work reveals the importance of ECH and electron acoustic waves in reshaping sub‐thermal electron distributions and improves our understanding on the potential effects of wave‐particle interactions in trapping ionospheric electron outflows and forming anisotropic (field‐aligned) electron distributions in the plasma sheet.

Reconstruction Analysis of Global Ionospheric Outflow Patterns

AbstractIonospheric outflow supplies nearly all of the heavy ions observed within the magnetosphere, as well as a significant fraction of the proton density. While much is known about upflow and outflow energization processes, the full global pattern of outflow and its evolution is only known statistically or through numerical modeling. Because of the dominant role of heavy ions in several key physical processes, this unknown nature of the full outflow pattern leads to significant uncertainty in understanding geospace dynamics, especially surrounding storm intervals. That is, global models risk not accurately reproducing the main features of intense space storms because the amount of ionospheric outflow is poorly specified and thus magnetospheric composition and mass loading could be ill‐defined. This study defines a potential mission to observe ionospheric outflow from several platforms, allowing for a reasonable and sufficient reconstruction of the full outflow pattern on an orbital cadence. An observing system simulation experiment is conducted, revealing that four well‐placed satellites are sufficient for reasonably accurate outflow reconstructions. The science scope of this mission could include the following: reveal the global structure of ionospheric outflow; relate outflow patterns to geomagnetic activity level; and determine the spatial and temporal nature of outflow composition. The science objectives could be focused to be achieved with minimal instrumentation (only a low‐energy ion spectrometer to obtain outflow reconstructions) or with a larger scientific scope by including contextual instrumentation.

Solar wind effect on the multi-fluid plasma expansion in the Venusian upper ionosphere

Inspired by the observations suggesting that at altitudes of about 1000 km the interaction between solar wind streams and Venus’ ionosphere plasma leads to ions acceleration and outflow, the influence of different solar wind physical parameters, such as densities, temperatures and initial streaming velocities, has been studied. The ionosphere plasma system consists of two positive ion populations O+, H+ and electrons along with the solar wind streaming protons and electrons. We calculated the generated oxygen and hydrogen ions flow velocities and the electric fields. In addition, we calculated rough estimates for the escaping flux of ion populations (O+, H+) from Venus’ ionosphere and compared them to observations. To a large extent, we found that the estimates match. We also discuss the relevance of ionospheric ion acceleration and outflow from Venus’ upper.

Temporal Evolution of O+ Population in the Near‐Earth Plasma Sheet During Geomagnetic Storms as Observed by the Magnetospheric Multiscale Mission

AbstractDuring geomagnetic storms, the increase in energy input into the ionosphere in the form of Poynting flux and electron precipitation leads to an enhanced ionospheric outflow that results in an increase of the O+ content in the magnetosphere. Using different missions and instrumentation, two main ionospheric sources have been identified for the oxygen ions reaching the inner magnetosphere during geomagnetic storms: the dayside cusp, and the night side auroral region. Evidence of both pathways have been presented in the literature. However, the relative contribution of each of these pathways to the enhancement of O+ observed in near‐Earth plasma sheet, as well as the dynamics involved during the development of geomagnetic storms remains an open question. Here, we present the first statistical study to date to address this question, in the form of a superposed epoch analysis of O+ and H+ moments obtained by the Magnetospheric Multiscale (MMS) mission throughout the main phase of 90 geomagnetic storms with a minimum SYM‐H of at least −50 nT. The results show a clear increase in the oxygen density in the near‐Earth plasma sheet, with values further from Earth remaining low. Temperature values for both species show an increase with the progress of the storms. These results combined suggest that, during the main phase of geomagnetic storms, most of the oxygen ions observed in the near‐Earth plasma sheet are traveling directly from the nightside auroral region.

Magnetospheric‐Ionospheric‐Atmospheric Implications From the Juno Flyby of Ganymede

AbstractJuno flew over the northern mid‐latitudes of Ganymede during orbit 34 of the Juno mission, reaching an altitude of 1,053 km (16:56:07.972 UTC) at a sub spacecraft latitude/longitude of 33.66N, 57.5W degrees on 7 June 2021. Between 16:43 and 17:02 UT, Juno pierced Ganymede's magnetosphere at a velocity relative to Ganymede of 18.57 km s−1. Juno's instrumentation provided a unique opportunity to sample the local environment of Ganymede and its magnetosphere. We present measurements of the composition of the polar ionospheric outflow and the energetic electrons that penetrate Ganymede's atmosphere and produce its aurora. When these new observations are combined with modeling, conclusions can be drawn that affect our understanding of the atmosphere of Ganymede. The measured JADE precipitating plasma electrons provide an energy flux beyond that needed to create the observed oxygen emissions measured by UVS, but the electron energy spectrum is optically thin to the sparse atmosphere and does not provide the observed oxygen ultraviolet emission unless the O2 column density is increased by over an order of magnitude compared to previous atmospheric models. More than 99% of the electron energy flux passes through the atmosphere into the ice, thereby increasing the H2 and O2 content of the atmosphere. The increased H2 and O2 production is largely responsible for increasing the oxygen column density to a level that produces within known uncertainties the OI135.6 and OI130.4 nm emissions when bombarded by the electron energy flux observed by JADE. This suggests that past modeling efforts have underestimated the density of the atmosphere by over an order of magnitude.

Kinetic Modeling of Ionospheric Outflows in Pressure Cooker Environments

AbstractPlasma escape from the high‐latitude ionosphere (ion outflow) serves as a significant source of heavy plasma to the magnetospheric plasma sheet and ring current regions. Outflows alter mass density and reconnection rates, hence global responses of the magnetosphere. A new fully kinetic and semi‐kinetic model, KAOS (Kinetic model of Auroral ion OutflowS), is constructed from first principles which traces large numbers of individual O+ ion macro‐particles along curved magnetic field lines, using a guiding‐center approximation, in order to facilitate calculation of ion distribution functions and moments. Particle forces include mirror and parallel electric field forces, a self‐consistent ambipolar electric field, and a parameterized source of ion cyclotron resonance wave heating, thought to be central to the transverse energization of ions. The model is initiated with a steady‐state ion density altitude profile and Maxwellian velocity distribution and particle trajectories are advanced via a direct simulation Monte Carlo scheme. This outlines the implementation of the kinetic outflow model, demonstrates the model's ability to achieve near‐hydrostatic equilibrium necessary for simulation spin‐up, and investigates L‐shell dependent wave heating and pressure cookers scenarios. This paper illustrates the model initialization process and numerical investigations of L‐shell dependent outflows and pressure cooker environments and serves to advance our understanding of the drivers and particle dynamics in the auroral ionosphere.

A Model of Ganymede's Magnetic and Plasma Environment During the Juno PJ34 Flyby

AbstractUsing a hybrid model (kinetic ions, fluid electrons), we provide context for plasma and magnetic field observations from Juno's PJ34 flyby of Ganymede on 07 June 2021. We consider five model configurations that successively increase the complexity of Ganymede's atmosphere and ionosphere by including additional particle species and ionization mechanisms. We examine the density and flow patterns of pick‐up ions with small , intermediate (H2O+), and large masses in Ganymede's interaction region. The results are validated by comparing the modeled magnetic field and ion densities against time series from Juno's magnetometer and plasma instruments. Our major findings are: (a) Ganymede's internal dipole dominated the magnetic field signature observed inside the moon's magnetosphere, while plasma currents shaped the field perturbations within the “wake” region detected along the Jupiter‐averted magnetopause. (b) Ganymede's pick‐up tail leaves a subtle, but clearly discernible imprint in the magnetic field downstream of the moon. (c) Heavy pick‐up ions dominate ionospheric outflow and form a tail with steep outer boundaries. (d) During the flyby, the position of Ganymede's Jupiter‐facing magnetopause varied in time due to Kelvin‐Helmholtz waves traveling along the boundary layer. As such, the location of the Jupiter‐facing magnetopause observed by Juno represents only a single snapshot of this time‐dependent process. (e) Ionospheric hydrogen ions are partially generated outside of Ganymede's magnetopause, forming a dilute corona that surrounds the moon's magnetosphere. (f) Most H2O+ ions are produced at low latitudes where field lines are closed, resulting in a very dilute pick‐up tail for this species.

Ionospheric Plasma Density Gradients Associated With Night‐Side Energetic Electron Precipitation

AbstractEnergetic electron precipitation from the equatorial magnetosphere into the atmosphere plays an important role in magnetosphere‐ionosphere coupling: precipitating electrons alter ionospheric properties, whereas ionospheric outflows modify equatorial plasma conditions affecting electromagnetic wave generation and energetic electron scattering. However, ionospheric measurements cannot be directly related to wave and energetic electron properties measured by high‐altitude, near‐equatorial spacecraft, due to large mapping uncertainties. We aim to resolve this by projecting low‐altitude measurements of energetic electron precipitation by ELFIN CubeSats onto total electron content (TEC) maps serving as a proxy for ionospheric density structures. We examine three types of precipitation on the nightside: precipitation of <200 keV electrons in the plasma sheet, bursty precipitation of <500 keV electrons by whistler‐mode waves, and relativistic (>500 keV) electron precipitation by EMIC waves. All three types of precipitation show distinct features in TEC horizontal gradients, and we discuss possible implications of these features.

The Variation of Ionospheric O+ and H+ Outflow on Storm Timescales

AbstractGeomagnetic storms are primarily driven by stream interaction regions (SIRs) and coronal mass ejections (CMEs). Since SIR and CME storms have different solar wind and magnetic field characteristics, the magnetospheric response may vary accordingly. Using FAST/TEAMS data, we investigate the variation of ionospheric O+ and H+ outflow as a function of geomagnetic storm phase during SIR and CME magnetic storms. The effects of storm size and solar EUV flux, including solar cycle and seasonal effects, on storm time ionospheric outflow, are also investigated. The results show that for both CME and SIR storms, the O+ and H+ fluences peak during the main phase, and then declines in the recovery phase. However, for CME storms, there is also significant increase during the initial phase. Because the outflow starts during the initial phase in CME storms, there is time for the O+ to reach the plasma sheet before the start of the main phase. Since plasma is convected into the ring current from the plasma sheet during the main phase, this may explain why more O+ is observed in the ring current during CME storms than during SIR storms. We also find that outflow fluence is higher for intense storms than moderate storms and is higher during solar maximum than solar minimum.

Science Case for a Global Ionospheric Outflow Mission

Science Case for a Global Ionospheric Outflow Mission

The Effect of Compression Induced Chorus Waves on 10–100 s eV Electron Precipitation

AbstractOn 7 January 2014, a solar storm erupted, which eventually compressed the Earth's magnetosphere leading to the generation of chorus waves. These waves enhanced local wave‐particle interactions and led to the precipitation of electrons from 10 s eV to 100 s keV. This paper shows observations of a low energy cutoff in the precipitation spectrum from Van Allen Probe B Helium Oxygen Proton Electron measurements. This low energy cutoff is well replicated by the predicted loss calculated from pitch angle diffusion coefficients from wave and plasma observations on Probe B. To our knowledge, this is the first time a single spacecraft has been used to demonstrate an accurate theoretical prediction for chorus wave‐induced precipitation and its low energy cutoff. The specific properties of the precipitating soft electron spectrum have implications for ionospheric activity, with the lowest energies mainly contributing to thermospheric and ionospheric upwelling, which influences satellite drag and ionospheric outflow.

A Multi‐Satellite Case Study of Low‐Energy H+ Asymmetric Field‐Aligned Distributions Observed by MMS in the Earth's Magnetosphere

AbstractField‐aligned distributions of H+ with energy less than 100 eV in the magnetosphere originate from the ionosphere through ionospheric outflow. Recently, a statistical interhemispheric asymmetry in this ion outflow was reported. In this study, we investigate a case study of asymmetric field‐aligned distributions (E < 100 eV) resulting from asymmetric outflow observed on 13 August 2019 by the Hot Plasma Composition Analyzer onboard the Magnetospheric Multiscale mission. During the reported event, the average phase space density in the parallel direction (pitch angle < 45°) is about 4.7 times higher than in the anti‐parallel direction (pitch angle > 135°), indicating greater outflow from the southern hemisphere. In this event, electron precipitation fluxes are also asymmetric between hemispheres, with more electrons traveling to the southern hemisphere than to the northern hemisphere. Global MHD simulation results (SWMF/BAT‐S‐RUS available from the CCMC) confirm that several ionospheric parameters (joule heating, convection velocities, and energy fluxes of precipitating electrons) are more enhanced in the southern hemisphere. We suggest that these interhemispheric asymmetries might cause a stronger outflow from the southern hemisphere than from the northern hemisphere, to produce the observed asymmetric field‐aligned distributions for the low‐energy H+.

Multiscale hybrid modeling of the impact response of the Earth’s magnetotail to ionospheric O+ outflow

Ionospheric outflow plays an important role in coupling the ionosphere with the solar wind-magnetosphere system. Previous multi-fluid MHD studies explored the global influence of oxygen ions of ionospheric origin (O+) on magnetospheric dynamics. A detailed exploration of the interaction of ionospheric ions with the magnetotail requires kinetic treatment for ions. We perform a self-consistent investigation of these processes with a three-dimensional space-time adaptive hybrid code, HYPERS, powered by an intelligent Event-driven Multi-Agent Planning System (EMAPS). By comparing simulations with and without outflow we conclude that oxygen ions, flowing from the ionosphere through the lobes into the tail or directly entering the inner magnetosphere, are able to significantly modify the magnetotail configuration and induce X-points and current sheet structures at locations where magnetic reconnection does not occur in a simulation without outflow, potentially very close to the Earth. This finding may have implications for interpreting substorms and magnetotail reconnection events observed for southward magnetic field simultaneously with significant contents of oxygen ions of ionospheric origin.

Cluster Observation of Ion Outflow in Middle Altitude LLBL/Cusp from Different Origins

The ionosphere is the ionized part of the upper atmosphere that is caused mainly by photoionization by solar extreme ultraviolet (EUV) emission and the atmospheric photochemistry process. The ionospheric ions escape from the ionosphere and populate the Earth’s magnetosphere. In this case study, ion outflows from two different origins were obtained by spacecraft Cluster C1 in the magnetospheric cusp region. One of the outflows was from the reflection of the dispersed solar wind particles. The other was the ionospheric outflow passing through the low latitude boundary layer of the cusp (LLBL/cusp), which was energized by downward Poynting flux. Similar to the reflected solar wind particles, outflowing ionospheric cold ions could also extend to the high-latitude region with magnetic field line convection, which mixed it up with solar wind particles. Based on the Cluster observation in the cusp region, two different origins of the outflowing particles were determined, and their unique mechanisms of formation were discussed. Results suggest that the strong electric field associated with solar wind particle precipitation may additionally accelerate the cold ionospheric ion flow in the LLBL/cusp.

Connecting Energy Input With Ionospheric Upflow and Outflow

AbstractThe connection between energy inputs and the generation of ion upflows and outflows is a topic of keen scientific interest and the subject of a number of empirical studies. Despite this interest, it remains uncertain how different ion species respond to energy input, what defines the upper and lower bounds of the ion flux, and what role solar illumination plays in regulating the relationship between energy input and ion upflows/outflows. This work simulates how ion flux scales with low and high altitude energization, and to a combination of both. Furthermore, we examine the influence of solar illumination on these relationships by considering how the scaling of ion flux with energy input changes over the solar cycle, comparing solar minimum and maximum, as well as how they change from day to night conditions. We find O+ flux tends to respond more strongly to energy inputs than H+ flux, with the O+ flux often exhibiting a lower activation energy and a greater dynamic range. The lower bound of the ion flux at 4,000 km is typically defined by the polar wind H+, although O+ upflows can dominate at low altitudes in the presence of significant frictional heating of the ion gas or soft electron precipitation. However, when significant soft electron precipitation and wave‐particle interactions are present simultaneously the lower bound of the ion flux at 4,000 km is defined by the O+. Finally, we find a difference between the steady state response of the outflow to energy input and the peak response.