Debris flows can be described as rapid gravity-induced mass movements controlled by topography that are usually triggered as a consequence of storm rainfalls. One of the problems when dealing with debris flow recognition is that the eroded surface is usually very shallow and it can be masked by vegetation or fast weathering as early as one-two years after a landslide has occurred. For this reason, even areas that are highly susceptible to debris flow might suffer of a lack of reliable landslide inventories. However, these inventories are necessary for susceptibility assessment. Model transferability, which is based on calibrating a susceptibility model in a training area in order to predict the distribution of debris flows in a target area, might provide an efficient solution to dealing with this limit. However, when applying a transferability procedure, a key point is the optimal selection of the predictors to be included for calibrating the model in the source area. In this paper, the issue of optimal factor selection is analysed by comparing the predictive performances obtained following three different factor selection criteria. The study includes: i) a test of the similarity between the source and the target areas; ii) the calibration of the susceptibility model in the (training) source area, using different criteria for the selection of the predictors; iii) the validation of the models, both at the source (self-validation, through random partition) and at the target (transferring, through spatial partition) areas. The debris flow susceptibility is evaluated here using binary logistic regression through a R-scripted based procedure.Two separate study areas were selected in the Messina province (southern Italy) in its Ionian (Itala catchment) and Tyrrhenian sides (Saponara catchment), each hit by a severe debris flow event (in 2009 and 2011, respectively).The investigation attested that the best fitting model in the calibration areas resulted poorly performing in predicting the landslides of the test target area. At the same time, the susceptibility models calibrated with an optimal set of covariates in the source area allowed us to produce a robust and accurate prediction image for the debris flows activated in the Saponara catchment in 2011, exploiting only the data known after the Itala-2009 event.