For an additional biodiversity layer, we used an additive benefit

function zonation (Belote et al. 2021). This method took 1697 suitable species ranges (CONUS GAP 2018), each coarsened to ~5km resolution, calculated biodiversity prioritization based on beta diversity. Then using the ABF method, which is a zonation complementarity-based model

that iteratively removes grid cells based on core-area zonation (more heavily rates species richness), it assigned priorities to grid cell locations that best represent species composition (maximizing gamma diversity by prioritizing sites with high complementarity). We reclassified the raster with the top third of values (greater than .794957) = 1, other values = 0, converted to apolygon layer, and extracted values of 1.