The heliospheric energetic neutral atoms (ENAs) are products of charge exchange between solar wind and pick-up ions and interstellar neutral atoms. They are created in different regions of the heliosphere and its boundary region with the interstellar medium, constituting different ENA populations, and they carry information about their parent populations and production processes. Thus, ENAs enable mapping of the global structure of the heliosphere and the processes within and at its edge. Three instruments have provided sky maps of the heliospheric ENAs from 200 eV up to 44 keV over the solar activity cycle. The IBEX-Lo and IBEX-Hi instruments on board the Interstellar Boundary Explorer (IBEX) have provided ENA sky maps from 200 eV to 4.3 keV (central energy) from 2009 throughout Solar Cycle 24. The Ion and Neutral Camera (INCA) on board the Cassini–Huygens mission provided sky maps of the ENAs from 8–44 keV (central energy) on the spacecraft route to and in orbit around Saturn through Solar Cycles 23 and 24. We compare large-scale structures of ENA enhancements present across the sky maps in a wide energy range based on IBEX-Lo, IBEX-Hi, and INCA. They include Ribbon, heliotail lobes, and upwind ENA enhancement. We report on similarities and differences observed, including the evolution of the Ribbon from low to higher energies, and the presence of confined north and south heliotail lobes up to 44 keV.

This paper presents a comprehensive exploration of the Interstellar Boundary Explorer energetic neutral atom (ENA) ribbon, focusing on its spatial and temporal variations over 14 yr. Methodological advancements, including a refined map modeling procedure and a new ribbon separation technique with appropriate error propagation, enable a detailed investigation of the ribbon’s features. Utilizing statistically robust metrics, this study reveals details of the ribbon across energy and time. Key findings include energy- and time-dependent variations in flux, angular radius, ribbon profile width, and higher moments. By applying these metrics, we reveal new complexity to the evolution of the ribbon over time, highlighting the nuanced relationship between it and the solar wind. Furthermore, the study examines for the first time the ribbon as it passes through the starboard/heliotail region (Lon _EC 120°–180°), revealing properties distinct from other portions of the ribbon. The analysis uncovers an anticorrelation between ribbon width and flux, which provides quantitative support for a multisource ribbon created by a combination of solar wind neutrals that generate a spatiall narrow ribbon component and heliosheath neutrals giving rise to a broad component. Finally, differences in the temporal evolution of the ENA flux at different energies provide additional support that the location of the ribbon source region is beyond the heliopause.

Strong Thermal Emission Velocity Enhancement (STEVE) is a latitudinally narrow, purple-band emission observed at subauroral latitudes. Stable Auroral Red (SAR) arcs characterized by major red emission, and red/green arcs with both red and green emissions also occur at subauroral latitudes. Characteristics of magnetospheric source plasma and electromagnetic fields of these three types of arcs have not been fully understood because of the limited conjugate observations between magnetosphere and the ground. In this study, we report 11 conjugate observations (2 STEVEs, 7 SAR arcs, and 2 red/green arcs), using all-sky images obtained at seven ground stations over more than four years from January 2017 to April 2021 and magnetospheric satellites (Arase and Van Allen Probes). We found that, in the inner magnetosphere, the source region of STEVEs and red/green arcs were located outside the plasmasphere, and that of the SAR arc was in the region of spatial overlap between the plasmasphere and ring current region. Electromagnetic waves at frequencies below 1 Hz were observed for STEVEs and red/green arcs. SuperDARN radar data showed a strong westward plasma flow in the ionosphere, especially during STEVE events, whereas the plasma flows associated with SAR arcs and red/green arcs were generally weaker and variable. The STEVE and SAR arc can appear simultaneously at slightly different latitudes and STEVEs and red/green arcs can transform into SAR arcs. These first comprehensive ground-satellite measurements of three types of subauroral-latitude auroras increase our understanding on similarlity, differences, and coupling of these auroras in the ionosphere and the magnetosphere.

This study analyzes energetic neutral atom (ENA) spectral properties across distinct regions of globally distributed flux (GDF) sky maps, using Interstellar Boundary Explorer data from a full solar cycle, corrected for time dispersion. By time-shifting the data to the heliosheath using GDF source distances from D. B. Reisenfeld et al., we achieve a more accurate representation of heliosheath GDF energy spectra. We quantify ENA spectral characteristics, heliosheath line-of-sight-integrated proton pressure, and heliosheath proton temperature, comparing these to solar wind properties at 1 au and interplanetary scintillation-derived solar wind data. Our findings show that the spectral index is generally anticorrelated with heliosheath proton temperature and pressure, except in the central tail, where a partial positive correlation is observed. The lowest spectral index values occur when high-latitude heliosheath regions are dominated by fast solar wind from polar coronal holes. The south pole exhibits the flattest energy spectra due to plasma heating from both fast solar wind and a late-2014 pressure pulse. The central tail shows shorter variability (5–6 yr) for spectral index and heliosheath proton temperature, while proton pressure follows the 11 yr solar cycle. Most spectral shapes exhibit a “knee” distribution, peaking during solar maximum, with an “ankle” shape observed only at the south pole during solar cycle transitions. Asymmetry in proton pressure in the lobes is driven by the draping effect of the local interstellar magnetic field. This study provides insights into the energetic properties of GDF across the heliosphere, enhancing our understanding of the heliospheric environment.

The Interstellar Boundary Explorer (IBEX) has been observing the outer heliosphere and its interactions with the very local interstellar medium (VLISM) via measurements of energetic neutral atoms (ENAs) for over 14 yr. We discovered the IBEX Ribbon—a structure completely unanticipated by any prior theory or model—that almost certainly resides beyond the heliopause in the VLISM. We also characterized the other major source of heliospheric ENAs, the globally distributed flux (GDF), produced largely in the heliosheath between the termination shock and heliopause. In this study, we make three major new contributions. First, we validate, provide, and analyze the most recent 3 yr of IBEX-Hi (0.5–6 keV FWHM) data (2020–2022) for the first time. Second, we link these observations to the prior 11 yr of observations, exploring long-term variations. Finally, we provide the first IBEX team-validated Ribbon/GDF separation scheme and separated maps. Because of the uncertainty in separating different line-of-sight integrated sources, we provide not just best guess (median) maps, but also maps with upper and lower reasonable values of Ribbon and GDF fluxes, along with bounding fluxes that add the uncertainties to the upper and lower values. This allows theories and models to be compared with a range of possible values that the IBEX team believes are consistent with data. These observations, along with the reanalysis of the prior 11 yr of IBEX-Hi data, provide new insights and even further develop our detailed understanding of the heliosphere’s interaction with the local interstellar medium unlocked by IBEX.

The present study compares a single-band chorus wave against a banded chorus wave observed by Van Allen Probes at adjacent times, and demonstrates that the single-band chorus wave is associated with an anisotropic electron population over a broad energy range, while the banded chorus wave is accompanied by an electron phase space density plateau and an electron anisotropy reduction around Landau resonant energies. We further compare banded chorus waves with different spectral gap widths, and show that a wider spectral gap is associated with electron isotropization extending to higher energies with respect to the equatorial Landau resonant energy. We suggest that early generated chorus waves isotropize electrons via Landau resonant acceleration, and the waves that propagate to higher latitudes isotropize electrons at higher energies. The isotropization extending to higher energies leads to a larger spectral gap of new chorus waves after electrons bounce back to the equator. Naturally occurring chorus waves in the Earth's magnetosphere typically consist of two frequency bands. The present study aims to explain what controls the bandwidth of the chorus frequency gap that separates chorus waves into two bands. We first compare a single-band chorus wave against a banded chorus wave observed by a Van Allen Probe satellite at adjacent times. The banded chorus wave is accompanied by an electron phase space density plateau and an electron anisotropy reduction due to Landau resonance, while this phenomenon is not clearly seen in association with the single-band chorus wave. We further compare banded chorus waves with different gap widths. Satellite observations indicate that a wider frequency gap is associated with electron isotropization extending to higher energies. We suggest that Landau resonant acceleration extending to high latitudes isotropizes electron distribution at high energies, leading to new chorus waves with a large frequency gap. In contrast, Landau resonance that stops at a relatively lower latitude (due to waves being damped) leads to new chorus waves with a smaller frequency gap. Freshly injected anisotropic electron population without a PSD plateau generates single band chorus waves Banded chorus waves are more common because electrons usually have already undergone isotropization at Landau resonant energies along the drift path Landau acceleration extending to higher energies occurring at higher latitudes leads to more pronounced electron isotropization and larger chorus spectral gaps

The present study uncovers the fine structures of magnetosonic waves by investigating the EFW waveforms measured by Van Allen Probes. We show that each harmonic of the magnetosonic wave may consist of a series of elementary rising-tone emissions, implying a nonlinear mechanism for the wave generation. By investigating an elementary rising-tone magnetosonic wave that spans a wide frequency range, we show that the frequency sweep rate is likely proportional to the wave frequency. Furthermore, we reveal that each elementary rising-tone magnetosonic waves consist of multiple mini-harmonics spaced at O+ gyrofrequency. We reveal that O+ ions can suppress the generation of magnetosonic waves at multiples of O+ gyrofrequency, resulting in the mini-harmonic structure. The commonly observed mini-harmonics indicate an energy transfer between different ion species.

Abstract Combined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 RE) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences <25%. In addition, the MHD model subsolar magnetopause is about 0.5 RE sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

The interaction of our protective heliosphere and the Very Local Interstellar Medium (VLISM) is the least explored and most rewarding frontier of space physics.New evidence amplifies the central role of the heliosphere in the evolution of the solar system along its 4.6billion-year journey around the galaxy.In addition to the dense clouds of plasma, gas and dust seeding the early proto solar nebula, recent supernovae have left the entire solar system exposed to extreme fluxes of interstellar material and cosmic radiation with far-reaching implications.Our current knowledge lacks the direct measurements necessary to understand how our star upholds its vast heliosphere and its potentially game-changing role in the evolution of our galactic home.Interstellar Probe provides new, required measurements over more than a solar cycle to uncover the physical processes starting near the Sun responsible for creating our dynamic heliosphere.In April 2022, the pragmatic Interstellar Probe Mission Concept Study was completed after four years, detailing a Large Strategic heliophysics mission that would transect the heliosphere from 1 au to the VLISM.Its journey provides rich science for generations across heliophysics and presents an opportunity to push the frontier of space exploration farther than ever done before.Modest crossdivisional investments enable high-value planetary science and astrophysics, deepening our understanding of the emergence of our habitable planetary system.A trajectory through the forward hemisphere of the heliosphere would be accomplished by a launch in the 2036-2042 timeframe using conventional chemical propulsion and a heavy-lift launch vehicle, such as the Space Launch System (SLS).A Jupiter Gravity Assist could propel an 860-kg spacecraft with an 87-kg payload of ten instruments delivering a unified view of the global heliosphere, reaching the VLISM after 16 years.The spacecraft is designed to a 50-year nominal lifetime using modern-day technology based on successful missions like New Horizons.Two next-generation Radioisotope Thermal Generators (RTGs) would ensure 300 We at end of nominal mission at 375 au and could enable exploration even beyond 500 au.

An explosion of largely independent discoveries over the last decades has revolutionized our understanding of solar interactions from the Sun's surface, through the solar wind, and into the outer heliosphere.However, the key linkages between the parts of this complex system remain unresolved -a system-wide understanding is both finally possible and urgently needed.The fundamental question remains unanswered: how are the Sun's polar magnetic fields related with photospheric flows, and connected out into the interplanetary magnetic field and 3D solar wind?The time-dependent changes, observed from the Sun through the outer reaches of the solar wind separating our solar system from the local interstellar medium, play a critical role in determining the levels of radiation in the space environment, and the hazards associated with abrupt solar events that threaten satellites in space as well as the future of human space exploration.As the Sun appears to be moving away from the longest solar grand maximum in 9,300 years, the questions concerning polar fields, their evolution and relationship with the solar dynamo, with the heliospheric magnetic field and solar wind, and with the solar system's boundaries in the interstellar medium, need to be resolved urgently.The fundamental discoveries of the Ulysses mission unveiled the 3D structure of the solar wind and initially explored how this structure is intrinsically linked to the solar corona.Together the 3D observations of our changing Sun and outer heliospheric science, which have emerged through global observations, place the science of the Sun-2-Stars (S2S) mission as perhaps the most exciting and fastest developing area in Heliophysics.The S2S mission utilizes all existing and flight proven technologies including a Jupiter Gravity assist to launch into the first nearly circular-polar orbit (65º inclination, ~5 au, 12.5 yr) enabling a remarkable suite of in situ observations, imaging of the polar-3D Sun to the most distant reaches of our solar system and out beyond into the realms of the interstellar medium.The S2S ( ® ☆) mission concept links the changing Sun, including never-viewed polar plasma structure, and the first imaging of solar ejecta on the ecliptic, with the first contemporaneous measurements of solar wind, pickup ions and suprathermal ions from beyond the ionization cavity and at high latitudes.S2S includes the first energetic neutral atom heliospheric mapping from high latitudes, and interstellar neutral atom and interstellar dust from high latitudes to discover how our home in space is interlinked from the Sun's interior into the solar wind and out beyond our interstellar boundaries.S2S is the next vital and urgent step in Heliophysics that fundamentally transforms and unifies our view of the Sun, connecting solar interactions through the 3D heliosphere and out into the realm of the local galactic environment between the stars.

Diffuse aurora at the Earth’s high latitude regions is mainly caused by the low-energy (0.1–30 keV) electron precipitation which carries the major energy flux into the nightside upper atmosphere. Previous studies have demonstrated that combined scattering by the upper- and lower- band chorus waves acts as the dominant cause of diffuse auroral precipitation, but that is not necessarily the case as these two types of waves do not always occur simultaneously, with the lower-band more often. Here we report that the lower-band chorus satisfying the preferred condition can generate their second harmonics so as to trigger the diffuse auroral electron precipitation. We find that the lower-band chorus alone can only cause the precipitation of electrons greater than 4 keV, while the self-consistently generated second harmonic is weak but still able to result in the electron precipitation below 4 keV. The combined effect of those modes results in the observed pancake electron distributions and the diffuse aurora. Our results clearly demonstrate an alternative but universal mechanism of chorus-driven diffuse aurora in the Earth, which may also apply to the auroral formation in other planetary magnetospheres.

Isolated proton auroras (IPAs) appearing at subauroral latitudes are generated by energetic protons precipitating from the magnetosphere through interaction with electromagnetic ion cyclotron (EMIC) waves. Thus, an IPA is the ionospheric projection of the spatial and temporal variation of wave‐particle interaction regions in the magnetosphere. In this study, we conducted unique multi‐event analysis of simultaneous observations of IPAs and their source regions on 22 April, 7 September, and 22 March 2018, using all‐sky imagers at subauroral latitudes and the Van Allen Probes. When the satellite footprint passed over the IPAs associated with ground Pc1 geomagnetic pulsations, locally generated He + ‐band EMIC waves with the same frequencies as the ground Pc1 pulsations were observed in all events. The IPAs and EMIC waves had comparable narrow widths in the latitudinal direction. The EMIC waves appeared during the rapid enhancement of the ring current proton flux at energy range of ∼10–50 keV, while they disappeared at the rapid decrease of the electron density. From these results, we conclude that the boundaries of the localized IPAs and EMIC waves were determined by the overlap region of energetic proton enhancement and the plasmasphere. This overlap of ring‐current protons and plasmasphere is a favorable condition for the pitch‐angle scattering of protons by the EMIC waves. Characteristic magnetic and electric field variations with the IPAs were not observed by the satellite, indicating that the IPAs were not accompanied by field‐aligned currents comparable to that of oval auroral arcs.

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.

We derive annual skymaps of the proton temperature in the inner heliosheath (IHS), and track their temporal evolution over the years from 2009 to 2016 of Interstellar Boundary Explorer observations. Other associated thermodynamic parameters also determined are the density, kappa, that is, the parameter that characterizes kappa distributions, temperature rate, polytropic index, and entropy. We exploit the theory of kappa distributions and their connection with polytropes, to (i) express a new polytropic quantity Π that remains invariant along streamlines where temperature and density may vary, (ii) parameterize the proton flux in terms of the Π invariant and kappa, and (iii) derive the temperature and density, respectively, from the slope and intercept of the linear relationship between kappa and logarithm of Π. We find the following thermodynamic characteristics: (1) Temperature sky-maps and histograms shifted to their lowest values in 2012 and their highest in 2015; (2) Temperature negatively correlated with density, reflecting the subisothermal polytropic behavior; (3) Temperature positively correlated with kappa, revealing characteristics of the mechanism responsible for generating kappa distributions; (4) Processes in IHS are sub-isothermal tending toward isobaric, consistent with previously published results; (5) Linear relationship between kappa and polytropic indices, revealing characteristics of the particle potential energy; and (6) Entropy positively correlated with polytropic index, aligned with the underlying theory that entropy increases towards the isothermal state where the kappa distribution reduces to the Maxwell Boltzmann description.

The excitation of electrostatic and/or electromagnetic waves in the plasma universe is often associated with anisotropic velocity distributions of charged particles. In Earth's inner magnetosphere, this anisotropy can gradually develop as particles injected from the magnetotail drift around the Earth at different speeds depending on their energy and pitch angle. Here, we show that the perpendicular‐moving and bouncing ions can be separated more abruptly near the injection front. These pitch‐angle filters are localized magnetic dip structures formed by the diamagnetic behavior of the injected particles, which can trap perpendicular‐moving ions and allow bouncing ions to overtake. The resulting ion anisotropy facilitates the rapid generation of electromagnetic ion cyclotron (EMIC) waves, which in turn can largely reshape the Van Allen radiation belts. This scenario is examined by case and statistical observations, together with numerical simulations that reproduce most of the observational signatures, to support the causal relationship between magnetic dips, anisotropic ion distributions, and localized excitation of EMIC waves. Our study highlights the important roles of magnetic dips in the inner magnetospheric dynamics, as pitch‐angle filters of the injected ions and traveling hotspots of EMIC wave activities.

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 ℓ 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 present study addresses two basic questions related to banded chorus waves in the Earth’s magnetosphere: 1) are chorus spectral gaps formed near the equatorial source region or during propagation away from the equator? and 2) why are chorus spectral gaps usually located below 0.5 f ce ( f ce : electron gyro‐frequency)? By analyzing Van Allen Probes data, we demonstrate that chorus spectral gaps are observed in the source region where chorus waves propagate both in the parallel and anti‐parallel directions to the magnetic field. Chorus spectral gaps below 0.5 f ce are associated with electron parallel acceleration at energies above the equatorial Landau resonant energies. We explain that initially generated chorus waves quickly isotropize the electron distribution through Landau resonant acceleration, and the isotropization occurs for higher energies at higher latitudes. The isotropized population, after returning to the magnetic equator, leads to a chorus gap typically below 0.5 f ce by suppressing wave excitation.

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.

Isolated proton auroras (IPAs) appearing at subauroral latitudes are generated by energetic protons precipitating from the magnetosphere through interaction with electromagnetic ion cyclotron (EMIC) waves. An IPA thus indicates the spatial scale and temporal variation of wave‐particle interactions in the magnetosphere. In this study, a unique event of simultaneous ground and magnetospheric satellite observations of two IPAs were conducted on March 16, 2015, using an all‐sky imager at Athabasca, Canada and Van Allen Probes. The Van Allen Probes observed two isolated EMIC waves with frequencies of ∼1 and 0.4 Hz at L ≈ 5.0 when the satellite footprint crossed over the two IPAs. This suggests that the IPAs were caused by localized EMIC waves. Proton flux at 5–20 keV increased locally when the EMIC waves appeared. Electron flux at energies below ∼500 eV also increased. Temperature anisotropy of the energetic protons was estimated at 1.5–2.5 over a wide L ‐value range of 3.0–5.2. Electron density gradually decreased from L = 3.5 to 5.4, suggesting that the EMIC wave at L ≈ 5.0 was located in the gradual plasmapause. From these observations, we conclude that the localized IPAs and associated EMIC waves took place because of localized enhancement of energetic proton flux and plasma density structure near the plasmapause. The magnetic field observed by the satellite showed small variation during the wave observation, indicating that the IPAs were accompanied by the weak field‐aligned current.

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.