Xianzhe Jia

MESSENGER observations of Mercury's dayside magnetosphere under extreme solar wind conditions

CLJ and RMW acknowledge support from the Natural Sciences and Engineering Research Council of Canada, and CLJ acknowledges support from MESSENGER Participating Scientist grant NNX11AB84G. The MESSENGER project is supported by the NASA Discovery Program under contracts NASW- 00002 to the Carnegie Institution of Washington and NAS5-97271 to The Johns Hopkins University Applied Physics Laboratory.

Evidence of a Global Magma Ocean in Io’s Interior

Magnetic field measurements made near Jupiter’s moon Io strengthen the evidence for a magma ocean in its interior. Extensive volcanism and high-temperature lavas hint at a global magma reservoir in Io, but no direct evidence has been available. We exploited Jupiter’s rotating magnetic field as a sounding signal and show that the magnetometer data collected by the Galileo spacecraft near Io provide evidence of electromagnetic induction from a global conducting layer. We demonstrate that a completely solid mantle provides insufficient response to explain the magnetometer observations, but a global subsurface magma layer with a thickness of over 50 kilometers and a rock melt fraction of 20% or more is fully consistent with the observations. We also place a stronger upper limit of about 110 nanoteslas (surface equatorial field) on the dynamo dipolar field generated inside Io.

The search for a subsurface ocean in Ganymede with Hubble Space Telescope observations of its auroral ovals

We present a new approach to search for a subsurface ocean within Ganymede through observations and modeling of the dynamics of its auroral ovals. The locations of the auroral ovals oscillate due to Jupiters time-varying magnetospheric field seen in the rest frame of Ganymede. If an electrically conductive ocean is present, the external time-varying magnetic field is reduced due to induction within the ocean and the oscillation amplitude of the ovals decreases. Hubble Space Telescope (HST) observations show that the locations of the ovals oscillate on average by 2.0 ◦ ± 1.3 ◦ . Our model calculations predict a significantly stronger oscillation by 5.8 ◦ ± 1.3 ◦ without ocean compared to 2.2 ◦ ± 1.3 ◦ if an ocean is present. Because the ocean and the no-ocean hypotheses cannot be separated by simple visual inspection of individual HST images, we apply a statistical analysis including a Monte Carlo test to also address the uncertainty caused by the patchiness of observed emissions. The observations require a minimum electrical conductivity of 0.09 S/m for an ocean assumed to be located between 150 km and 250 km depth or alternatively a maximum depth of the top of the ocean at 330 km. Our analysis implies that Ganymedes dynamo possesses an outstandingly low quadrupole-to-dipole moment ratio. The new technique applied here is suited to probe the interior of other planetary bodies by monitoring their auroral response to time-varying magnetic fields.

Comparison of 3D kinetic and hydrodynamic models to ROSINA-COPS measurements of the neutral coma of 67P/Churyumov-Gerasimenko

67P/Churyumov-Gerasimenko (hereafter 67P) is a Jupiter-family comet and the object of investigation of the European Space Agency mission Rosetta. This report presents the first full 3D simulation results of 67P’s neutral gas coma. In this study we include results from a direct simulation Monte Carlo method, a hydrodynamic code, and a purely geometric calculation which computes the total illuminated surface area on the nucleus. All models include the triangulated 3D shape model of 67P as well as realistic illumination and shadowing conditions. The basic concept is the assumption that these illumination conditions on the nucleus are the main driver for the gas activity of the comet. As a consequence, the total production rate of 67P varies as a function of solar insolation. The best agreement between the model and the data is achieved when gas fluxes on the night side are in the range of 7% to 10% of the maximum flux, accounting for contributions from the most volatile components. To validate the output of our numerical simulations we compare the results of all three models to in situ gas number density measurements from the ROSINA COPS instrument. We are able to reproduce the overall features of these local neutral number density measurements of ROSINA COPS for the time period between early August 2014 and January 1 2015 with all three models. Some details in the measurements are not reproduced and warrant further investigation and refinement of the models. However, the overall assumption that illumination conditions on the nucleus are at least an important driver of the gas activity is validated by the models. According to our simulation results we find the total production rate of 67P to be constant between August and November 2014 with a value of about 1 × 1026 molecules s−1.

Global MHD simulations of Mercury's magnetosphere with coupled planetary interior: Induction effect of the planetary conducting core on the global interaction

Mercurys comparatively weak intrinsic magnetic field and its close proximity to the Sun lead to a magnetosphere that undergoes more direct space-weathering interactions than other planets. A unique aspect of Mercurys interaction system arises from the large ratio of the scale of the planet to the scale of the magnetosphere and the presence of a large-size core composed of highly conducting material. Consequently, there is strong feedback between the planetary interior and the magnetosphere, especially under conditions of strong external forcing. Understanding the coupled solar wind-magnetosphere-interior interaction at Mercury requires not only analysis of observations but also a modeling framework that is both comprehensive and inclusive. We have developed a new global MHD model for Mercury in which the planetary interior is modeled as layers of different electrical conductivities that electromagnetically couple to the surrounding plasma environment. This new modeling capability allows us to characterize the dynamical response of Mercury to time-varying external conditions in a self-consistent manner. Comparison of our model results with observations by the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER) spacecraft shows that the model provides a reasonably good representation of the global magnetosphere. To demonstrate the capability to model induction effects, we have performed idealized simulations in which Mercurys magnetosphere is impacted by a solar wind pressure enhancement. Our results show that due to the induction effect, Mercurys core exerts strong global influences on the way Mercury responds to changes in the external environment, including modifying the global magnetospheric structure and affecting the extent to which the solar wind directly impacts the surface. The global MHD model presented here represents a crucial step toward establishing a modeling framework that enables self-consistent characterization of Mercurys tightly coupled planetary interior-magnetosphere system.

Cassini in situ observations of long-duration magnetic reconnection in Saturn’s magnetotail

Cassini’s encounter with Saturn’s magnetotail — the long magnetosphere region stretching into space — has revealed that plasma exits the magnetosphere through long-duration magnetic reconnection, which ejects ten times more mass than estimated. Magnetic reconnection is a fundamental process in solar system and astrophysical plasmas, through which stored magnetic energy associated with current sheets is converted into thermal, kinetic and wave energy1,2,3,4. Magnetic reconnection is also thought to be a key process involved in shedding internally produced plasma from the giant magnetospheres at Jupiter and Saturn through topological reconfiguration of the magnetic field5,6. The region where magnetic fields reconnect is known as the diffusion region and in this letter we report on the first encounter of the Cassini spacecraft with a diffusion region in Saturn’s magnetotail. The data also show evidence of magnetic reconnection over a period of 19 h revealing that reconnection can, in fact, act for prolonged intervals in a rapidly rotating magnetosphere. We show that reconnection can be a significant pathway for internal plasma loss at Saturn6. This counters the view of reconnection as a transient method of internal plasma loss at Saturn5,7. These results, although directly relating to the magnetosphere of Saturn, have applications in the understanding of other rapidly rotating magnetospheres, including that of Jupiter and other astrophysical bodies.

Extended magnetohydrodynamics with embedded particle‐in‐cell simulation of Ganymede's magnetosphere

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 Jupiters magnetospheric plasma and Ganymedes 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 Ganymedes 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.

Ion composition in interchange injection events in Saturn's magnetosphere

Interchange injection events are commonly observed by the Cassini spacecraft in the region between about 6 and 12 Rs (1 Rs = 60,268 km) and even frequently beyond. In this study, 13 examples of interchange injection events are identified in Cassini/Cassini Plasma Spectrometer data under special conditions such that time-of-flight (TOF) mass spectra could be obtained from entirely within the events. Using the TOF data to separate the main ion species H+, H2+, and W+, approximate densities of each species are calculated under the assumption that all distributions were isotropic. The light-ion density ratios, H2+/H+, in the injection events are not discernibly different from those ratios in control intervals from the ambient plasma. However, the water-group ratio, W+/H+, is significantly lower than ambient. The comparison of the measured density ratios with the range of values observed throughout Saturns magnetosphere indicates that the values of W+/H+ that are as low as those observed within the injection events are found primarily beyond L~14 (where L is the equatorial crossing distance, in Saturn radius, of a dipole field line), indicating that the injection events are delivering plasma from the outer magnetosphere at times traveling at least 6 Rs.

Self-consistent multifluid MHD simulations of Europa's exospheric interaction with Jupiter's magnetosphere

The Jovian moon, Europa, hosts a thin neutral gas atmosphere, which is tightly coupled to Jupiters magnetosphere. Magnetospheric ions impacting the surface sputter off neutral atoms, which, upon ionization, carry currents that modify the magnetic field around the moon. The magnetic field in the plasma is also affected by Europas induced magnetic field. In this paper we investigate the environment of Europa using our multifluid MHD model and focus on the effects introduced by both the magnetospheric and the pickup ion populations. The model self-consistently derives the electron temperature that governs the electron impact ionization process, which is the major source of ionization in this environment. The resulting magnetic field is compared to measurements performed by the Galileo magnetometer, the bulk properties of the modeled thermal plasma population is compared to the Galileo Plasma Subsystem observations, and the modeled surface precipitation fluxes are compared to Galileo Ultraviolet Spectrometer observations. The model shows good agreement with the measured magnetic field and reproduces the basic features of the plasma interaction observed at the moon for both the E4 and the E26 flybys of the Galileo spacecraft. The simulation also produces perturbations asymmetric about the flow direction that account for observed asymmetries.

Flux transfer event observation at Saturn's dayside magnetopause by the Cassini spacecraft

We present the first observation of a flux rope at Saturns dayside magnetopause. This is an important result because it shows that the Saturnian magnetopause is conducive to multiple X-line reconnection and flux rope generation. Minimum variance analysis shows that the magnetic signature is consistent with a flux rope. The magnetic observations were well fitted to a constant-α force-free flux rope model. The radius and magnetic flux content of the rope are estimated to be 4600–8300 km and 0.2–0.8 MWb, respectively. Cassini also observed five traveling compression regions (remote signatures of flux ropes), in the adjacent magnetosphere. The magnetic flux content is compared to other estimates of flux opening via reconnection at Saturn.