Dissertação apresentada ao Programa de Pós-Graduação em Ciências Ambientais da Universidade do Extremo Sul Catarinense - UNESC, como requisito parcial para a obtenção do título de Mestre em Ciências Ambientais. / A perda da biodiversidade é um dos maiores problemas globais. A polinização tem papel chave para a conservação da biodiversidade global, prestando um serviço vital para os ecossistemas. A polinização por abelhas é a mais relevante dentre todas, pois respondem pela polinização da maioria das angiospermas. As abelhas estão em processo acelerado de desaparecimento e, como consequência disso, está causando um grande desequilíbrio ecológico. As mudanças climáticas estão entre as principais pressões globais de declínio das abelhas visto que fornecem as condições gerais para a ocorrência e desempenho de espécies silvestres. As variáveis climáticas, como a temperatura e a precipitação, têm se mostrado como fatores relevantes na interação planta-polinizador. O conhecimento detalhado da ecologia de cada espécie é fundamental para determinar sua potencial distribuição geográfica futura ou até mesmo sua extinção frente às mudanças climáticas. Considerando o cenário de mudanças climáticas futuro, levantamos a hipótese que as duas subespécies de Melipona quadrifasciata Lepeletier, 1836 deslocarão seus nichos climáticos e haverá uma redução espacial na área de ocorrência, acoplada com uma maior sobreposição no futuro. Este trabalho teve por objetivo geral avaliar os efeitos das mudanças climáticas na distribuição das duas subespécies de Melipona quadrifasciata no Brasil nos cenários atual e futuro. Foram utilizadas as subespécies M. quadrifasciata anthidioides com 137 registros e M. quadrifasciata quadrifasciata com 68 registros. As variáveis climáticas atuais e futuras foram extraídas do WorldClim na resolução de 10 arc-min e as variáveis futuras do modelo de circulação acopladas a CCSM4 no cenário de gases de efeito estufa RCP85. Os dados foram modelados no software MaxEnt com inclusão de camada de polarização para redução dos vieses de amostragem, avaliados pela curva ROC e TSS e foi calculada a sobreposição de nicho. Segundo o modelo, no cenário futuro a área de ocorrência potencial da subespécie M. quadrifasciata anthidioides sofre uma redução na região do Cerrado e aumenta na região da Mata Atlântica, além de aumentar sua ocorrência para o sul em direção ao Pampa. A área de ocorrência da subespécie M. quadrifasciata quadrifasciata sofre uma redução em toda sua área de distribuição. Os valores de AUC obtidos excederam 0,7 e o valor de TSS foi igual a 0,85 para todos os modelos. As subespécies mostraram pouca sobreposição de nicho atual (39,9%) e um aumento no futuro (40,6%). Apesar que as duas subespécies terão modificações nas suas distribuições geográficas, para M. quadrifasciata anthidioides o nicho climático poderá ser mais favorável caso esta tenha condições geográficas e ambiente biótico adequadas para migração. Enquanto isso, a subespécie M. quadrifasciata quadrifasciata terá condições climáticas futuras menos favoráveis para manter sua atual distribuição geográfica, sendo sugerido pelos modelos uma maior perda de habitat adequado.