Values for predicted probabilities of avian species occupancy were determined using colonization-extinction models (MacKenzie and others, 2003) as implemented in R (Version 3.4.4; https://www.r-project.org/) via the ‘colext’ function of the Unmarked package (Version 0.12-0; Fiske and Chandler 2011). Performance of a null model (without covariates) and 153 additional models that assessed the effects of geographic coordinates and habitat context covariates were evaluated using Akaike information criteria (AIC; Burnham and Anderson, 2002). When more than one model had substantial support, their respective model weights were used to spatially predict occupancy relative to covariate influence. Predictive model covariates and weights are provided in Appendix 4 (Twedt and Mini, 2020).