Plasma parameters (velocity and magnetic field components) distributions along the upwind direction from the heliospheric model. The vertical dotted lines match the location of discontinuity surfaces along the upwind direction. Although the interstellar plasma flow is supersonic (M ≈ 2), the existence of interstellar bow shock has been questioned so far because of the influence of the interstellar magnetic field (see, e.g. Izmodenov et al. 2005; McComas et al. 2012). However, the heliospheric model used provides the weak shock wave in the LISM.