Evaluating the impacts of sustainable land management practices on water quality in an agricultural catchment in Lower Austria using SWAT.

Managing agricultural watersheds in an environmentally friendly manner necessitate the strategic implementation of well-targeted sustainable land management (SLM) practices that limit soil and nonpoint source pollution losses and translocation. Watershed-scale SLM-scenario modeling has the potential to identify efficient and effective management strategies from the field to the integrated landscape level. In a case study targeting a 66-hectare watershed in Petzenkirchen, Lower Austria, the Soil and Water Assessment Tool (SWAT) was utilized to evaluate a variety of locally adoptable SLM practices. SWAT was calibrated and validated (monthly) at the catchment outlet for flow, sediment, nitrate-nitrogen (NO3-N), ammonium nitrogen (NH4-N), and mineralized phosphorus (PO4-P) using SWATplusR. Considering the locally existing agricultural practices and socioeconomic and environmental factors of the research area, four conservation practices were evaluated: baseline scenario, contour farming (CF), winter cover crops (CC), and a combination of no-till and cover crops (NT + CC). The NT + CC SLM practice was found to be the most effective soil conservation practice in reducing soil loss by around 80%, whereas CF obtained the best results for decreasing the nutrient loads of NO3-N and PO4-P by 11% and 35%, respectively. The findings of this study imply that the setup SWAT model can serve the context-specific performance assessment and eventual promotion of SLM interventions that mitigate on-site land degradation and the consequential off-site environmental pollution resulting from agricultural nonpoint sources.

Effect of land use-land cover change on the regimes of surface runoff-the case of Lake Basaka catchment (Ethiopia).

In this study, the effect of land use-land cover change (LULCC) on surface (direct) runoff was estimated for Lake Basaka catchment using the soil conservation services-curve number model in the geospatial information system (ArcInfo), assisted by remote sensing. The result indicated that Lake Basaka catchment experienced a significant LULCC. About 86% of forest coverage and 46% of grasslands were lost over the study period (1973-2015), which were shifted to open bushy woodlands, farms, lake water and wetlands. The runoff responses were observed to be increasing since 1970s, especially after the inception of large-scale irrigation schemes to the region. The highest increase of surface runoff was observed to occur after mid-1980s, which is in line with the significant LULCC and the corresponding increment of lake level in that period. The reduction in vegetation cover has resulted in an increase of runoff coefficient (rc) from 0.07 in the 1960s to about 0.23 in 2000s. The sensitivity analysis result indicated that about 70% of the increase runoff rate in the lake catchment is attributed to LULCC, and the remaining proportion is due to rainfall. However, the effect of extreme rainfall on runoff process could not be underemphasized since it has significant impact especially during extreme events (observed rc of 0.33 in 2008). Overall, when predicting the runoff response of the lake catchment, it is importance to take into account possible future LULCC and evolution.

Long-term data from field erosion plot studies in eastern Austria.

Soil erosion and runoff data are collected at three sites in eastern Austria using field erosion plots. Observed treatments include 1) conventional tillage with plough (CT), 2) mulch tillage with winter cover crop (MT), and 3) no-till with winter cover crop (NT). Data cover a time span from 1994 to 2018. They include data about surface runoff, soil loss, nitrogen, phosphorus and soil organic carbon losses as well as cop yields associated with the erosion processes. Interpretation of the data will be found in "Long-term experience with conservation tillage practices in Austria: impacts on soil erosion processes" [1].

Comparison of three types of laser optical disdrometers under natural rainfall conditions.

Optical disdrometers can be used to estimate rainfall erosivity; however, the relative accuracy of different disdrometers is unclear. This study compared three types of optical laser-based disdrometers to quantify differences in measured rainfall characteristics and to develop correction factors for kinetic energy (KE). Two identical PWS100 (Campbell Scientific), one Laser Precipitation Monitor (Thies Clima) and a first-generation Parsivel (OTT) were collocated with a weighing rain gauge (OTT Pluvio2) at a site in Austria. All disdrometers underestimated total rainfall compared to the rain gauge with relative biases from 2% to 29%. Differences in drop size distribution and velocity resulted in different KE estimates. By applying a linear regression to the KE-intensity relationship of each disdrometer, a correction factor for KE between the disdrometers was developed. This factor ranged from 1.15 to 1.36 and allowed comparison of KE between different disdrometer types despite differences in measured drop size and velocity.

Mapping monthly rainfall erosivity in Europe.

Rainfall erosivity as a dynamic factor of soil loss by water erosion is modelled intra-annually for the first time at European scale. The development of Rainfall Erosivity Database at European Scale (REDES) and its 2015 update with the extension to monthly component allowed to develop monthly and seasonal R-factor maps and assess rainfall erosivity both spatially and temporally. During winter months, significant rainfall erosivity is present only in part of the Mediterranean countries. A sudden increase of erosivity occurs in major part of European Union (except Mediterranean basin, western part of Britain and Ireland) in May and the highest values are registered during summer months. Starting from September, R-factor has a decreasing trend. The mean rainfall erosivity in summer is almost 4 times higher (315MJmmha-1h-1) compared to winter (87MJmmha-1h-1). The Cubist model has been selected among various statistical models to perform the spatial interpolation due to its excellent performance, ability to model non-linearity and interpretability. The monthly prediction is an order more difficult than the annual one as it is limited by the number of covariates and, for consistency, the sum of all months has to be close to annual erosivity. The performance of the Cubist models proved to be generally high, resulting in R2 values between 0.40 and 0.64 in cross-validation. The obtained months show an increasing trend of erosivity occurring from winter to summer starting from western to Eastern Europe. The maps also show a clear delineation of areas with different erosivity seasonal patterns, whose spatial outline was evidenced by cluster analysis. The monthly erosivity maps can be used to develop composite indicators that map both intra-annual variability and concentration of erosive events. Consequently, spatio-temporal mapping of rainfall erosivity permits to identify the months and the areas with highest risk of soil loss where conservation measures should be applied in different seasons of the year.

Global rainfall erosivity assessment based on high-temporal resolution rainfall records.

The exposure of the Earth's surface to the energetic input of rainfall is one of the key factors controlling water erosion. While water erosion is identified as the most serious cause of soil degradation globally, global patterns of rainfall erosivity remain poorly quantified and estimates have large uncertainties. This hampers the implementation of effective soil degradation mitigation and restoration strategies. Quantifying rainfall erosivity is challenging as it requires high temporal resolution(<30 min) and high fidelity rainfall recordings. We present the results of an extensive global data collection effort whereby we estimated rainfall erosivity for 3,625 stations covering 63 countries. This first ever Global Rainfall Erosivity Database was used to develop a global erosivity map at 30 arc-seconds(~1 km) based on a Gaussian Process Regression(GPR). Globally, the mean rainfall erosivity was estimated to be 2,190 MJ mm ha-1 h-1 yr-1, with the highest values in South America and the Caribbean countries, Central east Africa and South east Asia. The lowest values are mainly found in Canada, the Russian Federation, Northern Europe, Northern Africa and the Middle East. The tropical climate zone has the highest mean rainfall erosivity followed by the temperate whereas the lowest mean was estimated in the cold climate zone.

Rainfall erosivity in Europe.

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.

Reply to the comment on "Rainfall erosivity in Europe" by Auerswald et al.

Recently, in the Auerswald et al. (2015) comment on "Rainfall erosivity in Europe", 5 criticisms were addressed: i) the neglect of seasonal erosion indices, ii) the neglect of published studies and data, iii) the low temporal resolution of the data, especially of the maximum rain intensity, iv) the use of precipitation data instead of rain data and the subsequent deviation of the R-factor in Germany and Austria compared with previous studies, and v) the differences in considered time periods between countries. We reply as follows: (i) An evaluation of the seasonal erosion index at the European scale is, to our knowledge, not achievable at present with the available data but would be a future goal. Synchronous publication of the seasonal erosion index is not mandatory, specifically because seasonal soil loss ratios are not available at this scale to date. We are looking forward to the appropriate study by the authors of the comment, who assert that they have access to the required data. (ii) We discuss and evaluate relevant studies in our original work and in this reply; however, we cannot consider what is not available to the scientific community. (iii) The third point of critique was based on a misunderstanding by Auerswald et al. (2015), as we did indeed calculate the maximum intensity with the highest resolution of data available. (iv) The low R-factor values in Germany and the higher values in Austria compared with previous studies are not due to the involvement of snow but are rather due to a Pan-European interpolation. We argue that an interpolation across the borders of Austria creates a more reliable data set. (v) We agree that the use of a short time series or time series from different periods is generally a problem in all large-scale studies and requires improvement in the future. However, because this affects countries with a rather low variability of the R-factor in our study, we are confident that the overall results of the map are not biased. In conclusion, the Pan-European rainfall data compilation (REDES) was a great success and yielded data from 1541 stations with an average length of 17.1years and a temporal resolution of <60min. However, a Pan-European data collection will never be complete without the help and supply of data from its users. Thus, we invite the authors of the comment to share their data in the open REDES to help build even better rainfall-erosivity maps at regional or European scales.