Abstract. In the present paper, the rainfall erosivity factor (R factor) for the area of the Czech Republic is assessed. Based on 10 min data for 96 stations and corresponding R factor estimates, a number of spatial interpolation methods are applied and cross-validated. These methods include inverse distance weighting, standard, ordinary, and regression kriging with parameters estimated by the method of moments and restricted maximum likelihood, and a generalized least-squares (GLS) model. For the regression-based methods, various statistics of monthly precipitation as well as geographical indices are considered as covariates. In addition to the uncertainty originating from spatial interpolation, the uncertainty due to estimation of the rainfall kinetic energy (needed for calculation of the R factor) as well as the effect of record length and spatial coverage are also addressed. Finally, the contribution of each source of uncertainty is quantified. The average R factor for the area of the Czech Republic is 640 MJ ha−1 mm h−1, with values for the individual stations ranging between 320 and 1520 MJ ha−1 mm h−1. Among various spatial interpolation methods, the GLS model relating the R factor to the altitude, longitude, mean precipitation, and mean fraction of precipitation above the 95th percentile of monthly precipitation performed best. Application of the GLS model also reduced the uncertainty due to the record length, which is substantial when the R factor is estimated for individual sites. Our results revealed that reasonable estimates of the R factor can be obtained even from relatively short records (15–20 years), provided sufficient spatial coverage and covariates are available.