Abstract. Quantifying time-series of sulphur dioxide (SO2) emissions during explosive eruptions provides insight into volcanic processes, assists in volcanic hazard mitigation, and permits quantification of the climatic impact of major eruptions. While volcanic SO2 is routinely detected from space during eruptions, the retrieval of plume injection height and SO2 flux time-series remains challenging. Here we present a new numerical method based on forward-and backwardtrajectory analyses which enable such time-series to be robustly determined. The method is applied to satellite images of 10 volcanic eruption clouds through the integration of the HYSPLIT software with custom-designed Python routines in a fully automated manner. Plume injection height and SO2 flux time-series are computed with a period of ~10 minutes with low computational cost.Using this technique, we investigated the SO2 emissions from two sub-Plinian eruptions of Calbuco, Chile, produced in April 2015. We found a mean injection height above the vent of ~15 km for the two eruptions, with 15 overshooting tops reaching ~20 km. We calculated a total of 300±46 kt of SO2 released almost equally during both events, with 160±30 kt produced by the first event and 140±35 kt by the second. The retrieved SO2 flux time-series show an intense gas release during the first eruption (average flux of 2560 kt day -1 ), while a lower SO2 flux profile was seen for the second (average flux 560 kt day -1 ), suggesting that the first eruption was richer in SO2. This result is exemplified by plotting SO2 flux against retrieved plume height above the vent, revealing distinct trends for the two events. We propose that a pre-20 erupted exsolved volatile phase was present prior to the first event, which could have led to the necessary overpressure to trigger the eruption. The second eruption, instead, was mainly driven by syneruptive degassing. This hypothesis is supported by melt inclusion measurements of sulfur concentrations in plagioclase phenocrysts and groundmass glass of tephra samples through electron microprobe analysis.This work demonstrates that detailed interpretations of sub-surface magmatic processes during eruptions are 25 possible using satellite SO2 data. Quantitative comparisons of high temporal resolution plume height and SO2 flux timeseries offer a powerful tool to examine processes triggering and controlling eruptions. These novel tools open a new frontier in space-based volcanological research, and will be of great value when applied to remote, poorly monitored volcanoes, and to major eruptions that can have regional and global climate implications through, for example, influencing ozone depletion in the stratosphere and light scattering from stratospheric aerosols. 30Solid Earth Discuss., https://doi