In 2021, more than 50 paroxysmal episodes occurred at the South-East Crater (SEC) of Mt Etna, Italy. The 23–24 March lava fountain was one of the longest episodes and began with weak Strombolian explosions, gradually transitioning to lava fountaining. The eruption intensity then dropped more slowly than in previous episodes, resulting in pulsating Strombolian explosions dominated by ash emission. Thirty-four tephra samples were used to reconstruct the fallout dispersal and estimate the total erupted mass. Grain size, textural, petrological and geochemical analyses indicate different features and were compared with the gas phase (SO2 and HCl) in the volcanic plume. By applying stochastic global optimization to simulations of the temporal evolution of the eruption column height and tephra dispersal and deposition, the total erupted mass retrieved (6.76 × 108 kg) matches well the total erupted mass estimation by the ground-based deposit (8.03 ± 2.38 × 108 kg), reducing the column height throughout the episode from 6.44 to 4.5 km above sea level and resulting in a mass eruption rate ranging from 1.96 × 105 to 8.18 × 103 kg/s. The unusual duration of the March episode and the characteristics of the erupted products point to the change in explosive style and magma fragmentation from fountaining to ash emission phases, associated with a slower magma supply inducing a change in magma rheology and a final, prolonged ash generation. Furthermore, this study showed that using observational data and the variation in eruption source parameters for numerical simulations can improve the accuracy of predicting the dispersal plume, thus mitigating the potential impact of longer paroxysmal episodes.