Gully erosion is one of the erosive processes that mostly change the shape of the earth surface and has severe environmental and economic damages. The aim of this research is modeling between geo-environmental parameters effective in gully erosion and gully occurrence in the study area and gully erosion susceptibility mapping using evidential belief function (EBF) data driven model in toroud watershed that has high susceptibility to gully erosion. At first, a gully erosion inventory map is prepared, using extensive field surveys and 80 gullies which have been identified, 70 percentage (56 gully location) randomly selected to modeling, while the remaining 30 percentage (24 gully location) are used to validation. In modeling, if there was high correlation among parameters, reduce accuracy of model, thus has done multi-collinearity test among independent variables. Tolerance and the variance inflation factor (VIF) are two important indexes for multi-collinearity diagnosis. Finally 15 parameters including geomorphological, geological, environmental and hydrological are selected for modeling. In evidential belief function model four relationships were calculated: belief (Bel), disbelief (Dis), uncertainty (Unc), and plausibility (Pls) and belief function are used for gully erosion susceptibility mapping. Area under the curve are used for model validation. According to results, EBF model with prediction rate (1) and success rate (0.959) had excellent accuracy and capability in identification of prone areas to gully erosion in study area. The results indicates that 21.79 percentage (90.84 km2) in study area located in high and very high susceptibility class.