Predicting the spatial distribution of wintering golden eagles to inform full annual cycle conservation in western North America.

Wildlife conservation strategies focused on one season or population segment may fail to adequately protect populations, especially when a species' habitat preferences vary among seasons, age-classes, geographic regions, or other factors. Conservation of golden eagles (Aquila chrysaetos) is an example of such a complex scenario, in which the distribution, habitat use, and migratory strategies of this species of conservation concern vary by age-class, reproductive status, region, and season. Nonetheless, research aimed at mapping priority use areas to inform management of golden eagles in western North America has typically focused on territory-holding adults during the breeding period, largely to the exclusion of other seasons and life-history groups. To support population-wide conservation planning across the full annual cycle for golden eagles, we developed a distribution model for individuals in a season not typically evaluated-winter-and in an area of the interior western U.S. that is a high priority for conservation of the species. We used a large GPS-telemetry dataset and library of environmental variables to develop a machine-learning model to predict spatial variation in the relative intensity of use by golden eagles during winter in Wyoming, USA, and surrounding ecoregions. Based on a rigorous series of evaluations including cross-validation, withheld and independent data, our winter-season model accurately predicted spatial variation in intensity of use by multiple age- and life-history groups of eagles not associated with nesting territories (i.e., all age classes of long-distance migrants, and resident non-adults and adult "floaters", and movements of adult territory holders and their offspring outside their breeding territories). Important predictors in the model were wind and uplift (40.2% contribution), vegetation and landcover (27.9%), topography (14%), climate and weather (9.4%), and ecoregion (8.7%). Predicted areas of high-use winter habitat had relatively low spatial overlap with nesting habitat, suggesting a conservation strategy targeting high-use areas for one season would capture as much as half and as little as one quarter of high-use areas for the other season. The majority of predicted high-use habitat (top 10% quantile) occurred on private lands (55%); lands managed by states and the Bureau of Land Management (BLM) had a lower amount (33%), but higher concentration of high-use habitat than expected for their area (1.5-1.6x). These results will enable those involved in conservation and management of golden eagles in our study region to incorporate spatial prioritization of wintering habitat into their existing regulatory processes, land-use planning tasks, and conservation actions.

Modeling spatial variation in density of golden eagle nest sites in the western United States.

In order to contribute to conservation planning efforts for golden eagles (Aquila chrysaetos) in the western U.S., we developed nest site models using >6,500 nest site locations throughout a >3,483,000 km2 area of the western U.S. We developed models for twelve discrete modeling regions, and estimated relative density of nest sites for each region. Cross-validation showed that, in general, models accurately estimated relative nest site densities within regions and sub-regions. Areas estimated to have the highest densities of breeding golden eagles had from 132-2,660 times greater densities compared to the lowest density areas. Observed nest site densities were very similar to those reported from published studies. Large extents of each modeling region consisted of low predicted nest site density, while a small percentage of each modeling region contained disproportionately high nest site density. For example, we estimated that areas with relative nest density values <0.3 represented from 62.8-97.8% ([Formula: see text] = 82.5%) of each modeling area, and those areas contained from 14.7-30.0% ([Formula: see text] = 22.1%) of the nest sites. In contrast, areas with relative nest density values >0.5 represented from 1.0-12.8% ([Formula: see text] = 6.3%) of modeling areas, and those areas contained from 47.7-66.9% ([Formula: see text] = 57.3%) of the nest sites. Our findings have direct application to: 1) large-scale conservation planning efforts, 2) risk analyses for land-use proposals such as recreational trails or wind power development, and 3) identifying mitigation areas to offset the impacts of human disturbance.

Conservation planning for species recovery under the Endangered Species Act: A case study with the Northern Spotted Owl.

The northern spotted owl (Strix occidentalis caurina) was listed as threatened under the U.S. Endangered Species Act (ESA) in 1990. We applied modern spatial conservation theory and models to evaluate several candidate critical habitat networks, and sought an efficient conservation solution that encompassed the highest value lands for spotted owl recovery rather than maximizing the total area of potential critical habitat. We created a map of relative habitat suitability, which served as input to the spatial conservation prioritization program Zonation. We used the spatially-explicit individual-based population model HexSim to estimate and compare simulated spotted owl population outcomes among a suite of candidate critical habitat networks that varied in size and spatial arrangement under alternative scenarios of future habitat suitability and barred owl (S. varia) effects. We evaluated simulated spotted owl population outcomes, including total population size, and extinction and quasi-extinction likelihoods for 108 combinations of candidate critical habitat networks by habitat change by barred owl scenarios, both range-wide and within 11 distinct portions of the owl's range. Barred owl encounter rates and the amount and suitability of habitat had substantial effects on simulated spotted owl populations. When barred owl encounter rates were high, changes in the amount and suitability of habitat had minimal impacts on population performance. Under lowered barred owl encounter rates, candidate critical habitat networks that included most existing high suitability habitat supported a high likelihood of long-term population persistence. Barred owls are currently the primary driving force behind poor population performance of NSOs; however, our models demonstrated that a sufficient area of high suitability habitat remains essential for recovery when effects of barred owls can be reduced. The modeling approach we employed is sufficiently flexible to incorporate new information about spotted owls as it becomes available and could likely be applied to conservation planning for other species.

Hierarchical Bayesian spatial models for multispecies conservation planning and monitoring.

Biologists who develop and apply habitat models are often familiar with the statistical challenges posed by their data's spatial structure but are unsure of whether the use of complex spatial models will increase the utility of model results in planning. We compared the relative performance of nonspatial and hierarchical Bayesian spatial models for three vertebrate and invertebrate taxa of conservation concern (Church's sideband snails [Monadenia churchi], red tree voles [Arborimus longicaudus], and Pacific fishers [Martes pennanti pacifica]) that provide examples of a range of distributional extents and dispersal abilities. We used presence-absence data derived from regional monitoring programs to develop models with both landscape and site-level environmental covariates. We used Markov chain Monte Carlo algorithms and a conditional autoregressive or intrinsic conditional autoregressive model framework to fit spatial models. The fit of Bayesian spatial models was between 35 and 55% better than the fit of nonspatial analogue models. Bayesian spatial models outperformed analogous models developed with maximum entropy (Maxent) methods. Although the best spatial and nonspatial models included similar environmental variables, spatial models provided estimates of residual spatial effects that suggested how ecological processes might structure distribution patterns. Spatial models built from presence-absence data improved fit most for localized endemic species with ranges constrained by poorly known biogeographic factors and for widely distributed species suspected to be strongly affected by unmeasured environmental variables or population processes. By treating spatial effects as a variable of interest rather than a nuisance, hierarchical Bayesian spatial models, especially when they are based on a common broad-scale spatial lattice (here the national Forest Inventory and Analysis grid of 24 km(2) hexagons), can increase the relevance of habitat models to multispecies conservation planning.

Using forest inventory data to assess fisher resting habitat suitability in California.

The fisher (Martes pennanti) is a forest-dwelling carnivore whose current distribution and association with late-seral forest conditions make it vulnerable to stand-altering human activities or natural disturbances. Fishers select a variety of structures for daily resting bouts. These habitat elements, together with foraging and reproductive (denning) habitat, constitute the habitat requirements of fishers. We develop a model capable of predicting the suitability of fisher resting habitat using standard forest vegetation inventory data. The inventory data were derived from Forest Inventory and Analysis (FIA), a nationwide probability-based sample used to estimate forest characteristics. We developed the model by comparing vegetation and topographic data at 75 randomly selected fisher resting structures in the southern Sierra Nevada with 232 forest inventory plots. We collected vegetation data at fisher resting locations using the FIA vegetation sampling protocol and centering the 1-ha FIA plot on the resting structure. To distinguish used and available inventory plots, we used nonparametric logistic regression to evaluate a set of a priori biological models. The top model represented a dominant portion of the Akaike weights (0.87), explained 31.5% of the deviance, and included the following variables: average canopy closure, basal area of trees <51 cm diameter breast height (dbh), average hardwood dbh, maximum tree dbh, percentage slope, and the dbh of the largest conifer snag. Our use of routinely collected forest inventory data allows the assessment and monitoring of change in fisher resting habitat suitability over large regions with no additional sampling effort. Although models were constrained to include only variables available from the list of those measured using the FIA protocol, we did not find this to be a shortcoming. The model makes it possible to compare average resting habitat suitability values before and after forest management treatments, among administrative units, across regions and over time. Considering hundreds of plot estimates as a sample of habitat conditions over large spatial scales can bring a broad perspective, at high resolution, and efficiency to the assessment and monitoring of wildlife habitat.