Rainfall is one the main drivers of soil erosion. The erosive force of rainfall is expressed as rainfall erosivity. Rainfall erosivity considers the rainfall amount and intensity, and is most commonly expressed as the R-factor in the USLE model and its revised version, RUSLE. At national and continental levels, the scarce availability of data obliges soil erosion modellers to estimate this factor based on rainfall data with only low temporal resolution (daily, monthly, annual averages). The purpose of this study is to assess rainfall erosivity in Europe in the form of the RUSLE R-factor, based on the best available datasets. Data have been collected from 1541 precipitation stations in all European Union (EU) Member States and Switzerland, with temporal resolutions of 5 to 60 min. The R-factor values calculated from precipitation data of different temporal resolutions were normalised to R-factor values with temporal resolutions of 30 min using linear regression functions. Precipitation time series ranged from a minimum of 5 years to a maximum of 40 years. The average time series per precipitation station is around 17.1 years, the most datasets including the first decade of the 21st century. Gaussian Process Regression (GPR) has been used to interpolate the R-factor station values to a European rainfall erosivity map at 1 km resolution. The covariates used for the R-factor interpolation were climatic data (total precipitation, seasonal precipitation, precipitation of driest/wettest months, average temperature), elevation and latitude/longitude. The mean R-factor for the EU plus Switzerland is 722 MJ mm ha(-1) h(-1) yr(-1), with the highest values (>1000 MJ mm ha(-1) h(-1) yr(-1)) in the Mediterranean and alpine regions and the lowest (<500 MJ mm ha(-1) h(-1) yr(-1)) in the Nordic countries. The erosivity density (erosivity normalised to annual precipitation amounts) was also the highest in Mediterranean regions which implies high risk for erosive events and floods.

In the region of the European Alps, national and regional meteorological services operate rain‐gauge networks, which together, constitute one of the densest in situ observation systems in a large‐scale high‐mountain region. Data from these networks are consistently analyzed, in this study, to develop a pan‐Alpine grid dataset and to describe the region's mesoscale precipitation climate, including the occurrence of heavy precipitation and long dry periods. The analyses are based on a collation of high‐resolution rain‐gauge data from seven Alpine countries, with 5500 measurements per day on average, spanning the period 1971–2008. The dataset is an update of an earlier version with improved data density and more thorough quality control. The grid dataset has a grid spacing of 5 km, daily time resolution, and was constructed with a distance‐angular weighting scheme that integrates climatological precipitation–topography relationships. Scales effectively resolved in the dataset are coarser than the grid spacing and vary in time and space, depending on station density. We quantify the uncertainty of the dataset by cross‐validation and in relation to topographic complexity, data density and season. Results indicate that grid point estimates are systematically underestimated (overestimated) at large (small) precipitation intensities, when they are interpreted as point estimates. Our climatological analyses highlight interesting variations in indicators of daily precipitation that deviate from the pattern and course of mean precipitation and illustrate the complex role of topography. The daily Alpine precipitation grid dataset was developed as part of the EU funded EURO4M project and is freely available for scientific use.

Combined Global Surface Summary of Day and European Climate Assessment and Dataset daily meteorological data sets (around 9000 stations) were used to build spatio‐temporal geostatistical models and predict daily air temperature at ground resolution of 1 km for the global land mass. Predictions in space and time were made for the mean, maximum, and minimum temperatures using spatio‐temporal regression‐kriging with a time series of Moderate Resolution Imaging Spectroradiometer (MODIS) 8 day images, topographic layers (digital elevation model and topographic wetness index), and a geometric temperature trend as covariates. The accuracy of predicting daily temperatures was assessed using leave‐one‐out cross validation. To account for geographical point clustering of station data and get a more representative cross‐validation accuracy, predicted values were aggregated to blocks of land of size 500×500 km. Results show that the average accuracy for predicting mean, maximum, and minimum daily temperatures is root‐mean‐square error (RMSE) =±2°C for areas densely covered with stations and between ±2°C and ±4°C for areas with lower station density. The lowest prediction accuracy was observed at high altitudes (>1000 m) and in Antarctica with an RMSE around 6°C. The model and predictions were built for the year 2011 only, but the same methodology could be extended for the whole range of the MODIS land surface temperature images (2001 to today), i.e., to produce global archives of daily temperatures (a next‐generation http://WorldClim.org repository) and to feed various global environmental models.

A computational framework to generate daily temperature maps using time-series of publicly available MODIS MOD11A2 product Land Surface Temperature (LST) images (1 km resolution; 8-day composites) is illustrated using temperature measurements from the national network of meteorological stations (159) in Croatia. The input data set contains 57,282 ground measurements of daily temperature for the year 2008. Temperature was modeled as a function of latitude, longitude, distance from the sea, elevation, time, insolation, and the MODIS LST images. The original rasters were first converted to principal components to reduce noise and filter missing pixels in the LST images. The residual were next analyzed for spatio-temporal auto-correlation; sum-metric separable variograms were fitted to account for zonal and geometric space-time anisotropy. The final predictions were generated for time-slices of a 3D space-time cube, constructed in the R environment for statistical computing. The results show that the space-time regression model can explain a significant part of the variation in station-data (84%). MODIS LST 8-day (cloud-free) images are unbiased estimator of the daily temperature, but with relatively low precision (±4.1°C); however their added value is that they systematically improve detection of local changes in land surface temperature due to local meteorological conditions and/or active heat sources (urban areas, land cover classes). The results of 10–fold cross-validation show that use of spatio-temporal regression-kriging and incorporation of time-series of remote sensing images leads to significantly more accurate maps of temperature than if plain spatial techniques were used. The average (global) accuracy of mapping temperature was ±2.4°C. The regression-kriging explained 91% of variability in daily temperatures, compared to 44% for ordinary kriging. Further software advancement—interactive space-time variogram exploration and automated retrieval, resampling and filtering of MODIS images—are anticipated.

As a follow up and an advancement of the recently published Rainfall Erosivity Database at European Scale (REDES) and the respective mean annual R-factor map, the monthly aspect of rainfall erosivity has been added to REDES. Rainfall erosivity is crucial to be considered at a monthly resolution, for the optimization of land management (seasonal variation of vegetation cover and agricultural support practices) as well as natural hazard protection (landslides and flood prediction). We expanded REDES by 140 rainfall stations, thus covering areas where monthly R-factor values were missing (Slovakia, Poland) or former data density was not satisfactory (Austria, France, and Spain). The different time resolutions (from 5 to 60 min) of high temporal data require a conversion of monthly R-factor based on a pool of stations with available data at all time resolutions. Because the conversion factors show smaller monthly variability in winter (January: 1.54) than in summer (August: 2.13), applying conversion factors on a monthly basis is suggested. The estimated monthly conversion factors allow transferring the R-factor to the desired time resolution at a European scale. The June to September period contributes to 53% of the annual rainfall erosivity in Europe, with different spatial and temporal patterns depending on the region. The study also investigated the heterogeneous seasonal patterns in different regions of Europe: on average, the Northern and Central European countries exhibit the largest R-factor values in summer, while the Southern European countries do so from October to January. In almost all countries (excluding Ireland, United Kingdom and North France), the seasonal variability of rainfall erosivity is high. Very few areas (mainly located in Spain and France) show the largest from February to April. The average monthly erosivity density is very large in August (1.67) and July (1.63), while very small in January and February (0.37). This study addresses the need to develop monthly calibration factors for seasonal estimation of rainfall erosivity and presents the spatial patterns of monthly rainfall erosivity in European Union and Switzerland. Moreover, the study presents the regions and seasons under threat of rainfall erosivity.