To advance the development of conservation planning for rare species with small geographic ranges, we determined habitat associations of Siskiyou Mountains salamanders (Plethodon stormi) and developed habitat suitability models at Wne (10 ha), medium (40 ha), and broad (202 ha) spatial scales using available Geographic Information Systems data and logistic regression analysis with an information theoretic approach. Across spatial scales, there was very little support for models with structural habitat features, such as tree canopy cover and conifer diameter. Model-averaged 95% conWdence intervals for regression coeYcients and associated odds ratios indicated that the occurrence of Siskiyou Mountains salamanders was positively associated with rocky soils and PaciWc madrone (Abutus menziesii) and negatively associated with elevation and white Wr (Abies concolor); these associations were consistent across 3 spatial scales. The occurrence of this species also was positively associated with hardwood density at the medium spatial scale. Odds ratios projected that a 10% decrease in white Wr abundance would increase the odds of salamander occurrence 3.02-4.47 times, depending on spatial scale. We selected the model with rocky soils, white Wr, and Oregon white oak (Quercus garryana) as the best model across 3 spatial scales and created habitat suitability maps for Siskiyou Mountains salamanders by projecting habitat suitability scores across the landscape. Our habitat suitability models and maps are applicable to selection of priority conservation areas for Siskiyou Mountains salamanders, and our approach can be easily adapted to conservation of other rare species in any geographical location.