A robust method for spatial prediction of landslide hazard in roaded and roadless areas of forest is described. The method is based on assigning digital terrain attributes into continuous landform classes. The continuous landform classification is achieved by applying a fuzzy k-means approach to a watershed scale area before the classification is extrapolated to a broader region. The extrapolated fuzzy landform classes and datasets of road-related and non road-related landslides are then combined in a geographic information system (GIS) for the exploration of predictive correlations and model development. In particular, a Bayesian probabilistic modeling approach is illustrated using a case study of the Clearwater National Forest (CNF) in central Idaho, which experienced significant and widespread landslide events in recent years. The computed landslide hazard potential is presented on probabilistic maps for roaded and roadless areas. The maps can be used as a decision support tool in forest planning involving the maintenance, obliteration or development of new forest roads in steep mountainous terrain.