Habitat and climatic associations of climate-sensitive species along a southern range boundary

Climate change and habitat loss are recognized as important drivers of shifts in wildlife species' geographic distributions. While often considered independently, there is considerable overlap between these drivers, and understanding how they contribute to range shifts can predict future species assemblages and inform effective management. Our objective was to evaluate the impacts of habitat, climatic, and anthropogenic effects on the distributions of climate‐sensitive vertebrates along a southern range boundary in Northern Michigan, USA. We combined multiple sources of occurrence data, including harvest and citizen‐science data, then used hierarchical Bayesian spatial models to determine habitat and climatic associations for four climate‐sensitive vertebrate species (American marten [Martes americana], snowshoe hare [Lepus americanus], ruffed grouse [Bonasa umbellus], and moose [Alces alces]). We used total basal area of at‐risk forest types to represent habitat, and temperature and winter habitat indices to represent climate. Marten associated with upland spruce‐fir and lowland riparian forest types, hares with lowland conifer and aspen‐birch, grouse with lowland riparian hardwoods, and moose with upland spruce‐fir. Species differed in climatic drivers with hares positively associated with cooler annual temperatures, moose with cooler summer temperatures, and grouse with colder winter temperatures. Contrary to expectations, temperature variables outperformed winter habitat indices. Model performance varied greatly among species, as did predicted distributions along the southern edge of the Northwoods region. As multiple species were associated with lowland riparian and upland spruce‐fir habitats, these results provide potential for efficient prioritization of habitat management. Both direct and indirect effects from climate change are likely to impact the distribution of climate‐sensitive species in the future and the use of multiple data types and sources in the modelling of species distributions can result in more accurate predictions resulting in improved management at policy‐relevant scales.