Gully erosion has significantly increased and severely threatened agricultural production and food security, especially in the black soil region of Northeast China. The present study is set out to select the importance factors in gully erosion occurrence and compare similarities and differences in importance factors among watersheds, validate the applicability of the transformer model, and analyze the spatial distribution characteristic of gully erosion susceptibility maps (GESMs) and relationship between sloping farmland and gully erosion susceptibility area. 25 geo-environmental factors (GEFs) affecting the occurrence of gully erosion were identified through various data resource platforms and Arcgis10.2, and gully inventory maps were generated from remote sensing image interpretation and field survey. The mathematical relationships between GEFs and erosion gullies were established using random forest (RF), convolutional neural network (CNN), and transformer models after multicollinearity test. The 10-fold cross-validation and 8 indicators were used to comprehensively compare the model performances. Results showed that the 10 factors played a key role in gully erosion occurrence in the black soil region of Northeast China which were convergence index (CI), distance from river, rainfall, terrain ruggedness index (TRI), normalized difference vegetation index (NDVI), topographic wetness index (TWI), elevation, distance from road, drainage density, and slope respectively. Except TWI, elevation, and slope, other 7 importance factors are shared among watersheds but with different degrees of importance. The transformer model has better applicability. According to the GESMs, the low susceptibility areas were still dominant and the occurrence of gully erosion was mostly in the very high susceptibility area. Our results demonstrate that the very high susceptibility area was related to the sloping farmland closely and exceeding 75% of very high susceptibility areas were located on sloping farmland.