Cutaneous leishmaniasis (CL) is caused by a protozoan of the genus Leishmania and is transmitted by sand flies. The state of Espírito Santo (ES), an endemic area in southeast Brazil, has shown a considerably high prevalence in recent decades. Environmental niche modelling (ENM) is a useful tool for predicting potential disease risk. In this study, ENM was applied to sand fly species and CL cases in ES to identify the principal vector and risk areas of the disease. Sand flies were collected in 466 rural localities between 1997 and 2013 using active and passive capture. Insects were identified to the species level, and the localities were georeferenced. Twenty-one bioclimatic variables were selected from WorldClim. Maxent was used to construct models projecting the potential distribution for five Lutzomyia species and CL cases. ENMTools was used to overlap the species and the CL case models. The Kruskal-Wallis test was performed, adopting a 5% significance level. Approximately 250,000 specimens were captured, belonging to 43 species. The area under the curve (AUC) was considered acceptable for all models. The slope was considered relevant to the construction of the models for all the species identified. The overlay test identified Lutzomyia intermedia as the main vector of CL in southeast Brazil. ENM tools enable an analysis of the association among environmental variables, vector distributions and CL cases, which can be used to support epidemiologic and entomological vigilance actions to control the expansion of CL in vulnerable areas.