Soil erosion by water is a serious threat for the Mediterranean region. Raindrop impacts and consequent runoff generation are the main driving forces of this geomorphic process of soil degradation. The potential ability for rainfall to cause soil loss is expressed as rainfall erosivity, a key parameter required by most soil loss prediction models. In Italy, rainfall erosivity measurements are limited to few locations, preventing researchers from effectively assessing the geography and magnitude of soil loss across the country. The objectives of this study were to investigate the spatio-temporal distribution of rainfall erosivity in Italy and to develop a national-scale grid-based map of rainfall erosivity. Thus, annual rainfall erosivity values were measured and subsequently interpolated using a geostatistical approach. Time series of pluviographic records (10-years) with high temporal resolution (mostly 30-min) for 386 meteorological stations were analysed. Regression-kriging was used to interpolate rainfall erosivity values of the meteorological stations to an Italian rainfall erosivity map (500-m). A set of 23 environmental covariates was tested, of which seven covariates were selected based on a stepwise approach (mostly significant at the 0.01 level). The interpolation method showed a good performance for both the cross-validation data set ( = 0.777) and the fitting data set (R2 = 0.779).