Observed non-gravitational acceleration in the cometary equatorial frame, magnitude (upper panel), and corresponding latitude and longitude of the direction of the sublimation force (lower panel). The latitude ϕ_D is correlated with the anti-solar latitude (− ϕ_sun) (black line in the lower panel), while the longitude correlates with the night terminator λ_D = 180°. The shaded band indicates the variation across the set of initial conditions. The red graph in the upper panel shows the ROSINA-derived water production of 67P from Läuter et al. (2019).