“…All environmental layers were downloaded at high resolution (30 arc‐second, approximately 1 × 1 km) with a WGS84 projection, then converted to projected coordinates using the Mexico ITRF2008 Lambert Conformal Conical projection (900 × 900 m resolution). We chose a resolution of 1 km due to the known limited dispersal in the Arremon species complex (Castaño et al, 2019), which may reflect the importance of landscape features impeding gene flow. Geoprocessing, manipulation, and visualization of spatial data were performed using ArcGIS v.10.4.1 (ESRI), QGIS v.3.16.9 (QGIS Development Team), and the following geospatial packages in R: ‘dismo’ v.1.3‐3 (Hijmans et al, 2020), ‘ggmap’ v.3.0.0 (Kahle & Wickham, 2013), ‘ggplot2’ v.3.3.5 (Wickham, 2016), ‘raster’ v3.4‐13 (Hijmans, 2021), ‘rgdal’ v.1.‐23 (Bivand et al, 2021), ‘rgeos’ v. 0.5‐5 (Bivand & Rundel, 2020), and ‘sf’ v.1.0‐2 (Pebesma, 2018).…”