Authors: E. J. Zirnstein (Princeton University, NJ), B. L. Shrestha (Princeton University, NJ), P. Swaczyna (Space Research Centre PAS (CBK PAN), PL), M. A. Dayeh (Southwest Research Institute, TX), J. Heerikhuisen (University of Waikato, NZ), J. R. Szalay (Princeton University, NJ)

A “heliospheric” termination shock (HTS) surrounds our solar system at approximately 100 astronomical units from the Sun, where the expanding solar wind (SW) is compressed and heated before entering the heliosheath (HS) and encountering the interstellar medium. HTS and HS-accelerated particles govern the pressure balance with the interstellar medium, but little is known about the HTS’s global properties beyond in situ measurements from Voyager in only two directions of the sky. We fill this gap in knowledge with a novel and complex methodology: particle-in-cell, test particle, and MHD simulations, combined with a global minimization scheme to derive global HTS compression ratio sky maps. The methods utilize Interstellar Boundary Explorer observations of energetic neutral atoms produced from HTS and HS-accelerated pickup ions. Our results reveal unique, three-dimensional characteristics of the HTS, such as higher compression near the poles during solar minimum, north-south asymmetries from the disparate polar coronal holes’ evolution, and minimum compression near the flanks likely from SW slowing by enhanced MeV particle pressure.