22Snakebite envenoming is an important public health problem in Iran, despite its risk not 23 being quantified. This study aims to use venomous snakes' habitat suitability as an 24 indicator of snakebite risk, to identify high-priority areas for snakebite management 25 across the country. Thus, an ensemble approach using five distribution modeling 26 methods: Generalized Boosted Models, Generalized Additive Models, Maximum Entropy 27 Modeling Generalized Linear Models, and Random Forest was applied to produce a 28 spatial snakebite risk model for Iran. To achieve this, four venomous snakes' habitat 29 suitability (Macrovipera lebetina, Echis carinatus, Pseudocerastes persicus and Naja 30 oxiana) were modeled and then multiplied. These medically important snakes are 31 responsible for the most snakebite incidents in Iran. Multiplying habitat suitability 32 models of the four snakes showed that the northeast of Iran (west of Khorasan-e-Razavi 33 province) has the highest snakebite risk in the country. In addition, villages that were at 34 risk of envenoming from the four snakes were identified. Results revealed that 51,112 35 villages are at risk of envenoming from M. lebetina, 30,339 from E. carinatus, 51,657 36 from P. persicus and 12,124 from N. oxiana. This paper demonstrates application of 37 species distribution modeling in public health research and identified potential snakebite 38 risk areas in Iran by using venomous snakes' habitat suitability models as an indicating 39 factor. Results of this study can be used in snakebite and human-snake conflict 40 management in Iran. We recommend increasing public awareness of snakebite 3 41 envenoming and education of local people in areas which identified with the highest 42 snakebite risk. 43 4 62 identifying environmental drivers of species distribution [42−45] and predicting impacts 63 of climate change on biodiversity [46−51]. SDMs are successfully used to identify 64 suitable habitats of species even in areas with no distribution records [52−54]. Thus, 65 these models can be used to identify suitable habitats of venomous snakes as proxies of 66 snakebite risk [12, 55−57] in data poor regions like Iran. 67 The main goal of this paper was to apply SDMs and produce a spatial risk model 68 for snakebite in Iran. To do this, firstly, five distribution modeling methods [38]: 69 Generalized Boosted Models, Generalized Additive Models, Maximum Entropy 70 Modeling Generalized Linear Models, and Random Forest and distribution data of M. 71 lebetina, E. carinatus, P. persicus and N. oxiana were to produce their habitat suitability 72 models. Secondly, the five habitat suitability models [38−58] of each species were 73combined by ensemble approach and finally the four species ensemble models were 74 multiplied to identify potential snakebite risk. We also in addition, number of villages 75 that are at risk of envenoming by these four snakes determined in Iran. 76 77 Materials and methods 78 Occurrence data 79 Distribution records of the M. lebetina, E. carinatus, P. per...