Abstract
Global climate change has been identified as a potential driver of observed insect declines, yet in many regions, there are critical data gaps that make it difficult to assess how communities are responding to climate change. Poleward regions are of particular interest because warming is most rapid while biodiversity data are most sparse. Building on recent advances in occupancy modeling of presence‐only data, we reconstructed 50 years (1970–2019) of butterfly occupancy trends in response to rising minimum temperatures in one of the most under‐sampled regions of North America. Among 90 modeled species, we found that cold‐adapted species are far more often in decline compared with their warm‐adapted, more southernly distributed counterparts. Furthermore, in a post hoc analysis using species' traits, we find that species' range‐wide average annual temperature is the only consistent predictor of occupancy changes. Species with warmer ranges were most likely to be increasing in occupancy. This trend results in the majority of butterflies increasing in occupancy probability over the last 50 years. Our results provide the first look at macroscale butterfly biodiversity shifts in high‐latitude North America. These results highlight the potential of leveraging the wealth of presence‐only data, the most abundant source of biodiversity data, for inferring changes in species distributions.