Jessica E. Shyvers Nathan D. Van Schmidt D. Joanne Saher Julie A. Heinrichs Michael S. O'Donnell Cameron L. Aldridge 20240826 raster digital data https://doi.org/10.5066/P9VBT1ER Jessica E. Shyvers Nathan D. Van Schmidt D. Joanne Saher Julie A. Heinrichs Michael S. O'Donnell Cameron L. Aldridge 2024 publication Habitat restoration efforts to conserve wildlife species are often conducted along a range of local site conditions, with limited information available to gauge relative outcomes for habitat suitability among sites and identify those that may lead to the greatest returns on restoration investment. We leveraged existing resource selection function models to generate heatmaps of spatially varying habitat suitability improvement potential for the Gunnison Sage-grouse (Centrocercus minimus) based on a suite of habitat restoration actions deployed across crucial habitats within six remaining satellite populations. We first simulated expected change in model covariates (habitat features) from a suite of restoration actions (increasing sagebrush, herbaceous, or litter cover, non-sagebrush shrub management, installation of mesic improvement structures, and removal of invasive plants) to generate modified input layers for each. We then reran the original models using these modified layers and calculated the predicted change in habitat suitability across space. The resulting heatmaps identify areas with the greatest improvement potential for each restoration action to help guide strategic restoration planning for the species. This data release, for the Poncha Pass satellite population, includes a set of 13 total raster files. These include: 6 uncategorized heatmaps illustrating predicted change in Gunnison Sage-grouse habitat suitability across space following habitat restoration actions (either single or combined), 6 categorized heatmaps additionally showing areas where 1) new habitat was created, 2) non-habitat remained non-habitat despite management interventions, or 3) negative changes in suitability were observed, and 1 heatmap illustrating predicted changes in suitability following new or worsening plant invasions (cheatgrass, represented by annual herbaceous). Habitat restorations vary by population depending on the reference model. We only ran management action simulations when the reference model had covariates suitable for the simulation (for example, pinyon juniper removal was only run when pinyon juniper was a covariate; See Saher and others (2022) for model details). Raster file names are coded as follows: PPb = Poncha Pass Breeding PPs = Poncha Pass Summer ahrb_rm = annual herbaceous removal (decrease in cover) ahrb_inv = annual herbaceous invasion (increase in cover) combo = combined actions mes_impr = mesic improvements (increase in area) pns_decr = decrease non-sagebrush shrub cover pns_incr = increase non-sagebrush shrub cover C = CATEGORIZED Maps V = UNCATEGORIZED Maps X = INVASION (categorized) Maps The heatmaps generated illustrate spatial variation in the potential for habitat restoration actions, aimed at benefiting Gunnison Sage-grouse, to improve habitat suitability across the Poncha Pass satellite population. They are intended to be used in combination with local knowledge and expertise to identify priority areas for deploying limited management resources, thereby serving as valuable decision-support tools to guide strategic restoration planning efforts for recovery of the species. 20220301 publication date Complete None planned Poncha Pass Gunnison Sage-grouse Population Crucial Habitats -106.07030 -105.97530 38.42300 38.32920 None Gunnison Sage-grouse Centrocercus minimus pinyon-juniper cheatgrass sagebrush mesic improvements USGS Metadata Identifier USGS:6658d316d34ef3137d35fbb0 None Colorado Gunnison Sage-grouse range Sagebrush biome None Centrocercus minimus Kingdom Animalia Subkingdom Bilateria Infrakingdom Deuterostomia Phylum Chordata Subphylum Vertebrata Infraphylum Gnathostomata Superclass Tetrapoda Class Aves Order Galliformes Family Phasianidae Subfamily Tetraoninae Genus Centrocercus Species Centrocercus minimus TSN: 677540 None. Please see 'Distribution Info' for details. None. Users are advised to read the dataset's metadata thoroughly to understand appropriate use and data limitations. Cameron L. Aldridge U.S. Geological Survey, ROCKY MOUNTAIN REGION Supervisory Research Ecologist mailing

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US (970) 266-9100 aldridgec@usgs.gov This dataset, "Maps of habitat suitability improvement potential for the Poncha Pass Gunnison Sage-grouse (Centrocercus minimus) satellite population" is 642KB in size and was built using ArcGIS Pro (version 2.8.0, copyright 2021 Esri Inc.). D. Joanne Saher Michael S. O’Donnell Cameron L. Aldridge Julie A. Heinrichs 2022 publication Global Ecology and Conservation vol. 35 n/a Elsevier BV ppg. e01935 https://doi.org/10.1016/j.gecco.2021.e01935 The values provided in this dataset represent relative predicted change in habitat suitability based on simulated changes in environmental covariates, therefore the relationship to true values is unknown. For details on the original resource selection function models used in development of these layers refer to Saher, D.J., M.S. O’Donnell, C.L. Aldridge, and J.A., Heinrichs. 2022. Balancing model generality and specificity in management-focused habitat selection models for Gunnison Sage-grouse. Global Ecology and Conservation. e01935 and reference maps from the associated Saher and others (2021) data release (https://doi.org/10.5066/P93WFW13). These reference models were generated from a resource selection function model describing seasonal Gunnison sage-grouse habitat. Saher and others (2021) verified the values occurred between 0 and 1, represented relative suitability of habitat. They used modeled results and use locations to identify threshold values where any value less than a threshold reflected a relative suitability of habitat and selected the threshold value as the resource selection function value that captured 95% of known telemetry locations for a given population model, season, and scale. These data are generated from breeding and summer patch resource selection function models developed by Saher and others (2021, 2022) describing seasonal Gunnison sage-grouse habitat. The values represent a relative rank of seasonal habitat suitability and occur between 0 and 1 and are not comparable across or within populations. The relative change in habitat suitability surfaces represent given seasons (life stages [breeding or summer]) at the patch scale (crucial habitats). The predicted surfaces are defined exclusively for the patch-scale ( crucial habitats) resource selection function maps (Saher and others, 2021) that captured 95% of use locations. All covariate data sources reflect 30-meter (m) by 30-m spatial resolution. Not applicable. D. Joanne Saher Michael S. O'Donnell Cameron L. Aldridge Julie A. Heinrichs 2021 raster digital data Saher, J., O'Donnell, M.S., Aldridge, C.L., and Heinrichs, J.A., 2021, Gunnison sage-grouse habitat suitability of six satellite populations in southwestern Colorado: San Miguel, Crawford, Piñon Mesa, Dove Creek, Cerro Summit-Cimarron-Sims, and Poncha Pass: U.S. Geological Survey data release, https://doi.org/10.5066/P93WFW13 https://doi.org/10.5066/P93WFW13 Digital and/or Hardcopy 2021 publication date Reference Habitat Suitability Maps These data represent the reference maps used to calculate change in suitability following modifications to the environmental (covariate) inputs to simulate management actions. They represent the unmodified suitability conditions for each patch model from our analysis. Matthew Rigge Collin Homer Lauren Cleeves Debra K. Meyer Brett Bunde Hua Shi George Xian Spencer Schell Matthew Bobo 20200128 publication Remote Sensing vol. 12, issue 3 n/a MDPI AG ppg. 412 https://doi.org/10.3390/rs12030412 https://www.mrlc.gov/data/type/rcmap-basemap-(2016) Digital and/or Hardcopy 2020 publication date Rangeland Vegetation These data represent environmental covariate inputs (litter cover, sagebrush cover, non-sagebrush shrub cover, shrub height) used in the original reference resource selection function models. We simulated change from habitat management actions by modifying values in these layers and rerunning the reference models to generate our change in habitat suitability maps. Michael J. Falkowski Jeffrey S. Evans David E. Naugle Christian A. Hagen Scott A. Carleton Jeremy D. Maestas Azad Henareh Khalyani Aaron J. Poznanovic Andrew J. Lawrence 201701 publication Rangeland Ecology & Management vol. 70, issue 1 n/a Elsevier BV ppg. 15-24 https://doi.org/10.1016/j.rama.2016.08.002 https://map.sagegrouseinitiative.com/ecosystem/tree-cover Digital and/or Hardcopy 2017 publication date Conifer PJ These data represent environmental covariate inputs for conifer pinyon-juniper (binary presence/absence of 1-4% or 4-10% cover/acre) used in the original reference resource selection function models. We simulated change from habitat management actions by modifying values in these layers and rerunning the reference models to generate our change in habitat suitability maps. U.S. Geological Survey 20170226 raster digital data https://earthexplorer.usgs.gov/ Digital and/or Hardcopy 20150601 20150626 ground condition Landsat (mesic) These data represent environmental covariate inputs for mesic habitats (binary presence/absence) used in the original reference resource selection function models. We simulated change from habitat management actions by modifying values in these layers and rerunning the reference models to generate our change in habitat suitability maps. The Nature Conservancy Teresa Chapman 2022 vector digital data Data either are not available or have limited availability owing to restrictions enforced by The Nature Conservancy. You may request access to the data by contacting Teresa Chapman, The Nature Conservancy. Digital and/or Hardcopy 2021 Date generated Mesic Potential The stream reach polygon dataset delineates areas of potential mesic habitat improvement or restoration. The data were converted to both binary and presence only (1 and no data) rasters and used to simulate potential mesic habitat creation or restoration within the Gunnison Sage-grouse satellite population crucial habitats. We used model equations and covariate inputs from Saher and others (2022), to recreate existing seasonal resource selection function maps representing habitat suitability for the six remaining Gunnison Sage-grouse satellite populations in Colorado: Crawford, Cerro Summit-Cimarron-Sims Mesa, Piñon Mesa, San Miguel, Poncha Pass, and Dove Creek (ArcGIS Pro© [v.2.8.0, copyright 2021 Esri Inc.). Our analysis was restricted to breeding patch-scale (crucial habitats) and summer patch-scale habitats (mesic modification only) maps. The recreated maps represented reference values for calculating relative change in habitat suitability following modifications to environmental covariate inputs (Steps 2-5) Reference Habitat Suitability Maps 20220301 Jessica E Shyvers U.S. Geological Survey, ROCKY MOUNTAIN REGION Ecologist mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9213 jshyvers@usgs.gov We simulated and mapped expected changes associated with each habitat management action by modifying the original 30-m x 30-m geospatial input layers (pre-moving window) and re-running the Saher and others (2022) reference models. Modifications to input layers, simulating change from specific actions, were made within a 1-kilometer (km) buffer of the reference map extents (to encompass all patch-scale moving window extents). Modified covariate layers were generated for three action categories: removal of invasive plants (both native and exotic), vegetative cover modifications (increases or decreases), and mesic habitat improvements and were made in the following manner for each type of action using ArcGIS Pro (v.2.8.0, copyright 2021 Esri Inc.): 1) For removal of invasive plants, we reclassified all pinyon-juniper conifer and annual herbaceous values from the original covariate layers to zero (to simulate total removal). 2) For other vegetation cover modifications (increase percent sagebrush cover, decrease mean non-sagebrush shrub cover, increase percent litter cover), we modified all pixel values from the original covariate layers by plus or minus 1 standard deviation of original geospatial layer values within the 1-km buffered crucial habitat extents and recalculated the moving windows for each model using the Focal Statistics tool in ArcGIS Pro. We used the following modifications (plus or minus 1 standard deviation) to simulate improvements to Gunnison sage-grouse habitats: non-sage shrub (Poncha Pass) = plus 9.6441775004 non-sage shrub (San Miguel) = plus 9.5884763452 litter cover (Piñon Mesa) = plus 9.0567769826 herbaceous cover (San Miguel) = plus 9.6414313335 percent other sagebrush (Dove Creek) = plus 4.1313071418 sagebrush cover (Crawford) = plus 7.0540571074 sagebrush cover (Cerro Summit-Cimarron- Sims Mesa) = plus 6.3300767588 sagebrush height (Piñon Mesa) = plus 13.5135271392 sagebrush height (San Miguel) = plus 17.5334076893 shrub height (Piñon Mesa) = minus 41.3606987755 non-sagebrush shrub (Poncha Pass) = minus 9.6441775004 non-sagebrush shrub (San Miguel) = minus 9.5884763452 3) To simulate mesic improvement actions, we mosaicked all potential mesic habitat expansions (from the "Mesic Potential" source data identified by The Nature Conservancy), to the original mesic input layers and and recalculated the moving windows for each model using the Focal Statistics tool in ArcGIS Pro (v.2.8.0, copyright 2021 Esri Inc.). We recalculated decay functions (euclidean distance tool & raster calculator [Exp(-"EucDist_PM_dec_test"/570)] tools) in ArcGIS Pro). 4) To simulate new or worsening annual herbaceous invasions, we modified all pixel values from the original covariate layers by plus 1 standard deviation of original geospatial layer values within the 1-km buffered crucial habitat extents and recalculated the moving windows for each model using the Focal Statistics tool in ArcGIS Pro. We used the following standard deviation modifications to simulate new or worsening invasions and capped values in the modified layers at 100% cover to remain within logical bounds.: Crawford = 0.5963902232 Poncha Pass = 0.5811664663 San Miguel = 2.8249578108 Cerro Summit-Cimarron-Sims Mesa = 1.9312548554 5) To simulate new or worsening conifer pinyon-juniper encroachment, we modified all pixel values from the original moving window input values used by Saher and others (2022) by plus 1 standard deviation within the 1-km buffered crucial habitat extents. We used the following standard deviation modifications to simulate new or worsening encroachment and capped values in the modified layers at 100% cover to remain within logical bounds: Piñon Mesa = 0.20182 San Miguel = 7.88449 Crawford = 10.46857 Cerro Summit-Cimarron-Sims Mesa = 11.40124 Reference Habitat Suitability Maps Rangeland Vegetation Conifer PJ Landsat (mesic) Mesic Potential 20220301 Jessica E Shyvers U.S. Geological Survey, ROCKY MOUNTAIN REGION Ecologist mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9213 jshyvers@usgs.gov We generated new resource selection function maps by substituting the modified input layers (generated in Step 2) into the reference model equations from Saher and others (2022) and reran the models using the ArcPy package in ArcGIS Pro (v.2.8.0, copyright 2021 Esri Inc.) to generate a comparative resource selection function layer for each action. Modified covariates were assessed one-at-a-time, except for models that included multiple covariates calculated from the same original geospatial input layer (for example, percent mesic cover [moving window value] and distance to mesic [decay function] in Dove Creek), in which case we recalculated and substituted both covariates simultaneously. For combination actions (two management actions conducted simultaneously) we substituted modified input layers for both actions in consideration. 20220301 Jessica E Shyvers U.S. Geological Survey, ROCKY MOUNTAIN REGION Ecologist mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9213 jshyvers@usgs.gov Linear predictor values, associated with each location used in the Saher and others (2022), were converted to ranked percentiles based on the range of potential values for each population. These percentiles provided the reference values by which the reference and modified resource selection function surfaces were compared. Using the percentiles was necessary because the reference models used the exponential form to estimate resource selection on the exponential scale, and therefore unit changes across the dimensionless exponential function did not directly correspond to a probability of selection as they would in a logistic model. We converted values for each reference habitat suitability map (Saher and others, 2021) and our modified maps (generated in Step 3) to ranked percentile values using the linear predictor tables (listing RSF values for each habitat use location) from the Saher and others (2022) models. All RSF values on the landscape were binned using the RSF values for each habitat use location from the table, and then reclassified to the binned percentile using the "ReclassByTable" function in the ArcPy package in ArcGIS Pro (v.2.8.0, copyright 2021 Esri Inc.). The resulting MapRSF values from the reference maps were then subtracted from the MapRSF values in our modified maps to generate uncategorized percentile change maps for each simulated management action (or combination of actions) in each satellite population using the equation: Change in habitat suitability = -1*(MapRSFmodified – MapRSFreference). Calculations were performed using the ArcPy package in ArcGIS Pro© (v.2.8.0, copyright 2021 Esri Inc.) and the resulting Uncategorized Percentile Change Maps. Uncategorized maps for single and combined management actions were visualized using the VALUE.lyrx symbology layers in ArcGIS Pro. These symbology layers are included in this data release for ease of use. Uncategorized maps are indicated by a the “V” suffix in the file name. For these surfaces a value of 0 indicates that there is no change in habitat suitability. 20220301 Jessica E Shyvers U.S. Geological Survey, ROCKY MOUNTAIN REGION Ecologist mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9213 jshyvers@usgs.gov Categorization - For all layers depicting predicted change from single or combined management actions we used nested conditional statements to identify and reclassify pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1, 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. We defined “habitat” as those pixels with habitat suitability ranks exceeding the value that captured 95% of use locations in the reference models; those falling below were considered “non-habitat” Categorized versions of the full value maps depicting these transitions were visualized using CATEGORIZED and INVASION.lyrx symbology layers in ArcGIS Pro, in addition to visualizing negative vs positive percentile change values where these transitions did not occur. These symbology layers are included in this data release for ease of use. Categorized surfaces are indicated by a “C” or “X” suffix in the file name. 20220301 Jessica E Shyvers U.S. Geological Survey, ROCKY MOUNTAIN REGION Ecologist mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9213 jshyvers@usgs.gov Raster Grid Cell 327 241 1 Albers Conical Equal Area 29.5 45.5 -96.0 23.0 0.0 0.0 row and column 30.0 30.0 meters WGS_1984 WGS 84 6378137.0 298.257223563 PPb_ahrb_invX.tif Raster geospatial data file. Producer Defined Value Categorical raster surface showing the change in habitat suitability (breeding season) as a result of increasing annual herbaceous cover (invasion). Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 -1.0 PPb_ahrb_rmC.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPb_ahrb_rmV.tif, showing the change in habitat suitability (breeding season) as a result of the removal of annual herbaceous cover. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPb_ahrb_rmV.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (breeding season) due to the removal of annual herbaceous cover. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.0 0.99194002151489 PPb_combo_C.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPb_combo_V.tif, showing the change in habitat suitability (breeding season) as a result of decreasing annual herbaceous cover and increasing non-sagebrush shrub cover. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPb_combo_V.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (breeding season) due to decreasing annual herbaceous cover and increasing non-sagebrush shrub cover. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.0 0.99194002151489 PPb_mes_imprC.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPb_mes_imprV.tif, showing the change in habitat suitability (breeding season) as a result of increasing the amount of mesic area. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPb_mes_imprV.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (breeding season) due to increasing the amount of mesic area. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.0 0.58064997196198 PPb_pns_decrC.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPb_pns_decrV.tif, showing the change in habitat suitability (breeding season) as a result of decreasing non-sagebrush shrub cover. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPb_pns_decrV.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (breeding season) due to decreasing non-sagebrush shrub cover. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.40323001146317 0.064510002732277 PPb_pns_incrC.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPb_pns_incrV.tif, showing the change in habitat suitability (breeding season) as a result of increasing non-sagebrush shrub cover. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPb_pns_incrV.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (breeding season) due to increasing non-sagebrush shrub cover. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.0 0.91128998994827 PPs_mes_imprC.tif Raster geospatial data file. Producer Defined Value Categorical raster surface, derived from PPs_mes_imprV.tif, showing the change in habitat suitability (summer) as a result of increasing the amount of mesic area. Pixels that transitioned from 1) non-habitat to habitat (habitat created) were coded as value = 1 , 2) non-habitat remaining non-habitat (despite management actions) were coded as value = -2, and 3) for layers depicting predicted change from new or worsening plant invasions (pinyon-juniper or annual herbaceous) we identified habitat to non-habitat transitions (habitat loss, shown as negative change) coded as value =-1. Values between 0 and 1 (habitat that remained habitat after management action) were not categorized, instead they are presented as simply changes in habitat suitability, where values close to 0 represent little change and those closer to 1 indicate the greatest change relative to the reference habitat suitability surface. This was done so that the impact of the management action can be visualized even when there was no change in category. A value of 0 implies no change within the habitat category. Producer Defined -2.0 1.0 PPs_mes_imprV.tif Raster geospatial data file. Producer Defined Value Change in habitat suitability (summer) due to increasing the amount of mesic area. When present, positive values indicate where management action is predicted to improve habitat suitability, negative values indicate areas where a decline in habitat suitability would be expected as a result of management action, relative to the reference habitat suitability surface. A value of 0 indicates that there is no change in habitat suitability. Producer Defined -0.74658000469208 0.56849002838135 Raster file names are coded as follows: PPb = Poncha Pass Breeding PPs = Poncha Pass Summer ahrb_rm = annual herbaceous removal (decrease in cover) ahrb_inv = annual herbaceous invasion (increase in cover) combo = combined actions mes_impr = mesic improvements (increase in area) pns_decr = decrease non-sagebrush shrub cover pns_incr = increase non-sagebrush shrub cover C = CATEGORIZED Maps V = UNCATEGORIZED Maps X = INVASION (categorized) Maps Shyvers, J.E., N.D. Van Schmidt, D.J., Saher, J.A. Heinrichs, M.S. O'Donnell, and C.L. Aldridge. Maps of habitat suitability improvement potential for the Gunnison Sage-grouse (Centrocercus minimus) satellite populations in Southwestern Colorado. U.S. Geological Survey, ScienceBase data release. GS ScienceBase U.S. Geological Survey mailing address

Denver Federal Center, Building 810, Mail Stop 30

Denver CO 80225 United States 1-888-275-8747 970-226-9230 sciencebase@usgs.gov Unless otherwise stated, all data, metadata and related materials are considered to satisfy the quality standards relative to the purpose for which the data were collected. Although these data and associated metadata have been reviewed for accuracy and completeness and approved for release by the U.S. Geological Survey (USGS), no warranty expressed or implied is made regarding the display or utility of the data on any other system or for general or scientific purposes, nor on all computer systems, nor shall the act of distribution constitute any such warranty. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. Digital Data https://doi.org/10.5066/P9VBT1ER None 20240826 Fort Data Management U.S. Geological Survey, Fort Collins Science Center Fort Data Management mailing address

2150 Centre Avenue Bldg C

Fort Collins CO 80526 US 970-226-9100 fortdatamanagement@usgs.gov FGDC Biological Data Profile of the Content Standard for Digital Geospatial Metadata FGDC-STD-001.1-1999