We define the spatial clustering of X-rays within Jupiter’s northern auroral regions by classifying their distributions into ‘X-ray auroral structures’. Using data from Chandra during Juno’s main mission observations (24 May 2016 – 8 September 2019), we define five X-ray structures based on their ionospheric location and calculate the distribution of auroral photons. The morphology and ionospheric location of these structures allow us to explore the possibility of numerous X-ray auroral magnetospheric drivers. We compare these distributions to Hubble Space Telescope (HST)and Juno (Waves and MAG) data, and a 1D solar wind propagation model to infer the state of Jupiter’s magnetosphere. Our results suggest that the five sub-classes of ‘X-ray structures’ fall under two broad morphologies: fully polar and low latitude emissions. Visibility modelling of each structure suggests the non-uniformity of the photon distributions across the Chandra intervals are likely associated with the switching on/off of magnetospheric drivers as opposed to geometrical effects. The combination of ultraviolet (UV) and X-ray morphological structures is a powerful tool to elucidate the behaviour of both electrons and ions and their link to solar wind/magnetospheric conditions in the absence of an upstream solar monitor.

To help understand and determine the driver of jovian auroral X-rays, we present the first statistical study to focus on the morphology and dynamics of the jovian northern hot spot (NHS) using Chandra data. The catalogue we explore dates from 18 December 2000 up to and including 8 September 2019. Using a numerical criterion, we characterize the typical and extreme behaviour of the concentrated NHS emissions across the catalogue. The mean power of the NHS is found to be 1.91 GW with a maximum brightness of 2.02 Rayleighs (R), representing by far the brightest parts of the jovian X-ray spectrum. We report a statistically significant region of emissions at the NHS center which is always present, the averaged hot spot nucleus (AHSNuc), with mean power of 0.57 GW and inferred average brightness of ∼ 1.2 R. We use a flux equivalence mapping model to link this distinct region of X-ray output to a likely source location and find that the majority of mappable NHS photons emanate from the pre-dusk to pre-midnight sector, coincident with the dusk flank boundary. A smaller cluster maps to the noon magnetopause boundary, dominated by the AHSNuc, suggesting that there may be multiple drivers of X-ray emissions. On application of timing analysis techniques (Rayleigh, Monte Carlo, Jackknife), we identify several instances of statistically significant quasi-periodic oscillations (QPOs) in the NHS photons ranging from ∼ 2.3-min to 36.4-min, suggesting possible links with ultra-low frequency activity on the magnetopause boundary (e.g. dayside reconnection, Kelvin-Helmholtz instabilities).

The Earth’s magnetosheath and cusps are the sources of soft X-rays. In the accompanying paper (Part 1) and this paper, we discuss the methods of finding the magnetopause position by analyzing the X-ray images. We use the software developed for the Soft X-ray Imager (SXI) on board the forthcoming Solar wind - Magnetosphere - Ionosphere Link Explorer (SMILE) mission. We show how to find the maximum SXI count rate in noisy count maps. We verify the assumption that the maximum of the X-ray emissivity integrated along the Line-of-Sight (Ix) is tangent to the magnetopause. We consider two cases using two MHD models and apply different methods of magnetospheric masking. Overall, the magnetopause is located close to the maximum Ix gradient or between the maximum Ix gradient and the maximum Ix depending on the method used. But since the angular distance between the maximum Ix gradient and the maximum Ix is relatively small (about 3{degree sign}), the maximum Ix might be used as an indicator of the outer boundary of a wide magnetopause layer usually obtained in MHD simulations.

The magnetopause standoff distance characterizes global magnetospheric compression and deformation in response to changes in the solar wind dynamic pressure and interplanetary magnetic field orientation. We cannot derive this parameter from in-situ spacecraft measurements. However, time-series of the magnetopause standoff distance can be obtained in the near future using observations by soft X-ray imagers. In two companion papers, we describe methods of finding the standoff distance from X-ray images. In Part 1, we present the results of MHD simulations which we use for the calculation of the X-ray emissivity in the magnetosheath and cusps. Some MHD models predict relatively high density in the magnetosphere, larger than observed in the data. Correcting this, we develop magnetospheric masking methods to separate the magnetosphere from the magnetosheath and cusps. We simulate the X-ray emissivity in the magnetosheath for different solar wind conditions and dipole tilts.

We model an interval of sustained northward interplanetary magnetic field, for which we have a comprehensive set of observational data. This interval is associated with the arrival of an interplanetary coronal mass ejection. The solar wind densities at the time are particularly high and the interplanetary magnetic field is primarily northward. This results in strong auroral emissions within the polar cap in a cusp spot, which we associate with lobe reconnection at the high-latitude magnetopause. We also observe areas of upwards field-aligned current within the summer Northern Hemisphere polar cap that exhibit large current magnitudes. The model is able to reproduce the spatial distribution of the field-aligned currents well, even under changing conditions in the incoming interplanetary magnetic field. Discrepancies exist between the modeled and observed current magnitudes. Notably, the winter Southern Hemisphere exhibits much lower current magnitudes overall. We also model a sharp transition of the location of magnetopause reconnection. This changes rapidly from a subsolar location at the low-latitude magnetopause under southward interplanetary magnetic field conditions, to a high-latitude lobe reconnection location when the field is northward. This occurs during a fast rotation of the IMF at the shock front of a magnetic cloud.

We present a statistical study of Jupiter’s disk X-ray emissions using 19 years of Chandra X-Ray Observatory (CXO) observations. Previous work has suggested that these emissions are consistent with solar X-rays elastically scattered from Jupiter’s upper atmosphere. We showcase a new Pulse Invariant (PI) filtering method that minimises instrumental effects which may produce unphysical trends in photon counts across the nearly-two-decade span of the observations. We compare the CXO results with solar X-ray flux data from the Geostationary Operational Environmental Satellites (GOES) X-ray Sensor (XRS) for the wavelength band 1-8 Å (long channel), to quantify the correlation between solar activity and jovian disk counts. We find a statistically significant Pearson’s Correlation Coefficient (PCC) of 0.9, which confirms that emitted jovian disk X-rays are predominantly governed by solar activity. We also utilise the high spatial resolution of the High Resolution Camera Instrument (HRC-I) on board the CXO to map the disk photons to their positions on Jupiter’s surface. Voronoi tessellation diagrams were constructed with the JRM09 (Juno Reference Model through Perijove 9) internal field model overlaid to identify any spatial preference of equatorial photons. After accounting for area and scattering across the curved surface of the planet, we find a preference of jovian disk emission at 2-3.5 Gauss surface magnetic field strength. This suggests that a portion of the disk X-rays may be linked to processes other than solar scattering: the spatial preference associated with magnetic field strength may imply increased precipitation from the radiation belts, as previously postulated.