In order to determine the susceptibility of gully erosion at a small watershed scale on the Loess Plateau of China, three hybrid models were developed. These models were based on the Multi-Attributive Border Approximation Area Comparison (MABAC), frequency ratio (FR), CatBoost (CB), LightGBM (LG), and extremely-randomized tree (ET). Based on the Unmanned Aerial Vehicles (UAV) photos, a total of 83 gullies with 12,150 gully pixels and 8 conditioning variables were extracted and used to create the gully inventory database. The correlations between the conditioning parameters and the pixels of gullies were then determined using FR, and the relative importance of these conditioning factors was quantified using machine learning. Then, for gully erosion susceptibility mapping (GESM), three hybrid gully erosion susceptibility models called MABAC-FR-CB, MABAC-FR-LG, and MABAC-FR-ET were developed. The performance of three hybrid models was assessed using the receiver operating characteristic curve (ROC) and the Kappa coefficient. The results claimed that slope steepness greatly influenced the erosion of the gully. The MABAC-FR-ET performed the most precisely, with area under curvature (AUC) of 0.998 and a Kappa of 0.952. As a result, it was determined that MABAC-FR-ET is the most exact and accurate method for predicting the susceptibility to gully erosion in the study watershed.