skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 1, 2026

Title: Combining past and contemporary species occurrences with ordinal species distribution modeling to investigate responses to climate change
Many organisms leave evidence of their former occurrence, such as scat, abandoned burrows, middens, ancient eDNA or fossils, which indicate areas from which a species has since disappeared. However, combining this evidence with contemporary occurrences within a single modeling framework remains challenging. Traditional binary species‐distribution modeling reduces occurrence to two temporally coarse states (present/absent), so thus cannot leverage the information inherent in temporal sequences of evidence of past occurrence. In contrast, ordinal modeling can use the natural time‐varying order of states (e.g. never occupied versus previously occupied versus currently occupied) to provide greater insights into range shifts. We demonstrate the power of ordinal modeling for identifying the major influences of biogeographic and climatic variables on current and past occupancy of the American pikaOchotona princeps, a climate‐sensitive mammal. Sampling over five years across the species' southernmost, warm‐edge range limit, we tested the effects of these variables at 570 habitat patches where occurrence was classified either as binary or ordinal. The two analyses produced different top models and predictors – ordinal modeling highlighted chronic cold as the most‐important predictor of occurrence, whereas binary modeling indicated primacy of average summer‐long temperatures. Colder wintertime temperatures were associated in ordinal models with higher likelihood of occurrence, which we hypothesize reflect longer retention of insulative and meltwater‐provisioning snowpacks. Our binary results mirrored those of other past pika investigations employing binary analysis, wherein warmer temperatures decrease likelihood of occurrence. Because both ordinal‐ and binary‐analysis top models included climatic and biogeographic factors, results constitute important considerations for climate‐adaptation planning. Cross‐time evidences of species occurrences remain underutilized for assessing responses to climate change. Compared to multi‐state occupancy modeling, which presumes all states occur in the same time period, ordinal models enable use of historical evidence of species' occurrence to identify factors driving species' distributions more finely across time.  more » « less
Award ID(s):
2022036
PAR ID:
10621587
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Ecography
Date Published:
Journal Name:
Ecography
Volume:
2025
Issue:
2
ISSN:
0906-7590
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Long-term monitoring of habitat occupancy can reveal patterns of habitat use, population dynamics, and factors controlling species distribution. The American pika (Ochotona princeps), a small mammal found in rocky habitats throughout western North America, has been targeted for occupancy studies due to its relatively conspicuous behavior and its unusual adaptations for surviving long, cold winters without hibernation. These adaptations include an unusually high resting metabolic rate and maintenance of body temperatures near the lethal maximum for this species, which would appear to compromise the pika's ability to survive warmer summers. Recent monitoring as well as projections based on future climate scenarios have suggested this species is experiencing a period of range retraction due to warming summers and/or loss of insulating winter snow cover. Niwot Ridge is situated ideally to test competing hypotheses about the trajectory and drivers of pika range shift. The pika is still common throughout the Colorado Rockies, but published models differ markedly regarding projections of the pika’s future distribution in this region. Niwot Ridge has experienced warmer summers as well as shorter periods of insulating snow cover in recent years, and there is evidence that pikas are now less common than they once were in at least one area on the ridge. This study is designed to provide robust data on pika population trends through long-term monitoring of occupancy in a spatially balanced random sample of pika habitat patches centered on Niwot Ridge. Survey plots (n = 72) were selected according to a Generalized Random-Tessellation Stratified (GRTS) algorithm, stratified dichotomously by elevation, average annual snow accumulation (SWE), and probabilities of pika occurrence based on previous data. Each plot extends 12 m in radius from a GRTS point. To ensure that each plot contains at least 10% cover of talus, plot coordinates were adjusted (usually less than 50 m) or replaced using the GRTS oversample to select the next available and suitable plot within the same categories of elevation, SWE and probability of occurrence (see "pika-survey-GRTS-plot-tracking-record.cr.data.csv" for plot strata, survey schedules, GRTS sequence, and records of plot replacement or location adjustments). Trained technicians survey plots for pikas and fresh pika sign (food caches and fecal pellets) as well as metrics of habitat quality. Each year, 48 of the 72 plots are surveyed in a rotating panel design (24 plots are surveyed annually, 24 in even years and 24 in odd years). Plots are surveyed in August when pikas are engaged in food caching and other conspicuous behaviors related to territory establishment and defense. Data collected at each plot are detailed in a survey manual ("pika_survey.cr.methods.docx"). Each plot is outfitted with a data logger (sensor) to record sub-surface temperature several times each day. Photos of plot and sensor locations are used in navigation and sensor retrieval. Each survey is completed during a brief (half-hour) visit to the plot to service the sensor and to record habitat and pika data. A subset of plots (n = 12) are selected for double surveys each year to allow estimation of pika detection probability. Estimates of detection probability are also informed by data on time to detection of pikas and pika sign recorded during each survey. Samples of fresh pika fecal pellets are collected from occupied plots and are stored as vouchers of pika presence and for use in studies of population genetics and physiology, including studies of physiological stress in relation to habitat quality and microclimate. 
    more » « less
  2. Long-term monitoring of habitat occupancy can reveal patterns of habitat use, population dynamics, and factors controlling species distribution. The American pika (Ochotona princeps), a small mammal found in rocky habitats throughout western North America, has been targeted for occupancy studies due to its relatively conspicuous behavior and its unusual adaptations for surviving long, cold winters without hibernation. These adaptations include an unusually high resting metabolic rate and maintenance of body temperatures near the lethal maximum for this species, which would appear to compromise the pika's ability to survive warmer summers. Recent monitoring as well as projections based on future climate scenarios have suggested this species is experiencing a period of range retraction due to warming summers and/or loss of insulating winter snow cover. Niwot Ridge is situated ideally to test competing hypotheses about the trajectory and drivers of pika range shift. The pika is still common throughout the Colorado Rockies, but published models differ markedly regarding projections of the pika’s future distribution in this region. Niwot Ridge has experienced warmer summers as well as shorter periods of insulating snow cover in recent years, and there is evidence that pikas are now less common than they once were in at least one area on the ridge. This study is designed to provide robust data on pika population trends through long-term monitoring of occupancy in a spatially balanced random sample of pika habitat patches centered on Niwot Ridge. Survey plots (n = 72) were selected according to a Generalized Random-Tessellation Stratified (GRTS) algorithm, stratified dichotomously by elevation, average annual snow accumulation (SWE), and probabilities of pika occurrence based on previous data. Each plot extends 12 m in radius from a GRTS point. To ensure that each plot contains at least 10% cover of talus, plot coordinates were adjusted (usually less than 50 m) or replaced using the GRTS oversample to select the next available and suitable plot within the same categories of elevation, SWE and probability of occurrence (see "pika-survey-GRTS-plot-tracking-record.cr.data.csv" for plot strata, survey schedules, GRTS sequence, and records of plot replacement or location adjustments). Trained technicians survey plots for pikas and fresh pika sign (food caches and fecal pellets) as well as metrics of habitat quality. Each year, 48 of the 72 plots are surveyed in a rotating panel design (24 plots are surveyed annually, 24 in even years and 24 in odd years). Plots are surveyed in August when pikas are engaged in food caching and other conspicuous behaviors related to territory establishment and defense. Data collected at each plot are detailed in a survey manual ("pika_survey.cr.methods.docx"). Each plot is outfitted with a data logger (sensor) to record sub-surface temperature several times each day. Photos of plot and sensor locations are used in navigation and sensor retrieval. Each survey is completed during a brief (half-hour) visit to the plot to service the sensor and to record habitat and pika data. A subset of plots (n = 12) are selected for double surveys each year to allow estimation of pika detection probability. Estimates of detection probability are also informed by data on time to detection of pikas and pika sign recorded during each survey. Samples of fresh pika fecal pellets are collected from occupied plots and are stored as vouchers of pika presence and for use in studies of population genetics and physiology, including studies of physiological stress in relation to habitat quality and microclimate. 
    more » « less
  3. Abstract Rapid climate change is threatening biodiversity via habitat loss, range shifts, increases in invasive species, novel species interactions, and other unforeseen changes. Coastal and estuarine species are especially vulnerable to the impacts of climate change due to sea level rise and may be severely impacted in the next several decades. Species distribution modeling can project the potential future distributions of species under scenarios of climate change using bioclimatic data and georeferenced occurrence data. However, models projecting suitable habitat into the future are impossible to ground truth. One solution is to develop species distribution models for the present and project them to periods in the recent past where distributions are known to test model performance before making projections into the future. Here, we develop models using abiotic environmental variables to quantify the current suitable habitat available to eight Neotropical coastal species: four mangrove species and four salt marsh species. Using a novel model validation approach that leverages newly available monthly climatic data from 1960 to 2018, we project these niche models into two time periods in the recent past (i.e., within the past half century) when either mangrove or salt marsh dominance was documented via other data sources. Models were hindcast‐validated and then used to project the suitable habitat of all species at four time periods in the future under a model of climate change. For all future time periods, the projected suitable habitat of mangrove species decreased, and suitable habitat declined more severely in salt marsh species. 
    more » « less
  4. Abstract AimEcological niches shape species commonness and rarity, yet, the relative importance of different niche mechanisms within and across ecosystems remains unresolved. We tested the influence of niche breadth (range of environmental conditions where species occur) and niche position (marginality of a species’ environmental distribution relative to the mean environmental conditions of a region) on tree‐species abundance and occupancy across three biogeographic regions. LocationArgentinian Andes; Bolivian Amazon; Missouri Ozarks. Time period2002–2010. Major taxa studiedTrees. MethodsWe calculated abiotic‐niche breadths and abiotic‐niche positions using 16 climate, soil and topographic variables. For each region, we used model selection to test the relative influence of niche breadth and niche position on local abundance and occupancy in regional‐scale networks of 0.1‐ha forest plots. To account for species–environment associations caused by other mechanisms (e.g., dispersal), we used null models that randomized associations between species occurrences and environmental variables. ResultsWe found strong support for the niche‐position hypothesis. In all regions, species with higher local abundance and occupancy occurred in non‐marginal environments. Observed relationships between occupancy and niche position also differed from random species–environment associations in all regions. Surprisingly, we found little support for the niche‐breadth hypothesis. Observed relationships between both local abundance and niche breadth, and occupancy and niche breadth, did not differ from random species–environment associations. Main conclusionNiche position was more important than niche breadth in shaping species commonness and rarity across temperate, sub‐tropical and tropical forests. In all forests, tree species with widespread geographic distributions were associated with environmental conditions commonly found throughout the region, suggesting that niche position has similar effects on species occupancy across contrasting biogeographic regions. Our findings imply that conservation efforts aimed at protecting populations of common and rare tree species should prioritize conservation of both common and rare habitats. 
    more » « less
  5. Environmental conditions are dynamic, and plants respond to those dynamics on multiple time scales. Disequilibrium occurs when a response occurs more slowly than the driving environmental changes. We review evidence regarding disequilibrium in plant distributions, including their responses to paleoclimate changes, recent climate change and new species introductions. There is strong evidence that plant species distributions are often in some disequilibrium with their environmental conditions.This disequilibrium poses a challenge when projecting future species distributions using species distribution models (SDMs). Classically, SDMs assume that the set of species occurrences is an unbiased sample of the suitable environmental conditions. However, a species in disequilibrium with the environment may have higher‐than‐expected occurrence probabilities (e.g. due to extinction debts) or lower‐than‐expected occurrence probabilities (e.g. due to dispersal limitation) in different areas. If unaccounted for, this will lead to biased estimates of the environmental suitability.We review methods for avoiding such biases in SDMs, ranging from simple thinning of the occurrence dataset to complex dynamic and process‐based models. Such models require large data inputs, natural history knowledge and technical expertise, so implementing them can be challenging. Despite this, we advocate for their increased use, since process‐based models provide the best potential to account for biases in model training data and to then represent the dynamics of species occupancy as ranges shift.Synthesis. Occurrence records for a species are often in disequilibrium with climate. SDMs trained on such data will produce biased estimates of a species' niche unless this disequilibrium is addressed in the modelling. A range of tools, spanning a wide gradient of complexity and realism, can resolve this bias. 
    more » « less