Dreiss et al. (2022) combined eight refugia datasets (current

climate diversity, ecotypic diversity, land facet diversity, landscape diversity, bird macrorefugia, climate dissimilarity, climate velocity, and tree macrorefugia) using a weighted principal components analysis. We took the ecoregion scale composite layer, reclassified to a 0/1 binary using the top third of values (6.7694) and converted 1’s to polygons. Because this layer already represents several facets of ecosystem diversity, we chose to use it as the sole indicator of ecosystem function for our first draft intactness model.