This study updates knowledge on historical geographic distribution of sand fly species through identifying altitudinal and bioclimatic patterns in leishmaniasis endemic areas in Mexico. We analyze and identify sand fly specimens obtained through national efforts by the Instituto de Diagnóstico y Referencia Epidemiológicos (InDRE, Mexico), collected between 1995 and 2001, 2008–2012, and 2017–2023, and add bibliographic information (427 additional records). After a principal components analysis of WorldClim variables and altitudinal ranges, variables that better explain the distribution of sand fly species were chosen (BIO6, BIO12, and BIO16, explaining 72 % of variation). A total of 1,187 specimens of 22 species were retrieved from eight states, providing 29 new municipalities and 48 new localities, being Lutzomyia cruciata, Micropygomyia chiapanensis, and Psathyromyia shannoni the most common species. We presented new historical records of distribution for sand fly species from Morelos (3), Oaxaca (7) and Tabasco (1). The 82.7 % of sand fly species analyzed were distributed in areas with altitudinal ranges below 420 m. The anthropophilic species Psathyromyia shannoni, and Lutzomyia cruciata showed the greatest variability regarding altitudinal range, and climatic preferences, while several wild species showed abiotic preferences. It is likely that the effect of urbanization and climate change generate new beneficial biotopes for the proliferation of the vector sand fly species. Complementary studies that consider seasonality, vegetation types, and change in land use could provide new information to better understand the spread of vector-borne diseases.