Pickup Ions Research Articles

Interstellar Neutral He Parameters from Crossing Parameter Tubes with the Interstellar Mapping and Acceleration Probe Informed by 10 yr of Interstellar Boundary Explorer Observations

Abstract 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 ISN ∞ = 25.99 ± 0.18 km s−1, λ ISN ∞ = 75 .° 28 ± 0 .° 13 , β ISN ∞ = −5 .° 200 ± 0 .° 075 , and T ISN ∞ = 7496 ± 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 λ ISN ∞ = 75 .° 28 − 2.21 + 2.27 and the corresponding range of other ISN parameters is V ISN ∞ = 25.99 − 1.76 + 1.86 km s−1, β ISN ∞ = −5 .° 200 − 0.085 + 0.093 , and T ISN ∞ = 7496 − 1528 + 1274 K. This enhances the full χ 2 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.

Energetic pickup proton population downstream of the termination shock as revealed by IBEX-Hi data

ABSTRACT Pickup protons originate as a result of the ionization of hydrogen atoms in the supersonic solar wind, forming the suprathermal component of protons in the heliosphere. While they are being picked by the heliospheric magnetic field and convected into the heliosheath, the pickup protons may suffer stochastic acceleration from the solar wind turbulence in the region from the Sun up to the heliospheric termination shock, where they can also experience shock-drift acceleration or reflection from the cross-shock potential. These processes create a high-energy tail in the pickup ion energy distribution. The properties of this energetic pickup proton population are still not well defined, in spite of the fact that they are vital for models that simulate energetic neutral atom fluxes. We consider two scenarios for the pickup proton velocity distribution downstream of the heliospheric termination shock (a filled shell with an energetic power-law tail, and bi-Maxwellian). Based on a numerical kinetic model and observations of the energetic neutral atom fluxes from the inner heliosheath by the IBEX-Hi instrument, we characterize the pickup proton distribution and provide estimations of the properties of the energetic pickup proton population downstream of the termination shock.

Simulation of the Scattering of Continuously Injected Pickup Ions outside the Heliopause

Abstract Hybrid simulations in 2D space and 3D velocity dimensions with continuous injection of pickup ions (PUIs) provide insight into the plasma processes that are responsible for the pitch angle scattering of PUIs outside the heliopause. The present investigation includes for the first time continuous injection of PUIs and shows how the scattering depends on the energy of the PUIs and the strength of the background magnetic field as well as the dependence on the injection rate of the time for the isotropization of the pitch angle distribution. The results demonstrate that, with the gradual injection of PUIs of a narrow ring velocity distribution perpendicular to the background magnetic field, oblique mirror mode waves develop first, followed by the growth of quasiparallel propagating ion cyclotron waves. Subsequently, the PUIs are scattered by the excited waves and gradually approach an isotropic distribution. A time for isotropization is defined to be the time at which T ∣∣/T ⊥, i.e., the ratio of the parallel to perpendicular PUI thermal energy changes from ≈0 to ≈0.15. By varying the PUI injection rate, estimates of the time for the PUI distribution to be isotropized are presented. The isotropization time obtained is shorter, ≈ months, than the time, ≈ years, required by the conventional secondary ENA mechanism to explain the IBEX ENA ribbon.

Low-frequency Waves due to Newborn Interstellar Pickup He+ Observed by the Ulysses Spacecraft

Abstract We have surveyed magnetic field data from the Ulysses spacecraft and found examples of magnetic waves with the expected characteristics that point to excitation by newborn pickup He+. With interstellar neutrals as the likely source for the pickup ions, we have modeled the ion production rates and used them to produce wave excitation rates that we compare to the background turbulence rates. The source ions are thought to be always present, but the waves are seen when growth rates are comparable to or exceed the turbulence rates. With the exception of the fast latitude scans, and unlike the waves excited by newborn interstellar pickup H+, the waves are seen throughout the Ulysses orbit.

Solar Orbiter’s encounter with the tail of comet C/2019 Y4 (ATLAS): Magnetic field draping and cometary pick-up ion waves

Context.Solar Orbiter is expected to have flown close to the tail of comet C/2019 Y4 (ATLAS) during the spacecraft’s first perihelion in June 2020. Models predict a possible crossing of the comet tails by the spacecraft at a distance from the Sun of approximately 0.5 AU.Aims.This study is aimed at identifying possible signatures of the interaction of the solar wind plasma with material released by comet ATLAS, including the detection of draped magnetic field as well as the presence of cometary pick-up ions and of ion-scale waves excited by associated instabilities. This encounter provides us with the first opportunity of addressing such dynamics in the inner Heliosphere and improving our understanding of the plasma interaction between comets and the solar wind.Methods.We analysed data from all in situ instruments on board Solar Orbiter and compared their independent measurements in order to identify and characterize the nature of structures and waves observed in the plasma when the encounter was predicted.Results.We identified a magnetic field structure observed at the start of 4 June, associated with a full magnetic reversal, a local deceleration of the flow and large plasma density, and enhanced dust and energetic ions events. The cross-comparison of all these observations support a possible cometary origin for this structure and suggests the presence of magnetic field draping around some low-field and high-density object. Inside and around this large scale structure, several ion-scale wave-forms are detected that are consistent with small-scale waves and structures generated by cometary pick-up ion instabilities.Conclusions.Solar Orbiter measurements are consistent with the crossing through a magnetic and plasma structure of cometary origin embedded in the ambient solar wind. We suggest that this corresponds to the magnetotail of one of the fragments of comet ATLAS or to a portion of the tail that was previously disconnected and advected past the spacecraft by the solar wind.

Induced Magnetic Fields and Plasma Motions in the Inner Part of the Martian Magnetosphere

AbstractAnalysis of Mars Atmosphere and Volatile Evolution (MAVEN)/Supra‐Thermal And Thermal Ion Composition observations in the Martian upper atmosphere, bounded at higher altitudes by the shocked solar wind, shows that the draping of interplanetary magnetic field penetrates down to low altitudes (∼200−250 km) and governs dynamics of the ionosphere. The upper ionospheric plasma is driven into motion flowing around Mars similar to the shocked solar wind in the adjacent magnetosheath. Such a fluid‐like motion is accompanied by ion acceleration caused by the bending of the magnetic field, leading to ion extraction and finally to ion pickup. Extraction of ions and their acceleration produces a recoil effect of the bulk ionosphere in the opposite direction. This provides a strong asymmetry in ion dynamics in two different hemispheres, accompanied by wrapping of the magnetic field lines around Mars and respective reconnection.

Hybrid simulations of the ring-beam instabilities driven by the pickup ions in the outer heliosheath

ABSTRACT One-dimensional hybrid simulations are performed to investigate the plasma instabilities driven by the pickup ions in the outer heliosheath, which are believed to have a ring-beam velocity distribution in general. The modes propagating parallel and antiparallel to the background magnetic field are studied for pickup ion ring-beam distributions of different pickup angles, following the linear kinetic instability analysis presented in our previous work. The simulation results first confirm the unstable modes predicted by the linear analysis and further reveal how the pickup ions are scattered by the excited waves. While the maximum growth rate of the low-frequency, left-helicity waves occurs at the pickup angle of 82° (consistent with the linear analysis results), the saturation level of the fluctuating magnetic field is the highest around 45° pickup angle. For the pickup ion parameters examined, the pitch-angle spread of the pickup ions at the end of the simulations decreases with increasing pickup angle. Our study goes beyond previous studies of the pickup-ion-driven instabilities in the outer heliosheath which mainly focused on ring-like distributions of around 90° pickup angle. With the caveats that the pickup ion parameters are in the Alfvén cyclotron stability gap and that the one-dimensional simulations only include parallel and antiparallel-propagating modes, our results show that the scattering of pickup ions at small pickup angles is limited to one hemisphere so they can drift away from the ribbon direction as required by the spatial retention scenario of the energetic neutral atom ribbon observed by the Interstellar Boundary EXplorer.

Signature of a Heliotail Organized by the Solar Magnetic Field and the Role of Nonideal Processes in Modeled IBEX ENA Maps: A Comparison of the BU and Moscow MHD Models

Abstract Energetic neutral atom (ENA) models typically require post-processing routines to convert the distributions of plasma and H atoms into ENA maps. Here we investigate how two kinetic-MHD models of the heliosphere (the BU and Moscow models) manifest in modeled ENA maps using the same prescription and how they compare with Interstellar Boundary Explorer (IBEX) observations. Both MHD models treat the solar wind as a single-ion plasma for protons, which include thermal solar wind ions, pick-up ions (PUIs), and electrons. Our ENA prescription partitions the plasma into three distinct ion populations (thermal solar wind, PUIs transmitted and ones energized at the termination shock) and models the populations with Maxwellian distributions. Both kinetic-MHD heliospheric models produce a heliotail with heliosheath plasma that is organized by the solar magnetic field into two distinct north and south columns that become lobes of high mass flux flowing down the heliotail; however, in the BU model, the ISM flows between the two lobes at distances in the heliotail larger than 300 au. While our prescription produces similar ENA maps for the two different plasma and H atom solutions at the IBEX-Hi energy range (0.5–6 keV), the modeled ENA maps require a scaling factor of ∼2 to be in agreement with the data. This problem is present in other ENA models with the Maxwellian approximation of multiple ion species and indicates that either a higher neutral density or some acceleration of PUIs in the heliosheath is required.

The Transport and Evolution of MHD Turbulence throughout the Heliosphere: Models and Observations

A detailed study of solar wind turbulence throughout the heliosphere in both the upwind and downwind directions is presented. We use an incompressible magnetohydrodynamic (MHD) turbulence model that includes the effects of electrons, the separation of turbulence energy into proton and electron heating, the electron heat flux, and Coulomb collisions between protons and electrons. We derive expressions for the turbulence cascade rate corresponding to the energy in forward and backward propagating modes, the fluctuating kinetic and magnetic energy, the normalized cross-helicity, and the normalized residual energy, and calculate the turbulence cascade rate from 0.17 to 75 au in the upwind and downwind directions. Finally, we use the turbulence transport models to derive cosmic ray (CR) parallel and perpendicular mean free paths (mfps) in the upwind and downwind heliocentric directions. We find that turbulence in the upwind and downwind directions is different, in part because of the asymmetric distribution of new born pickup ions in the two directions, which results in the CR mfps being different in the two directions. This is important for models that describe the modulation of cosmic rays by the solar wind.

Future Venus missions and flybys: A collection of possible measurements with mass spectrometers and plasma instruments

This study contains predictions for mass spectrometry and plasma instrument measurements during upcoming Venus flybys of BepiColombo and Solar Orbiter and discusses the possibility of a phosphine detection with mass spectrometry in Venus’ upper atmosphere. The results extend the ones published previously in Gruchola et al. (2019) , where predictions for the proposed ESA mission EnVision and the Venus flyby of JUICE were included. Both the Venus flyby of BepiColombo and the one of Solar Orbiter will take place in August 2021, only 1 day apart. BepiColombo, carrying the neutral mass gas spectrometer STROFIO, could probe the atmosphere around closest approach and obtain data on the thermal and hot neutral particle populations in the upper atmosphere. According to this study, the thermal hydrogen population as well as the hot species H, C, N and O, including isotopes, should be visible to STROFIO. Especially data on the abundance of the hot species would yield important insight into the planetary escape processes. The Solar Orbiter on the other hand carries the plasma instrument SWA-HIS, designed to measure the energetic solar wind ions. During its second Venus gravity assist it will traverse the magnetosheath, the region between bow shock and ion composition boundary, where in addition to the solar wind ions energetic planetary ions are present. The planetary pickup ions can be measured by SWA-HIS, providing important information on the outflow of planetary ions and the ionization processes in the upper atmosphere itself. The recent reports of phosphine in Venus’ cloud decks with an abundance of 20 ppb at 80 km probably overestimate the actual PH 3 abundance, as the data are currently being reanalyzed (Greaves et al., 2021). However, even with this upper limit of 20 ppb a phosphine detection with a mass spectrometer e.g. on-board ESA’s proposed EnVision mission, seems unlikely. To resolve the PH 3 peak and the 16 O 18 O peak a mass resolution of almost 10’000 is required, and the PH 2 D peak is masked by the Cl fragment peak of HCl. Furthermore, NGMS on-board Pioneer Venus with a mass resolution of around 440 did most likely not detect phosphine, as it scanned only a few masses in the mass range of interest where more abundant species than phosphine are present.

Collisional Evolution of the Inner Zodiacal Cloud

Abstract The zodiacal cloud is one of the largest structures in the solar system and strongly governed by meteoroid collisions near the Sun. Collisional erosion occurs throughout the zodiacal cloud, yet it is historically difficult to directly measure and has never been observed for discrete meteoroid streams. After six orbits with Parker Solar Probe (PSP), its dust impact rates are consistent with at least three distinct populations: bound zodiacal dust grains on elliptic orbits (α-meteoroids), unbound β-meteoroids on hyperbolic orbits, and a third population of impactors that may be either direct observations of discrete meteoroid streams or their collisional by-products (“β-streams”). The β-stream from the Geminids meteoroid stream is a favorable candidate for the third impactor population. β-streams of varying intensities are expected to be produced by all meteoroid streams, particularly in the inner solar system, and are a universal phenomenon in all exozodiacal disks. We find the majority of collisional erosion of the zodiacal cloud occurs in the range of 10–20 solar radii and expect this region to also produce the majority of pickup ions due to dust in the inner solar system. A zodiacal erosion rate of at least ∼100 kg s−1 and flux of β-meteoroids at 1 au of (0.4–0.8) × 10−4 m−2 s−1 are found to be consistent with the observed impact rates. The β-meteoroids investigated here are not found to be primarily responsible for the inner source of pickup ions, suggesting nanograins susceptible to electromagnetic forces with radii below ∼50 nm are the inner source of pickup ions. We expect the peak deposited energy flux to PSP due to dust to increase in subsequent orbits, up to 7 times that experienced during its sixth orbit.

Turbulent Acceleration of Interstellar Pickup Ions at the Heliospheric Termination Shock Forms the Global ENA Spectrum

Abstract The heliospheric energetic neutral atom spectrum observed by the Interstellar Boundary Explorer (IBEX) reveals that the heliosheath proton distribution is consistent with a power law. The origin of the spectrum is likely from interstellar pickup ions (PUIs) accelerated at the heliospheric termination shock (HTS). We present an explanation of the proton spectrum origin using a test particle simulation of PUIs accelerated at the HTS. PUIs experience preferential heating by the motional electric field in the shock foot, but do not develop a power-law tail without the presence of turbulence at wavenumbers (k) close to the PUI gyroradius scale (R g). Voyager 2 observations of the magnetic field downstream of the HTS indicate a moderate amount of turbulence at kR g ≅ 1, δ B / B 0 2 ≅ 0.01 , which we find to be sufficient for producing a downstream suprathermal PUI tail but not at intensities observed by IBEX. Within the shock ramp, however, Voyager observed the turbulence power at much smaller scales to be nearly 100 times stronger, suggesting the possibility of strong turbulence at the PUI gyroradius scale. We show that a proton distribution can develop a power law downstream of the HTS consistent with IBEX observations if δ B / B 0 2 ≳ 0.1 at kR g ≅ 1 in the shock foot. Shock drift acceleration of PUIs by the motional electric field is aided by interactions with turbulence upstream of the shock overshoot. Steepening of the IBEX proton spectrum in directions farther from the heliospheric nose suggests the HTS compression ratio and/or turbulence power weakens near the heliotail.

Remote sensing of cometary bow shocks: modelled asymmetric outgassing and pickup ion observations

ABSTRACT Despite the long escort by the ESA Rosetta mission, direct observations of a fully developed bow shock around 67P/Churyumov–Gerasimenko have not been reported. Expanding on our previous work on indirect observations of a shock, we model the large-scale features in cometary pickup ions, and compare the results with the ESA Rosetta Plasma Consortium Ion Composition Analyser ion spectrometer measurements over the pre-perihelion portion of the escort phase. Using our hybrid plasma simulation, an empirical, asymmetric outgassing model for 67P, and varied interplanetary magnetic field (IMF) clock angles, we model the evolution of the large-scale plasma environment. We find that the subsolar bow shock standoff distance is enhanced by asymmetric outgassing with a factor of 2 to 3, reaching up to $18\,000\, \rm {km}$ approaching perihelion. We find that distinct spectral features in simulated pickup ion distributions are present for simulations with shock-like structures, with the details of the spectral features depending on shock standoff distance, heliocentric distance, and IMF configuration. Asymmetric outgassing along with IMF clock angle is found to have a strong effect on the location of the spectral features, while the IMF clock angle causes no significant effect on the bow shock standoff distance. These dependences further complicate the interpretation of the ion observations made by Rosetta. Our data-model comparison shows that the large-scale cometary plasma environment can be probed by remote sensing the pickup ions, at least when the comet’s activity is comparable to that of 67P, and the solar wind parameters are known.

Stability analysis of the pickup ion ring-beam distributions in the outer heliosheath

ABSTRACT The electromagnetic instabilities driven by the pickup ions in the outer heliosheath play a crucial role in interpreting the energetic neutral atom (ENA) ribbon observed by the Interstellar Boundary EXplorer (IBEX). Previous studies on the stability of the outer heliosheath pickup ions focused on ring-like velocity distributions with pickup angles equal or close to 90°. In this study, we investigate the instabilities driven by a ring-beam velocity distribution of the pickup ions in the outer heliosheath using linear instability analysis. The modes propagating along (both parallel and antiparallel to) the background magnetic field are examined, and the emphasis is on how the instabilities change with the pickup angle, α, varying from zero to 90°. The parallel thermal spread of the pickup ions is chosen to make the plasma parameters lie in the Alfvén cyclotron ‘stability gap’ to exclude the well-studied Alfvén cyclotron instability. Our linear instability analysis reveals that the unstable modes all have left helicity but may propagate parallel or antiparallel to the background magnetic field. At small α, the parallel-propagating modes are unstable in two separate wavenumber (and frequency) ranges. The two unstable ranges merge into one when α increases to 82°. In addition, the antiparallel-propagating modes become unstable when α is above 5°, but the parallel-propagating modes have larger growth rates until α becomes larger than αc = 82.8°. The growth rates of the parallel-propagating modes quickly drop below zero when α goes beyond αc, while the antiparallel-propagating modes remain unstable until α approaches 90°.

Multipoint Density Measurements of Geocoronal Pickup Ions

AbstractObservations of the plasma environment outside of the Earth's magnetosphere by the Magnetospheric Multiscale Mission Hot Plasma Composition Analyzers (MMS‐HPCAs) allow density measurements of geocoronal hydrogen pickup ions (PUIs). A study of this environment from December 4, 2019 to December 6, 2019, part of a single orbital period of the MMS spacecraft, revealed an inverse square dependence on PUI density. This dependence corresponds well with an inverse‐cube dependence on the neutral hydrogen density within the geocorona.

Variations in the Pickup Ion Density Structure in Response to the Growth of the Kelvin–Helmholtz Instability along the Heliopause

Abstract Features of the response of pickup ions (PUIs) to the Kelvin–Helmholtz instability (KHI) on the heliopause (HP) are examined by means of two-dimensional hybrid simulations. We assume the supersonic neutral solar wind as the source of PUIs gyrating about the magnetic field in the outer heliosheath. These PUIs become energetic neutral atoms (ENAs) via charge exchange with interstellar hydrogen, and a portion of these ENAs are detected by spacecraft such as the Interstellar Boundary Explorer (IBEX). To evaluate the possibility of identifying the KHI on the HP from ENA observations, we assume that an imprint of the KHI may be displayed in spatial and temporal variations in the observed ENA profile. As an alternative to ENA, the column density of PUIs integrated across the HP is calculated. The KH-inducing vortex forces not only background protons but also PUIs to roll up deep in the inner heliosheath. The KH vortex also results in the emission of magnetosonic pulses that sweep PUIs in the outer heliosheath and lead to their local confinement. These effects elongate the PUIs’ spatial distribution in the direction normal to the HP. The appearance of the local confined structure in the PUIs’ column density is consequently confirmed, and this feature can be confirmed as the KHI evolution. Although the simulation cannot be quantitatively compared with the observations currently available because its resolution is too low, we expect that the derived properties will be useful for diagnosing the nature of HP fluctuation in future missions.

Boundary Conditions at the Heliospheric Termination Shock with Pickup Ions

Abstract In a collisionless shock the directed flow energy is converted mainly in thermal energy of the plasma species. At the termination shock (TS) a substantial portion of energy goes into heating of pickup ions (PUIs), while heating of the solar wind protons (SW) is weaker than it would be without PUI. Heating of both species is nonadiabatic. Downstream pressure of the mixture is determined by conservation laws for the whole mixture. SW heating is sensitive to the details of the shock front, while heating of PUI is not. The profile is an analytical approximation of the observed TS. Here the downstream temperature and pressure of PUIs are obtained for the first time using test particle analysis in a model shock profile for various magnetic compression ratios and shock angles. The profile is an analytical approximation of the observed TS. The results of the analysis are used in the pressure balance equation and the corresponding SW heating is estimated. The analysis is supported by full particle simulations, except for the SW heating that was not studied using test particle data due to its apparent dependence on fine structure of the shock front.

Energetic Neutral Atom Fluxes from the Heliosheath: Constraints from in situ Measurements and Models

Abstract Voyager 2 observations throughout the heliosheath from the termination shock to the heliopause are used to normalize and constrain model pickup ion (PUI) fluxes. Integrating normalized PUI fluxes along the Voyager 2 trajectory through the heliosheath, and combining these integral fluxes with the energy-dependent charge-exchange cross section and the neutral hydrogen density, produces semi-empirical estimates of the energetic neutral atom (ENA) fluxes from the heliosheath. These estimated ENA fluxes are compared with observed ENA fluxes from the Interstellar Boundary Explorer (IBEX) to determine what percentage of the observed fluxes at each IBEX energy are from the heliosheath. These percentages are a maximum of ∼10% for most energies and depend strongly on termination shock properties, plasma density, bulk plasma flow characteristics, the shape of the heliopause, and turbulent energy diffusion in the heliosheath.

MMS Observations of Energized He+ Pickup Ions at Quasiperpendicular Shocks

Abstract Accelerated He+ pickup ions (PUIs) downstream of quasiperpendicular shocks are studied as a function of the fast-mode Mach number (Mf) and shock obliquity (θ Bn). We analyze 10 quasiperpendicular shocks with Mach numbers in the range [1, 7] observed by the Magnetospheric MultiScale (MMS) mission, and compare upstream and downstream He+ velocity distribution functions. For each shock event, we characterize the upstream PUI distribution and derive reduced 1D velocity distributions for the selected upstream and downstream intervals. We also compare the upstream-to-downstream ratio of spectral indices, computed from the He+ perpendicular distributions, to Mf and θ Bn. We find a positive correlation of this spectral index ratio and Mf, which suggests that perpendicular energization of He+ PUIs is enhanced as the shock becomes stronger. These results inform modeling efforts of PUIs and shock-acceleration processes, particularly those taking place at the termination shock.

Strength of the Termination Shock Inferred from the Globally Distributed Energetic Neutral Atom Flux from IBEX

Abstract 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.