This study demonstrated the potential of methods derived from geomorphometry and regression models to evaluate landslide susceptibility in a study area located in southern Colombia. From a morphometric stance, the first step was to evaluate the quality of DEM sources by comparison to control points obtained by static-mode GPS. The PALSAR_RTC_hi data was selected for having the best accuracy of heights and was used to derivate terrain parameters at SAGA software. Then, the Principal Component Analysis selected variables with low collinearity, and we classified twelve landforms using fuzzy k-means algorithm, which was compared to a geomorphological map by using the multinomial logistic regression method in R software. We got a Kappa coincidence index of about 30%. The resulting landslide susceptibility mapping took dependent (a mask with unstable-stable cells from an existing landslide inventory) and independent variables (selected morphometric ones). The binary logistic regression showed the propensity of the area to be adversely affected by landslides. This model's performance was tested with a ROC curve over a sample, with 20% of landslide database resulting in an Area Under the Curve of 0,55. This result was contrasted with a spatial prediction model of debris flow, explaining the high frequency of avalanches.