We propose a probabilistic model to determine landslide hazard at the basin scale. The model predicts where landslides will occur, how frequently they will occur, and how large they will be. We test the model in the Staffora River basin, in the northern Apennines, Italy. For the study area, we prepare a multi-temporal inventory map through the interpretation of multiple sets of aerial photographs taken between 1955 and 1999. We partition the basin into 2243 geo-morpho-hydrological units, and obtain the probability of spatial occurrence of landslides by discriminant analysis of thematic variables, including morphological, lithological, structural and land use. For each mapping unit, we obtain the landslide recurrence by dividing the total number of landslide events inventoried in the unit by the time span of the investigated period. Assuming that landslide recurrence will remain the same in the future, and adopting a Poisson probability model, we determine the exceedance probability of having one or more landslides in each mapping unit, for different periods. We obtain the probability of landslide size by analysing the frequency-area statistics of landslides, obtained from the multi-temporal inventory map. Assuming independence, we obtain a quantitative estimate of landslide hazard for each mapping unit as the joint probability of landslide size, of landslide temporal occurrence and of landslide spatial occurrence.