EGUsphere

Copernicus Publications

Göttingen, Germany

10.5194/egusphere-2026-4819

Return periods of erosive rain in Germany

Auerswald

Karl

https://orcid.org/0000-0001-5275-4320

¹ ² Fischer

Franziska

³ Winterrath

Tanja

⁴

School of Life Science, Technical University of Munich TUM, Freising, 85354, Germany

Bayerische Landesanstalt für Landwirtschaft, Freising, 85354, Germany

Norwegian Institute of Bioeconomic Research, Ås, 1430, Norway

Deutscher Wetterdienst, Abteilung Hydrometeorologie, Offenbach, 63067, Germany

19 08 2026

2026 1 19

2026

This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/

This article is available from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4819/

The full text article is available as a PDF file from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-4819/egusphere-2026-4819.pdf

Return periods (RPs) of erosive rainfall events are required for erosion risk assessment and for implementing the German Soil Protection Law, which defines farmers' responsibility for erosion damage expected at least once within 10 years. However, no established method exists to quantify the RP of individual erosive rainfall events. Here, we develop and evaluate a framework for estimating the conditional distribution of rainfall erosivities, called erosion indices (EI), to quantify RPs of erosive rainfall events in Germany. The analysis combines 26,566 erosive events from 115 rain gauge (RG) stations and 1,007,437 from weather radar (RR) measurements covering Germany between 2001 and 2017. Generalised Extreme Value distributions were fitted to compare measurement methods, while empirical cumulative distribution functions were used to calculate RPs directly because RP estimates were highly sensitive to small deviations in exceedance probabilities. The number of erosive events increased systematically with the long-term mean annual erosivity factor R. Radar-based event frequencies followed a third-order polynomial relationship with R, whereas RG data showed an approximately linear relationship within the smaller range of observed R factors covering only the lower half of the RR data. Differences between RG and RR event frequencies mainly affected small erosive events (EI < 10 N h‑1), whereas larger events relevant to RP estimation were captured similarly by both methods after applying published scaling corrections. Across Germany, the return level EI for a given RP could be described by a unified equation relating EIRP to both RP and the local R factor. For a 10-year RP, the resulting relationship was EI10yr = 3.46 × R0.59. The EI10yr increased from 40 N h‑1 to 100 N h‑1 when the R factor increased from 40 N h‑1 a‑1 to 500 N h‑1 a‑1. Validation against RG-derived estimates yielded a root mean squared error of 8.3 N h‑1 despite the datasets' contrasting spatial and temporal resolutions and measurement methods. The results demonstrate that RPs for erosive rainfall events can be quantified consistently using radar-derived erosivity data, adjusted for local climatic conditions via the R factor. The R factor also captured the large interannual variability and thus can likely predict climate change effects on event erosivities. The derived equations provide a practical basis for erosion risk assessment, near-real-time erosion forecasting, and evaluation of legal responsibility for erosion damage under changing climatic conditions. The study further shows that conventional rainfall depth–duration–frequency relationships cannot substitute for erosivity-based RPs because rainfall erosivity exhibits a fundamentally different dependence on event duration and rainfall intensity.