Register      Login
Wildlife Research Wildlife Research Society
Ecology, management and conservation in natural and modified habitats
RESEARCH ARTICLE (Open Access)

The importance of seasonal resource selection when managing a threatened species: targeting conservation actions within critical habitat designations for the Gunnison sage-grouse

M. B. Rice A D E , A. D. Apa B and L. A. Wiechman C D
+ Author Affiliations
- Author Affiliations

A Colorado Parks and Wildlife, 317 W Prospect Avenue, Fort Collins, CO 80526, USA.

B Colorado Parks and Wildlife, 711 Independent Avenue, Grand Junction, CO 81505, USA.

C Department of Fish, Wildlife, and Conservation Biology, Colorado State University, Colorado, Fort Collins, CO 80526, USA.

D Present address: US Fish and Wildlife Service, 1201 Oakridge Drive, Fort Collins, CO 80525, USA.

E Corresponding author. Email: rice1min@gmail.com

Wildlife Research 44(5) 407-417 https://doi.org/10.1071/WR17027
Submitted: 28 July 2016  Accepted: 16 June 2017   Published: 25 September 2017

Journal Compilation © CSIRO 2017 Open Access CC BY-NC-ND

Abstract

Context: The ability to identify priority habitat is critical for species of conservation concern. The designation of critical habitat under the US Endangered Species Act 1973 identifies areas occupied by the species that are important for conservation and may need special management or protection. However, relatively few species’ critical habitats designations incorporate habitat suitability models or seasonal specificity, even when that information exists. Gunnison sage-grouse (GUSG) have declined substantially from their historical range and were listed as threatened by the US Fish and Wildlife Service (USFWS) in November 2014. GUSG are distributed into eight isolated populations in Colorado and Utah, and one population, the Gunnison Basin (GB), has been the focus of much research.

Aims: To provide season-specific resource selection models to improve targeted conservation actions within the designated critical habitat in the GB.

Methods: We utilised radio-telemetry data from GUSG captured and monitored from 2004 to 2010. We were able to estimate resource selection models for the breeding (1 April–15 July) and summer (16 July–30 September) seasons in the GB using vegetation, topographical and anthropogenic variables. We compared the seasonal models with the existing critical habitat to investigate whether the more specific seasonal models helped identify priority habitat for GUSG.

Key results: The predictive surface for the breeding model indicated higher use of large areas of sagebrush, whereas the predictive surface for the summer model predicted use of more diverse habitats. The breeding and summer models (combined) matched the current critical habitat designation 68.5% of the time. We found that although the overall habitat was similar between the critical habitat designation and our combined models, the pattern and configuration of the habitat were very different.

Conclusions: These models highlight areas with favourable environmental variables and spatial juxtaposition to establish priority habitat within the critical habitat designated by USFWS. More seasonally specific resource selection models will assist in identifying specific areas within the critical habitat designation to concentrate habitat improvements, conservation and restoration within the GB.

Implications: This information can be used to provide insight into the patterns of seasonal habitat selection and can identify priority GUSG habitat to incorporate into critical habitat designation for targeted management actions.

Additional keywords: Centrocercus minimus, Colorado, critical habitat, Gunnison sage-grouse, resource selection, species distribution.

Introduction

The conservation of species-at-risk and the habitats they rely upon is a concern for resource managers. The ability to identify these priority habitats is an important and urgent first step in most conservation strategies (Johnson et al. 2004; Fedy et al. 2014). Identifying high quality habitat can be used to support conservation decisions regarding invasive species risk assessment, critical habitat designation, property acquisition and translocation of threatened or endangered species, but the use of these decisions in guiding conservation policies is still scarce (Guisan et al. 2013).

When species are listed under the US Endangered Species Act (ESA) 1973, there is a requirement to evaluate ‘critical habitat’ (Greenwald et al. 2012). Critical habitat is defined by the US Fish and Wildlife Service (USFWS) as the ‘specific areas within the geographic area, occupied by the species at the time it was listed, that contain the physical or biological features that are essential to the conservation of endangered and threatened species and that may need special management or protection’ (USFWS 2015). The designation of critical habitat requires the use of the best scientific data available (Kalen 2014; Murphy and Weiland 2016), but most critical habitat designations are identified based on known locations or general habitat features (Camaclang et al. 2015). Relatively few critical habitat designations are identified using habitat suitability models or knowledge of spatial structure, because it requires additional data (Camaclang et al. 2015). As a result, a critical habitat designation without using habitat suitability models, when the data is available, can lead to broad and sweeping geographic designations that have little utility in siting on-the-ground conservation efforts because they are neither ‘specific’ to the species nor occupied.

Gunnison sage-grouse (Centrocercus minimus; hereafter GUSG) are considered threatened under the ESA and critical habitat has been designated (USFWS 2014a). The distribution of GUSG has declined to an estimated 10.3% of its historical range (Schroeder et al. 2004). Studies have linked the grouse’s decline to habitat loss, degradation and fragmentation of sagebrush landscapes (Oyler-McCance et al. 2001; GSRSC 2005; Bukowski and Baker 2013). Current GUSG were recognised as a separate species from Greater sage-grouse (C. urophasianus), with the primary differences including size, plumage, courtship display and genetics (Young et al. 2000). GUSG distribution is limited to seven isolated populations in south-west Colorado and one in south-eastern Utah (Schroeder et al. 2004). Most of these populations are considered satellite except for one population, the Gunnison Basin (GB), which contains 85–90% of all GUSG rangewide, ~63% of occupied habitat and ~60% of leks (Fig. 1; GSRSC 2005; USFWS 2014b).


Fig. 1.  The Gunnison Basin population in south-west Colorado, USA along with the satellite populations (Crawford, Pinon Mesa, Cerro/Cimmarron/Sims, San Miguel, Dove Creek, and Poncha Pass) and the area in the Gunnison Basin covered by each dataset including the original data (2004–10), validation dataset A (2002), validation dataset B (2010–11), and active lek locations. Validation dataset B is the only dataset that includes the western portion of the Gunnison Basin.
Click to zoom

The current critical habitat designation outlined by the USFWS uses the occupied habitat (called the species area mapping) defined by Colorado Parks and Wildlife (CPW) for the GUSG rangewide conservation plan (GSRSC 2005). Although this was the best available science in 2005, these maps were a rudimentary outline of the geographic range of GUSG based on vegetation cover types that GUSG could occupy, according to professional judgement and observation, and were mostly based on breeding season (i.e. lek locations) occupancy. There is a basic assumption in the range delineation that within the GUSG range all the plant communities are equally valuable; this is a potentially misleading oversimplification that ignores seasonal use differences and the landscape diversity that can influence GUSG movements and habitat use (Connelly et al. 2011). The GSRSC (2005) acknowledged these limitations and developed specific conservation strategies encouraging seasonal habitat mapping efforts to better identify areas to prioritise for protection, because the current occupied range may not reflect species behaviour or movement patterns (GSRSC 2005). Ultimately, the USFWS used this occupied range boundary as the basis for critical habitat designation and refined critical habitat further by excluding existing candidate conservation agreements with assurances, conservation easements and other management agreements with conservation measures applicable to GUSG. These exclusions reduced critical habitat by ~15% (USFWS 2014b), but many of these excluded areas may or may not be important GUSG habitat. Their exclusion from designation is one reason why the current critical habitat designation excludes biologically relevant habitat. In addition, it is recognised that to appropriately manage GUSG populations and habitat, managers must have a better understanding of seasonal habitat use to conserve the important segments of the landscape associated with distinct seasonal habitats throughout their annual cycle, regardless of land ownership (Connelly et al. 2000; Connelly et al. 2011; Fedy et al. 2014). The critical habitat designation by the USFWS acknowledges local scale seasonal habitat structure for GUSG, but does not (and cannot) use that information in designating critical habitat (USFWS 2014b), thereby ignoring the seasonal movement patterns of GUSG on the landscape.

The USFWS designates ~84% of the GB as critical habitat, but there is no differentiation of spatial or temporal variability in quality and use. For over a decade, fine scale models based on detailed patterns of animal use have shown potential to identify crucial habitat not obvious at larger scales (Klar et al. 2008). These more detailed seasonal models are needed for identifying restoration and rehabilitation of areas important to the recovery and viability of GUSG (Wisdom et al. 2011). For example, previous models (Aldridge et al. 2012) have identified nesting habitat for GUSG and suggest that ~50% of the GB is important nesting habitat – less than the current designated critical habitat (USFWS 2014b). In an effort to define seasonal habitat selection patterns, we collected radio-telemetry data from over 200 GUSG from 2004 to 2010, along with two independent validation datasets that could be used to model GUSG habitat selection in a more robust fashion specific to the GB during two seasons (breeding, 1 April–15 July; and summer, 16 July–30 September) and covering those 6 months of the year.

We had two main objectives for this study. First, we developed habitat selection models for the breeding and summer seasons using GUSG radio-telemetry data collected in the GB population from 2004 to 2010. Second, we compared the specific seasonal models with the critical habitat designated during the ESA threatened listing process. Because of the broad generalisation of critical habitat designation by the USFWS, we believe that our analysis will assist the USFWS in the development of a recovery plan for GUSG to more specifically focus on conservation actions in the most important areas within the GB.


Materials and methods

Study area

Our study area was the Gunnison Basin population of GUSG in Gunnison and Saguache counties, Colorado, USA. This area comprises 85–90% of all GUSG (GSRSC 2005; Fig. 1). The Gunnison Basin is a 200 km2 intermontane basin ranging in elevation from 2300 to 2900 m (Hupp and Braun 1989). Mountainous terrain areas border the north, east and south-east sides of this population and are not commonly used by GUSG.

Data

We captured GUSG from March to early May from 2004 to 2010, using spotlighting techniques (Giesen et al. 1982; Wakkinen et al. 1992). We fit GUSG with a 16 or 17 g necklace-style VHF radio-transmitter (model A4050 by Advanced Telemetry systems, Isanti, MN, USA or model R12B by Holohil Systems, Ontario, Canada) equipped with a 4-h mortality sensor, and also fit each bird with a uniquely numbered leg band (National Band and Tag Co., Newport, KY USA). The transmitter was <2% of an average GUSG female (1270 g, s.d. 90 g) or male (2110 g, s.d. 190 g). Trapping and handling protocols were approved by the Colorado Parks and Wildlife Animal Care and Use Committee (permit # 02-2005).

Following release, we located radio-marked individuals on the ground using hand-held Yagi antennas once every 1–3 days (from date of capture through September) to monitor status (dead or alive) and movement patterns. Monthly flights were conducted during the winter to assess survival, but due to logistical restraints, on the ground telemetry locations were not collected. We used triangulation to estimate locations using maximum likelihood estimates generated in program LOCATE II (Nams VO (1990) Locate II. Pacer, Truro: Nova Scotia, Canada).

Locations were assigned to one of two seasons: breeding or summer. The cutoff dates were based on the Gunnison sage-grouse Conservation Plan, which guides GUSG management in Colorado (GSRSC 2005). We used a resource selection framework comparing used vs available locations, so we randomly generated a sample of ‘available’ locations within the same geographical extent of the GB (Stokland et al. 2011). We conducted a sensitivity analysis of the available sample size for each season (Northrup et al. 2013) and found that coefficients converged at n = 9000 available samples in both seasons. Average daily movement distances were estimated across all the birds in both seasons and used to buffer all presence and available locations within that season (180.5 m in the breeding season and 223.0 m in the summer season). Due to terrain constraints and access issues, there were a few telemetry errors that were larger than this average daily movement. Therefore, we used this buffer to remove those locations and summarised all habitat variables at a biologically relevant scale.

We classified vegetation type into eight biologically meaningful categories (available as Supplementary Material 1 on the journal website) using the Colorado basinwide vegetation layer (i.e. sagebrush (Artemisia spp.), forest, etc.). This land cover layer was constructed in 2005 from 25-m resolution landsat imagery as part of the Colorado Vegetation Classification Project administered by CPW in collaboration with the Bureau of Land Management and the USDA Forest Service. We believe that the vegetation cover type was consistent and reflected current conditions in the GB as there were few relevant environmental perturbations (e.g. wild or managed fire) before or during our study. From the eight categories, we excluded those that consisted of <0.1% of the cover within the buffers (e.g. urban, riparian and water). Thus, the vegetation categories we used in model development were irrigated agriculture (pastureland), sagebrush, grassland, bare and forest, which have all been shown to be important predictors of habitat use in previous studies on both Greater sage-grouse and GUSG (Rice et al. 2013; Stanley et al. 2015; Walker et al. 2016). In addition, these vegetation categories have been cited as influential in multiple seasons for GUSG habitat and are often the focus of management actions in Colorado (GSRSC 2005).

We obtained elevation data from the USA Geological Survey (USGS) digital elevation model and used the national hydrography dataset to measure density of water bodies (Dzialak et al. 2012). We used the National Wetlands Inventory layer from USFWS to measure the distance to wetlands. Development variables included roads, easements and address points (hereafter referred to as residential) collected by Gunnison County. Using these layers, we calculated distance to highways (paved roads only), distance to residential, distance to easements, residential density and road density (includes paved and unpaved roads). We also measured the distance to sagebrush, because observations of GRSG feeding on the edges of vegetation cover near the ecotone with sagebrush have been recorded, particularly in the summer season (Connelly et al. 2011). A list of the variables summarised by seasonal buffer is in Table 1.


Table 1.  Variables used in Gunnison sage-grouse habitat models with the mean value and standard errors for presence and available buffers in the breeding (1 April–15 July) and summer (16 July–30 September) seasons in Gunnison Basin, Colorado, USA (2004–10)
Click to zoom

Model building by season

We first calculated Pearson correlation coefficients (r) for all the habitat variables (r > 0.70; McGarigal et al. 2000) to remove highly correlated variables that may cause multicollinearity issues in the model and removed the variable with a larger P-value. Table 1 indicates the final variables used in each season after removing correlated variables. We fit our used and available location data with a binomial generalised linear model using the package ‘lme4’ in program R (R Core Team 2013; Bates et al. 2013). We used a random intercept for each individual grouse within each season to account for unbalanced sampling among animals (Gillies et al. 2006).

We constructed a set of alternative models from all linear combinations of the habitat variables in each season (McAlpine et al. 2008). We generated predictions from each of the best-fitting models within the 95% confidence set and averaged the predictions into a final model for each season to strengthen inference (Burnham and Anderson 2002). The unstandardised coefficients for each of the seasonal models are equivalent to selection ratios (Manly et al. 2002) and exp(βi) can be interpreted directly as the odds ratios. We also used standardised coefficients (averaged across the same set of models) to assess the relative effects of different covariates measured at different scales (Schielzeth 2010). To do this we centred and scaled all variables in each season (mean = 0, s.d. = 1;McAlpine et al. 2008).

The model weights in the 95% set were recalculated to sum to 1 and each model within the set was used to create a model prediction surface (Anderson 2008). Each prediction surface within the model set was multiplied by its weight and then added together to produce a final model averaged prediction surface (Aldridge et al. 2012). This process allows all plausible models in a set to be used in multimodel inferences for spatial predictions (Anderson 2008).

We used the final model averaged prediction surface for each season to create a prediction surface in ArcMap 10.1 based on the associated habitat variables (ArcGIS 10.1; Environmental Research Systems Institute, Redlands, CA). We applied the logistic equation in the following form to create the relative probability of GUSG presence across the Gunnison Basin:

E1

where xn are covariates at location i for bird j with fixed effect regression coefficients βn, and β0 is the mean intercept. The logistic function was used to create a relative probability of presence surface with values between 0 and 1 across the GB for each season (1 = high, 0 = low).

In order to compare our models with the critical habitat, we completed an error matrix for the breeding model, the summer model and the combination of both seasonal models to critical habitat. We used prediction values of >50% relative probability of presence to define occupied habitat and all relative probability of presence values of <49% as unoccupied habitat. We used overall accuracy and the true skill statistic (TSS) as measures of comparison accuracy (Allouche et al. 2006). The TSS ranges from –1 to 1, where 1 equals perfect agreement and values of 0 or less indicate a performance no better than random (Allouche et al. 2006).

Model validation

The gold standard for any model is to test model predictions against independently collected data (Wiens et al. 2008). We used two independent telemetry datasets for validation of the breeding and summer seasonal models as well as lek location data in the breeding season. We captured and marked GUSG for both independent datasets using the same techniques described earlier. The first dataset was collected in 2002 (validation set A) and the second dataset was collected during 2010–2011 (validation set B). In both studies, we circled each bird at a <30 m radius and manually corrected the Universal Transverse Mercator (UTM) location to reduce observation error. In addition, because it was important to test predictions in areas where data had not been collected for the original model development, and to verify that the model results could extrapolate across the GB, validation set B was collected in the western portion of GB where our seasonal models were extrapolated (Fig. 1). We separated each validation dataset into seasons based on the original models (breeding: 1 April–July 15; summer: 16 July–30 September). For the validation datasets, we calculated the proportion of locations or leks for which the model predicted >50% relative probability of presence for each seasonal prediction model (Sawyer et al. 2007).


Results

We collected 7643 locations between April 2004 and September 2010, from 210 radio-marked GUSG. There were 188 GUSG and 3922 locations for the breeding season and 171 GUSG and 3721 locations for the summer season. Movement buffers for GUSG in the GB tended to be smaller during the breeding season (180.5 m) than during the summer (223.0 m). GUSG locations in both the breeding and summer seasons occurred in areas with higher proportions of sagebrush, higher water density, proximity to sagebrush, lower road density, lower residential density and further from residential than the available locations (Table 1).

Our breeding model suggests that GUSG locations were 21 times more likely to be within sagebrush and 16 times more likely to be within grassland than available locations (Table 2, Fig. 2a). Additionally, GUSG locations were five times more likely to be further from wetlands than available locations, but two times as likely to be in areas of high water density. We also detected smaller effects where GUSG were closer to easements, further from residential development and highways, occupying slightly lower elevations and within lower residential densities (Table 2).


Table 2.  Top model standardised and unstandardised coefficients for Gunnison sage-grouse resource selection in Gunnison Basin, Colorado, USA (2004–10) in the breeding season (1 April–15 July) including the standard error (s.e.), lower (LCI) and upper (UCI) 95% confidence intervals and odds ratios
Click to zoom


Fig. 2.  (a) Relative probability of Gunnison sage-grouse use in Gunnison Basin, Colorado, USA (2004–12) during the breeding season (1 April–15 July). (b) Relative probability of Gunnison sage-grouse use in the Gunnison Basin, Colorado, USA (2004–12) during the summer season (16 July–30 September).
Click to zoom

Our summer model indicated that GUSG locations were 20 times as likely to be in sagebrush, 57 times as likely to be in irrigated agriculture, and 12 times as likely to be in bare cover types compared with available locations, indicating a diverse use of habitat types (Table 3, Fig. 2b). Smaller effects included GUSG being further from wetlands and residential development, closer to highways and in areas with low road density and high water density during the summer season.


Table 3.  Top model standardised and unstandardised coefficients for Gunnison sage-grouse resource selection in Gunnison Basin, Colorado, USA (2004–10) in the summer season (16 July–30 September) including the standard error (s.e.), lower (LCI) and upper (UCI) 95% confidence intervals and odds ratios
Click to zoom

Both of our models validated well with the independent datasets with slightly better results for the breeding model. Our validation set A had 507 locations (breeding, n = 262; summer, n = 245) and validation set B had 1600 locations (breeding, n = 762; summer, n = 838). The percentages of locations with a relative probability >0.50 for the breeding model were 83.0% for validation set A, 84.0% for validation set B and 85% for leks. The percentages of locations with a relative probability >0.50 for the summer model were 75.0% for validation set A and 79% for validation set B.

The current critical habitat designation without exclusions was approximately the same as the rangewide map provided by Colorado Parks and Wildlife, but was reduced by ~15% with exclusions removed (Table 4). For comparison, we estimated the area of occupied habitat for the breeding and summer models by summing the predictive surface with values >0.50 relative probability of presence threshold, assuming that values above 0.50 represent the habitat most likely to support GUSG (Table 4). Although the critical habitat is not divided by seasons, the area encompassed by the breeding and summer models was 27.9% and 17.9% smaller, respectively, compared with critical habitat (Table 4). Because we created two models based on habitat use during the breeding and summer seasons, whereas critical habitat represents overall habitat (regardless of seasonal use), direct comparison between the current critical habitat and our models is difficult. However, we intersected the breeding and summer seasonal models and used this as a conservative way to look at overall potential differences compared with critical habitat. By intersecting the breeding and summer models, the area was 4.8% smaller compared with critical habitat (Table 4), however, this approach devalues the benefit and importance of using the seasonally specific models, given the pattern and configuration of the occupied and unoccupied habitat is quite different to the critical habitat. Overall, the occupied and unoccupied habitat in the critical habitat map agreed with the combined breeding and summer models 68.5% of the time (Supplementary Material 2). Therefore, 31.5% of the cells in the critical habitat designation did not align with our combined models (Supplementary Material 2). The breeding and summer models agreed with the critical habitat in 56.4% and 58.4% of the cells, respectively. We found that our combined model predicted unoccupied habitat in areas the critical habitat predicted occupied in 20.9% of the cells (Supplementary Material 2). In contrast, we found our combined model predicted occupied habitat where the critical habitat predicated unoccupied habitat in 10.6% of the cells. Therefore, our models identified areas in the unoccupied critical habitat designation that were predicted to be relatively better quality habitat and areas in the occupied critical habitat predicted to be relatively poorer in quality (Fig. 3). This illustrates that although the overall amount of habitat was similar between the critical habitat designation and our models, the pattern and configuration of that habitat were very different.


Table 4.  Approximate area in hectares estimated by the seasonal models using the 50% relatively probability cutoff compared with both the Colorado Parks and Wildlife (CPW) occupied range map for Gunnison sage-grouse and the critical habitat designated by USFWS in 2014 (raster values have been rounded so totals may not be exactly the same)
T4


Fig. 3.  Occupied and potential critical habitat designated by US Fish and Wildlife Service in 2014 compared with the breeding (1 April–15 July), summer (16 July–30 September) and combined seasonal model predictions (occupied habitat is based on >0.50 relative probability cutoff value).
Click to zoom


Discussion

Our breeding season model confirms similar habitat patterns found in the only other landscape modelling effort for GUSG that predicted GUSG use in higher proportions of sagebrush, further from roads and lower residential densities (Aldridge et al. 2012). Our analysis provides the first predictive summer season model for GUSG. Our results show that GUSG are dependent on large patches of sagebrush during both the breeding and summer seasons, although they tend to have more diverse habitat selection during the summer when they utilise numerous habitats, including irrigated agriculture and bare ground. Use of smaller sagebrush openings in a diverse suite of cover types has been found within Colorado with GRSG (Hausleitner 2003; Thompson 2012). In addition, use of agriculture is not uncommon during the summer season and has been documented in the literature (Connelly et al. 2011) and by many residents in the GB (Knapp et al. 2013).

The USFWS (2014b) recognises that local scale habitat structure and quality vary spatially, with some areas providing habitat for one or more seasons, but to date, critical habitat has not been defined on a seasonal basis. Although the USFWS (2014b) identified fine scale habitat structure attributes, they acknowledged that the data to delineate seasonal habitats at a landscape level did not exist at the time of designation (USFWS 2014b). We suggest that our seasonal models provide an analysis of a long-term dataset that the USFWS can use in addition to the estimated critical habitat within the GB to address these information gaps while establishing recovery goals and siting conservation actions. Designation of critical habitat does not lend itself to seasonal models, because any given site is usually classified as either critical habitat or not. We are not proposing to change how critical habitat is designated, but rather to include more data-driven models when they are available as a second tier of designation. In addition, the spatial arrangement and relative probability values provide a more accurate, data-driven assessment of the environmental variables that influence GUSG habitat use within the GB. This information provides a much stronger ability to identify on the ground conservation measures for species conservation when the information is available and can assist in targeted recovery of the species.

There are multiple reasons why multi-scale spatially explicit predictive surfaces could help conservation-reliant species such as GUSG (Scott et al. 2010). First, the critical habitat designation does not explicitly address the environmental variables important to GUSG. The incorporation of data-driven relationships with environmental variables allows management agencies to direct their management efforts on these variables in a recovery effort. For example, the summer model indicates that multiple habitat types are important to GUSG; this might guide management actions. If restoration and rehabilitation of habitat is a key component to GUSG recovery (Wisdom et al. 2011), management needs to identify finer scale priorities for habitat management. Accomplishing meaningful habitat protection requires more than just a large spatial area, but rather broad-scale land use planning (Kalen 2014), which our models can support.

Second, our analysis provides the most current and best science available on GUSG resource selection in the GB. We agree with the recommendation by Murphy and Weiland (2016) to use contemporary data that reflect spatial and temporal patterns in resource use by the species as one way to meet the best available science directive in the ESA (section 4(b) (2)) as well as ensuring data is current, quantitative and documented (Clark et al. 2002). This analysis attempts to fulfil these recommendations because it directly incorporates telemetry data over a 6-year period throughout two seasons of the GUSG lifecycle.

Finally, there is value in evaluating critical habitat at multiple spatial scales (Shirk et al. 2014) because species may respond differently at larger or smaller scales (Cunningham et al. 2014). Consideration of scale is necessary for deciding how habitat data should be applied in resource management (Boyce et al. 2003). Thus, we suggest incorporating both the general critical habitat outlined by the USFWS and also our finer scale models to focus and prioritise habitat restoration and recovery efforts. We caution managers that drawing conclusions about habitat selection based on observations at any one scale may misconstrue the importance of variables driving system behaviour overall (Doherty et al. 2010). Using a hierarchical approach to conservation planning is not novel (Johnson et al. 2004), but we are unaware of any use of the application for critical habitat designation, nor the use in policy actions related to ESA in the USA.

Our models provide important insights into habitat use patterns of GUSG in the breeding and summer seasons, but information illustrating winter habitat use patterns is still lacking. We recommend allocating future resources and research on collecting data during the winter season. The GB population has been the focus of much research (Aldridge et al. 2012; Davis et al. 2014; Davis et al. 2015), so we caution applying our models to the other seven satellite populations (<10–15% rangewide population; Fig. 1), where data is more scarce and local seasonal use patterns may differ. More monitoring and research could be directed towards the satellite populations to improve these data gaps.

Gunnison Basin stakeholders are keenly aware of the seasonal habitat use patterns of GUSG (Knapp et al. 2013) and have the perception, real or assumed, about the restrictions that can be imposed by critical habitat designation on the public and private sectors (Salzman 1990). Because of these concerns, stakeholders assisted us in selecting variables for our model development; they are invested and engaged in the recovery of GUSG. Gunnison County has filed an intent to sue the federal government for numerous issues, such as failing to use the best available science and issues related to the critical habitat designation (Gunnison County 2014), including failure to demonstrate that critical habitat designation is essential to conserving the species, designating habitat that was not beneficial to the species and designating critical habitat based on flawed and unsubstantiated conjecture regarding historic range (Gunnison County 2014). Our analysis provides an opportunity for the USFWS to use a more detailed, data-driven analysis when identifying where to target conservation actions and inform recovery goals. We do not believe that our models can replace the current critical habitat designation because the intent of critical habitat incorporates economic valuation of the landscape and pre-existing conservation easements by law. However, by utilising these models in the recovery planning process, there’s an opportunity to do so with a secondary product that has stakeholder support, which could lead to more successful, collaborative, strategically sited conservation efforts. We recommend that the USFWS take an adaptive approach to conserving this population by using our models as new information combined with stakeholder input to guide conservation efforts going forward. Using a hierarchical approach to GUSG recovery – by including our models to strategically identify seasonal areas to concentrate habitat improvements, conservation and restoration within the GB – may assist in gaining local support in GUSG conservation and recovery to accelerate the delisting process. This hierarchical approach for GUSG could then be used by USFWS on other species that have more detailed data available to inform critical habitat designations.


Conflicts of interest

The authors declare no conflicts of interest.



Acknowledgements

We would like to thank Colorado Parks and Wildlife for funding this study as well as CPW personnel including M. Phillips, J. Wenum, N. Seward and P. Street. We thank J. Cochran, M. Pelletier and the Gunnison sage-grouse Strategic Committee for their support and input. We thank the technicians who assisted on this project as well as the landowners who allowed access to their properties. We also thank J. Ivan and D. Bilyeu for comments on earlier drafts of this manuscript. The findings and conclusions in this manuscript are those of the authors and do not necessarily represent the views of the US Fish and Wildlife Service.


References

Aldridge, C. L., Saher, D. J., Childers, T. M., Stahlnecker, K. E., and Bowen, Z. H. (2012). Crucial nesting habitat for Gunnison sage-grouse: a spatially explicit hierarchical approach. The Journal of Wildlife Management 76, 391–406.
Crucial nesting habitat for Gunnison sage-grouse: a spatially explicit hierarchical approach.Crossref | GoogleScholarGoogle Scholar |

Allouche, O., Tsoar, A., and Kadmon, R. (2006). Assessing the accuracy of species distribution models: prevalence, kappa, and the true skill statistic (TSS). Journal of Applied Ecology 43, 1223–1232.
Assessing the accuracy of species distribution models: prevalence, kappa, and the true skill statistic (TSS).Crossref | GoogleScholarGoogle Scholar |

Anderson, D. R. (2008). ‘Model Based Inference in the Life Sciences: a Primer on Evidence.’ (Springer Science+Business Media, LLC: New York.)

Bates, D., Maechler, M., Bolker, B., and Walker, S. (2013). Lme4: linear mixed effects models using Eigen and S4. R package version 1.0–4. Available at http://CRAN.R-project.org/package=lme4 [Accessed 01 June 2014].

Boyce, M. S., Turner, M. G., Fryxell, J., and Turchin, P. (2003). Scale and heterogeneity in habitat selection by elk in Yellowstone National Park. Ecoscience 10, 421–431.
Scale and heterogeneity in habitat selection by elk in Yellowstone National Park.Crossref | GoogleScholarGoogle Scholar |

Bukowski, B. E., and Baker, W. L. (2013). Historical fire in sagebrush landscapes of Gunnison sage-grouse range from land survey records. Journal of Arid Environments 98, 1–9.
Historical fire in sagebrush landscapes of Gunnison sage-grouse range from land survey records.Crossref | GoogleScholarGoogle Scholar |

Burnham, K. P., and Anderson, D. R. (2002). ‘Model Selection and Multimodel Inference: Practical Information-Theoretic Approach.’ (Springer-Verlag: New York).

Camaclang, A. E., Maron, M., Martin, T. G., and Possingham, H. P. (2015). Current practices in the identification of critical habitat for threatened species. Conservation Biology 29, 482–492.
Current practices in the identification of critical habitat for threatened species.Crossref | GoogleScholarGoogle Scholar |

Clark, J. A., Hoekstra, J. M., Boersma, P. D., and Kareiva, P. (2002). Improving U.S. endangered species act recovery plans: key findings and recommendations of the SCB recovery plan project. Conservation Biology 16, 1510–1519.
Improving U.S. endangered species act recovery plans: key findings and recommendations of the SCB recovery plan project.Crossref | GoogleScholarGoogle Scholar |

Connelly, J. W., Schroeder, M. A., Sands, A. R., and Braun, C. E. (2000). Guidelines to manage sage grouse populations and their habitats. Wildlife Society Bulletin 28, 967–985.

Connelly, J. W., Rinkes, E. T., and Braun, C. E. (2011). Characteristics of greater sage-grouse habitats: a landscape species at micro- and macro-scales. In ‘Studies in Avian Biology; Greater Sage-grouse: Ecology and Conservation of a Landscape Species and its Habitats’. (Eds S. T. Knick and J. W. Connelly.) pp. 69–83. (University of California: Berkeley, CA.)

Cunningham, R., Lindenmayer, D., Barton, P., Ikin, K., Crane, M., Michael, D., Okada, S., Gibbons, P., and Stein, J. (2014). Cross-sectional and temporal relationships between bird occupancy and vegetation cover at multiple spatial scales. Ecological Applications 24, 1275–1288.
Cross-sectional and temporal relationships between bird occupancy and vegetation cover at multiple spatial scales.Crossref | GoogleScholarGoogle Scholar |

Davis, A. J., Hooten, M. B., Phillips, M. L., and Doherty, P. F. (2014). An integrated modeling approach to estimating Gunnison sage-grouse population dynamics: combining index and demographic data. Ecology and Evolution 4, 4247–4257.

Davis, A. J., Phillips, M. L., and Doherty, P. F. (2015). Nest success of Gunnison sage-grouse in Colorado, USA. PLoS One 10, e0136310.

Doherty, K. E., Naugle, D. E., and Walker, B. L. (2010). Greater sage-grouse nesting habitat: the importance of managing at multiple scales. The Journal of Wildlife Management 74, 1544–1553.
Greater sage-grouse nesting habitat: the importance of managing at multiple scales.Crossref | GoogleScholarGoogle Scholar |

Dzialak, M. R., Olson, C. V., Harju, S. M., Webb, S. L., and Winstead, J. B. (2012). Temporal and hierarchical spatial components of animal occurrence: conserving seasonal habitat for greater sage-grouse. Ecosphere 3, 30.
Temporal and hierarchical spatial components of animal occurrence: conserving seasonal habitat for greater sage-grouse.Crossref | GoogleScholarGoogle Scholar |

Fedy, B. C., Doherty, K. E., Aldridge, C. L., O’Donnell, M., Beck, J. L., Bedrosian, B., Gummer, D., Holloran, M. J., Johnson, G. D., Kaczor, N. W., Kirol, C. P., Mandich, C. A., Marshall, D., McKee, G., Olson, C., Pratt, A. C., Swanson, C. C., and Walker, B. L. (2014). Habitat prioritization across large landscapes, multiple seasons, and novel areas: an example using Greater sage-grouse in Wyoming. Wildlife Monographs 190, 1–39.
Habitat prioritization across large landscapes, multiple seasons, and novel areas: an example using Greater sage-grouse in Wyoming.Crossref | GoogleScholarGoogle Scholar |

Giesen, K. M., Schoenberg, T. J., and Braun, C. E. (1982). Methods for trapping sage grouse in Colorado. Wildlife Society Bulletin 10, 224–231.

Gillies, C. S., Hebblewhite, M., Nielsen, S. E., Krawchuk, M. A., Aldridge, C. L., Frair, J. L., Saher, D. J., Stevens, C. E., and Jerde, C. L. (2006). Application of random effects to the study of resource selection by animals. Journal of Animal Ecology 75, 887–898.
Application of random effects to the study of resource selection by animals.Crossref | GoogleScholarGoogle Scholar |

Greenwald, N. D., Suckling, K. F., and Pimm, S. L. (2012). Critical habitat and the role of peer review in government decisions. Bioscience 62, 686–690.
Critical habitat and the role of peer review in government decisions.Crossref | GoogleScholarGoogle Scholar |

Guisan, A., Tingley, R., Baumgartner, J. B., Naujokaitis-Lwis, I., Sutcliffe, P. R., Tulloch, A. I. T., Regan, T. J., Brotons, L., McDonald-Madden, E., Mantyka-Pringle, C., Martin, T. G., Rhodes, J. R., Maggini, R., Setterfield, S. A., Elith, J., Schwartz, M. W., Wintle, B. A., Broennimann, O., Austin, M., Ferrier, S., Kearney, M. R., Possingham, H. P., and Buckley, Y. M. (2013). Predicting species distributions for conservation decisions. Ecology Letters 16, 1424–1435.
Predicting species distributions for conservation decisions.Crossref | GoogleScholarGoogle Scholar |

Gunnison County (2014). Gunnison County notice of intent to sue. Available at http://www.gunnisoncounty.org/docuentcenter/view/3623 [Accessed 5 January 2017].

Gunnison Sage-grouse Rangewide Steering Committee (GSRSC) (2005). ‘Gunnison Sage-grouse Rangewide Conservation Plan.’ (Colorado Division of Wildlife: Denver, CO.)

Hausleitner, D. (2003). Population dynamics, habitat use and movements of greater sage-grouse in Moffat county, Colorado. M.S. Thesis, University of Idaho, Moscow, ID.

Hupp, J. W., and Braun, C. E. (1989). Topographic distribution of sage grouse foraging in winter. The Journal of Wildlife Management 53, 823–829.
Topographic distribution of sage grouse foraging in winter.Crossref | GoogleScholarGoogle Scholar |

Johnson, C. J., Seip, D. R., and Boyce, M. S. (2004). A quantitative approach to conservation planning: using resource selection functions to map the distribution of mountain caribou at multiple spatial scales. Journal of Applied Ecology 41, 238–251.
A quantitative approach to conservation planning: using resource selection functions to map the distribution of mountain caribou at multiple spatial scales.Crossref | GoogleScholarGoogle Scholar |

Kalen, S. (2014). Landscape shifting paradigm for the endangered species act: an integrated critical habitat recovery program. Natural Resources Journal 55, 47–104.

Klar, N., Fernández, N., Kramer-Schadt, S., Herrmann, M., Trinzen, M., Büttner, I., and Niemitz, C. (2008). Habitat selection models for the European wildcat conservation. Biological Conservation 141, 308–319.
Habitat selection models for the European wildcat conservation.Crossref | GoogleScholarGoogle Scholar |

Knapp, C. N., Cochran, J., Chapin, S. F., Kofinas, G., and Sayre, N. (2013). Putting local knowledge and context to work for Gunnison sage-grouse conservation. Human–Wildlife Interactions 7, 195–213.

Manly, B. F. J., McDonald, L. L., Thomas, D. L., McDonald, T. L., and Erickson, W. P. (2002). ‘Resource Selection by Animals: Statistical Design and Analysis for Field Studies.’ (Kluwer Academic Publishers: Dordrecht, The Netherlands.)

McAlpine, C. A., Rhodes, J. R., Bowen, M. E., Lunney, D., Callaghan, J. G., Mitchell, D. L., and Possingham, H. P. (2008). Can multiscale models of species’ distribution be generalized from region to region? A case study of koala. Journal of Applied Ecology 45, 558–567.
Can multiscale models of species’ distribution be generalized from region to region? A case study of koala.Crossref | GoogleScholarGoogle Scholar |

McGarigal, K., Cushman, S. A., and Stafford, S. G. (2000). ‘Multivariate Statistics for Wildlife and Ecology Research.’ (Springer-Verlag: New York.)

Murphy, D. D., and Weiland, P. S. (2016). Guidance on the use of best available science under the U.S. Endangered Species Act. Environmental Management 58, 1–15.
Guidance on the use of best available science under the U.S. Endangered Species Act.Crossref | GoogleScholarGoogle Scholar |

Northrup, J. M., Hooten, M. B., Anderson, C. R., and Wittemyer, G. (2013). Practical guidance on characterizing availability in resource selection functions under a use-availability design. Ecology 94, 1456–1463.
Practical guidance on characterizing availability in resource selection functions under a use-availability design.Crossref | GoogleScholarGoogle Scholar |

Oyler-McCance, S. J., Burnham, K. P., and Braun, C. E. (2001). Influence of changes in sagebrush on Gunnison sage grouse in southwestern Colorado. The Southwestern Naturalist 46, 323–331.
Influence of changes in sagebrush on Gunnison sage grouse in southwestern Colorado.Crossref | GoogleScholarGoogle Scholar |

R Core Team (2013). ‘R: a language and environment for statistical computing.’ (R Foundation for Statistical Computing: Vienna.) Available at http://www.R-project.org/.

Rice, M. B., Apa, A. D., Phillips, M. L., Gammonly, J. H., Petch, B. F., and Eichhoff, K. (2013). Analysis of regional species models based on radio-telemetry datasets from multiple small-scale studies. The Journal of Wildlife Management 77, 821–831.
Analysis of regional species models based on radio-telemetry datasets from multiple small-scale studies.Crossref | GoogleScholarGoogle Scholar |

Salzman, J. (1990). Evolution and application of critical habitat under the Endangered Species Act. The Harvard Environmental Law Review 14, 311–342.

Sawyer, H., Nielson, R. M., Lindzey, F. G., Keith, L., Powell, J. H., and Abraham, A. A. (2007). Habitat selection of Rocky Mountain elk in a nonforested environment. The Journal of Wildlife Management 71, 868–874.
Habitat selection of Rocky Mountain elk in a nonforested environment.Crossref | GoogleScholarGoogle Scholar |

Schielzeth, H. (2010). Simple means to improve the interpretability of regression coefficients. Methods in Ecology and Evolution 1, 103–113.
Simple means to improve the interpretability of regression coefficients.Crossref | GoogleScholarGoogle Scholar |

Schroeder, M. A., Aldridge, C. L., Apa, A. D., Bohne, J. R., Braun, C. E., Bunnell, S. D., Connelly, J. W., Deibert, P. A., Gardner, S. C., Hilliard, M. A., Kobriger, G. D., McAdam, S. M., and McCArthy, Schroeder, M. A., Aldridge, C. L., Apa, A. D., Bohne, J. R., Braun, C. E., Bunnell, S. D., Connelly, J. W., Deibert, P. A., Gardner, S. C., Hilliard, M. A., Kobriger, G. D., McAdam, S. M., and McCArthy, (2004). Distribution of sage-grouse in North America. The Condor 106, 363–376.
Distribution of sage-grouse in North America.Crossref | GoogleScholarGoogle Scholar |

Scott, J. M., Goble, D. D., and Wiens, J. A. (2010). Conservation-reliant species and the future of conservation. Conservation Letters 3, 91–97.
Conservation-reliant species and the future of conservation.Crossref | GoogleScholarGoogle Scholar |

Shirk, A. J., Raphael, M. G., and Cushman, S. A. (2014). Spatiotemporal variation in resource selection: insights from the American marten (Martes americana). Ecological Applications 24, 1434–1444.
Spatiotemporal variation in resource selection: insights from the American marten (Martes americana).Crossref | GoogleScholarGoogle Scholar |

Stanley, T. R., Aldridge, C. L., Saher, D. J., and Childers, T. M. (2015). Daily nest survival rates of Gunnison sage-grouse (Centrocercus minimus): assessing local- and landscape-scale drivers. The Wilson Journal of Ornithology 127, 59–71.
Daily nest survival rates of Gunnison sage-grouse (Centrocercus minimus): assessing local- and landscape-scale drivers.Crossref | GoogleScholarGoogle Scholar |

Stokland, J. N., Halvorsen, R., and Stoa, B. (2011). Species distribution modeling- effect of design and sample size of pseudo-absence observations. Ecological Modelling 222, 1800–1809.
Species distribution modeling- effect of design and sample size of pseudo-absence observations.Crossref | GoogleScholarGoogle Scholar |

Thompson, T. R. (2012). Dispersal ecology of greater sage-grouse in northwestern Colorado: evidence from demographic and genetics methods. Ph.D. Dissertation, University of Idaho, Moscow, ID.

US Fish and Wildlife Service (USFWS) (2014a). Endangered and Threatened Wildlife and Plants; Determination for the Gunnison Sage-grouse as a Threatened or Endangered Species. Final Rule. 50 CFR Part 17. Federal Register 79, 69192–69310.

USFWS (2014b). Endangered and Threatened Wildlife and Plants; Designation of Critical Habitat for Gunnison Sage-grouse. 50 CFR Part 17. Federal Register 79, 69312–69363.

USFWS (2015). Critical habitat: what is it? USFWS Endangered Species Program, Falls Church, VA. Available at https://www.fws.gov/midwest/endangered/saving/CriticalHabitatFactSheet.html [accessed May 15, 2015.]

Wakkinen, W. L., Reese, K. P., Connelly, J. W., and Fischer, R. A. (1992). An improved spotlighting technique for capturing sage-grouse. Wildlife Society Bulletin 20, 425–426.

Walker, B. L., Apa, A. D., and Eichhoff, K. (2016). Mapping and prioritizing seasonal habitats for greater sage-grouse in northwestern Colorado. The Journal of Wildlife Management 80, 63–77.
Mapping and prioritizing seasonal habitats for greater sage-grouse in northwestern Colorado.Crossref | GoogleScholarGoogle Scholar |

Wiens, T. S., Dale, B. C., Boyce, M. S., and Kershaw, G. P. (2008). Three way k-fold cross validation of resource selection functions. Ecological Modelling 212, 244–255.
Three way k-fold cross validation of resource selection functions.Crossref | GoogleScholarGoogle Scholar |

Wisdom, M. J., Meinke, C. W., Knick, S. T., and Schroeder, M. A. (2011). Factors associated with extirpation of sage-grouse. In ‘Studies in Avian Biology; Greater Sage-grouse: Ecology and Conservation of a Landscape Species and its Habitats’. (Eds S. T. Knick and J. W. Connelly.) pp. 69–83. (University of California: Berkeley, CA.)

Young, J. R., Braun, C. E., Oyler-McCance, S. J., Hupp, J. W., and Quinn, T. W. (2000). A new species of sage-grouse (Phasianidae: centrocercus) from southwestern Colorado. The Wilson Bulletin 112, 445–453.
A new species of sage-grouse (Phasianidae: centrocercus) from southwestern Colorado.Crossref | GoogleScholarGoogle Scholar |