Climate Tolerances and Habitat Requirements Jointly Shape the Elevational Distribution of the American Pika (Ochotona princeps), with Implications for Climate Change Effects

An observational analysis of Canada Jay (Perisoreus canadensis) foraging and caching ecology in Denali National Park and Preserve, Alaska, USA

Arctic and subarctic wildlife are among the most vulnerable species to climate change. Canada Jays (Perisoreus canadensis (Linnaeus, 1776)) are generalist residents of northern boreal forests and scatter-hoard food to insulate against food scarcity during winter. Unlike most scatter-hoarders, however, Canada Jays primarily cache perishable food, rendering their caches more susceptible to climate change induced degradation and loss. Here we use a mostly noninvasive approach to document Canada Jay foraging ecology among a population in interior Alaska, USA, including the types of food acquired, foraging and caching rates, and cache longevity and loss. We also tested for associations between foraging and caching rates with reproductive metrics to assess possible relationships among food and productivity. We found that Canada Jays have a varied diet that changed seasonally, and responded to a record-setting warm spring by directing foraging efforts away from cache recovery and towards the emergence of fresh food. We did not find evidence for relationships between foraging and caching rate with reproductive output, possibly owing to small sample sizes. We found that caches were recovered quickly (<4 weeks) and frequently lost to conspecific and heterospecific competitors. Our study suggests that Canada Jays may be better poised to respond to changes in cache integrity and food availability than has been previously recognized.

Genome-wide analysis reveals associations between climate and regional patterns of adaptive divergence and dispersal in American pikas

Understanding the role of adaptation in species' responses to climate change is important for evaluating the evolutionary potential of populations and informing conservation efforts. Population genomics provides a useful approach for identifying putative signatures of selection and the underlying environmental factors or biological processes that may be involved. Here, we employed a population genomic approach within a space-for-time study design to investigate the genetic basis of local adaptation and reconstruct patterns of movement across rapidly changing environments in a thermally sensitive mammal, the American pika (Ochotona princeps). Using genotypic data at 49,074 single-nucleotide polymorphisms (SNPs), we analyzed patterns of genome-wide diversity, structure, and migration along three independent elevational transects located at the northern extent (Tweedsmuir South Provincial Park, British Columbia, Canada) and core (North Cascades National Park, Washington, USA) of the Cascades lineage. We identified 899 robust outlier SNPs within- and among-transects. Of those annotated to genes with known function, many were linked with cellular processes related to climate stress including ATP-binding, ATP citrate synthase activity, ATPase activity, hormone activity, metal ion-binding, and protein-binding. Moreover, we detected evidence for contrasting patterns of directional migration along transects across geographic regions that suggest an increased propensity for American pikas to disperse among lower elevation populations at higher latitudes where environments are generally cooler. Ultimately, our data indicate that fine-scale demographic patterns and adaptive processes may vary among populations of American pikas, providing an important context for evaluating biotic responses to climate change in this species and other alpine-adapted mammals.

Factors influencing distributional shifts and abundance at the range core of a climate-sensitive mammal

Species are frequently responding to contemporary climate change by shifting to higher elevations and poleward to track suitable climate space. However, depending on local conditions and species' sensitivity, the nature of these shifts can be highly variable and difficult to predict. Here, we examine how the American pika (Ochotona princeps), a philopatric, montane lagomorph, responds to climatic gradients at three spatial scales. Using mixed-effects modeling in an information-theoretic approach, we evaluated a priori model suites regarding predictors of site occupancy, relative abundance, and elevational-range retraction across 760 talus patches, nested within 64 watersheds across the Northern Rocky Mountains of North America, during 2017-2020. The top environmental predictors differed across these response metrics. Warmer temperatures in summer and winter were associated with lower occupancy, lower relative abundances, and greater elevational retraction across watersheds. Occupancy was also strongly influenced by habitat patch size, but only when combined with climate metrics such as actual evapotranspiration. Using a second analytical approach, acute heat stress and summer precipitation best explained retraction residuals (i.e., the relative extent of retraction given the original elevational range of occupancy). Despite the study domain occurring near the species' geographic-range center, where populations might have higher abundances and be at lower risk of climate-related stress, 33.9% of patches showed evidence of recent extirpations. Pika-extirpated sites averaged 1.44℃ warmer in summer than did occupied sites. Additionally, the minimum elevation of pika occupancy has retracted upslope in 69% of watersheds (mean: 281 m). Our results emphasize the nuance associated with evaluating species' range dynamics in response to climate gradients, variability, and temperature exceedances, especially in regions where species occupy gradients of conditions that may constitute multiple range edges. Furthermore, this study highlights the importance of evaluating diverse drivers across response metrics to improve the predictive accuracy of widely used, correlative models.

Temporal vs. spatial variation in stress-associated metabolites within a population of climate-sensitive small mammals

Temporal variation in stress might signify changes in an animal's internal or external environment, while spatial variation in stress might signify variation in the quality of the habitats that individual animals experience. Habitat-induced variations in stress might be easiest to detect in highly territorial animals, and especially in species that do not take advantage of common strategies for modulating habitat-induced stress, such as migration (escape in space) or hibernation (escape in time). Spatial and temporal variation in response to potential stressors has received little study in wild animals, especially at scales appropriate for relating stress to specific habitat characteristics. Here, we use the American pika (Ochotona princeps), a territorial small mammal, to investigate stress response within and among territories. For individually territorial animals such as pikas, differences in habitat quality should lead to differences in stress exhibited by territory owners. We indexed stress using stress-associated hormone metabolites in feces collected non-invasively from pika territories every 2 weeks from June to September 2018. We hypothesized that differences in territory quality would lead to spatial differences in mean stress and that seasonal variation in physiology or the physical environment would lead to synchronous variation across territories through time. We used linear mixed-effects models to explore spatiotemporal variation in stress using fixed effects of day-of-year and broad habitat characteristics (elevation, aspect, site), as well as local variation in habitat characteristics hypothesized to affect territory quality for this saxicolous species (talus depth, clast size, available forage types). We found that temporal variation within territories was greater than spatial variation among territories, suggesting that shared seasonal stressors are more influential than differences in individual habitat quality. This approach could be used in other wildlife studies to refine our understanding of habitat quality and its effect on individual stress levels as a driver of population decline.

Conservation status of American pikas (Ochotona princeps)

The American pika (Ochotona princeps) is commonly perceived as a species that is at high risk of extinction due to climate change. The purpose of this review is two-fold: to evaluate the claim that climate change is threatening pikas with extinction, and to summarize the conservation status of the American pika. Most American pikas inhabit major cordilleras, such as the Rocky Mountain, Sierra Nevada, and Cascade ranges. Occupancy of potential pika habitat in these ranges is uniformly high and no discernible climate signal has been found that discriminates between the many occupied and relatively few unoccupied sites that have been recently surveyed. Pikas therefore are thriving across most of their range. The story differs in more marginal parts of the species range, primarily across the Great Basin, where a higher percentage of available habitat is unoccupied. A comprehensive review of Great Basin pikas revealed that occupied sites, sites of recent extirpation, and old sites, were regularly found within the same geographic and climatic space as extant sites, and suggested that pikas in the Great Basin tolerated a broader set of habitat and climatic conditions than previously understood. Studies of a small subset of extirpated sites in the Great Basin and in California found that climate variables (most notably measures of hot temperature) were associated more often with extirpated sites than occupied sites. Importantly, upward contraction of the lower elevation boundary also was found at some sites. However, models that incorporated variables other than climate (such as availability of upslope talus habitat) often were better predictors of site persistence. Many extirpations occurred on small habitat patches, which were subject to stochastic extinction, as informed by a long-term pika metapopulation study in Bodie, California. In addition, several sites may have been compromised by cattle grazing or other anthropogenic factors. In contrast, several low, hot sites (Bodie, Mono Craters, Craters of the Moon National Monument and Preserve, Lava Beds National Monument, Columbia River Gorge) retain active pika populations, demonstrating the adaptive capacity and resilience of pikas in response to adverse environmental conditions. Pikas cope with warm temperatures by retreating into cool interstices of their talus habitat and augment their restricted daytime foraging with nocturnal activity. Pikas exhibit significant flexibility in their foraging tactics and are highly selective in their choice of available vegetation. The trait that places pikas at greatest risk from climate change is their poor dispersal capability. Dispersal is more restricted in hotter environments, and isolated low-elevation sites that become extirpated are unlikely to be recolonized in a warming climate. The narrative that American pikas are going extinct appears to be an overreach. Pikas are doing well across most of their range, but there are limited, low-elevation losses that are likely to be permanent in what is currently marginal pika habitat. The resilience of pikas in the face of climate change, and their ability or inability to persist in marginal, hot environments, will continue to contribute to our understanding of the impact of climate change on individual species.

Elevation and landscape change drive the distribution of a montane, endemic grassland bird

CONTEXT

Tropical montane habitats support high biodiversity and are hotspots of endemism, with grasslands being integral components of many such landscapes. The montane grasslands of the Western Ghats have seen extensive land-use change over anthropogenic timescales. The factors influencing the ability of grassland-dependent species to persist in habitats experiencing loss and fragmentation, particularly in montane grasslands, are poorly known.

OBJECTIVES

We studied the relationship between the Nilgiri pipit Anthus nilghiriensis, a threatened endemic bird that typifies these montane grasslands, and its habitat, across most of its global distribution. We examined what habitat features make remnant grasslands viable, which is necessary for their effective management.

METHODS

We conducted 663 surveys in 170 sites and used both single-season occupancy modeling and N-mixture modeling to account for processes influencing detection, presence, and abundance.

RESULTS

Elevation had a positive influence on species presence, patch size had a moderate positive influence, and patch isolation had a moderate negative influence. Species abundance was positively influenced by elevation and characteristics related to habitat structure, and negatively influenced by the presence of invasive woody vegetation.

CONCLUSIONS

The strong effect of elevation on the highly range-restricted Nilgiri pipit is likely to make it vulnerable to climate change. This highly range-restricted species is locally extinct at several locations, and persists at low densities in remnants of its habitat left by recent fragmentation. Our findings indicate a need to control and reverse the spread of exotic woody invasives to preserve the grasslands themselves and the specialist species dependent upon them.

Investigating the role of environment in pika (Ochotona) body size patterns across taxonomic levels, space, and time

Body size is an important trait in animals because it influences a multitude of additional life history traits. The causal mechanisms underlying body size patterns across spatial, temporal, and taxonomic hierarchies are debated, and of renewed interest in this era of climate change. Here, we tested multiple hypotheses regarding body mass patterns at the intraspecific and interspecific levels. We investigated body size patterns within a climate-sensitive small mammal species, Ochotona princeps (n = 2,873 individuals), across their range with local environmental variables. We also examined body mass of populations over time to determine if body size has evolved in situ in response to environmental change. At the interspecific level we compared the mean mass of 26 pika species (genus Ochotona) to determine if environmental temperatures, food availability, habitat variability, or range area influence body size. We found correlations between temperature, vegetation, and particularly precipitation variables, with body mass within O. princeps, but no linear relationship between body size and any climate or habitat variable for Ochotona species. Body size trends in relation to climate were stronger at the intraspecific than the interspecific level. Our results suggest that body size within O. princeps likely is related to food availability, and that body size evolution is not always a viable response to temperature change. Different mechanisms may be driving body size at the interspecific and intraspecific levels and factors other than environment, such as biotic interactions, may also be influential in determining body size over space and time.

Ecological consequences of anomalies in atmospheric moisture and snowpack

Although increased frequency of extreme-weather events is one of the most secure predictions associated with contemporary climate change, effects of such events on distribution and abundance of climate-sensitive species remain poorly understood. Montane ecosystems may be especially sensitive to extreme weather because of complex abiotic and biotic interactions that propagate from climate-driven reductions in snowpack. Snowpack not only protects subnivean biotas from extreme cold, but also influences forage availability through timing of melt-off and water availability. We related relative abundances of an alpine mammal, the American pika (Ochotona princeps), to measures of weather and snowpack dynamics over an 8-yr period that included before and after a year of record-low snowpack in Washington, USA. We sought to (1) quantify any change in pika abundance associated with the snowpack anomaly and (2) identify aspects of weather and snowpack that influenced abundance of pikas. Pikas showed a 1-yr lag response to the snowpack anomaly and exhibited marked declines in abundance at elevations below 1,400 m simultaneous with increased abundances at higher elevations. Atmospheric moisture, indexed by vapor pressure deficit (VPD), was especially important, evidenced by strong support for the top-ranked model that included the interaction of VPD with snowpack duration. Notably, our novel application of VPD from gridded climate data for analyses of animal abundances shows strong potential for improving species distribution models because VPD represents an important aspect of weather that influences the physiology and habitat of biota. Pikas were apparently affected by cold stress without snowpack at mid elevations, whereas changes to forage associated with snowpack and VPD were influential at high and low elevations. Our results reveal context dependency in pika responses to weather and illustrate how snow drought can lead to rapid change in the abundance of subnivean animals.

Evolution for extreme living: variation in mitochondrial cytochrome c oxidase genes correlated with elevation in pikas (genus Ochotona)

The genus Ochotona (pikas) is a clade of cold-tolerant lagomorphs that includes many high-elevation species. Pikas offer a unique opportunity to study adaptations and potential limitations of an ecologically important mammal to high-elevation hypoxia. We analyzed the evolution of 3 mitochondrial genes encoding the catalytic core of cytochrome c oxidase (COX) in 10 pika species occupying elevations from sea level to 5000 m. COX is an enzyme highly reliant on oxygen and essential for cell function. One amino acid property, the equilibrium constant (ionization of COOH), was found to be under selection in the overall protein complex. We observed a strong relationship between the net value change in this property and the elevation each species occupies, with higher-elevation species having potentially more efficient proteins. We also found evidence of selection in low-elevation species for potentially less efficient COX, perhaps trading efficiency for heat production in the absence of hypoxia. Our results suggest that different pika species may have evolved elevation-specific COX proteins, specialization that may indicate limitations in their ability to shift their elevational ranges in response to future climate change.

Alpine glacial relict species losing out to climate change: The case of the fragmented mountain hare population (Lepus timidus) in the Alps

Alpine and Arctic species are considered to be particularly vulnerable to climate change, which is expected to cause habitat loss, fragmentation and-ultimately-extinction of cold-adapted species. However, the impact of climate change on glacial relict populations is not well understood, and specific recommendations for adaptive conservation management are lacking. We focused on the mountain hare (Lepus timidus) as a model species and modelled species distribution in combination with patch and landscape-based connectivity metrics. They were derived from graph-theory models to quantify changes in species distribution and to estimate the current and future importance of habitat patches for overall population connectivity. Models were calibrated based on 1,046 locations of species presence distributed across three biogeographic regions in the Swiss Alps and extrapolated according to two IPCC scenarios of climate change (RCP 4.5 & 8.5), each represented by three downscaled global climate models. The models predicted an average habitat loss of 35% (22%-55%) by 2100, mainly due to an increase in temperature during the reproductive season. An increase in habitat fragmentation was reflected in a 43% decrease in patch size, a 17% increase in the number of habitat patches and a 34% increase in inter-patch distance. However, the predicted changes in habitat availability and connectivity varied considerably between biogeographic regions: Whereas the greatest habitat losses with an increase in inter-patch distance were predicted at the southern and northern edges of the species' Alpine distribution, the greatest increase in patch number and decrease in patch size is expected in the central Swiss Alps. Finally, both the number of isolated habitat patches and the number of patches crucial for maintaining the habitat network increased under the different variants of climate change. Focusing conservation action on the central Swiss Alps may help mitigate the predicted effects of climate change on population connectivity.

Science and Sentiment: Grinnell's Fact-Based Philosophy of Biodiversity Conservation

At the beginning of the twentieth century, the biologist Joseph Grinnell made a distinction between science and sentiment for producing fact-based generalizations on how to conserve biodiversity. We are inspired by Grinnellian science, which successfully produced a century-long impact on studying and conserving biodiversity that runs orthogonal to some familiar philosophical distinctions such as fact versus value, emotion versus reason and basic versus applied science. According to Grinnell, unlike sentiment-based generalizations, a fact-based generalization traces its diverse commitments and thus becomes tractable for its audience. We argue that foregrounding tractability better explains Grinnell's practice in the context of his time as well as in the context of current discourse among scientists over the political "biases" of biodiversity research and its problem of "reproducibility."

Mechanistic variables can enhance predictive models of endotherm distributions: the American pika under current, past, and future climates

How climate constrains species' distributions through time and space is an important question in the context of conservation planning for climate change. Despite increasing awareness of the need to incorporate mechanism into species distribution models (SDMs), mechanistic modeling of endotherm distributions remains limited in this literature. Using the American pika (Ochotona princeps) as an example, we present a framework whereby mechanism can be incorporated into endotherm SDMs. Pika distribution has repeatedly been found to be constrained by warm temperatures, so we used Niche Mapper, a mechanistic heat-balance model, to convert macroclimate data to pika-specific surface activity time in summer across the western United States. We then explored the difference between using a macroclimate predictor (summer temperature) and using a mechanistic predictor (predicted surface activity time) in SDMs. Both approaches accurately predicted pika presences in current and past climate regimes. However, the activity models predicted 8-19% less habitat loss in response to annual temperature increases of ~3-5 °C predicted in the region by 2070, suggesting that pikas may be able to buffer some climate change effects through behavioral thermoregulation that can be captured by mechanistic modeling. Incorporating mechanism added value to the modeling by providing increased confidence in areas where different modeling approaches agreed and providing a range of outcomes in areas of disagreement. It also provided a more proximate variable relating animal distribution to climate, allowing investigations into how unique habitat characteristics and intraspecific phenotypic variation may allow pikas to exist in areas outside those predicted by generic SDMs. Only a small number of easily obtainable data are required to parameterize this mechanistic model for any endotherm, and its use can improve SDM predictions by explicitly modeling a widely applicable direct physiological effect: climate-imposed restrictions on activity. This more complete understanding is necessary to inform climate adaptation actions, management strategies, and conservation plans.

Predictors of Current and Longer-Term Patterns of Abundance of American Pikas (Ochotona princeps) across a Leading-Edge Protected Area

American pikas (Ochotona princeps) have been heralded as indicators of montane-mammal response to contemporary climate change. Pikas no longer occupy the driest and lowest-elevation sites in numerous parts of their geographic range. Conversely, pikas have exhibited higher rates of occupancy and persistence in Rocky Mountain and Sierra Nevada montane 'mainlands'. Research and monitoring efforts on pikas across the western USA have collectively shown the nuance and complexity with which climate will often act on species in diverse topographic and climatic contexts. However, to date no studies have investigated habitat, distribution, and abundance of pikas across hundreds of sites within a remote wilderness area. Additionally, relatively little is known about whether climate acts most strongly on pikas through direct or indirect (e.g., vegetation-mediated) mechanisms. During 2007-2009, we collectively hiked >16,000 km throughout the 410,077-ha Glacier National Park, Montana, USA, in an effort to identify topographic, microrefugial, and vegetative characteristics predictive of pika abundance. We identified 411 apparently pika-suitable habitat patches with binoculars (in situ), and surveyed 314 of them for pika signs. Ranking of alternative logistic-regression models based on AICc scores revealed that short-term pika abundances were positively associated with intermediate elevations, greater cover of mosses, and taller forbs, and decreased each year, for a total decline of 68% during the three-year study; whereas longer-term abundances were associated only with static variables (longitude, elevation, gradient) and were lower on north-facing slopes. Earlier Julian date and time of day of the survey (i.e., midday vs. not) were associated with lower observed pika abundance. We recommend that wildlife monitoring account for this seasonal and diel variation when surveying pikas. Broad-scale information on status and abundance determinants of montane mammals, especially for remote protected areas, is crucial for land and wildlife-resource managers trying to anticipate mammalian responses to climate change.