A Full Solar Cycle of Interstellar Boundary Explorer (IBEX) Observations

NASA's Interstellar Boundary Explorer (IBEX) mission has operated in space for a full solar activity cycle (Solar Cycle 24), and IBEX observations have exposed the global three-dimensional structure of the heliosphere and its interaction with the very local interstellar medium for the first time. Here, we extend the prior IBEX observations of energetic neutral atoms (ENAs) by adding a comprehensive analysis of four additional years (2016 through 2019). We document several improvements and rerelease the entire 11 yr, IBEX-Hi data set. The new observations track the continuing expansion of the outer heliosphere's response to the large solar wind pressure increase in late 2014. We find that the intensification of ENAs from the heliosheath continued to expand progressively over time to directions farther from the initial, closest direction to the heliospheric boundaries, ∼20° south of the upwind direction. This expansion extended beyond the south pole in 2018 and the north pole in 2019, demonstrating that the termination shock and heliopause are closer in the south. The heliotail has not yet responded, indicating that the boundaries are significantly farther away in the downwind direction. Finally, the slow solar wind (∼1 keV) ENAs just started to intensify from the closest regions of the IBEX Ribbon. This is about two and a half years after the initial response from heliosheath ENAs and about four and a half years after the increase in solar wind output, both clearly implicating a "secondary ENA" source in the draped interstellar magnetic field, just beyond the heliopause.

Properties of the inner heliosheath (IHS) plasma are inferred from energetic neutral atom (ENA) observations by ∼1 au spacecraft. However, the Compton–Getting effect due to the plasma velocity relative to the spacecraft is rarely taken into account, even though the plasma speed is a significant fraction of the ENA speed. In this study, we transform Interstellar Boundary Explorer (IBEX) ENA spectra to the IHS plasma frame using flow profiles from a 3D heliosphere simulation. We find that proton spectra in the plasma frame are steeper by ∼30% to 5% at ∼0.5 to 6 keV, respectively, compared to ENAs in the spacecraft frame. While radial plasma flows contribute most to the Compton–Getting effect, transverse flows at mid/high latitudes and the heliosphere flanks account for up to ∼30% of the frame transformation for IBEX-Hi at ∼0.7 keV and up to ∼60% for IBEX-Lo at ∼0.1 keV. We determine that the majority of IHS proton fluxes derived from IBEX-Hi measurements in 2009–2016 are statistically consistent with power-law distributions, with mean proton index ∼2.1 and standard deviation ∼0.4. We find significantly fewer spectral breaks in IBEX observations compared to early analyses, which we determine were a product of the "ion gun" background prevalent in ∼2009–2012 before corrections made by McComas et al. in subsequent data releases. We recommend that future analyses of the IHS plasma utilizing ENA measurements take into account the Compton–Getting effect including radial and transverse flows, particularly IBEX and Interstellar Mapping and Acceleration Probe measurements below ∼10 keV.

In this study, we estimate the heliospheric termination shock (HTS) compression ratio at multiple directions in the sky from a quantitative comparison of the observed and simulated inner heliosheath (IHS) energetic neutral atom (ENA) fluxes. We use a 3D steady-state simulation of the heliosphere to simulate the ENA fluxes by postprocessing the MHD plasma using a multi-Maxwellian distribution for protons in the IHS. The simulated ENA fluxes are compared with time exposure–averaged IBEX-Hi data for the first 3 yr of the mission. The quantitative comparison is performed by calculating the fractional difference in the spectral slope between the observed and simulated ENA fluxes for a range of compression ratios, where the simulated ENA spectrum is varied as a function of downstream pickup ion temperature as a function of compression ratio. The estimated compression ratio in a particular direction is determined by the minimum value of the fractional difference in spectral slope. Our study shows that the compression ratio estimated by this method is in close agreement with the large-scale compression ratio observed by Voyager 2 in its travel direction. Also, the compression ratio in other directions near the ecliptic plane is similar to the compression ratio at the Voyager 2 direction. The weakest shock compression is found to be on the port side of the heliosphere at direction (27°, 15°). This is the first study to estimate the HTS compression ratio at multiple directions in the sky from IBEX data.

The Interstellar Boundary Explorer (IBEX) mission has shown that variations in the energetic neutral atom (ENA) flux from the outer heliosphere are associated with the solar cycle and longer-term variations in the solar wind (SW). In particular, there is a good correlation between the dynamic pressure of the outbound SW and variations in the later-observed IBEX ENA flux. The time difference between observations of the outbound SW and the heliospheric ENAs with which they correlate ranges from approximately 2 to 6 yr or more, depending on ENA energy and look direction. This time difference can be used as a means of "sounding" the heliosheath, that is, finding the average distance to the ENA source region in a particular direction. We apply this method to build a 3D map of the heliosphere. We use IBEX ENA data collected over a complete solar cycle, from 2009 through 2019, corrected for survival probability to the inner heliosphere. Here we divide the data into 56 "macropixels" covering the entire sky. As each point in the sky is sampled once every 6 months, this gives us a time series of 22 points macropixel^–1 on which to time-correlate. Consistent with prior studies and heliospheric models, we find that the shortest distance to the heliopause, d_HP, is slightly south of the nose direction (d_HP ∼ 110–120 au), with a flaring toward the flanks and poles (d_HP ∼ 160–180 au). The heliosphere extends at least ∼350 au tailward, which is the distance limit of the technique.

Remote imaging of plasmas in the heliosphere and very local interstellar medium is possible with energetic neutral atoms (ENAs), created through the charge exchange of protons with interstellar neutral atoms. ENA observations collected by the Interstellar Boundary Explorer (IBEX) revealed two distinctive sources. One source is the globally distributed flux (GDF), which extends over the entire sky and varies over large spatial scales. The other source encompasses only a narrow circular band in the sky and is called the IBEX ribbon. Here, we utilize the observed difference in spatial scales of these two ENA sources to separate them. We find that linear combinations of spherical harmonics up to degree ${{\ell }}_{\max }=3$ can reproduce most of the ENA fluxes observed outside the ribbon region. We use these combinations to model the GDF and the difference between the observed fluxes and the GDF yields estimation of the ribbon emission. The separated ribbon responds with a longer time delay to the solar wind changes than the GDF, suggesting a more distant source of the ribbon ENAs. Moreover, we locate the direction of the maximum plasma pressure based on the GDF. This direction is 17°.2 ± 0°.5 away from the upwind direction within the plane containing the interstellar flow and interstellar magnetic field vectors. This deflection is consistent with the expected position of the maximum external pressure at the heliopause. The maps with separated ribbon and GDF are posted concurrently with this paper and can be used to further study these two sources.

The Sun's motion through the interstellar medium leads to an interstellar neutral (ISN) wind through the heliosphere. Several ISN species, including He, moderately depleted by ionization are observed with pickup ions and directly imaged. Since 2009, analyzed Interstellar Boundary Explorer (IBEX) observations returned a precise 4D parameter tube associated with the bulk velocity vector and the temperature of ISN flow distribution. This 4D parameter tube is typically expressed in terms of the ISN speed, the inflow latitudinal direction, and the temperature as a function of the inflow longitudinal direction and the local flow Mach number. We have used IBEX observations and those from other spacecraft to reduce statistical parameter uncertainties: ${V}_{\mathrm{ISN}\infty }=25.99\pm 0.18$ km s⁻¹, ${\lambda }_{\mathrm{ISN}\infty }=75\buildrel{\circ}\over{.} 28\pm 0\buildrel{\circ}\over{.} 13$, ${\beta }_{\mathrm{ISN}\infty }={\rm{-5}}\buildrel{\circ}\over{.} 200\pm 0\buildrel{\circ}\over{.} 075$, and ${T}_{\mathrm{ISN}\infty }=7496\pm 172$ K. IBEX ISN viewing is restricted almost perpendicular to the Earth–Sun line, which limits observations in ecliptic longitude to ∼130° ± 30° and results in relatively small uncertainties across the IBEX parameter tube but large uncertainties along it. Operations over the last three years enabled the IBEX spin axis to drift to the maximum operational offset (7°) west of the Sun, helping to break the ISN parameter degeneracy by weakly crossing the IBEX parameter tubes: the range of possible inflow longitudes extends over the range ${\lambda }_{\mathrm{ISN}\infty }=75\buildrel{\circ}\over{.} {28}_{-2.21}^{+2.27}$ and the corresponding range of other ISN parameters is ${V}_{\mathrm{ISN}\infty }={25.99}_{-1.76}^{+1.86}$ km s⁻¹, ${\beta }_{\mathrm{ISN}\infty }={\rm{-5}}\buildrel{\circ}\over{.} {200}_{-0.085}^{+0.093}$, and ${T}_{\mathrm{ISN}\infty }={7496}_{-1528}^{+1274}$ K. This enhances the full χ² analysis of ISN parameters through comparison with detailed models. The next-generation IBEX-Lo sensor on IMAP will be mounted on a pivot platform, enabling IMAP-Lo to follow the ISN flow over almost the entire spacecraft orbit around the Sun. A near-continuous set of 4D parameter tube orientations on IMAP will be observed for He and for O, Ne, and H that cross at varying angles to substantially reduce the ISN flow parameter uncertainties and mitigate systematic uncertainties (e.g., from ionization effects and the presence of secondary components) to derive the precise parameters of the primary and secondary local interstellar plasma flows.

The IBEX-Lo instrument on board the Interstellar Boundary Explorer (IBEX) mission samples interstellar neutral (ISN) helium atoms penetrating the heliosphere from the very local interstellar medium (VLISM). In this study, we analyze the IBEX-Lo ISN helium observations covering a complete solar cycle, from 2009 through 2020 using a comprehensive uncertainty analysis including statistical and systematic sources. We employ the Warsaw Test Particle Model to simulate ISN helium fluxes at IBEX, which are subsequently compared with the observed count rate in the three lowest energy steps of IBEX-Lo. The χ² analysis shows that the ISN helium flows from ecliptic $\left(\lambda ,\beta \right)=(255\buildrel{\circ}\over{.} 59\pm 0\buildrel{\circ}\over{.} 23,5\buildrel{\circ}\over{.} 14\pm 0\buildrel{\circ}\over{.} 08)$, with speed v_HP = 25.86 ± 0.21 km s⁻¹ and temperature T_HP = 7450 ±140 K at the heliopause. Accounting for gravitational attraction and elastic collisions, the ISN helium speed and temperature in the pristine VLISM far from the heliopause are v_VLISM = 25.9 km s⁻¹ and T_VLISM = 6150 K, respectively. The time evolution of the ISN helium fluxes at 1 au over 12 yr suggests significant changes in the IBEX-Lo detection efficiency, higher ionization rates of ISN helium atoms in the heliosphere than assumed in the model, or an additional unaccounted for signal source in the analyzed observations. Nevertheless, we do not find any indication of the evolution of the derived parameters of ISN helium over the period analyzed. Finally, we argue that the continued operation of IBEX-Lo to overlap with the Interstellar Mapping and Acceleration Probe will be pivotal in tracking possible physical changes in the VLISM.

Measurements of starlight polarized by aligned interstellar dust grains are used to probe the relation between the orientation of the ambient interstellar magnetic field (ISMF) and the ISMF traced by the ribbons of energetic neutral atoms discovered by the Interstellar Boundary Explorer spacecraft. We utilize polarization data, many acquired specifically for this study, to trace the configuration of the ISMF within 40 pc. A statistical analysis yields a best-fit ISMF orientation, B_magpol, aligned with Galactic coordinates ℓ = 42°, b = 49°. Further analysis shows the ISMF is more orderly for "downfield" stars located over 90° from B_magpol. The data subset of downfield stars yields an orientation for the nearby ISMF at ecliptic coordinates λ, β ≈ 219° ± 15°, 43° ± 9° (Galactic coordinates l, b ≈ 40°, 56°, ±17°). This best-fit ISMF orientation from polarization data is close to the field direction obtained from ribbon models. This agreement suggests that the ISMF shaping the heliosphere belongs to an extended ordered magnetic field. Extended filamentary structures are found throughout the sky. A previously discovered filament traversing the heliosphere nose region, "Filament A," extends over 300° of the sky, and crosses the upwind direction of interstellar dust flowing into the heliosphere. Filament A overlaps the locations of the Voyager kilohertz emissions, three quasar intraday variables, cosmic microwave background (CMB) components, and the inflow direction of interstellar grains sampled by Ulysses and Galileo. These features are likely located in the upstream outer heliosheath where ISMF drapes over the heliosphere, suggesting Filament A coincides with a dusty magnetized plasma. A filament 55° long is aligned with a possible shock interface between local interstellar clouds. A dark spot in the CMB is seen within 5° of the filament and within 10° of the downfield ISMF direction. Two large magnetic arcs are centered on the directions of the heliotail. The overlap between CMB components and the aligned dust grains forming Filament A indicates the configuration of dust entrained in the ISMF interacting with the heliosphere provides a measurable foreground to the CMB.

The Interstellar Boundary Explorer (IBEX) is a NASA satellite in Earth orbit, dedicated to observing both interstellar neutral atoms entering the heliosphere and energetic neutral atoms (ENAs) from the interstellar boundaries from roughly 10 eV to 6 keV. This work presents the averaged maps, energy spectra, and temporal variability of heliospheric ENA intensities measured with the IBEX-Lo instrument at 1 au at energies between 10 eV and 2 keV, covering one entire solar cycle from 2009 through 2019. These results expand the range in time and energy for studying the globally distributed ENA flux and the IBEX Ribbon. The observed ENA intensities exceed model predictions, in particular below 500 eV. Moreover, the ENA intensities between 50–200 eV energy show an unexpected rise and fall around the year 2015 in most sky regions.

Energetic neutral atom (ENA) measurements by IBEX reveal that the heliotail comprises an energy-dependent multilobe structure. We examine the heliotail evolution over 11 yr of IBEX observations covering a full solar cycle (SC). We find the following: (1) The heliotail structure persists over the entire SC, comprising three ENA-enhanced and two ENA-suppressed lobes. (2) Lobe sizes and locations are generally stable but exhibit variations in ENA fluxes driven by the SC. (3) Lobe centers follow a cyclic behavior over multiple SC phases, indicating direct signatures of slow and fast solar wind (SW) interactions in the inner heliosheath (IHS). (4) The tilted plane passing through the port–starboard lobes' centers oscillates in latitude but maintains its tilt from the ecliptic plane, likely a consequence of the interstellar magnetic field draping around the heliosphere. (5) The transition of the central heliotail from a single lobe at ∼1.1 keV to two lobes above ∼2 keV is SC-dependent and directly reflects the IHS plasma properties, i.e., when ENA fluxes from fast SW from the polar coronal holes change over time. (6) The central lobe exhibits a substructure that is enhanced and offset from the downwind direction, possibly indicating an asymmetric ENA emission or an asymmetry in the parent plasma distribution. These results reveal the general stability of the heliotail structure over time and distinct variations in individual lobes' properties in relation to the SC phases. Furthermore, results show the effects of multiple SC phases in the tail, reflecting different ENA travel times and source histories.