Abstract: The investigation of the spatiotemporal heterogeneity of rainfall erosivity is significant in regional soil erosion risk assessment as well as water and soil conservation planning. Based on long-period(1961-2016) daily rainfall data from 63 meteorological stations located in Qin-Ba mountains region and its surrounding area, this study calculated values of rainfall erosivity and analyzed the spatial and temporal variation of rainfall erosivity using Kriging spatial interpolation method, Mann-kendall trend test, Pettitt test and other methods. Results showed that the regional average annual rainfall erosivity was 3 696 MJ·mm/(hm²·h·a). The rainfall erosivity showed a unimodal intra-annual variation with a maximum value in July accounting for 26.6 % of the total year and also had a strong seasonal pattern with maximum rainfall erosivity in summer and minimum in winter. The rainfall erosivity reached its maximum in the 1980s and minimum in the 1990s. The annual rainfall erosivity had distinct phases, whereas had no significant trends and mutations. It turned out that the spatial distribution of the average annual rainfall erosivity in the region increased from north to south varying from 787 to 8 858 MJ·mm/(hm² ·h·a). Furthermore, the average annual rainfall erosivity of Qin-Ba mountains region generally decreased with the increasing of latitude(P<0.01) and altitude(P<0.01).