[Objective] Rainfall erosivity is an important indicator for measuring the potential ability of rainfall to cause soil erosion, which is crucial for understanding surface soil erosion, soil degradation, and ecological protection. It is of great significance for predicting soil erosion processes, assessing regional risks, and preventing soil erosion. [Method] This study adopted a simple algorithm model based on daily rainfall, and combined Mann Kendall method and wavelet analysis method to conduct long-term data analysis on 67 stations in the hilly and gully area of the Loess Plateau, exploring their trends and periodic changes. At the same time, spatial interpolation method was used to reveal the spatial distribution characteristics of rainfall erosivity. [Results] The average annual rainfall erosivity range in this region is 739.97-1749.03MJ . mm . hm-2 . h-1, decreasing spatially from southeast to northwest. In terms of interannual variation, the overall erosive power of rainfall shows a weak increasing trend, with the highest annual value being 1915.09MJ mm hm_-2 h^_-1 in 2013 and the lowest annual value being 741.18 MJ mm hm_-2 h^_-1 in 1974. The intensity of rainfall erosivity in the warm season is significantly stronger than that in the cold season, with an average annual rainfall erosivity of 599.01-1742.95 MJ mm hm_-2 h^_-1 in the warm season, while it is only 27.61-171.11 MJ mm hm_-2 h_-1 ^{in the cold season. [}Conclusion] Trend testing shows that there has been no significant change in rainfall erosivity, but there are periodic variations of 43 years and 8 years.