Although climate acts as a fundamental constraint on the distribution of organisms, understanding how this relationship between climate and distribution varies over a species' range is critical for addressing the potential impacts of accelerated climate change on biodiversity. Bioclimatic niche models provide compelling evidence that many species will experience range shifts under scenarios of global change, yet these broad, macroecological perspectives lack specificity at local scales, where unique combinations of environment, biota, and history conspire against generalizations. We explored how these idiosyncrasies of place affect the climate-distribution relationship of the American pika (Ochotona princeps) by replicating intensive field surveys across bioclimatic gradients in eight U.S. national parks. At macroecological scales, the importance of climate as a constraint on pika distribution appears unequivocal; forecasts suggest that the species' range will contract sharply in coming decades. However, the species persists outside of its modeled bioclimatic envelope in many locations, fueling uncertainty and debate over its conservation status. Using a Bayesian hierarchical approach, we modeled variation in local patterns of pika distribution along topographic position, vegetation cover, elevation, temperature, and precipitation gradients in each park landscape. We also accounted for annual turnover in site occupancy probabilities. Topographic position and vegetation cover influenced occurrence in all parks. After accounting for these factors, pika occurrence varied widely among parks along bioclimatic gradients. Precipitation by itself was not a particularly influential predictor. However, measures of heat stress appeared most influential in the driest parks, suggesting an interaction between the strength of climate effects and the position of parks along precipitation gradients. The combination of high elevation, cold temperatures, and high precipitation lowered occurrence probabilities in some parks, suggesting an upper elevational limit for pikas in some environments. Our results demonstrate that the idiosyncrasies of place influence both the nature and strength of the climate-distribution relationship for the American pika. Fine-grained, but geographically extensive, studies replicated across multiple landscapes offer insights important to assessing the impacts of climate change that otherwise may be masked at macroecological scales. The hierarchical approach to modeling provides a coherent conceptual and technical framework for gaining these insights.

BioOne Complete (complete.BioOne.org) is a full-text database of 200 subscribed and open-access titles in the biological, ecological, and environmental sciences published by nonprofit societies, associations, museums, institutions, and presses.

No abstract

Ecological niche theory holds that species distributions are shaped by a large and complex suite of interacting factors. Species distribution models (SDMs) are increasingly used to describe species' niches and predict the effects of future environmental change, including climate change. Currently, SDMs often fail to capture the complexity of species' niches, resulting in predictions that are generally limited to climate-occupancy interactions. Here, we explore the potential impact of climate change on the American pika using a replicated place-based approach that incorporates climate, gene flow, habitat configuration, and microhabitat complexity into SDMs. Using contemporary presence-absence data from occupancy surveys, genetic data to infer connectivity between habitat patches, and 21 environmental niche variables, we built separate SDMs for pika populations inhabiting eight US National Park Service units representing the habitat and climatic breadth of the species across the western United States. We then predicted occurrence probability under current (1981-2010) and three future time periods (out to 2100). Occurrence probabilities and the relative importance of predictor variables varied widely among study areas, revealing important local-scale differences in the realized niche of the American pika. This variation resulted in diverse and - in some cases - highly divergent future potential occupancy patterns for pikas, ranging from complete extirpation in some study areas to stable occupancy patterns in others. Habitat composition and connectivity, which are rarely incorporated in SDM projections, were influential in predicting pika occupancy in all study areas and frequently outranked climate variables. Our findings illustrate the importance of a place-based approach to species distribution modeling that includes fine-scale factors when assessing current and future climate impacts on species' distributions, especially when predictions are intended to manage and conserve species of concern within individual protected areas.

Strategic conservation efforts for cryptic species, especially bats, are hindered by limited understanding of distribution and population trends. Integrating long‐term encounter surveys with multi‐season occupancy models provides a solution whereby inferences about changing occupancy probabilities and latent changes in abundance can be supported. When harnessed to a Bayesian inferential paradigm, this modeling framework offers flexibility for conservation programs that need to update prior model‐based understanding about at‐risk species with new data. This scenario is exemplified by a bat monitoring program in the Pacific Northwestern United States in which results from 8 years of surveys from 2003 to 2010 require updating with new data from 2016 to 2018. The new data were collected after the arrival of bat white‐nose syndrome and expansion of wind power generation, stressors expected to cause population declines in at least two vulnerable species, little brown bat (Myotis lucifugus) and the hoary bat (Lasiurus cinereus). We used multi‐season occupancy models with empirically informed prior distributions drawn from previous occupancy results (2003–2010) to assess evidence of contemporary decline in these two species. Empirically informed priors provided the bridge across the two monitoring periods and increased precision of parameter posterior distributions, but did not alter inferences relative to use of vague priors. We found evidence of region‐wide summertime decline for the hoary bat (trueλ^ = 0.86 ± 0.10) since 2010, but no evidence of decline for the little brown bat (trueλ^ = 1.1 ± 0.10). White‐nose syndrome was documented in the region in 2016 and may not yet have caused regional impact to the little brown bat. However, our discovery of hoary bat decline is consistent with the hypothesis that the longer duration and greater geographic extent of the wind energy stressor (collision and barotrauma) have impacted the species. These hypotheses can be evaluated and updated over time within our framework of pre–post impact monitoring and modeling. Our approach provides the foundation for a strategic evidence‐based conservation system and contributes to a growing preponderance of evidence from multiple lines of inquiry that bat species are declining.

Aim Bat mortality rates from white-nose syndrome and wind power development are unprecedented. Cryptic and wide-ranging behaviours of bats make them difficult to survey, and population estimation is often intractable. We advance a model-based framework for making spatially explicit predictions about summertime distributions of bats from capture and acoustic surveys. Motivated by species-energy and life-history theory, our models describe hypotheses about spatio-temporal variation in bat distributions along environmental gradients and life-history attributes, providing a statistical basis for conservation decision-making.Location Oregon and Washington, USA.Methods We developed Bayesian hierarchical models for 14 bat species from an 8-year monitoring dataset across a~430,000 km 2 study area. Models accounted for imperfect detection and were temporally dynamic. We mapped predicted occurrence probabilities and prediction uncertainties as baselines for assessing future declines.Results Forest cover, snag abundance and cliffs were important predictors for most species. Species occurrence patterns varied along elevation and precipitation gradients, suggesting a potential hump-shaped diversity-productivity relationship. Annual turnover in occurrence was generally low, and occurrence probabilities were stable among most species. We found modest evidence that turnover covaried with the relative riskiness of bat roosting and migration. The fringed myotis (Myotis thysanodes), canyon bat (Parastrellus hesperus) and pallid bat (Antrozous pallidus) were rare; fringed myotis occurrence probabilities declined over the study period. We simulated anticipated declines to demonstrate that mapped occurrence probabilities, updated over time, provide an intuitive way to assess bat conservation status for a broad audience.Main conclusions Landscape keystone structures associated with roosting habitat emerged as regionally important predictors of bat distributions. The challenges of bat monitoring have constrained previous species distribution modelling efforts to temporally static presence-only approaches. Our approach extends to broader spatial and temporal scales than has been possible in the past for bats, making a substantial increase in capacity for bat conservation.

Bats face unprecedented threats from habitat loss, climate change, disease, and wind power development, and populations of many species are in decline. A better ability to quantify bat population status and trend is urgently needed in order to develop effective conservation strategies. We used a Bayesian autoregressive approach to develop dynamic distribution models for Myotis lucifugus, the little brown bat, across a large portion of northwestern USA, using a four-year detection history matrix obtained from a regional monitoring program. This widespread and abundant species has experienced precipitous local population declines in northeastern USA resulting from the novel disease white-nose syndrome, and is facing likely range-wide declines. Our models were temporally dynamic and accounted for imperfect detection. Drawing on species-energy theory, we included measures of net primary productivity (NPP) and forest cover in models, predicting that M. lucifugus occurrence probabilities would covary positively along those gradients. Despite its common status, M. lucifugus was only detected during -50% of the surveys in occupied sample units. The overall naive estimate for the proportion of the study region occupied by the species was 0.69, but after accounting for imperfect detection, this increased to -0.90. Our models provide evidence of an association between NPP and forest cover and M. lucifugus distribution, with implications for the projected effects of accelerated climate change in the region, which include net aridification as snowpack and stream flows decline. Annual turnover, the probability that an occupied sample unit was a newly occupied one, was estimated to be low (-0.04-0.14), resulting in flat trend estimated with relatively high precision (SD = 0.04). We mapped the variation in predicted occurrence probabilities and corresponding prediction uncertainty along the productivity gradient. Our results provide a much needed baseline against which future anticipated declines in M. lucifugus occurrence can be measured. The dynamic distribution modeling approach has broad applicability to regional bat monitoring efforts now underway in several countries and we suggest ways to improve and expand our grid-based monitoring program to gain robust insights into bat population status and trend across large portions of North America.

Summary1. Acoustic surveys have become a common survey method for bats and other vocal taxa. Previous work shows that bat echolocation may be misidentified, but common analytic methods, such as occupancy models, assume that misidentifications do not occur. Unless rare, such misidentifications could lead to incorrect inferences with significant management implications. 2. We fit a false-positive occupancy model to data from paired bat detector and mist-net surveys to estimate probability of presence when survey data may include false positives. We compared estimated occupancy and detection rates to those obtained from a standard occupancy model. We also derived a formula to estimate the probability that bats were present at a site given its detection history. As an example, we analysed survey data for little brown bats Myotis lucifugus from 135 sites in Washington and Oregon, USA. 3. We estimated that at an unoccupied site, acoustic surveys had a 14% chance per night of producing spurious M. lucifugus detections. Estimated detection rates were higher and occupancy rates were lower under the false-positive model, relative to a standard occupancy model. Un-modelled false positives also affected inferences about occupancy at individual sites. For example, probability of occupancy at individual sites with acoustic detections but no captures ranged from 2% to 100% under the false-positive occupancy model, but was always 100% under a standard occupancy model. 4. Synthesis and applications. Our results suggest that false positives sufficient to affect inferences may be common in acoustic surveys for bats. We demonstrate an approach that can estimate occupancy, regardless of the false-positive rate, when acoustic surveys are paired with capture surveys. Applications of this approach include monitoring the spread of WhiteNose Syndrome, estimating the impact of climate change and informing conservation listing decisions. We calculate a site-specific probability of occupancy, conditional on survey results, which could inform local permitting decisions, such as for wind energy projects. More generally, the magnitude of false positives suggests that false-positive occupancy models can improve accuracy in research and monitoring of bats and provide wildlife managers with more reliable information.