Susceptibility of intrusion-related landslides in an active volcano was evaluated coupling the landslide susceptibility estimation by random forest (RF), and the probabilistic volcanic vent opening distribution, as proxy for magma injection, using the QVAST tool. In order to develop and test the method proposed here, the RF/QVAST approach was adopted for Stromboli volcano (Southern Italy) since it experienced moderate to huge instability events, it is geomorphologically prone to instability events, and it is affected by active intense volcanic activity that can produce slope instability. The main destabilizing factors of the volcanic flanks are the slope, the aspect, the terrain roughness, the land cover and the litho-technical features of the outcropping rocks. Estimation of volcanic susceptibility shows that the areas with high probability of new vent opening are located in the north-western unstable volcano flank (Sciara del Fuoco), in the volcano summit and the north-eastern volcano flank coherent with the possible reactivation of the eruptive fissures related to the regional tectonic setting. The areas with higher probability of intrusion-related landslides are located in the upper part of the Sciara del Fuoco, while the rest of the island show moderate to low probability of intrusion-related landslide occurrence.