Towards the Reproducibility in Soil Erosion Modeling: a New Pan- European Soil Erosion Map

Soil erosion by water is a widespread phenomenon throughout Europe and has the potentiality, with his on-site and off-site effects, to affect water quality, food security and floods. Despite the implementation of numerous and different models to estimate soil erosion by water in Europe, there is still a lack of harmonization of assessment methodologies. Very often, different approaches give soil erosion rates significantly different and even when the same model is applied to the same region the results may be different. This can be due to the way the model is implemented (i.e. with the selection of different algorithms when available) and/or to the use of datasets having different resolution or accuracy. Scientific computation is emerging as one of the central topic of the scientific method, to overcome these problems there is thus the necessity to develop reproducible computational method where codes and data are available. The present study is an illustration of such an approach. Using only public available datasets, we applied the Revised Universal Soil loss Equation (RUSLE) to locate the most sensitive areas to soil erosion in Europe. A significant effort was made to select the better simplified equations to be used when a strict application of the RUSLE model is not possible. In particular for the computation of the Rainfall Erosivity factor (R) the reproducible research paradigm was applied. The calculation of the R factor was implemented using public datasets and the GNU R language and an easily reproducible validation procedure based on measured precipitation time series was applied using Matlab language. Designing the computational modeling architecture with the aim to ease as much as possible the future reuse of the model in analyzing climate change scenarios is also a challenging goal of the research.