Scale size of cometary bow shocks

With the upcoming Comet Interceptor mission aiming for a flyby of an hitherto unknown long-period comet, we investigate the expected scale size of the plasma environment to be encountered during this mission. As the target comet is not known, and may not be known before the launch of Comet Interceptor in 2029, we do not know the expected outgassing rate from the nucleus. Therefore, we have no knowledge of the expected scale size of the plasma environment, which can vary by orders of magnitude. Taking the bow shock size as a characteristic size of the plasma environment, we are interested in knowing how this grows with increasing outgassing rate. Previous cometary flyby missions have generated a small statistical dataset of outgassing rates vs. bow shock distances, while computer simulations of the solar wind interaction with various comets have yielded additional datapoints of this. We combine the measured values with a large fraction of these simulations to build up a dataset that spans over four orders of magnitude in both outgassing rate and bow shock distance. The bow shock distances are normalized to the solar wind conditions (400 km/s, 5 cm-3) and ionisation rate (7e-7 s-1) at 1 AU, and also to a flow velocity of 1 km/s of the outgassing neutrals. We then compare this dataset with the gas-dynamic model of Biermann et al., (1967) which was later expanded by Koenders et al., (2013) and find a good model-data agreement. Furthermore, assuming that the bow shock takes the shape of a conic section (as has been found empirically to be the case for most planetary bow shocks) we provide an outgassing rate-dependent bow shock model. This might be useful when planning the operation time-line of Comet Interceptor, or for any other future cometary flyby mission.