In this study, we focus on the eruption of Mount Etna on Christmas 2018, which emitted great amounts of SO2 from 24th to 30th December into the free troposphere. Simulations based on two different estimations of SO2 emission fluxes are conducted with the chemistry-transport model MOCAGE in order to study the impact of these estimations on the volcanic plume modeling. The two flux emissions used are retrieved (1) from the ground-based network FLAME, located on the flank of the volcano, and (2) from the spaceborne instrument SEVIRI onboard the geostationary satellite MSG. Multiple spaceborne observations, in the infrared and ultraviolet bands, are used to evaluate the model results. Overall, the model results match well with the plume location over the period of the eruption showing the good transport of the volcanic plume by the model, which is linked to the use of a realistic estimation of the altitude of injection of the emissions. However, there are some discrepancies in the plume concentrations of SO2 between the two simulations, which are due to the differences between the two emission flux estimations used that are large on some of the days. These differences are linked to uncertainties in the retrieval methods and observations used to derive SO2 volcanic fluxes. We find that the uncertainties in the satellite-retrieved column of SO2 used for the evaluation of the simulations, linked to the instrument sensitivity and/or the retrieval algorithm, are sometimes nearly as large as the differences between the two simulations. This shows a limitation of the use of satellite retrievals of SO2 concentrations to quantitatively validate modeled volcanic plumes. In the paper, we also discuss approaches to improve the simulation of SO2 concentrations in volcanic plumes through model improvements and also via more advanced methods to more effectively use satellite-derived products.