Abstract The decline of wild bee populations causes the decline of bee-pollinated plant populations through the deterioration of pollination services. Since high bee species richness generally involves high functional group diversity, protecting areas of high bee species richness will help to maintain pollination services for plants. However, those areas do not always include the habitats of bee species with specialized functions that expand the range of plants being pollinated. To map important areas for protecting native bee species and their functions, we estimated the distributions and functional range of 13 bumble bee species and 1 honey bee species in Japan. The distributions were estimated from an ensemble of six species distribution models using bee occurrence data and environmental data. The functional range of bee species was estimated by combining the estimated distributions and proboscis length, which frequently corresponds to the floral shape of the plant species they pollinate. The estimated species richness was high in western Hokkaido and the estimated functional range was wide in central Honshu. Our method is useful to see whether areas important for high species richness of pollinators differ from those for rare species or their functions.