Predicting abundance of desert riparian birds: validation and calibration of the Effective Area Model.

Reliable prediction of the effects of landscape change on species abundance is critical to land managers who must make frequent, rapid decisions with long-term consequences. However, due to inherent temporal and spatial variability in ecological systems, previous attempts to predict species abundance in novel locations and/or time frames have been largely unsuccessful. The Effective Area Model (EAM) uses change in habitat composition and geometry coupled with response of animals to habitat edges to predict change in species abundance at a landscape scale. Our research goals were to validate EAM abundance predictions in new locations and to develop a calibration framework that enables absolute abundance predictions in novel regions or time frames. For model validation, we compared the EAM to a null model excluding edge effects in terms of accurate prediction of species abundance. The EAM outperformed the null model for 83.3% of species (N=12) for which it was possible to discern a difference when considering 50 validation sites. Likewise, the EAM outperformed the null model when considering subsets of validation sites categorized on the basis of four variables (isolation, presence of water, region, and focal habitat). Additionally, we explored a framework for producing calibrated models to decrease prediction error given inherent temporal and spatial variability in abundance. We calibrated the EAM to new locations using linear regression between observed and predicted abundance with and without additional habitat covariates. We found that model adjustments for unexplained variability in time and space, as well as variability that can be explained by incorporating additional covariates, improved EAM predictions. Calibrated EAM abundance estimates with additional site-level variables explained a significant amount of variability (P < 0.05) in observed abundance for 17 of 20 species, with R2 values >25% for 12 species, >48% for six species, and >60% for four species when considering all predictive models. The calibration framework described in this paper can be used to predict absolute abundance in sites different from those in which data were collected if the target population of sites to which one would like to statistically infer is sampled in a probabilistic way.     

Nest selection by snow petrels Pagodroma nivea in East Antarctica: Validating predictive habitat selection models at the continental scale

Little is known on the factors controlling distribution and abundance of snow petrels in Antarctica. Studying habitat selection through modeling may provide useful information on the relationships between this species and its environment, especially relevant in a climate change context, where habitat availability may change. Validating the predictive capability of habitat selection models with independent data is a vital step in assessing the performance of such models and their potential for predicting species’ distribution in poorly documented areas.From the results of ground surveys conducted in the Casey region (2002–2003, Wilkes Land, East Antarctica), habitat selection models based on a dataset of 4000 nests were created to predict the nesting distribution of snow petrels as a function of topography and substrate. In this study, the Casey models were tested at Mawson, 3800 km away from Casey. The location and characteristics of approximately 7700 snow petrel nests were collected during ground surveys (Summer 2004–2005). Using GIS, predictive maps of nest distribution were produced for the Mawson region with the models derived from the Casey datasets and predictions were compared to the observed data. Models performance was assessed using classification matrixes and Receiver operating characteristic (ROC) curves. Overall correct classification rates for the Casey models varied from 57% to 90%. However, two geomorphologically different sub-regions (coastal islands and inland mountains) were clearly distinguished in terms of habitat selection by Casey model predictions but also by the specific variations in coefficients of terms in new models, derived from the Mawson data sets. Observed variations in the snow petrel aggregations were found to be related to local habitat availability.We discuss the applicability of various types of models (GLM, CT) and investigate the effect of scale on the prediction of snow petrel habitats. While the Casey models created with data collected at the nest scale did not perform well at Mawson due to regional variations in nest micro-characteristics, the predictive performance of models created with data compiled at a coarser scale (habitat units) was satisfactory. Substrate type was the most robust predictor of nest presence between Casey and Mawson. This study demonstrate that it is possible to predict at the large scale the presence of snow petrel nests based on simple predictors such as topography and substrate, which can be obtained from aerial photography. Such methodologies have valuable applications in the management and conservation of this top predator and associated resources and may be applied to other Antarctic, Sub-Antarctic and lower latitudes species and in a variety of habitats.     

Using habitat suitability index and particle dispersion models for early detection of marine invaders.

Eradication and control of invasive species are often possible only if populations are detected when they are small and localized. To be efficient, detection surveys should be targeted at locations where there is the greatest risk of incursions. We examine the utility of habitat suitability index (HSI) and particle dispersion models for targeting sampling for marine pests. Habitat suitability index models are a simple way to identify suitable habitat when species distribution data are lacking. We compared the performance of HSI models with statistical models derived from independent data from New Zealand on the distribution of two nonindigenous bivalves: Theora lubrica and Musculista senhousia. Logistic regression models developed using the HSI scores as predictors of the presence/absence of Theora and Musculista explained 26.7% and 6.2% of the deviance in the data, respectively. Odds ratios for the HSI scores were greater than unity, indicating that they were genuine predictors of the presence/ absence of each species. The fit and predictive accuracy of each logistic model were improved when simulated patterns of dispersion from the nearest port were added as a predictor variable. Nevertheless, the combined model explained, at best, 46.5% of the deviance in the distribution of Theora and correctly predicted 56% of true presences and 50% of all cases. Omission errors were between 6% and 16%. Although statistical distribution models built directly from environmental predictors always outperformed the equivalent HSI models, the gain in model fit and accuracy was modest. High residual deviance in both types of model suggests that the distributions realized by Theora and Musculista in the field data were influenced by factors not explicitly modeled as explanatory variables and by error in the environmental data used to project suitable habitat for the species. Our results highlight the difficulty of accurately predicting the distribution of invasive marine species that exhibit low habitat occupancy and patchy distributions in time and space. Although the HSI and statistical models had utility as predictors of the likely distribution of nonindigenous marine species, the level of spatial accuracy achieved with them may be well below expectations for sensitive surveillance programs.     

Mapping risk for nest predation on a barrier island

Barrier islands and coastal beach systems provide nesting habitat for marine and estuarine turtles. Densely settled coastal areas may subsidize nest predators. Our purpose was to inform conservation by providing a greater understanding of habitat-based risk factors for nest predation, for an estuarine turtle. We expected that habitat conditions at predated nests would differ from random locations at two spatial extents. We developed and validated an island-wide model for the distribution of predated Diamondback terrapin nests using locations of 198 predated nests collected during exhaustive searches at Fisherman Island National Wildlife Refuge, USA. We used aerial photographs to identify all areas of possible nesting habitat and searched each and surrounding environments for nests, collecting location and random-point microhabitat data. We built models for the probability of finding a predated nest using an equal number of random points and validated them with a reserve set (N?=?67). Five variables in 9 a priori models were used and the best selected model (AIC weight 0.98) reflected positive associations with sand patches near marshes and roadways. Model validation had an average capture rate of predated nests of 84.14 % (26.17–97.38 %, Q1 77.53 %, median 88.07 %, Q3 95.08 %). Microhabitat selection results suggest that nests placed at the edges of sand patches adjacent to upland shrub/forest and marsh systems are vulnerable to predation. Forests and marshes provide cover and alternative resources for predators and roadways provide access; a suggestion is to focus nest protection efforts on the edges of dunes, near dense vegetation and roads.     

Predicting the Range of Chinese Mitten Crabs in Europe

Abstract:  Ecological niche modeling provides a means for predicting the potential future distribution of a nonindigenous species based on environmental characteristics of the species' native range. We applied this method to the Chinese mitten crab (Eriocheir sinensis) , a catadromous crustacean with a long history of invasion in Europe. We used genetic algorithm for rule-set prediction to predict the potential European distribution of mitten crab based on its distribution in 42 locations in its native Asia. The climatic variables, air temperature, number of days, amount of precipitation, and wetness index, contributed significantly to predictions of native distribution limits. Although the genetic algorithm for rule-set prediction model was developed for the native range, the species' extensive distribution in Europe ( n = 434) allowed independent validation of the predictions. Application of the model to Europe was successful, with 84% of occurrences in regions predicted to be suitable by >80% of the models and <4% of occurrences in areas predicted suitable by <50% of the models (mainly along the northern range). At the watershed scale, areas with established mitten crab populations had significantly higher habitat matching than sites that were not invaded. The independent validation of the Asian-based model by the European distribution revealed that predictions were highly accurate. The model also identified large areas of Europe, particularly along the Mediterranean coast, as vulnerable to future invasion. These predictions can be used to develop strategies to control the spread of mitten crab by preventing introductions into vulnerable areas.     

Modelling the productivity of naturalised pasture in the North Island,New Zealand: a decision tree approach

Decision tree, one of the data mining methods, has been widely used as a modelling approach and has shown better predictive ability than traditional approaches (e.g. regression). However, very little is known from the literature about how the decision tree performs in predicting pasture productivity. In this study, decision tree models were developed to investigate and predict the annual and seasonal productivity of naturalised hill-pasture in the North Island, New Zealand, and were compared with regression models with respect to model fit, validation and predictive accuracy. The results indicated that the decision tree models for annual and seasonal pasture productivity all had a smaller average squared error (ASE) and a higher percentage of correctly predicted cases than the corresponding regression models. The decision tree model for annual pasture productivity had an ASE which was only half of that of the regression model, and correctly predicted 90% of the cases in the model validation which was 10.8 percentage points higher than that of the regression model. Furthermore, the decision tree models for annual and seasonal pasture productivity also clearly revealed the relative importance of environmental and management variables in influencing pasture productivity, and the interaction among these variables. Spring rainfall was the most significant factor influencing annual pasture productivity, while hill slope was the most significant factor influencing spring and winter pasture productivity, and annual P fertiliser input and autumn rainfall were the most significant factors influencing summer and autumn pasture productivity. One limitation of using the decision tree to predict pasture productivity was that it did not generate a continuous prediction, and thus could not detect the influence of small changes in environmental and management variables on pasture productivity.     

Habitat classification modeling with incomplete data: pushing the habitat envelope.

Habitat classification models (HCMs) are invaluable tools for species conservation, land-use planning, reserve design, and metapopulation assessments, particularly at broad spatial scales. However, species occurrence data are often lacking and typically limited to presence points at broad scales. This lack of absence data precludes the use of many statistical techniques for HCMs. One option is to generate pseudo-absence points so that the many available statistical modeling tools can bb used. Traditional techniques generate pseudo-absence points at random across broadly defined species ranges, often failing to include biological knowledge concerning the species-habitat relationship. We incorporated biological knowledge of the species-habitat relationship into pseudo-absence points by creating habitat envelopes that constrain the region from which points were randomly selected. We define a habitat envelope as an ecological representation of a species, or species feature's (e.g., nest) observed distribution (i.e., realized niche) based on a single attribute, or the spatial intersection of multiple attributes. We created HCMs for Northern Goshawk (Accipiter gentilis atricapillus) nest habitat during the breeding season across Utah forests with extant nest presence points and ecologically based pseudo-absence points using logistic regression. Predictor variables were derived from 30-m USDA Landfire and 250-m Forest Inventory and Analysis (FIA) map products. These habitat-envelope-based models were then compared to null envelope models which use traditional practices for generating pseudo-absences. Models were assessed for fit and predictive capability using metrics such as kappa, threshold-independent receiver operating characteristic (ROC) plots, adjusted deviance (D(adj)2), and cross-validation, and were also assessed for ecological relevance. For all cases, habitat envelope-based models outperformed null envelope models and were more ecologically relevant, suggesting that incorporating biological knowledge into pseudo-absence point generation is a powerful tool for species habitat assessments. Furthermore, given some a priori knowledge of the species-habitat relationship, ecologically based pseudo-absence points can be applied to any species, ecosystem, data resolution, and spatial extent.     

Modeling and Predicting Species Occurrence Using Broad-Scale Environmental Variables: an Example with Butterflies of the Great Basin

Abstract: If occurrence of individual species can be modeled as a function of easily quantified environmental variables (e.g., derived from a geographic information system [GIS]) and the predictions of these models are demonstrably successful, then the scientific foundation for management planning will be strengthened. We used Bayesian logistic regression to develop predictive models for resident butterflies in the central Great Basin of western North America. Species inventory data and values for 14 environmental variables from 49 locations (segments of canyons) in the Toquima Range ( Nevada, U.S.A.) were used to build the models. Squares of the environmental variables were also used to accommodate possibly nonmonotonic responses. We obtained statistically significant models for 36 of 56 (64%) resident species of butterflies. The models explained 8–72% of the deviance in occurrence of those species. Each of the independent variables was significant in at least one model, and squared versions of five variables contributed to models. Elevation was included in more than half of the models. Models included one to four variables; only one variable was significant in about half the models. We conducted preliminary tests of two of our models by using an existing set of data on the occurrence of butterflies in the neighboring Toiyabe Range. We compared conventional logistic classification with posterior probability distributions derived from Bayesian modeling. For the latter, we restricted our predictions to locations with a high ( 70%) probability of predicted presence or absence. We will perform further tests after conducting inventories at new locations in the Toquima Range and nearby Shoshone Mountains, for which we have computed environmental variables by using remotely acquired topographic data, digital-terrain and microclimatic models, and GIS computation.     

Crustacean habitat potential mapping in a tidal flat using remote sensing and GIS

To make a macrofaunal (crustacean) habitat potential map, the spatial distribution of ecological variables in the Hwangdo tidal flat, Korea, was explored. Spatial variables were mapped using remote sensing and a geographic information system (GIS) combined with field observations. A frequency ratio (FR) and logistic regression (LR) model were employed to map the macrofauna potential area for the Ilyoplax dentimerosa, a crustacean species. Spatial variables affecting the tidal macrofauna distribution were selected based on abundance and biomass and used within a spatial database derived from remotely sensed data of various types of sensors. The spatial variables included the intertidal digital elevation model (DEM), slope, distance from a tidal channel, tidal channel density, surface sediment facies, spectral reflectance of the near infrared (NIR) bands and the tidal exposure duration. The relation between the I. dentimerosa and each spatial variable was calculated using the FR and LR. The species was randomly divided into a training set (70%) to analyse habitat potential using FR and LR and a test set (30%) to validate the predicted habitat potential map. The relations were overlaid to produce a habitat potential map with the species potential index (SPI) value for each pixel. The potential habitat maps were compared with the surveyed habitat locations such as validation data set. The comparison results showed that the LR model (accuracy is 85.28%) is better in prediction than the FR (accuracy is 78.96%) model. The performance of models gave satisfactory accuracies. The LR provides the quantitative influence of variables on a potential habitat of species; otherwise, the FR shows the quantitative influence of a class in each variable. The combination of a GIS-based frequency ratio and logistic regression models and remote sensing with field observations is an effective method to determine locations favorable for macrofaunal species occurrences in a tidal flat.     

Predicting barrier passage and habitat suitability for migratory fish species

Fish migrate to spawn, feed, seek refuge from predators, and escape harmful environmental conditions. The success of upstream migration is limited by the presence of barriers that can impede the passage of fish. We used a spatially explicit modeling strategy to examine the effects of barriers on passage for 21 native and non-native migratory fish species and the amount of suitable habitat blocked for each species. Spatially derived physical parameter estimates and literature based fish capabilities and tolerances were used to predict fish passage success and habitat suitability. Both the fish passage and the habitat suitability models accurately predicted fish presence above barriers for most common, non-stocked species. The fish passage model predicted that barriers greater than or equal to 6 m block all migratory species. Chinook salmon (Oncorhynchus tshawytscha) was expected to be blocked the least. The habitat suitability model predicted that low gradient streams with intact habitat quality were likely to support the highest number of fish species. The fish passage and habitat suitability models were intended to be used by environmental managers as strategy development tools to prioritize candidate dams for field assessment and make decisions regarding the management of migratory fish populations.     

Chinook salmon use of spawning patches: relative roles of habitat quality, size, and connectivity.

Declines in many native fish populations have led to reassessments of management goals and shifted priorities from consumptive uses to species preservation. As management has shifted, relevant environmental characteristics have evolved from traditional metrics that described local habitat quality to characterizations of habitat size and connectivity. Despite the implications this shift has for how habitats may be prioritized for conservation, it has been rare to assess the relative importance of these habitat components. We used an information-theoretic approach to select the best models from sets of logistic regressions that linked habitat quality, size, and connectivity to the occurrence of chinook salmon (Oncorhynchus tshawytscha) nests. Spawning distributions were censused annually from 1995 to 2004, and data were complemented with field measurements that described habitat quality in 43 suitable spawning patches across a stream network that drained 1150 km2 in central Idaho. Results indicated that the most plausible models were dominated by measures of habitat size and connectivity, whereas habitat quality was of minor importance. Connectivity was the strongest predictor of nest occurrence, but connectivity interacted with habitat size, which became relatively more important when populations were reduced. Comparison of observed nest distributions to null model predictions confirmed that the habitat size association was driven by a biological mechanism when populations were small, but this association may have been an area-related sampling artifact at higher abundances. The implications for habitat management are that the size and connectivity of existing habitat networks should be maintained whenever possible. In situations where habitat restoration is occurring, expansion of existing areas or creation of new habitats in key areas that increase connectivity may be beneficial. Information about habitat size and connectivity also could be used to strategically prioritize areas for improvement of local habitat quality, with areas not meeting minimum thresholds being deemed inappropriate for pursuit of restoration activities.     

Contrasting ecological information content in whaling archives with modern cetacean surveys for conservation planning and identification of historical distribution changes

Many species are restricted to a marginal or suboptimal fraction of their historical range due to anthropogenic impacts, making it hard to interpret their ecological preferences from modern-day data alone. However, inferring past ecological states is limited by the availability of robust data and biases in historical archives, posing a challenge for policy makers . To highlight how historical records can be used to understand the ecological requirements of threatened species and inform conservation, we investigated sperm whale (Physeter macrocephalus) distribution in the Western Indian Ocean. We assessed differences in information content and habitat suitability predictions based on whale occurrence data from Yankee whaling logs (1792–1912) and from modern cetacean surveys (1995–2020). We built maximum entropy habitat suitability models containing static (bathymetry-derived) variables to compare models comprising historical-only and modern-only data. Using both historical and modern habitat suitability predictions  we assessed marine protected area (MPA) placement by contrasting suitability in- and outside MPAs. The historical model predicted high habitat suitability in shelf and coastal regions near continents and islands, whereas the modern model predicted a less coastal distribution with high habitat suitability more restricted to areas of steep topography. The proportion of high habitat suitability inside versus outside MPAs was higher when applying the historical predictions than the modern predictions, suggesting that different marine spatial planning optimums can be reached from either data sources. Moreover, differences in relative habitat suitability predictions between eras were consistent with the historical depletion of sperm whales from coastal regions, which were easily accessed and targeted by whalers, resulting in a modern distribution limited more to steep continental margins and remote oceanic ridges. The use of historical data can provide important new insights and, through cautious interpretation, inform conservation planning and policy, for example, by identifying refugee species and regions of anticipated population recovery.     

Metacommunity dynamics over 16 years in a pyrogenic shrubland

Metacommunity theory allows predictions about the dynamics of potentially interacting species' assemblages that are linked by dispersal, but strong empirical tests of the theory are rare. We analyzed the metacommunity dynamics of Florida rosemary scrub, a patchily distributed pyrogenic community, to test predictions about turnover rates, community nestedness, and responses to patch size, arrangement, and quality. We collected occurrence data for 45 plant species from 88 rosemary scrub patches in 1989 and 2005 and used growth form, mechanism of regeneration after fire, and degree of habitat specialization to categorize species by life history. We tested whether patch size, fire history, and structural connectivity (a measure of proximity and size of surrounding patches) could be used to predict apparent extinctions and colonizations. In addition, we tested the accuracy of incidence-function models built with the patch survey data from 1989. After fire local extinction rates were higher for herbs than woody plants, higher for species that regenerated only from seed than species able to resprout, and higher for generalist than specialist species. Fewer rosemary specialists and a higher proportion of habitat generalists were extirpated on recently burned patches than on patches not burned between 1989 and 2005. Nestedness was highest for specialists among all life-history groups. Estimated model parameters from 1989 predicted the observed (1989-2005) extinction rates and the number of patches with persistent populations of individual species. These results indicate that species with different life-history strategies within the same metacommunity can have substantially different responses to patch configuration and quality. Real metacommunities may not conform to certain assumptions of simple models, but incidence-function models that consider only patch size, configuration, and quality can have significant predictive accuracy.     

The roles of climate, phylogenetic relatedness, introduction effort, and reproductive traits in the establishment of non-native reptiles and amphibians

We developed a method to predict the potential of non-native reptiles and amphibians (herpetofauna) to establish populations. This method may inform efforts to prevent the introduction of invasive non-native species. We used boosted regression trees to determine whether nine variables influence establishment success of introduced herpetofauna in California and Florida. We used an independent data set to assess model performance. Propagule pressure was the variable most strongly associated with establishment success. Species with short juvenile periods and species with phylogenetically more distant relatives in regional biotas were more likely to establish than species that start breeding later and those that have close relatives. Average climate match (the similarity of climate between native and non-native range) and life form were also important. Frogs and lizards were the taxonomic groups most likely to establish, whereas a much lower proportion of snakes and turtles established. We used results from our best model to compile a spreadsheet-based model for easy use and interpretation. Probability scores obtained from the spreadsheet model were strongly correlated with establishment success as were probabilities predicted for independent data by the boosted regression tree model. However, the error rate for predictions made with independent data was much higher than with cross validation using training data. This difference in predictive power does not preclude use of the model to assess the probability of establishment of herpetofauna because (1) the independent data had no information for two variables (meaning the full predictive capacity of the model could not be realized) and (2) the model structure is consistent with the recent literature on the primary determinants of establishment success for herpetofauna. It may still be difficult to predict the establishment probability of poorly studied taxa, but it is clear that non-native species (especially lizards and frogs) that mature early and come from environments similar to that of the introduction region have the highest probability of establishment.     

Improving effectiveness of systematic conservation planning with density data

Systematic conservation planning aims to design networks of protected areas that meet conservation goals across large landscapes. The optimal design of these conservation networks is most frequently based on the modeled habitat suitability or probability of occurrence of species, despite evidence that model predictions may not be highly correlated with species density. We hypothesized that conservation networks designed using species density distributions more efficiently conserve populations of all species considered than networks designed using probability of occurrence models. To test this hypothesis, we used the Zonation conservation prioritization algorithm to evaluate conservation network designs based on probability of occurrence versus density models for 26 land bird species in the U.S. Pacific Northwest. We assessed the efficacy of each conservation network based on predicted species densities and predicted species diversity. High‐density model Zonation rankings protected more individuals per species when networks protected the highest priority 10‐40% of the landscape. Compared with density‐based models, the occurrence‐based models protected more individuals in the lowest 50% priority areas of the landscape. The 2 approaches conserved species diversity in similar ways: predicted diversity was higher in higher priority locations in both conservation networks. We conclude that both density and probability of occurrence models can be useful for setting conservation priorities but that density‐based models are best suited for identifying the highest priority areas. Developing methods to aggregate species count data from unrelated monitoring efforts and making these data widely available through ecoinformatics portals such as the Avian Knowledge Network will enable species count data to be more widely incorporated into systematic conservation planning efforts.     

Linking occurrence and fitness to persistence: habitat-based approach for endangered greater sage-grouse.

Detailed empirical models predicting both species occurrence and fitness across a landscape are necessary to understand processes related to population persistence. Failure to consider both occurrence and fitness may result in incorrect assessments of habitat importance leading to inappropriate management strategies. We took a two-stage approach to identifying critical nesting and brood-rearing habitat for the endangered Greater Sage-Grouse (Centrocercus urophasianus) in Alberta at a landscape scale. First, we used logistic regression to develop spatial models predicting the relative probability of use (occurrence) for Sage-Grouse nests and broods. Secondly, we used Cox proportional hazards survival models to identify the most risky habitats across the landscape. We combined these two approaches to identify Sage-Grouse habitats that pose minimal risk of failure (source habitats) and attractive sink habitats that pose increased risk (ecological traps). Our models showed that Sage-Grouse select for heterogeneous patches of moderate sagebrush cover (quadratic relationship) and avoid anthropogenic edge habitat for nesting. Nests were more successful in heterogeneous habitats, but nest success was independent of anthropogenic features. Similarly, broods selected heterogeneous high-productivity habitats with sagebrush while avoiding human developments, cultivated cropland, and high densities of oil wells. Chick mortalities tended to occur in proximity to oil and gas developments and along riparian habitats. For nests and broods, respectively, approximately 10% and 5% of the study area was considered source habitat, whereas 19% and 15% of habitat was attractive sink habitat. Limited source habitats appear to be the main reason for poor nest success (39%) and low chick survival (12%). Our habitat models identify areas of protection priority and areas that require immediate management attention to enhance recruitment to secure the viability of this population. This novel approach to habitat-based population viability modeling has merit for many species of concern.     

Predicting species distributions from samples collected along roadsides

Predictive models of species distributions are typically developed with data collected along roads. Roadside sampling may provide a biased (nonrandom) sample; however, it is currently unknown whether roadside sampling limits the accuracy of predictions generated by species distribution models. We tested whether roadside sampling affects the accuracy of predictions generated by species distribution models by using a prospective sampling strategy designed specifically to address this issue. We built models from roadside data and validated model predictions at paired locations on unpaved roads and 200 m away from roads (off road), spatially and temporally independent from the data used for model building. We predicted species distributions of 15 bird species on the basis of point-count data from a landbird monitoring program in Montana and Idaho (U.S.A.). We used hierarchical occupancy models to account for imperfect detection. We expected predictions of species distributions derived from roadside-sampling data would be less accurate when validated with data from off-road sampling than when it was validated with data from roadside sampling and that model accuracy would be differentially affected by whether species were generalists, associated with edges, or associated with interior forest. Model performance measures (kappa, area under the curve of a receiver operating characteristic plot, and true skill statistic) did not differ between model predictions of roadside and off-road distributions of species. Furthermore, performance measures did not differ among edge, generalist, and interior species, despite a difference in vegetation structure along roadsides and off road and that 2 of the 15 species were more likely to occur along roadsides. If the range of environmental gradients is surveyed in roadside-sampling efforts, our results suggest that surveys along unpaved roads can be a valuable, unbiased source of information for species distribution models.