Search

Spacecraft formation flying in co-planar identical orbits with oppositely directed apse lines (eccentricity vectors) allows simultaneous in-situ measurements in different altitudes and time-sequential measurements in a limited altitude range (rapid revisit); a key enabling factor in understanding many poorly understood phenomena in our atmosphere and its interaction with the magnetosphere. The dynamics and control problem for flying formations in mirrored orbits are investigated using the Mechanisms of Energetic Mass Ejection - Explorer (MEME-X) mission concept as a case study. It is determined that the proposed two-craft and four-craft configurations require routine maintenance to correct their relative motion. The formation maintenance costs are quantified for both configurations and determined to be within the capabilities of small satellites. A high-fidelity end-to-end simulation is developed to model the entire two-craft mission concept, from orbit insertion to decommissioning, providing an approximate baseline for the flight dynamics costs of similar missions using multiple spacecraft flying in mirrored orbits.

The solar wind has been seen as the major source of hot magnetospheric plasma since the early 1960’s. More recent theoretical and observational studies have shown that the cold (few eV) polar wind and warmer polar cusp plasma that flow continuously upward from the ionosphere can be a very significant source of ions in the magnetosphere and can become accelerated to the energies characteristic of the plasma sheet, ring current, and warm plasma cloak. Previous studies have also shown the presence of solar wind ions in these magnetospheric regions. These studies are based principally on proxy measurements of the ratios of He++/H+ and the high charge states of O+/H+. The resultant admixture of ionospheric ions and solar wind ions that results has been difficult to quantify, since the dominant H+ ions originating in the ionosphere and solar wind are indistinguishable. The ionospheric ions are already inside the magnetosphere and are filling it from the inside out with direct access from the ionosphere to the center of the magnetotail. The solar wind ions on the other hand must gain access through the outer boundaries of the magnetosphere, filling the magnetosphere from the outside in. These solar wind particles must then diffuse or drift from the flanks of the magnetosphere to the near-midnight reconnection region of the tail which takes more time to reach (hours) than the continuously large outflowing ionospheric polar wind (10’s of min). In this paper we examine the magnetospheric filling using the trajectories of the different ion sources to unravel the intermixing process rather than trying to interpret only the proxy ratios. We compare the timing of the access of the ionospheric and solar wind sources and we use new merged ionosphere-magnetosphere multi-fluid MHD modeling to separate and compare the ionospheric and solar wind H+ source strengths. The rapid access of the initially cold polar wind and warm polar cusp ions flowing down-tail in the lobes into the mid-plane of the magnetotail, suggests that, coupled with a southward turning of the IMF Bz, these ions can play a key triggering role in the onset of substorms and subsequent large storms.

The roles of heavy ions have long been an important subject in the magnetospheric physics since the first discovery of O+ ions in the magnetosphere as it hinted to the connection between the ionosp...

Ionospheric ions (mainly H+, He+, and O+) escape from the ionosphere and populate the Earth's magnetosphere. Their thermal energies are usually low when they first escape the ionosphere, typically a few electron volt to tens of electron volt, but they are energized in their journey through the magnetosphere. The ionospheric population is variable, and it makes significant contributions to the magnetospheric mass density in key regions where magnetic reconnection is at work. Solar wind—magnetosphere coupling occurs primarily via magnetic reconnection, a key plasma process that enables transfer of mass and energy into the near-Earth space environment. Reconnection leads to the triggering of magnetospheric storms, auroras, energetic particle precipitation and a host of other magnetospheric phenomena. Several works in the last decades have attempted to statistically quantify the amount of ionospheric plasma supplied to the magnetosphere, including the two key regions where magnetic reconnection occurs: the dayside magnetopause and the magnetotail. Recent in situ observations by the Magnetospheric Multiscale spacecraft and associated modeling have advanced our current understanding of how ionospheric ions alter the magnetic reconnection process, including its onset and efficiency. This article compiles the current understanding of the ionospheric plasma supply to the magnetosphere. It reviews both the quantification of these sources and their effects on the process of magnetic reconnection. It also provides a global description of how the ionospheric ion contribution modifies the way the solar wind couples to the Earth's magnetosphere and how these ions modify the global dynamics of the near-Earth space environment. Plain Language Summary Above the neutral atmosphere, space is filled with charged particles, which are tied to the Earth's magnetic field. The particles come from two sources, the solar wind and the Earth's upper atmosphere. Most of the solar wind particles are deflected by the Earth´s magnetic field, but some can penetrate into near-Earth space. The ionized layer of the upper atmosphere is continuously ejecting particles into space, which have low energies and are difficult to measure. We investigate the relative importance of the two charged particle sources for the dynamics of plasma processes in near-Earth space. In particular, we consider the effects of these sources in magnetic reconnection. Magnetic reconnection allows initially separated plasma regions to become magnetically connected and mix, and converts magnetic energy to kinetic energy of charged particles. Magnetic reconnection is the main driver of geomagnetic activity in the near-Earth space, and is responsible for the release of energy that drives a variety of space weather effects. We highlight the fact that plasma from the ionized upper atmosphere contributes a significant part of the density in the key regions where magnetic reconnection is at work, and that this contribution is larger when the geomagnetic activity is high.

The escape of heavy ions from the Earth atmosphere is consequences of energization and transport mechanisms, including photoionization, electron precipitation, ion-electron-neutral chemistry and co...

We test the hypothesis that their dominant driver of a planetary ambipolar electric field is the ionospheric electron pressure gradient (∇P(e)). The ionospheres of Venus and Mars are mapped using Langmuir probe measurements from NASA’s Pioneer Venus Orbiter (PVO) and Mars Atmosphere and Volatile Evolution (MAVEN) missions. We then determine the component of the ionospheric potential drop that can be explained by the electron pressure gradient drop along a simple draped field line. At Mars, this calculation is consistent with the mean potential drops measured statistically by MAVEN. However, at Venus, contrary to our current understanding, the thermal electron pressure gradient alone cannot explain Venus’ strong ambipolar field. These results strongly motivate a return to Venus with a comprehensive plasmas and fields package, similar to that on MAVEN, to investigate the physics of atmospheric escape at Earth’s closest analog.

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths’ ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

On 27 June 1996, the NASA Galileo spacecraft made humanity’s first flyby of Jupiter’s largest moon, Ganymede, discovering that it is the only moon known to possess an internally generated magnetic field. Resurrecting the original Galileo Plasma Subsystem (PLS) data analysis software, we processed the raw PLS data from G01 and for the first time present the properties of plasmas encountered. Entry into the magnetosphere of Ganymede occurred near the confluence of the magnetopause and plasma sheet. Reconnection-driven plasma flows were observed (consistent with an Earth-like Dungey cycle), which may be a result of reconnection in the plasma sheet, magnetopause, or might be Ganymede’s equivalent of a Low-Latitude Boundary Layer. Dropouts in plasma density combined with velocity perturbations afterward suggest that Galileo briefly crossed the cusps into closed magnetic field lines. Galileo then crossed the cusps, where field-aligned precipitating ions were observed flowing down into the surface, at a location consistent with observations by the Hubble Space Telescope. The density of plasma outflowing from Ganymede jumped an order of magnitude around closest approach over the north polar cap. The abrupt increase may be a result of crossing the cusp or may represent an altitude-dependent boundary such as an ionopause. More diffuse, warmer field-aligned outflows were observed in the lobes. Fluxes of particles near the moon on the nightside were significantly lower than on the dayside, possibly resulting from a diurnal cycle of the ionosphere and/or neutral atmosphere.

We present a new expansion of the Polar Wind Outflow Model (PWOM) to include kinetic ions using the Particle-in-Cell (PIC) approach with Monte Carlo collisions. This implementation uses the original hydrodynamic solution at low altitudes for efficiency, and couples to the kinetic solution at higher altitudes to account for kinetic effects important for ionospheric outflow. The modeling approach also includes wave-particle interactions, suprathermal electrons, and an hybrid parallel computing approach combining shared and distributed memory paralellization. The resulting model is thus a comprehensive, global, model of ionospheric outflow that can be run efficiently on large supercomputing clusters. We demonstrate the model’s capability to study a range of problems starting with the comparison of kinetic and hydrodynamic solutions along a single field line in the sunlit polar cap, and progressing to the altitude evolution of the ion conic distribution in the cusp region. The interplay between convection and the cusp on the global outflow solution is also examined. Finally, we demonstrate the impact of these new model features on the magnetosphere by presenting the first 2-way coupled ionospheric outflow-magnetosphere calculation including kinetic ion effects.

Recent studies suggest that the hydrogen (H) density in the exosphere and geocorona might differ from previously assumed values by factors as large as 2. We use the SAMI3 (Sami3 is Also a Model of the Ionosphere) and Comprehensive Inner Magnetosphere-Ionosphere models to evaluate scenarios where the hydrogen density is reduced or enhanced, by a factor of 2, relative to values given by commonly used empirical models. We show that the rate of plasmasphere refilling following a geomagnetic storm varies nearly linearly with the hydrogen density. We also show that the ring current associated with a geomagnetic storm decays more rapidly when H is increased. With respect to these two space weather effects, increased exosphere hydrogen density is associated with reduced threats to space assets during and following a geomagnetic storm.

We examine the mechanisms responsible for the dropout of energetic electron flux during 31 May to 1 June 2013 using Van Allen Probe (Radiation Belt Storm Probes (RBSP)) electron flux data and simulations with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. During the storm main phase, L-shells at RBSP locations are greater than ~8, which are connected to open drift shells. Consequently, diminished electron fluxes were observed over a wide range of energies. The combination of drift shell splitting, magnetopause shadowing, and drift loss all results in butterfly electron pitch angle distributions (PADs) at the nightside. During storm sudden commencement, RBSP observations display electron butterfly PADs over a wide range of energies. However, it is difficult to determine whether there are butterfly PADs during the storm main phase since the maximum observable equatorial pitch angle from RBSP is not larger than ~40° during this period. To investigate the causes of the dropout, the CIMI model is used as a global 4-D kinetic inner magnetosphere model. The CIMI model reproduces the dropout with very similar timing and flux levels and PADs along the RBSP trajectory for 593 keV. Furthermore, the CIMI simulation shows butterfly PADs for 593 keV during the storm main phase. Based on comparison of observations and simulations, we suggest that the dropout during this event mainly results from magnetopause shadowing.

This study presents a newly coupled model capable of treating the superthermal electron population in the global polar wind solution. The model combines the hydrodynamic Polar Wind Outflow Model (PWOM) with the kinetic SuperThermal Electron Transport (STET) code. The resulting PWOM-STET coupled model is described and then used to investigate the role of photoelectrons in the polar wind. We present polar wind results along single stationary field lines under dayside and nightside conditions, as well as the global solution reconstructed from nearly 1000 moving field lines. The model results show significant day-night asymmetries in the polar wind solution owing to the higher ionization and photoelectron fluxes on the dayside compared to the nightside. Field line motion is found to modify this dependence and create global structure by transporting field lines through different conditions of illumination and through the localized effects of Joule heating.

The growth and saturation of electromagnetic ion cyclotron (EMIC) waves is essential to the magnetospheric dynamics. Determining and isolating the effects of multiple ion parameters such as temperatures, anisotropies, and relative abundances is important for quantifying these processes in the magnetospheric plasma. In order to study these process, we utilize a 2.5-D hybrid model (where ions are modeled with the particle-in-cell (PIC) method, and electrons are modeled as background neutralizing fluid) to study the nonlinear electromagnetic wave-particle interactions of hot H+, cold H+, cold He+, and cold or hot O+ ions for a broad range of typical magnetospheric parameters. The excitation of EMIC waves is driven by the temperature anisotropy of hot H+ in our model. As a result, we quantify the parametric dependence of the linear growth, the nonlinear saturation level of perpendicular magnetic fluctuations, and the temporal evolution of the ion temperature anisotropies. We establish the relation between key plasma parameters and the saturated EMIC wave power, using either power law fits or a nonlinear regression method. We construct the dispersion relation of the waves using the results of the model and investigate the energy content in the various branches of the dispersion (k∥−ω space), showing that the different modes can generate wave power in different regions of k space. We find that large O+ concentration reduces the growth and saturated amplitude of the waves; but the waves are less sensitive to the temperature of the O+ in the temperature range relevant to the magnetosphere.

Changes in the heavy ion composition in the terrestrial ionosphere and magnetosphere can have significant impact on particle dynamics in the Earth’s magnetosphere-ionosphere system. Most instrument...

We investigate the space environment conditions near the Earth-size planet TOI~700~d using a set of numerical models for the stellar corona and wind, the planetary magnetosphere, and the planetary ionosphere. We drive our simulations using a scaled-down stellar input and a scaled-up solar input in order to obtain two independent solutions. We find that for the particular parameters used in our study, the stellar wind conditions near the planet are not very extreme -- slightly stronger than that near the Earth in terms of the stellar wind ram pressure and the intensity of the interplanetary magnetic field. Thus, the space environment near TOI700-d may not be extremely harmful to the planetary atmosphere, assuming the planet resembles the Earth. Nevertheless, we stress that the stellar input parameters and the actual planetary parameters are unconstrained, and different parameters may result in a much greater effect on the atmosphere of TOI700-d. Finally, we compare our results to solar wind measurements in the solar system and stress that modest stellar wind conditions may not guarantee atmospheric retention of exoplanets., Comment: accepted to ApJ

Electromagnetic ion cyclotron (EMIC) waves have long been considered to be a significant loss mechanism for relativistic electrons. This has most often been attributed to resonant interactions with the highest amplitude waves. But recent observations have suggested that the dominant energy of electrons precipitated to the atmosphere may often be relatively low, less than 1 MeV, whereas the minimum resonant energy of the highest amplitude waves is often greater than 2 MeV. Here we use relativistic electron test particle simulations in the wavefields of a hybrid code simulation of EMIC waves in dipole geometry in order to show that significant pitch angle scattering can occur due to interaction with low-amplitude short-wavelength EMIC waves. In the case we examined, these waves are in the H band (at frequencies above the He+ gyrofrequency), even though the highest amplitude waves were in the He band frequency range (below the He+ gyrofrequency). We also present wave power distributions for 29 EMIC simulations in straight magnetic field line geometry that show that the high wave number portion of the spectrum is in every case mostly due to the H band waves. Though He band waves are often associated with relativistic electron precipitation, it is possible that the He band waves do not directly scatter the sub-megaelectron volts (sub-MeV) electrons, but that the presence of He band waves is associated with high plasma density which lowers the minimum resonant energy so that these electrons can more easily resonate with the H band waves.

International audience; When the interplanetary magnetic field is northward for a period of time, O + from the high-latitude ionosphere escapes along reconnected magnetic field lines into the dayside magnetopause boundary layer. Dual-lobe reconnection closes these field lines, which traps O + and mass loads the boundary layer. This O + is an additional source of magnetospheric plasma that interacts with magnetosheath plasma through magnetic reconnection. This mass loading and interaction is illustrated through analysis of a magnetopause crossing by the Magnetospheric Multiscale spacecraft. While in the O +-rich boundary layer, the interplanetary magnetic field turns southward. As the Magnetospheric Multiscale spacecraft cross the high-shear magnetopause, reconnection signatures are observed. While the reconnection rate is likely reduced by the mass loading, reconnection is not suppressed at the magnetopause. The high-latitude dayside ionosphere is therefore a source of magnetospheric ions that contributes often to transient reduction in the reconnection rate at the dayside magnetopause.

International audience; Parallel electric fields and their associated electric potential structures play a crucial role in ionospheric-magnetospheric interactions at any planet. Although there is abundant evidence that parallel electric fields play key roles in Martian ionospheric outflow and auroral electron acceleration, the fields themselves are challenging to directly measure due to their relatively weak nature. Using measurements by the Solar Wind Electron Analyzer instrument aboard the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) Mars Scout, we present the discovery and measurement of a substantial (ΦMars=7.7 ± 0.6 V) parallel electric potential drop on closed magnetic field lines spanning the terminator from day to night above the great impact basin of Utopia Planitia, a region largely free of crustal magnetic fields. A survey of the previous 26 orbits passing over a range of longitudes revealed similar signatures on seven orbits, with a mean potential drop (ΦMars) of 10.9 ± 0.8 V, suggestive that although trans-terminator electric fields of comparable strength are not ubiquitous, they may be common, at least at these northerly latitudes.

Both electron cyclotron harmonic (ECH) waves and whistler mode chorus waves resonate with electrons of the Earths plasma sheet in the energy range from tens of eV to several keV and produce the electron diffuse aurora at ionospheric altitudes. Interaction of these superthermal electrons with the neutral atmosphere leads to the production of secondary electrons (E500600 eV) and, as a result, leads to the activation of lower energy superthermal electron spectra that can escape back to the magnetosphere and contribute to the thermal electron energy deposition processes in the magnetospheric plasma. The ECH and whistler mode chorus waves, however, can also interact with the secondary electrons that are coming from both of the magnetically conjugated ionospheres after they have been produced by initially precipitated high-energy electrons that came from the plasma sheet. After their degradation and subsequent reflection in magnetically conjugate atmospheric regions, both the secondary electrons and the precipitating electrons with high (E600 eV) initial energies will travel back through the loss cone, become trapped in the magnetosphere, and redistribute the energy content of the magnetosphere-ionosphere system. Thus, scattering of the secondary electrons by ECH and whistler mode chorus waves leads to an increase of the fraction of superthermal electron energy deposited into the core magnetospheric plasma.

Understanding what processes govern atmospheric escape and the loss of planetary water is of paramount importance for understanding how life in the universe can exist. One mechanism thought to be important at all planets is an “ambipolar” electric field that helps ions overcome gravity. We report the discovery and first quantitative extraterrestrial measurements of such a field at the planet Venus. Unexpectedly, despite comparable gravity, we show the field to be five times stronger than in Earth's similar ionosphere. Contrary to our understanding, Venus would still lose heavy ions (including oxygen and all water-group species) to space, even if there were no stripping by the solar wind. We therefore find that it is possible for planets to lose heavy ions to space entirely through electric forces in their ionospheres and such an “electric wind” must be considered when studying the evolution and potential habitability of any planet in any star system.

An energetic process is needed to convert N2 into compounds essential for life. Simulations show that interactions between powerful solar flares and Earth’s magnetic field could have facilitated nitrogen fixation in the early atmosphere.

A flux dropout is a sudden and sizable decrease in the energetic electron population of the outer radiation belt on the time scale of a few hours. We simulated a flux dropout of highly relativistic 2.5 MeV electrons using the Radiation Belt Environment model, incorporating the pitch angle diffusion coefficients caused by electromagnetic ion cyclotron (EMIC) waves for the geomagnetic storm events of 23-26 October 2002. This simulation showed a remarkable decrease in the 2.5 MeV electron flux during main phase of the storm, compared to those without EMIC waves. This decrease was independent of magnetopause shadowing or drift loss to the magnetopause. We suggest that the flux decrease was likely to be primarily due to pitch angle scattering to the loss cone by EMIC waves. Furthermore, the 2.5 MeV electron flux calculated with EMIC waves correspond very well with that observed from Solar Anomalous and Magnetospheric Particle EXplorer spacecraft. EMIC wave scattering is therefore likely one of the key mechanisms to understand flux dropouts. We modeled EMIC wave intensities by the Kp index. However, the calculated dropout is a several hours earlier than the observed one. We propose that Kp is not the best parameter to predict EMIC waves.

The Earth inner magnetosphere is a donut-shape large reservoir of magnetized plasma, energetic particles (keV – MeV), and various plasma waves. Plasma waves play important roles in the Earth’s magnetosphere. Plasma waves with very low frequency (VLF) of $\sim$ kHz and ultra low frequency (ULF) of $\sim$ Hz, can interact with energetic electrons and ions trapped in the magnetic field in the magnetosphere at approximately 3 – 10 Earth radii. Thus, these waves can change pitch-angle (angle to the magnetic field line in momentum space) and result in charged particle precipitation into the ionosphere and upper atmosphere. This can in turn cause charged-particle-impact ionization and change ionospheric densities and conductivities [e.g., 1]. Here we comprehensively study charged particle precipitation from the inner magnetosphere into ionosphere and their effects. We simulate ion and electron dynamics in the inner magnetosphere during moderate and strong storms, using the CIMI/BATS-R-US model [2, 3], which solve the bounce-averaged Boltzmann’ transport and magnetohydrodynamics equations including quasi-linear wave-particle interactions. We also simulate ion and electron precipitation into the mid-latitude ionosphere and estimate their impact on the ionosphere calculating ionospheric chemistry and conductivity, which affects magnetosphere-ionsphere coupled current system and Joule heating in the ionsphere.

During the 14 November 2012 geomagnetic storm, the Van Allen Probes spacecraft observed a number of sharp decreases (“dropouts”) in particle fluxes for ions and electrons of different energies. In this paper, we investigate the global magnetosphere dynamics and magnetosphere–ionosphere (M–I) coupling during the dropout events using multipoint measurements by Van Allen Probes, TWINS, and AMPERE together with the output of the two-way coupled global BATS-R-US–CRCM model. We find different behavior for two pairs of dropouts. For one pair, the same pattern was repeated: (1) weak nightside Region 1 and 2 Birkeland currents before and during the dropout; (2) intensification of Region 2 currents after the dropout; and (3) a particle injection detected by TWINS after the dropout. The model predicted similar behavior of Birkeland currents. TWINS low-altitude emissions demonstrated high variability during these intervals, indicating high geomagnetic activity in the near-Earth tail region. For the second pair of dropouts, the structure of both Birkeland currents and ENA emissions was relatively stable. The model also showed quasi-stationary behavior of Birkeland currents and simulated ENA emissions with gradual ring current buildup. We confirm that the first pair of dropouts was caused by large-scale motions of the OCB (open–closed boundary) during substorm activity. We show the new result that this OCB motion was associated with global changes in Birkeland (M–I coupling) currents and strong modulation of low-altitude ion precipitation. The second pair of dropouts is the result of smaller OCB disturbances not related to magnetospheric substorms. The local observations of the first pair of dropouts result from a global magnetospheric reconfiguration, which is manifested by ion injections and enhanced ion precipitation detected by TWINS and changes in the structure of Birkeland currents detected by AMPERE. This study demonstrates that multipoint measurements along with the global model results enable the reconstruction of a more complete system-level picture of the dropout events and provides insight into M–I coupling aspects that have not previously been investigated. Keywords. Magnetospheric physics (magnetosphere–ionosphere interactions; magnetospheric configuration and dynamics); space plasma physics (numerical simulation studies)

The warm plasma cloak is a source of magnetospheric plasma that contain significant O+. When the O+ density in the magnetosphere near the magnetopause is >0.2 cm‐3 and the H+ density is 20% due to mass‐loading only about 2% to 4% of the time. However, during geomagnetic storms, O+ dominates the mass density of the warm plasma cloak and these mass densities are very high. Therefore, a separate study is conducted to determine the effect of the warm plasma cloak on magnetopause reconnection during geomagnetically disturbed times. This study shows that the warm plasma cloak reduces the reconnection rate significantly about 25% of the time during disturbed conditions. publishedVersion

We have recently developed a new modeling capability to embed the implicit particle-in-cell (PIC) model iPIC3D into the Block-Adaptive-Tree-Solarwind-Roe-Upwind-Scheme magnetohydrodynamic (MHD) model. The MHD with embedded PIC domains (MHD-EPIC) algorithm is a two-way coupled kinetic-fluid model. As one of the very first applications of the MHD-EPIC algorithm, we simulate the interaction between Jupiter's magnetospheric plasma and Ganymede's magnetosphere. We compare the MHD-EPIC simulations with pure Hall MHD simulations and compare both model results with Galileo observations to assess the importance of kinetic effects in controlling the configuration and dynamics of Ganymede's magnetosphere. We find that the Hall MHD and MHD-EPIC solutions are qualitatively similar, but there are significant quantitative differences. In particular, the density and pressure inside the magnetosphere show different distributions. For our baseline grid resolution the PIC solution is more dynamic than the Hall MHD simulation and it compares significantly better with the Galileo magnetic measurements than the Hall MHD solution. The power spectra of the observed and simulated magnetic field fluctuations agree extremely well for the MHD-EPIC model. The MHD-EPIC simulation also produced a few flux transfer events (FTEs) that have magnetic signatures very similar to an observed event. The simulation shows that the FTEs often exhibit complex 3-D structures with their orientations changing substantially between the equatorial plane and the Galileo trajectory, which explains the magnetic signatures observed during the magnetopause crossings. The computational cost of the MHD-EPIC simulation was only about 4 times more than that of the Hall MHD simulation.

In this work, we demonstrate how to track magnetic separators in three-dimensional simulated magnetic fields with or without magnetic nulls, apply these techniques to enhance our understanding of reconnection at the magnetopause. We present three methods for locating magnetic separators and apply them to 3-D resistive MHD simulations of the Earth's magnetosphere using the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme code. The techniques for finding separators and determining the reconnection rate are insensitive to interplanetary magnetic field (IMF) clock angle and can in principle be applied to any magnetospheric model. Moreover, the techniques have a number of advantages over prior separator finding techniques applied to the magnetosphere. The present work examines cases of high and low resistivity for two clock angles. We go beyond previous work examine the separator during Flux Transfer Events (FTEs). Our analysis of reconnection on the magnetopause yields a number of interesting conclusions: Reconnection occurs all along the separator even during predominately northward IMF cases. Multiple separators form in low-resistivity conditions, and in the region of an FTE the separator splits into distinct branches. Moreover, the local contribution to the reconnection rate, as determined by the local parallel electric field, drops in the vicinity of the FTE with respect to the value when there are none.

The precipitation of high-energy magnetospheric electrons (E greater than 500-600 electronvolts) in the diffuse aurora contributes significant energy flux into Earth's ionosphere. In the diffuse aurora, precipitating electrons initially injected from the plasmasheet via wave-particle interaction processes degrade in the atmosphere toward lower energies and produce secondary electrons via impact ionization of the neutral atmosphere. These initially precipitating electrons of magnetospheric origin can be additionally reflected back into the magnetosphere by the two magnetically conjugated atmospheres, leading to a series of multiple reflections that can greatly influence the initially precipitating flux at the upper ionospheric boundary (700-800 kilometers) and the resultant population of secondary electrons and electrons cascading toward lower energies. We present the solution of the Boltzmann.Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E is less than or equal to 600 electronvolts) and their energy interplay in the magnetosphere and two conjugated ionospheres. This solution takes into account the role of multiple atmospheric reflections of the precipitated electrons that were initially moved into the loss cone via wave.particle interaction processes in Earth's plasmasheet.

An important mechanism in the generation of polar wind outflow is the ambipolar electric potential which assists ions in overcoming gravity and is a key mechanism for Terrestrial ionospheric escape. At Mars, open field lines are not confined to the poles, and outflow of ionospheric electrons is observed far into the tail. It has thus been hypothesized that a similar electric potential may be present at Mars, contributing to global ionospheric loss. However, no direct measurements of this potential have been made. In this pilot study, we examine photoelectron spectra measured by the Solar Wind Electron Analyzer instrument on the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) Mars Scout to put an initial upper bound on the total potential drop in the ionosphere of Mars of Φ♂≾⊥2V , with the possibility of a further ≾4.5 V potential drop above this in the magnetotail. If the total potential drop was close to the upper limit, then strong outflows of major ionospheric species (H+, O+, and O2+) would be expected. However, if most of the potential drop is confined below the spacecraft, as expected by current theory, then such a potential would not be sufficient on its own to accelerate O2+ to escape velocities, but would be sufficient for lighter ions. However, any potential would contribute to atmospheric loss through the enhancement of Jeans escape.

The twin Van Allen Probes spacecraft witnessed a series of lobe encounters between 0200 and 0515 UT on 14 November 2012. Although lobe entry had been observed previously by other spacecraft, the two Van Allen Probe spacecraft allow us to observe the motion of the boundary for the first time. Moreover, this event is unique in that it consists of a series of six quasi-periodic lobe entries. The events occurred on the dawn flank between 4 and 6.6 local time and at altitudes between 5.6 and 6.2 RE. During the events Dst dropped to less than −100nT with the IMF being strongly southward (Bz = −15nT) and eastward (By = 20 nT). Observations by LANL-GEO spacecraft at geosynchronous orbit also show lobe encounters on the dawn and dusk flanks. The two spacecraft configuration provides strong evidence that these periodic entries into the lobe are the result of local expansions of the OCB propagating from the tail and passing over the Van Allen Probes. Examination of pitch angle binned data from the HOPE instrument shows spatially large, accelerated ion structures occurring near simultaneously at both spacecraft, with the presence of oxygen indicating that they have an ionospheric source. The outflows are dispersed in energy and are detected when the spacecraft are on both open and closed field lines. These events provide a chance to examine the global magnetic field topology in detail, as well as smaller-scale spatial and temporal characteristics of the OCB, allowing us to constrain the position of the open/closed field line boundary and compare it to a global MHD model using a novel method. This technique shows that the model can reproduce a periodic approach and retreat of the OCB from the spacecraft but can overestimate its distance by as much as 3 RE. The model appears to simulate the dynamic processes that cause the spacecraft to encounter the lobe but incorrectly maps the overall topology of the magnetosphere during these extreme conditions.

It is now well established that the ionosphere, because it acts as a significant source of plasma, plays a critical role in ring current dynamics. However, because the ring current deposits energy into the ionosphere, the inverse may also be true: the ring current can play a critical role in the dynamics of ionospheric outflow. This study uses a set of coupled, first-principles-based numerical models to test the dependence of ionospheric outflow on ring current-driven region 2 field-aligned currents (FACs). A moderate magnetospheric storm event is modeled with the Space Weather Modeling Framework using a global MHD code (Block Adaptive Tree Solar wind Roe-type Upwind Scheme, BATS-R-US), a polar wind model (Polar Wind Outflow Model), and a bounce-averaged kinetic ring current model (ring current atmosphere interaction model with self-consistent magnetic field, RAM-SCB). Initially, each code is two-way coupled to all others except for RAM-SCB, which receives inputs from the other models but is not allowed to feed back pressure into the MHD model. The simulation is repeated with pressure coupling activated, which drives strong pressure gradients and region 2 FACs in BATS-R-US. It is found that the region 2 FACs increase heavy ion outflow by up to 6 times overmore » the non-coupled results. The additional outflow further energizes the ring current, establishing an ionosphere-magnetosphere mass feedback loop. This study further demonstrates that ionospheric outflow is not merely a plasma source for the magnetosphere but an integral part in the nonlinear ionosphere-magnetosphere-ring current system.« less

We present a method to quantify the upper-limit of the energy transmitted from the intense stellar wind to the upper atmospheres of three of the Trappist-1 planets (e, f, and g). We use a formalism that treats the system as two electromagnetic regions, where the efficiency of the energy transmission between one region (the stellar wind at the planetary orbits) to the other (the planetary ionospheres) depends on the relation between the conductances and impedances of the two regions. Since the energy flux of the stellar wind is very high at these planetary orbits, we find that for the case of high transmission efficiency (when the conductances and impedances are close in magnitude), the energy dissipation in the upper planetary atmospheres is also very large. On average, the Ohmic energy can reach $0.5-1~W/m^2$, about 1\% of the stellar irradiance and 5-15 times the EUV irradiance. Here, using constant values for the ionospheric conductance, we demonstrate that the stellar wind energy could potentially drive large atmospheric heating in terrestrial planets, as well as in hot jupiters. More detailed calculations are needed to assess the ionospheric conductance and to determine more accurately the amount of heating the stellar wind can drive in close-orbit planets., 6 pages, 1 table, 2 figures, accepted to ApJ Letters

Many space weather effects intensify during geomagnetic storms potentially impacting Earth-based and spaceborne technological systems. In this chapter, we concentrate on geomagnetic storms’ effects in the magnetosphere and provide descriptions of currently available 3D global modeling tools designed to reproduce the near-Earth space environment during geomagnetic storms. The outputs of these models could be used as inputs for ionosphere/thermosphere models (2D maps of ionospheric electric fields and precipitation), geomagnetically induced current models (geomagnetic variations on the Earth's surface), or environment models for spacecraft charging (ring current or radiation belt electron fluxes at spacecraft locations). We briefly review previously published results on modeling of large/extreme geomagnetic storms. We also provide an example simulation of the June 22–23, 2015 magnetic storm with minimum Dst = −204 nT, showing the distribution of pressure and current density within the magnetosphere, provide comparisons with spacecraft measurements, and present a data-model comparison of modeled magnetograms with ground-based observations. We analyze sources of observed dB/dt variations, try to identify the responsible magnetospheric processes, and emphasize the role of the ring current at low latitudes. We also discuss challenges for the future improvement of two-way coupled models.

From 2 to 5 UT on 14 November 2012, the Van Allen Probes observed repeated particle flux dropouts during the main phase of a geomagnetic storm as the satellites traversed the post-midnight to dawnside inner magnetosphere. Each flux dropout corresponded to an abrupt change in the magnetic topology, i.e., from a more dipolar configuration to a configuration with magnetic field lines stretched in the dawn-dusk direction. Geosynchronous GOES spacecraft located in the dusk and near-midnight sectors and the LANL constellation with wide local time coverage also observed repeated flux dropouts and stretched field lines with similar occurrence patterns to those of the Van Allen Probe events. THEMIS recorded multiple transient abrupt expansions of the evening-side magnetopause ∼20–30 min prior to the sequential Van Allen Probes observations. Ground-based magnetograms and all sky images demonstrate repeatable features in conjunction with the dropouts. We combine the various in situ and ground-based measurements to define and understand the global spatiotemporal features associated with the dropouts observed by the Van Allen Probes. We discuss various proposed hypotheses for the mechanism that plausibly caused this storm-time dropout event as well as formulate a new hypothesis that explains the combined in situ and ground-based observations: the earthward motion of magnetic flux ropes containing lobe plasmas that form along an extended magnetotail reconnection line in the near-Earth plasma sheet.

There are two main theories for the origin of the diffuse auroral electron precipitation: first, pitch angle scattering by electrostatic electron cyclotron harmonic (ECH) waves, and second, by whistler mode waves. Precipitating electrons initially injected from the plasma sheet to the loss cone via wave-particle interaction processes degrade in the atmosphere toward lower energies and produce secondary electrons via impact ionization of the neutral atmosphere. These secondary electrons can escape back to the magnetosphere, become trapped on closed magnetic field lines, and deposit their energy back to the inner magnetosphere. ECH and whistler mode waves can also move electrons in the opposite direction, from the loss cone into the trap zone, if the source of such electrons exists in conjugate ionospheres located at the same field lines as the trapped magnetospheric electron population. Such a situation exists in the simulation scenario of superthermal electron energy interplay in the region of diffuse aurora presented and discussed by Khazanov et al. (2014) and will be quantified in this paper by taking into account the interaction of secondary electrons with ECH waves.

Relativistic electrons in the outer radiation belt are highly dynamic and respond to interplanetary solar wind structures interacting with the Earth's magnetic field. A known mechanism dictating electron dynamics is the drift-resonant interaction with ultra-low frequency (ULF) waves. The present work simulates the ring current and radiation belt electron populations in the bounce-averaged, kinetic Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model coupled with the Block Adaptive Tree Solar Wind Roe-type Upwind Scheme (BATS-R-US) global magnetospheric magnetohydrodynamic (MHD) code using an idealized ULF wave solar wind density driver. ULF waves generated with 10 minute periods (at 1.67 mHz frequencies) in the MHD model are characterized and the corresponding energization of electrons and radial transport of electron phase space density is presented. The drift-resonant electron energy is determined in the simulation and is consistent with the electron resonance conditions in dipolar magnetic fields. The present results will be an important component of understanding inner magnetospheric dynamics and how these inner magnetospheric populations interact with ULF waves resulting from interplanetary solar wind structures.

Simulation studies of the Earth's radiation belts and ring current are very useful in understanding the acceleration, transport, and loss of energetic particles. Recently, the Comprehensive Ring Current Model (CRCM) and the Radiation Belt Environment (RBE) model were merged to form a Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. CIMI solves for many essential quantities in the inner magnetosphere, including ion and electron distributions in the ring current and radiation belts, plasmaspheric density, Region 2 currents, convection potential, and precipitation in the ionosphere. It incorporates whistler mode chorus and hiss wave diffusion of energetic electrons in energy, pitch angle, and cross terms. CIMI thus represents a comprehensive model that considers the effects of the ring current and plasmasphere on the radiation belts. We have performed a CIMI simulation for the storm on 5–9 April 2010 and then compared our results with data from the Two Wide-angle Imaging Neutral-atom Spectrometers and Akebono satellites. We identify the dominant energization and loss processes for the ring current and radiation belts. We find that the interactions with the whistler mode chorus waves are the main cause of the flux increase of MeV electrons during the recovery phase of this particular storm. When a self-consistent electric field from the CRCM is used, the enhancement of MeV electrons is higher than when an empirical convection model is applied. We also demonstrate how CIMI can be a powerful tool for analyzing and interpreting data from the new Van Allen Probes mission.

There is a growing concern over possible severe societal consequences related to adverse space weather impacts on man-made technological infrastructure. In the last two decades, significant progress has been made toward the first-principles modeling of space weather events, and three-dimensional (3-D) global magnetohydrodynamics (MHD) models have been at the forefront of this transition, thereby playing a critical role in advancing our understanding of space weather. However, the modeling of extreme space weather events is still a major challenge even for the modern global MHD models. In this study, we introduce a specially adapted University of Michigan 3-D global MHD model for simulating extreme space weather events with a Dst footprint comparable to the Carrington superstorm of September 1859 based on the estimate by Tsurutani et. al. (2003). Results are presented for a simulation run with “very extreme” constructed/idealized solar wind boundary conditions driving the magnetosphere. In particular, we describe the reaction of the magnetosphere-ionosphere system and the associated induced geoelectric field on the ground to such extreme driving conditions. The model setup is further tested using input data for an observed space weather event of Halloween storm October 2003 to verify the MHD model consistence and to draw additional guidance for future work. This extreme space weather MHD model setup is designed specifically for practical application to the modeling of extreme geomagnetically induced electric fields, which can drive large currents in ground-based conductor systems such as power transmission grids. Therefore, our ultimate goal is to explore the level of geoelectric fields that can be induced from an assumed storm of the reported magnitude, i.e., Dst∼=−1600 nT.

[1] The diffuse aurora has recently been shown to be a major contributor of energy flux into the Earth's ionosphere. Therefore, a comprehensive theoretical analysis is required to understand its role in energy redistribution in the coupled ionosphere-magnetosphere system. In previous theoretical descriptions of precipitated magnetospheric electrons (E∼ 1 keV), the major focus has been the ionization and excitation rates of the neutral atmosphere and the energy deposition rate to thermal ionospheric electrons. However, these precipitating electrons will also produce secondary electrons via impact ionization of the neutral atmosphere. This paper presents the solution of the Boltzman-Landau kinetic equation that uniformly describes the entire electron distribution function in the diffuse aurora, including the affiliated production of secondary electrons (E< 600 eV) and their ionosphere-magnetosphere coupling processes. In this article, we discuss for the first time how diffuse electron precipitation into the atmosphere and the associated secondary electron production participate in ionosphere-magnetosphere energy redistribution.

[1] Plasma dynamics in the inner magnetosphere are greatly affected by variations in the ion composition. The ratio of hydrogen to oxygen has been shown to be highly dependent on geomagnetic activity. To investigate this dependence, we examine the timing of the injection and subsequent evolution of O+ in the ring current during the storm of 6 August 2011 as observed by Two Wide-angle Imaging Neutral-atom Spectrometers (TWINS) instruments. To help interpret magnetospheric evolution using the global O Energetic Neutral Atom (ENA) emission measured by TWINS, we have employed a multidomain modeling of the global magnetosphere using the Space Weather Modeling Framework. TWINS Energetic Neutral Atom (ENA) imagers have the capability to distinguish between H and O emission and thus the major ion constituents of the ring current. Global composition measurements from TWINS spacecraft show intensifications of the oxygen ENA emission and thus an increase in the transport of ionospheric oxygen into the ring current that occur during the main phase of the storm. Both the observations and the simulation suggests that the peak in O ENA emission is correlated with the substorm injections that occurred during this time. The model also shows a very dynamic magnetosphere that allows for loss of oxygen from the Earth-magnetosphere system through plasmoids capable of transporting oxygen down the tail throughout the magnetic storm. This can possibly be a predominant pathway for loss of oxygen from the magnetosphere.

[1] We have recently extended the global magnetohydrodynamic (MHD) model BATS-R-US to account for pressure anisotropy. Since the inner magnetosphere dynamics cannot be fully described even by anisotropic MHD, we coupled our anisotropic MHD model with two inner magnetospheric models: the Rice Convection Model (RCM) and the Comprehensive Ring Current Model (CRCM). The coupled models provide better representations of the near-Earth plasma, especially during geomagnetic storms. In this paper, we present the two-way coupling algorithms with both ring current models. The major difference between these two couplings is that the RCM assumes isotropic and constant pressures along closed field lines, while the CRCM resolves pitch angle anisotropy. For model validation, we report global magnetosphere simulations performed by the coupled models. The simulation results are compared to the results given by the coupled isotropic MHD and ring current models. We find that in the global MHD simulations coupled with ring current models, pressure anisotropy results in a thinner magnetosheath, a shorter tail, a much smaller Earthward plasma jet from the tail reconnection site, and is also important in controlling the magnetic field configuration. The comparisons with satellite data for the magnetospheric event simulations show improvements on reproducing the measured tail magnetic field and inner magnetospheric flow velocity when including pressure anisotropy in the ring current model coupled global MHD model.

[1] A new, efficient, and highly accurate method for tracing magnetic separators in global magnetospheric simulations with arbitrary clock angle is presented. The technique is to begin at a magnetic null and iteratively march along the separator by finding where four magnetic topologies meet on a spherical surface. The technique is verified using exact solutions for separators resulting from an analytic magnetic field model that superposes dipolar and uniform magnetic fields. Global resistive magnetohydrodynamic simulations are performed using the three-dimensional Block-Adaptive Tree Solar-wind Roe-type Upwind Scheme code with a uniform resistivity, in eight distinct simulations with interplanetary magnetic field (IMF) clock angles ranging from 0° (parallel) to 180° (antiparallel). Magnetic nulls and separators are found in the simulations, and it is shown that separators traced here are accurate for any clock angle, unlike the last closed field line on the Sun-Earth line that fails for southward IMF. Trends in magnetic null locations and the structure of magnetic separators as a function of clock angle are presented and compared with those from the analytic field model. There are many qualitative similarities between the two models, but quantitative differences are also noted. Dependence on solar wind density is briefly investigated.

[1] This paper reports the metrics-based results of the Dst index part of the 2008–2009 GEM Metrics Challenge. The 2008–2009 GEM Metrics Challenge asked modelers to submit results for four geomagnetic storm events and five different types of observations that can be modeled by statistical, climatological or physics-based models of the magnetosphere-ionosphere system. We present the results of 30 model settings that were run at the Community Coordinated Modeling Center and at the institutions of various modelers for these events. To measure the performance of each of the models against the observations, we use comparisons of 1 hour averaged model data with the Dst index issued by the World Data Center for Geomagnetism, Kyoto, Japan, and direct comparison of 1 minute model data with the 1 minute Dst index calculated by the United States Geological Survey. The latter index can be used to calculate spectral variability of model outputs in comparison to the index. We find that model rankings vary widely by skill score used. None of the models consistently perform best for all events. We find that empirical models perform well in general. Magnetohydrodynamics-based models of the global magnetosphere with inner magnetosphere physics (ring current model) included and stand-alone ring current models with properly defined boundary conditions perform well and are able to match or surpass results from empirical models. Unlike in similar studies, the statistical models used in this study found their challenge in the weakest events rather than the strongest events.

[1] We present the coupling methodology and validation of a fully coupled inner and global magnetosphere code using the infrastructure provided by the Space Weather Modeling Framework (SWMF). In this model, the Comprehensive Ring Current Model (CRCM) represents the inner magnetosphere, while the Block‒Adaptive‒Tree Solar‒Wind Roe‒Type Upwind Scheme (BATS‒R‒US) represents the global magnetosphere. The combined model is a global magnetospheric code with a realistic ring current and consistent electric and magnetic fields. The computational performance of the model was improved to surpass real‒time execution by the use of the Message Passing Interface (MPI) to parallelize the CRCM. Initial simulations under steady driving found that the coupled model resulted in a higher pressure in the inner magnetosphere and an inflated closed field‒line region as compared to simulations without inner‒magnetosphere coupling. Our validation effort was split into two studies. The first study examined the ability of the model to reproduce Dst for a range of events from the Geospace Environment Modeling (GEM) Dst Challenge. It also investigated the possibility of a baseline shift and compared two approaches to calculating Dst from the model. We found that the model did a reasonable job predicting Dst and Sym‒H according to our two metrics of prediction efficiency and predicted yield. The second study focused on the specific case of the 22 July 2009 moderate geomagnetic storm. In this study, we directly compare model predictions and observations for Dst, THEMIS energy spectragrams, TWINS ENA images, and GOES 11 and 12 magnetometer data. The model did an adequate job reproducing trends in the data. Moreover, we found that composition can have a large effect on the result.

A self-consistent approach to superthermal electron (SE) transport along closed field lines in the inner magnetosphere is used to examine the concept of plasmaspheric transparency, magnetospheric trapping, and SE energy deposition to the thermal electrons. The dayside SE population is generated both by photoionization of the thermosphere and by secondary electron production from impact ionization when the photoelectrons collide with upper atmospheric neutral particles. It is shown that a self-consistent approach to this problem produces significant changes, in comparison with other approaches, in the SE energy exchange between the plasmasphere and the two magnetically conjugate ionospheres. In particular, plasmaspheric transparency can vary by a factor of two depending on the thermal plasma content along the field line and the illumination conditions of the two conjugate ionospheres. This variation in plasmaspheric transparency as a function of thermal plasma and ionospheric conditions increases with L-shell, as the field line gets longer and the equatorial pitch angle extent of the fly-through zone gets smaller. The inference drawn from these results is that such a self-consistent approach to SE transport and energy deposition should be included to ensure robustness in ionosphere-magnetosphere modeling networks.

The fluxes of energetic particles in the radiation belts are found to be strongly controlled by the solar wind conditions. In order to understand and predict the radiation particle intensities, we have developed a physics-based Radiation Belt Environment (RBE) model that considers the influences from the solar wind, ring current and plasmasphere. Recently, an improved calculation of wave-particle interactions has been incorporated. In particular, the model now includes cross diffusion in energy and pitch-angle. We find that the exclusion of cross diffusion could cause significant overestimation of electron flux enhancement during storm recovery. The RBE model is also connected to MHD fields so that the response of the radiation belts to fast variations in the global magnetosphere can be studied. We are able to reproduce the rapid flux increase during a substorm dipolarization on 4 September 2008. The timing is much shorter than the time scale of wave associated acceleration. Published by Elsevier Ltd.

We have integrated the Fok radiation belt environment (RBE) model into the space weather modeling framework (SWMF). RBE is coupled to the global magnetohydrodynamics component (represented by the Block-Adaptive-Tree Solar-wind Roe-type Upwind Scheme, BATS-R-US, code) and the Ionosphere Electrodynamics component of the SWMF, following initial results using the Weimer empirical model for the ionospheric potential. The radiation belt (RB) model solves the convection-diffusion equation of the plasma in the energy range of 10 keV to a few MeV. In stand-alone mode RBE uses Tsyganenko's empirical models for the magnetic field, and Weimer's empirical model for the ionospheric potential. In the SWMF the BATS-R-US model provides the time dependent magnetic field by efficiently tracing the closed magnetic field-lines and passing the geometrical and field strength information to RBE at a regular cadence. The ionosphere electrodynamics component uses a two-dimensional vertical potential solver to provide new potential maps to the RBE model at regular intervals. We discuss the coupling algorithm and show some preliminary results with the coupled code. We run our newly coupled model for periods of steady solar wind conditions and compare our results to the RB model using an empirical magnetic field and potential model. We also simulate the RB for an active time period and find that there are substantial differences in the RB model results when changing either the magnetic field or the electric field, including the creation of an outer belt enhancement via rapid inward transport on the time scale of tens of minutes.

It is well established that the development of pitch angle anisotropy in ring current ions provides a favorable condition for electromagnetic ion cyclotron (EMIC) wave growth. Many studies have calculated the EMIC growth rate assuming a bi-Maxwellian or ring distribution of energetic ions in velocity space. To test the appropriateness of these assumptions, we compute the EMIC growth rate by using the realistic ion distributions output by our Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. The CIMI model calculates energetic ion and electron distributions and many other essential quantities in the inner magnetosphere and ionosphere. We found that the ring current distributions are often very different from a bi-Maxwellian or a ring distribution. The corresponding EMIC growth rates are likewise quite distinctive from those estimated by functional fits. In the simulation of an idealized magnetic storm, we found that much stronger EMIC growth near the dayside magnetopause was attained with realistic ring current distributions. On the other hand, with the assumption of a ring distribution, EMIC growth at the peak of the ring current pressure was overestimated. Our calculation using the realistic ion distribution function reproduces the observed occurrence of EMIC waves and demonstrates the need to use realistic ring current ion distributions when modeling EMIC wave excitation.

We present high resolution Hall MHD simulations of Ganymede's magnetosphere demonstrating that Hall electric fields in ion-scale magnetic reconnection layers have significant global effects not captured in resistive MHD simulations. Consistent with local kinetic simulations of magnetic reconnection, our global simulations show the development of intense field-aligned currents along the magnetic separatrices. These currents extend all the way down to the moon's surface, where they may contribute to Ganymede's aurora. Within the magnetopause and magnetotail current sheets, Hall currents in the reconnection plane accelerate ions to the local Alfv\'en speed in the out-of-plane direction, producing a global system of ion drift belts that circulates Jovian magnetospheric plasma throughout Ganymede's magnetosphere. We discuss some observable consequences of these Hall-induced currents and ion drifts: the appearance of a sub-Jovian "double magnetopause" structure, an Alfv\'enic ion jet extending across the upstream magnetopause and an asymmetric pattern of magnetopause Kelvin-Helmholtz waves., Comment: 14 pages, 12 figures; presented at Geospace Environment Modeling (GEM) workshop (June, 2014) and Fall American Geophysical Union (AGU) meeting (December, 2014); submitted to Journal of Geophysical Research, December 2014

We compare estimates for the ion fluxes of twelve expected constituents of the lunar exosphere with estimates for the ion fluxes ejected from the lunar surface by solar wind ions and electrons. Our estimates demonstrate that measurements of lunar ions will help constrain the abundances of many undetected species in the lunar exosphere, particularly AI and Si, because the expected ion flux levels from the exosphere exceed those from the surface. To correctly infer the relative abundances of exospheric ions and neutrals from Kaguya Ion Mass Analyzer (IMA) measurements, we must take into account the velocity distributions of local ions. The predicted spectrum underestimates the measured levels of 0+ relative to other lunar ion species, a result that may suggest contributions by molecular ions to the measured 0+ rates.