Habitat suitability models to make conservation decisions based on areas of high species richness and endemism

Biodiversity positively relates with the provisioning of ecosystem services and preserving areas with elevated diversity of highly-functional species could help to ensure human well-being. Most studies addressed to make these decisions use maps relying on species occurrences, where sites containing several species are proposed as priority conservation areas. These maps, however, may underestimate species richness because of the incompleteness of occurrence data. To improve this methodology, we propose using habitat suitability models to estimate the potential distribution of species from occurrence data, and later shaping richness maps by overlapping these predicted distribution ranges. We tested this proposal with Mexican oaks because they provide several ecosystem services and habitat suitability models of species were calibrated with MaxEnt. We used linear regressions to compare the outputs of these predictive maps with those of maps based on species occurrences only and, for both mapping methods, we assessed how much surface of sites with elevated richness and endemism of oaks is currently included within nature reserves. Both mapping methods indicated that oak species are concentrated in mountain regions of Mexico, but predictive maps based on habitat suitability models indicated higher oak richness and endemism that maps based on species occurrences only. Our results also indicated that nature reserves cover a small fraction of areas harboring elevated richness and endemism of oaks. These results suggest that estimating richness across extensive geographic regions using habitat suitability models quickly provides accurate information to make conservation decisions for highly-functional species groups.