Numerical modeling of tephra dispersal and deposition is essential for evaluation of volcanic hazards. Many models consider reasonable physical approximations in order to reduce computational times, but this may introduce a certain degree of uncertainty in the simulation outputs. The important step of uncertainty quantification is dealt in this paper with respect to a coupled version of a plume model (PLUME‐MoM) and a tephra dispersal model (HYSPLIT). The performances of this model are evaluated through simulations of four past eruptions of different magnitudes and styles from three Andean volcanoes, and the uncertainty is quantified by evaluating the differences between modeled and observed data of plume height (at different time steps above the vent) as well as mass loading and grain size at given stratigraphic sections. Different meteorological data sets were also tested and had a sensible influence on the model outputs. Other results highlight that the model tends to underestimate plume heights while overestimating mass loading values, especially for higher‐magnitude eruptions. Moreover, the advective part of HYSPLIT seems to work more efficiently than the diffusive part. Finally, though the coupled PLUME‐MoM/HYSPLIT model generally is less efficient in reproducing deposit grain sizes, we propose that it may be used for hazard map production for higher‐magnitude eruptions (sub‐Plinian or Plinian) for what concern mass loading.

Volcanic ash clouds are common, often unpredictable, phenomena generated during explosive eruptions. Mainly composed of very fine ash particles, they can be transported in the atmosphere at great distances from the source, having detrimental socio-economic impacts. However, proximal settling processes controlling the proportion (ε) of the very fine ash fraction distally transported in the atmosphere are still poorly understood. Yet, for the past two decades, some operational meteorological agencies have used a default value of ε = 5% as input for forecast models of atmospheric ash cloud concentration. Here we show from combined satellite and field data of sustained eruptions that ε actually varies by two orders of magnitude with respect to the mass eruption rate. Unexpectedly, we demonstrate that the most intense eruptions are in fact the least efficient (with ε = 0.1%) in transporting very fine ash through the atmosphere. This implies that the amount of very fine ash distally transported in the atmosphere is up to 50 times lower than previously anticipated. We explain this finding by the efficiency of collective particle settling in ash-rich clouds which enhance early and en masse fallout of very fine ash. This suggests that proximal sedimentation during powerful eruptions is more controlled by the concentration of ash than by the grain size. This has major consequences for decision-makers in charge of air traffic safety regulation, as well as for the understanding of proximal settling processes. Finally, we propose a new statistical model for predicting the source mass eruption rate with an unprecedentedly low level of uncertainty.