The prediction of forest productivity is essential for sustainable forest management, particularly in countries, like Italy, where forest is an important part of many protected areas. A spatial predictive probability model for forest productivity rates in Italy was developed over the period 1961-1990, based on 135 annually-resolved records of site productivity and auxiliary variables measured at 219 stations. Our analysis shows that the probability of finding high (> 7.3 m 3 ha −1 yr −1) and low (< 5.8 m 3 ha −1 yr −1) productivity rates changes across different regions of Italy. The generated spatial patterns contribute to a better understanding of the factors structuring the distribution of forest productivity in Italy because they reflect the dependence of temperature and water availability conditions on the latitudinal and altitudinal location of the study areas. We observed that the temperature control dominates forest productivity at high elevations and latitudes, whereas low-elevation sites in central and southern Italy are more sensitive to water availability. The proposed spatial probability modelling should be further assessed for its possible incorporation into forest management plans.