Abstract
The largest currently compiled database of plot runoff and soil loss data in Europe and the Mediterranean was analysed to investigate effects of land use on annual soil loss (SL), annual runoff (R) and annual runoff coefficient (RC). This database comprises 227 plot-measuring sites in Europe and the Mediterranean, with SL for 1056 plots (PL) representing 7024 plot-years (PY) and R for 804 PL representing 5327 PY. Despite large data variability, continental-wide trends are observed. Construction sites have the highest mean annual RC (57%) and SL (325 Mg.ha-1.yr-1). Bare soil, vineyards and tree crops have high mean annual RC (5–10%) and SL (10–20 Mg.ha-1.yr-1). Cropland and fallow show similar mean annual RC (8.0 and 7.3%), but lower SL (6.5 and 5.8 Mg.ha-1.yr-1). Plots with (semi-)natural vegetation cover show lowest mean annual RC (<5%) and SL (<1 Mg.ha-1.yr-1). Plot length and slope gradient correlations with R and SL depend on land-use type and are not concurrent for R and SL. Most land-use types show positive correlations between annual R and SL. Plots in cold climates have higher annual RC than plots in temperate and pan-Mediterranean climates. Annual SL in the pan-Mediterranean is less than in temperate zones, due to stony or clayey soils having a low erodibility. Annual RC in the pan-Mediterranean was higher than in temperate zones. Annual R increases strongly with increasing annual precipitation (P) above 500 mm.yr-1, while annual SL was found to stabilize at P > 500 mm.yr-1. For shrubland, annual SL was found to decrease for P > 250–500 mm.yr-1, which is attributed to an accompanying increase in vegetation cover. However, no such trend was found for R. The results allow a rapid assessment of the impact of land-use changes on annual R, RC and SL, based on field-measured plot data.
Keywords
I Introduction
Runoff and soil loss due to interrill and rill erosion are important processes of soil degradation that cause significant on-site and off-site problems (e.g. Boardman and Poesen, 2006; Montgomery, 2007; Poesen and Hooke, 1997). An integrated approach to these problems at subcontinental scale requires runoff and soil loss to be assessed for a wide range of representative environmental conditions. An extensive assessment and mapping approach for large areas (e.g. Evans, 2002; Le Gouée et al., 2010; Oldeman et al., 1991) may provide an overview of the scale of the problem and locate erosion hotspots, but to gain insight into these processes and to develop strategies to mitigate their impacts, more detailed field-measured experimental data that accurately quantify soil loss are needed. Recently, such quantitative assessment of annual soil loss (SL, Mg.ha-1.yr-1) rates at a pan-continental scale for Europe has seen several applications like risk assessment mapping through modelling (Kirkby et al., 2004), exploration of spatial variability and controlling factors of annual SL rates (Cerdan et al., 2006, 2010), assessment of scale effects on sediment production (Vanmaercke et al., 2012) and the development of indicator systems to identify and monitor problem areas (e.g. Gobin et al., 2004). All these continental-wide applications either directly make use of available field-measured soil loss data or conclude they would benefit from a validation with such data.
While aforementioned studies place strong emphasis on the assessment of annual SL rates, runoff plays an important role as a causal factor of annual SL and the relations between annual runoff (R, mm.yr-1), annual runoff coefficients (RC, %) and annual SL are not yet fully understood quantitatively. Nevertheless, several process-based erosion models use runoff in order to estimate annual SL rates (Merritt et al., 2003) and good knowledge on these relations is an important part of soil loss modelling (e.g. Jetten and Favis-Mortlock, 2006; Kinnell and Risse, 1998). Furthermore, annual R and RC also directly relate to on-site problems like agricultural productivity (Rockström et al., 2010) and off-site problems like export of nutrients and pesticides (Rossi Pisa et al., 1999), flash floods (e.g. García-Ruiz et al., 2010; Poesen and Hooke, 1997) and the potential activation of other sediment sources (such as river banks and gullies) further downstream (Vanmaercke et al., 2011a). Hence, it is important to assess both annual R and SL rates in conjunction, as well as the effect of key controlling factors on annual R and SL rates.
Over the last 60 years, numerous quantitative experimental studies on annual R and SL have been conducted throughout Europe and the Mediterranean region, using different experimental methods (e.g. runoff plots, rainfall simulations, rill volume measurements, tracer methods). From these, bounded runoff plot studies under natural rainfall conditions can be considered as the most used and standardized experimental method (e.g. Boix-Fayos et al., 2006; Cammeraat, 1993; Cerdan et al., 2006; Evans, 1995; Hudson, 1993). Studies on runoff and soil loss plots have been extensively used in large-scale coordinated research projects in the United States, leading to the development of the (R)USLE(2) equation (Renard et al., 1997; Wischmeier and Smith, 1960, 1978). However, projects of such an extent have not taken place in Europe where many individual runoff and soil loss plot studies have been reported. As a result, the findings of runoff plot studies in Europe and the Mediterranean region are dispersed over numerous scientific papers, reports and theses. They are mostly designed to analyse effects of erosion controlling factors on annual R and SL in a particular area. These individual studies have provided better insights into runoff and soil loss processes at local scales, but the diversity and natural variability of runoff and soil loss plot studies limit the potential to extrapolate these findings to other environmental conditions. Boardman (1998) showed that published soil erosion rates for large areas, based on a small number of observed data, should be interpreted with care or may not be relevant at all. Furthermore, an overall assessment of runoff and soil loss rates is also hampered by a high temporal variability (e.g. Bagarello et al., 2011; Martínez-Casasnovas et al., 2002; Ollesch and Vacca, 2002).
From this lack of an overview and difficulties in extrapolating local R and SL data to larger areas arises the need for a pan-European compilation of all available R and SL data. As a response to this need, national-scale data sets on soil erosion have recently been assembled for most countries in Europe, although the methodology used and the erosion processes considered differ (Baade and Rekolainen, 2006; Boardman and Poesen, 2006). With respect to annual R and SL, several recent studies have compiled field-measured plot runoff and/or soil loss data at the regional, national or subcontinental scale (Boardman and Poesen, 2006; Table 1).
Overview of studies reporting measured runoff and soil loss plot data without soil and water conservation techniques, collected at country, regional or subcontinental scale for Europe and the Mediterranean. Key: – = not reported in study, + = reported in study, NA = not available, R = annual runoff, SL = annual soil loss, #PY = number of plot-years. Ba = bare, Cr = cropland, Fa = fallow, Fo = forest, Gr = grassland, Pf = post-fire, Sh = shrubland, Tc = tree crops, Vi = vineyard.
From these compilations, a better insight into some key factors determining rates and variability of annual R and SL at (sub)continental and regional scales has been gained. The dominant control of land-use type on annual SL was illustrated by Cerdan et al. (2006, 2010), where also soil type, plot length and slope gradient were used to account for further variability. Kosmas et al. (1997) found that the relation between annual precipitation (P, mm.yr-1) and annual R and SL at plot scale at eight different sites in the northern Mediterranean is mainly influenced by land use and hence temporal and spatial patterns of vegetation cover. For shrubland, these authors also observed a vegetation feedback mechanism whereby with increasing annual P (up to P = 200–300 mm.yr-1) SL first increases, then decreases with increasing annual P. This effect is similar to the one described by Langbein and Schumm (1958), who demonstrated for catchments in the United States that the relation between annual P and catchment sediment yield does not only reflect the increasing erosion potential of higher annual P but also includes feedback effects from a larger vegetation biomass with increasing annual precipitation, effectively reducing sediment yield above an annual P threshold of about 254 to 381 mm.yr-1 (10–15 inch.yr-1). For the other land uses studied by Kosmas et al. (1997), annual R and SL were all positively related to annual P.
However, several key elements to obtain a comprehensive understanding of annual R, RC and SL rates and controls in Europe and the Mediterranean region are not fully considered in these studies. First and foremost, the existing overviews mainly consider annual SL, while annual R and RC are studied to a lesser extent (Table 1). While Europe-wide assessments of plot-scale annual SL exist (Cerdan et al., 2010; Vanmaercke et al., 2012), this is not the case for annual R. In addition, most of these studies assess the effects of one or more controlling factors on annual R and SL, but recognize that they lack information on other important controlling factors. As such, Cerdan et al. (2006, 2010) acknowledge the importance of annual precipitation but did not include this in their analysis. Similarly, Kosmas et al. (1997) do not assess the effect of plot length or slope gradient but nevertheless cite its importance. Which controlling factors are included in the analysis, and how they are assessed, may have a significant impact on the discussion of annual R and SL rates. For instance, Fleskens and Stroosnijder (2007) argued that average annual SL rates in olive groves are unlikely to exceed 10 Mg.ha-1.yr-1, which was contested by Gómez et al. (2008) on the basis that several scale and environmental factors like plot length were not sufficiently taken into account.
Hence, there are still unresolved questions with respect to the relationships between precipitation, runoff (coefficients) and soil loss, and the effect of different land uses on these relations. While several local-scale studies address part of these questions, it is not known whether the findings of these studies also apply to other regions. At a continental-wide scale, there are no studies that comprehensively and quantitatively explore all of these relations. Nevertheless, such analysis is of great use to support model assumptions and contribute to an integrated approach towards soil degradation. Therefore, this study aims (1) to provide an overview of both R and SL rates by interrill and rill erosion, measured at the plot scale for Europe and the Mediterranean region; (2) to assess the variability in both R and SL rates for different land uses and different climatic regions in the study area; and (3) to analyse the relationship between observed R and SL rates, and their relationship to annual precipitation (P, mm.yr-1).
For a list of abbrevations and symbols used in this paper, see the Appendix.
II Materials and methods
1 Runoff plot selection criteria and database compilation
A database was constructed with annual R and SL data, measured on bounded runoff plots under natural rainfall conditions in Europe and the Mediterranean region (Figure 1). Data were collected from scientific papers, books (Boardman and Poesen, 2006), project reports and PhD theses, and through personal communication with various researchers. Only runoff and/or soil loss measurements conducted on bounded runoff plots under natural rainfall conditions with a known land use, a minimum plot length of 5 m, and for a measuring period (MP) that is representative (see section II(2) below for more details) for at least one year were considered. Results from runoff plots that were treated with soil and water conservation techniques were not included in this analysis since they do not represent prevailing field conditions. A list of the excluded soil and water conservation techniques is given in Table 2a. While plots without any soil cover throughout the year, i.e. bare plots (Table 2b), are not a common land use practice, they are often used in soil erosion studies as reference plots, representing maximum potential SL for the study conditions. However, traditional agricultural practices on cropland often leave the soil bare after tillage during part of the year. Hence, plots with bare soil were included in the database for reference purposes.

Geographical distribution of plot runoff and soil loss measuring sites over Europe and the Mediterranean with indication of the climatic zones derived from the LANMAP2 classification (Metzger et al., 2005; Mücher et al., 2010). Inset (a): Canary Islands. #PL = number of plots, n = number of plot-measuring sites.
List of soil and water conservation techniques (SWCT) tested on runoff plots and thus excluded from the plot database used for the analysis of the effects of land use on annual runoff and annual soil loss in Europe and the Mediterranean. Ba = bare, Cr = cropland, Fa = fallow, Fo = forest, Sh = shrubland, Ra = rangeland, Pf = post-fire, Gr = grassland, Vi = vineyards, Tc = tree crops, Cs = construction sites.
Land-use types considered in the plot database for Europe and the Mediterranean
While most plots use collection tanks or flow samplers to determine the total runoff and soil loss by interrill and rill erosion, in a small number of studies, soil loss is determined by measuring the rill volume (Govers and Poesen, 1988; Jankauskas and Fullen, 2002; Jankauskas and Jankauskiene, 2003a; Jankauskas et al., 2004). Based on a literature survey by Govers and Poesen (1988), total soil loss for these studies was calculated by adding 25% to the measured rill soil loss to account for interrill soil loss.
Each runoff and soil loss plot (PL) in the database represents measurements of annual runoff and/or soil loss at a particular measuring site for a specific combination of land use, soil type, plot length and slope gradient. For each of these PL, the corresponding number of plot-years (PY) was also recorded, which indicates the number of years represented by the data for that PL. PY = 1 corresponds to a measuring period of one year on a single runoff plot. When measurements were conducted on several replicate plots with identical experimental setup and results were reported individually for the different replicates, they were included as different PL in the database. If the average annual R and SL value for the replicate plots was reported (e.g. Bagarello and Ferro, 2010; Bagarello et al., 2010a, 2010b; Lopes et al., 2002; Mohammad and Adam, 2010) the average values for all replicates were counted as one PL in the database, while the number of PY was considered to be the sum of all PY of the replicates.
Based on the description given by the authors, all PL were assigned to a land-use type (Table 2b). When different land-use types were present on a single runoff plot (e.g. cropland-fallow rotation or cropland-grassland rotation), the years having the same land use were grouped together as one PL. Hence, data from a runoff plot with a rotation of cropland and fallow is entered in the database as 2 PL, one for cropland and one for fallow. This approach does not take into account possible effects of crop cultivation and soil treatment prior to the measuring period, which may persist during the two following years (Wischmeier and Smith, 1978) and hence may explain part of the observed variability in R and SL rates. Furthermore, the location of each of the plot-measuring sites was determined, either from coordinates given by the authors or from maps and descriptions in the publications. Some studies report data for plots that were installed at different sites close to each other (e.g. different slope aspects or gradients in the same valley). Whenever it was possible to accurately distinguish between these different sites, they were incorporated as separate plot-measuring sites in the database but if this was impossible, the location of the study area where plots were located was included as one plot-measuring site. Subsequently, the climatic zone (CZ) of each of the plot-measuring sites was determined according to the LANMAP2 classification (Metzger et al., 2005; Mücher et al., 2010). Plot-measuring sites in the Near East and northern Africa that fall outside the LANMAP2 cover area were classified as Mediterranean (Figure 1). While soil properties like texture, rock fragment cover, soil organic matter content and soil erodibility are recognized as important determinants of runoff and soil loss (e.g. Cerdan, 2006, 2010; Poesen and Lavee, 1994; Poesen et al., 1994; Sanchis et al., 2008; Torri et al., 1997), quantitative data on these properties were not systematically reported in the literature from which plot data were extracted for the database. Hence, a quantitative analysis of the effect of these soil properties on annual R and SL could not be made.
2 Annual data and data extrapolation procedures
All data used in this analysis are annual data. More than 84% of #PL (corresponding to >90% of #PY) reported annual R and SL that were obtained during a measuring period of one or more years and are reported as annual values by the author. In some studies, however, measurements were not carried out for full years, but the authors indicated that the data were representative for full years because no or negligible runoff and soil loss occurred during the period that the runoff plots were not in operation. This was the case for some studies during the dry season in the Mediterranean region (e.g. Mohammad and Adam, 2010; Roxo et al., 1996) or during permanent snow cover and frozen soil in colder climates (e.g. Fulajtár and Janský, 2001). In these cases, the measuring period was considered to be full years when calculating the corresponding PY for that PL. Other studies where the MP is shorter than 1 year were only included in the plot database when the authors explicitly report data as annual (extrapolation by the authors), or a reasonable extrapolation could be made. This was only done if measurements were conducted for a period during which at least two-thirds of the annual rainfall depth was recorded and rainfall is distributed uniformly throughout the year. In these cases annual R and SL were estimated by linear extrapolation according to the corresponding annual P. If no annual P data were available, data were extrapolated linearly to annual values according to the corresponding number of days if measurements continued for a period of at least 80% of the year and rainfall is distributed uniformly throughout the year. Uniformity of rainfall was checked visually using the long-term average rainfall distribution for the plot-measuring site, as given by the New_LocClim program (FAO, 2006).
3 Database analyses
Plot length, slope gradient and soil characteristics need to be taken into account to explain variability in annual R, RC and SL. Previous studies indicated that the relationship between plot length or slope gradient and annual SL or R is non-linear (Cerdan et al., 2010; Nearing, 1997; Poesen and Bryan, 1989; Wischmeier, 1966; Wischmeier and Smith, 1978).
The effect of plot length and slope gradient on annual R and SL was assessed by calculating the non-parametric Spearman’s rank correlation coefficients (rs). For annual SL, also the correlations with the RUSLE length factor (LPL, equation 1; Renard et al., 1997), a slope factor (SPL, equation 2; Nearing, 1997) and the product of the LPL and SPL (LS-factor, equation 3) were calculated. For all PL with a land-use type for which annual SL was significantly correlated to the LS-factor, the annual unit plot soil loss (SLu) was calculated using equation 4 (Bagarello et al., 2010b).
where LPL= plot length factor, λ = plot length (m), SPL = plot slope gradient factor, Ө = plot slope angle (°), LSPL= plot LS-factor, SLu = annual unit plot soil loss (Mg.ha-1.yr-1).
As no standard procedure for the correction of annual R and RC values to a unit plot exists, no correction factor could be calculated for annual R and RC. Regardless of other factors, it can be expected that the reliability of average annual R and SL measurements increases with increasing PY, as they capture more of the occurring natural variability. According to the central limit theorem, the standard error of a sample is inversely related to the square root of the number of observations (Tijms, 2004). Therefore, average annual R, RC and SL values were calculated by weighting the reported mean value of each PL with the square root of the number of plot-years. While this is a very basic approach and does not take into account the complex temporal variation in annual SL and often non-normal distributions of annual SL time series (Maetens et al., 2011), it can be expected that this weighting results in more reliable estimates of average annual R, RC and SL rates (Cerdan et al., 2010; Vanmaercke et al., 2011b). A k-sample Kolmogorov-Smirnov test (KS-test) was used to test for significant differences in the distribution of annual R, RC, SL and SLu data between each combination of two land-use types, whereby the significance level of the test was adapted using Bonferroni’s Inequality to account for familywise error in multiple corrections (Brittain, 1987). The same procedure was applied to test for significant differences between CZ for each land-use type. To further explore the runoff-soil loss relationship, regression equations of the form:
were calculated, where the dependent variable Y was either annual SL or SLu, the independent variable X was either annual R or RC and a and b are empirical constants. Both unweighted and regressions weighted according to the square root of the number of plot-years of each plot were calculated.
Precipitation has been studied intensively as the causal factor of runoff and soil loss. The rainfall erosivity index (EI30) has been proposed as a good measure for relating precipitation to (potential) annual SL (Wischmeier and Smith, 1958). While annual P was reported in most of the studies and is generally widely available with detailed spatial and temporal coverage, EI30 values were not consistently reported for the PL included in the database and are not available from other data sources (Gabriëls, 2006). As an alternative, the Modified Fournier Index (Arnoldus, 1980; equation 6) has been used as a measure of climatic erosivity (Gabriëls, 2006):
where MFI = Modified Fournier Index (mm2.mm-1), p = average monthly precipitation (mm.month-1), P = average annual precipitation (mm.yr-1). MFI was calculated for each of the plot measurement sites, using monthly precipitation data obtained from the CRU CL 2.0 data set (New et al., 2002).
III Results
1 Description of the plot database
The plot database contains data from 227 plot-measuring sites throughout Europe and the Mediterranean region (Figure 1), compiled from 213 individual publications. Annual SL data are available for a total of 1056 PL, corresponding to a total of 7204 PY. Annual R data were available for 804 PL (5327 PY). For 766 of these PL, representing 5013 PY, both annual R and SL data are available. For 673 PL (corresponding to 4583 PY), both R and P are reported, allowing the calculation of annual RC. The distribution of PL and PY over the different countries in the study area and the references to the data sources are given in Table 3.
Overview by country in Europe and the Mediterranean of the number of plots (#PL), number of plot-years (#PY) and sources included in the plot database
The first recorded soil loss measurements in the database started in 1950 at Cean-Turda, Romania (Motoc et al., 1998) and later on the number of plots increased until 1994, after which the number of plots started to decline. The earliest publication discussing plot measurements that could be found dated from 1968 (Dubber, 1968), although most of the publications date from 1986 onwards (Figure 2). The average MP of all runoff and soil loss plots is 6.0 yrs (median: 4 yrs, mode: 1 yr) with a minimum of 1 and a maximum of 42 yrs at Podu-Iloaiei , Romania (Bucur et al., 2007) (Figure 3). For most PL (>84% of #PL and >90% of #PY), measurements continued throughout the year for at least one full year. Annual R and SL data calculated for a MP less than 1 yr which were extrapolated to a full year account for less than 3% of all PL and les than 1% of PY (Table 4). The distribution of the number of plots and plot-years over the different land-use types and CZ is given in Table 5, while the distributions of the number of plots and plot-years according to annual P, plot length and plot slope gradients are given in Figure 4. The different land-use types in the database also show different frequency distributions (Figure 5), with cropland and fallow mainly occurring on gentler slopes (<20%), while forest, vineyards, construction sites and tree crop plots are generally situated on steeper slopes (>20%).

Evolution of the total number of plots (# plots) in operation per year for which annual soil loss and/or runoff were recorded in Europe and the Mediterranean as well as the total number of publications from which plot data were extracted in this study (# publications) per publication year

Frequency distribution of the number of plots (# plots) as a function of the measuring period (MP) for which annual runoff or annual soil loss were measured continuously in Europe and the Mediterranean

Frequency distribution of the number of plots (# PL) and the number of plot-years (# PY) in Europe and the Mediterranean for which annual runoff and/or annual runoff coefficient, and/or annual soil loss data are available with respect to: (a) plot length; (b) plot slope gradient; and (c) annual precipitation (P). Total #PL = 1 096, total #PY = 7 533. NA = not available (i.e. plot length, slope gradient or P not reported).

Frequency distribution of the % of plots for the different land-use types in Europe and the Mediterranean according to plot slope gradient class. n = total number of plots for a given land-use type.
Data collection period of plot data in the plot database. R = annual runoff, RC = annual runoff coefficient, SL = annual soil loss. Number of plots (#PL) and number of plot-years (#PY) for which plot data were collected (1) during full years (Full Year), (2) during a representative part of the year during which almost all of the annual rain was recorded and the authors considered the data to be representative for a full year (Repr. for Full Year), (3) during a measuring period less than one year during which at least 67% of the annual rain was recorded and for which the data were linearly extrapolated to 100% of the annual rainfall (Extrapol. to Full Year).
Overview of the number of plots (#PL) and the number of plot-years (#PY) for which annual runoff (R), annual runoff coefficient (RC) and/or annual soil loss (SL) data for Europe and the Mediterranean are available. NA = not available.
2 Effects of plot length and slope gradient on annual runoff, runoff coefficient and soil loss
The range of annual R and SL values recorded in the plot database varies over almost four orders of magnitude (Figure 6). To allow a comparison of the results obtained on plots with different plot lengths and slope gradients, the importance of these topographic variables was examined using the correlation coefficient between plot length and slope gradient and R and RC (Table 6). Annual R was found to be positively correlated with plot length for plots having bare soil, cropland, grassland and tree crops. A significant negative relationship between annual R and plot length was found for plots under forest and plots recently affected by fire. With respect to the relation between slope gradient and annual R only a significant positive correlation was found for post-fire conditions. For grassland and shrubland the relation between R and plot slope gradient was found to be significantly negative. With the exception of vineyards, where a significant negative correlation between slope gradient and annual RC was observed, the same trends, albeit with slightly different rs values, were found for annual RC.

Frequency distribution of annual runoff (R), annual runoff coefficients (RC), annual soil loss (SL) and annual unit plot soil loss (SLu, equation 4) rates in the plot database for Europe and the Mediterranean. #PL = total number of plots.
Spearman’s rank correlation coefficient (rs) and p-values for the correlation between annual runoff (R), annual runoff coefficients (RC) and plot length and slope gradient on the one hand and annual soil loss (SL) and USLE plot length-factor (Lpl), slope gradient-factor (Spl) and LS-factor (LSpl) on the other hand for different land-use types in Europe and the Mediterranean. Values in bold indicate significance at α=0.05. NA = not available.
For annual SL, a significant positive correlation with the LS-factor (equation 3) was found for bare plots, cropland, fallow, shrubland and tree crops (Table 6). For cropland, the correlation was significant with both the LPL and SPL, while for bare and tree crops only a significant correlation with LPL was found. Nevertheless, the correlation between annual SL and SPL is also relatively strong. For shrubland, only SPL was significantly correlated with annual SL, while for vineyards there was only a significant correlation with LPL. Therefore, unit plot SL (SLu, equation 4) was calculated for plots where the land use was bare, cropland, fallow, shrubland or tree crops. Although the correlation between annual SL and LSPL was not significant for construction sites, this is likely due to the limited number of plots for this land-use type. Nevertheless, annual SLu values were also calculated for PL on construction sites since these PL all have a bare soil surface.
3 Characteristics of the frequency distributions of annual runoff and soil loss for various land uses
Weighted mean values and box-plots indicating the range of annual R, RC, SL and SLu by land-use type are shown in Figure 7. Weighted mean values are always higher than the median value for annual R, RC and SL, which indicates that annual R, RC and SL have a positively skewed distribution, such as the log-normal distribution (e.g. Bagarello et al., 2010a). Construction sites have consistently the highest annual R, RC and SL. After correction for the plot length and slope gradient, annual SLu values for construction sites are smaller, which is attributed to the steep slopes associated with the construction site plots (Figure 5). Nevertheless, annual SLu for construction sites remains considerably higher than for other land-use types. However, the number of plots on construction sites is low, so the corresponding mean values for annual R, RC and SL rates should be interpreted with caution, although they do indicate a high vulnerability to erosion of this land-use type. This is attributed to a presence of bare, disturbed soil with a low structural stability on relatively steep slopes causing very high soil loss rates (Borselli et al., 2006). The remainder of the land-use types can be divided into two groups. A first group consisting of bare plots or plots with some type of crop cultivation (i.e. cropland, fallow, tree crops and vineyards) shows mean annual R rates between 30 and 60 mm.yr-1, annual RC rates between 5 and 15% and annual SL rates between 1 and 20 Mg.ha-1.yr-1. The second group of land-use types consists of PL with a (semi-)natural vegetation cover (i.e. forest, post-fire, shrubland, rangeland and grassland), with annual R, RC and SL rates less than 30 mm.yr-1, 5% and 1 Mg.ha-1.yr-1, respectively (Table 7). While bare plots and plots with crop cultivation have consistently higher mean annual R, RC and SL rates than plots with (semi-)natural vegetation, the ranking of the weighted mean values for the different land-use types within these two groups varies for annual R, RC and SL (Figure 7). The difference between these two groups was confirmed by the application of the Kolmogorov-Smirnov test, whereby significant differences between land-use types were mostly observed between land-use types where crop cultivation was applied and land-use types with (semi-)natural vegetation. Furthermore, for 25 out of 55 pairwise combinations of all the 11 land-use types, the distribution of annual SL data was found to be significantly different, while for annual R and RC only 12 and 8 combinations were found to be significantly different, respectively.

Frequency distribution of mean annual runoff (R), mean annual runoff coefficient (RC), mean annual soil loss (SL), and mean annual unit plot soil loss (SLu, equation 4) of all plots in Europe and the Mediterranean, grouped by land-use type. Box-plots in grey indicate mean annual soil losses scaled to a unit plot (SLu: plot length = 22.1 m, plot slope gradient = 9%). For the other land-use types, the original, measured SL was used. Black dots indicate the mean R, RC, SL and SLu, weighted by number of plot-years (#PY). Box plots are ordered according to descending weighted mean.
Weighted mean, median and coefficient of variation (CV) for annual runoff (R), annual runoff coefficient (RC), annual soil loss (SL) and annual unit plot soil loss (SLu) by land-use type for all data and grouped by climatic zone in Europe and the Mediterranean. Weighting was done according to the square root of the number of plot-years. For number of plots and plot-years corresponding to each land use, see Table 5. NC = not calculated because the relation between LS-factor and soil loss for this land-use type was insignificant (see Table 6). NA = no data available.
4 Relationships of annual runoff coefficients and soil loss for various land uses
The median and distributions of all available annual RC and SL data for different land-use types are shown in Figure 8. This figure clearly indicates that there is a general trend towards higher median annual SL rates with higher median annual RC. For land uses under crop cultivation, tree crops have the highest median annual RC and SL. Vineyards have comparable median annual RC and 75th percentile annual SL values. However, median annual SL is markedly lower for vineyards, indicating that high annual SL rates in the database are more rare for vineyards than for tree crops. The interquartile ranges for cropland and fallow plots are similar for annual RC and annual SL. Median annual RC is lower for fallow plots, however, which indicates that low annual RC are more prevalent on fallow plots. Although annual RC rates on bare and cropland plots are similar, bare plots show higher annual SL rates than cropland plots. Plots with (semi-)natural vegetation are characterized by consistently low median annual SL (0.08–0.3 Mg.ha-1.yr-1), but still show somewhat more variation in median annual RC rates (0.5–1.1%).

Relation between median annual runoff coefficient (RC), median annual soil loss (SL) and median annual unit plot soil loss (SLu) of all plots in Europe and the Mediterranean by land use with indication of the 25th and 75th percentile for RC and SL, based on the plot database (for details, see Tables 5 and 7). n = number of land-use types.
Using all plot data for which respectively pairs of annual R-SL, RC-SL, R-SLu and RC-SLu were available, the best regression correlations (equation 5) are generally observed between annual R and SL (Table 8). Only for plots under grassland and post-fire better correlations between RC and SL were found. As a lack of data on annual P, plot length or plot slope gradient data does not allow the calculation of annual RC or SLu, respectively, a smaller number of plots was available to calculate annual RC-SL, R-SLu and RC-SLu regressions. Nevertheless, the results of the regression analyses did not change using a subdata set of 611 PL where all four variables (annual R, RC, SL and SLu) are available. In general, regression with unweighted annual R, RC, SL and SLu data yields slightly better correlation than weighted regression, but weighted regressions were nevertheless considered to be more appropriate as the number of plot-years is taken into account (see section II(3)).
Coefficients of determination (r2) and number of plots (#PL) for the unweighted and weighted annual runoff–annual soil loss (R versus SL), annual runoff coefficient–annual soil loss (RC versus SL), annual runoff–annual unit plot soil loss (R versus SLu) and annual runoff coefficient–annual unit plot soil loss (RC versus SLu) regressions. Weighting was done using the square root of the number of plot-years. For each land use, the regression with the highest r2 value is indicated in bold. n.s. = regression is not significant at α=0.05. NA = no regression was calculated as no SLu values were calculated for these land-use types.
Regressions between annual R and SL were found not to be significant for fallow, vineyards and construction sites and hence these land-use types were not included in Figures 9a and 9b. As a large number of plots is available for cropland and crop type is expected to have an important effect on annual SL rates in this land-use type (e.g. Auerswald et al., 2009; Gabriëls et al., 2003), a further subdivision according to crop type was made (Table 2b). The regression between annual R and annual SL was only significant for cereals (Figure 10), most likely due to a lack of sufficient data for the other crop types. Nevertheless, for maize, sugar beet and potatoes the available data show high annual RC and SL rates on nearly all of the plots where these crops were planted (Figure 10). Significant (α = 0.05) relationships between annual runoff (R) and annual soil loss (SL) for different land-use types in Europe and the Mediterranean. Regressions are weighted according to the square root of number of plot-years for each plot. n = number of plots. Comparison of the significant (α = 0.05) weighted regressions (SL = aRb) between annual runoff (R) and annual soil loss (SL) for different land-use types in Europe and the Mediterranean. Regressions are weighted according to the square root of the number of plot-years for each plot. For regression parameters, see Figure 9a. Relationship between annual runoff (R) and annual soil loss (SL) for cereal, maize, potato and sugar beet in Europe and the Mediterranean. Only the regression for plots with cereals was found to be significant (α = 0.05) and is given here.


5 Effects of land use on runoff coefficient and soil loss for different climatic zones
In order to investigate possible impacts of the CZ on the relation between annual RC and SL for different land-use types, plot data were grouped according to CZ. After initial data analysis, the division of all plot data into seven different CZ, some of which contain only a few PL, was found to be too detailed. Therefore, the Mediterranean and Anatolian CZ were grouped in a new CZ hereafter referred to as ‘pan-Mediterranean’ (M). Likewise, the Atlantic, Continental and Steppic CZ were grouped in a ‘temperate’ CZ (T) and the Boreal and Alpine CZ in a ‘cold’ CZ (C). Weigthed mean values, median values and coefficients of variation for R, RC, SL and SLu for each land use and CZ are given in Table 7. Shrubland, rangeland and post-fire occur mainly in the pan-Mediterranean zone and the few PL in the temperate zone are located relatively close to the pan-Mediterranean zone. Therefore, these land-use types were not further considered in this analysis. Also construction sites were disregarded due to insufficient plot data. For the other land-use types the differences in median value and the distribution of R and SL data by CZ are indicated in Figure 11.

Relation between median annual runoff coefficient (RC) and median annual soil loss (SL) for all plots in Europe and the Mediterranean by land-use type and climatic zone, with indication of the 25th and 75th percentiles for R and SL. M = pan-Mediterranean, T = temperate, C = cold. For the number of plots (#PL) and number of plot-years (#PY), see Table 5, for the division in climatic zones, see Figure 1.
Mean annual SL values are smaller in the pan-Mediterranean than in the other two CZ for all land-use types (Table 7). With respect to the median, median annual SL rates for bare, fallow, tree crops and vineyards are significantly lower in the pan-Mediterranean zone than in the other CZ, except the differences between the pan-Mediterranean and cold CZ for bare PL, and the difference between pan-Mediterranean and temperate CZ for tree crops. These non-significant differences in annual SL are likely due to insufficient data for bare plots in the cold CZ and tree crops in the temperate CZ. Median annual SL in the pan-Mediterranean zone is not smaller for cropland, grassland and forests. However, mean annual SL for all these land-use types is smaller in the pan-Mediterranean zone compared to the temperate zone, indicating that the highest annual SL on cropland occur more frequently in the temperate zone than in the pan-Mediterranean (Figure 11). Only differences between the pan-Mediterranean and temperate CZ for cropland and between pan-Mediterranean and cold CZ for grassland were found to be significant.
With respect to annual RC, median annual RC is significantly higher in the cold CZ than in the other CZ for cropland, fallow and grassland. Also, mean annual RC is the highest in the cold CZ for all land-use types. For grassland, there is no significant difference in median annual SL between the different CZ, but differences in median annual RC between different CZ are significantly different, indicating that while grassland has consistently low annual SL rates throughout the study area, annual RC rates differ depending on CZ.
6 Effects of annual precipitation on annual runoff and soil loss for different land uses
Figure 12 displays the relation between P and the weighted mean annual R, RC, SL and SLu for the different land-use types. As all PL on construction sites fall within the 250–500 mm annual P-class, they were not included in this figure. Separate graphs displaying the annual R, RC and SL by precipitation class for the different land-use types are presented in Figures 13a, 13b and 13c. It should be noted that for some land-use types only a few PL are available and hence some of the observed trends are uncertain. Nevertheless, clear trends can be observed. For all land-use types, annual R generally increases with increasing annual P (Figure 13a). For plots under cropland, there is already a substantial increase in annual R for 500–750 mm annual P, while for plots with a (semi-)natural vegetation cover, the most substantial increases in annual R occur in the 750–1000 mm and >1000 mm annual P classes. Application of the Kolmogorov-Smirnov test between subsequent precipitation classes indicated that the distribution of R data was only significantly different between the 250–500 mm annual P classes for cropland, fallow, post-fire forest and shrubland plots. For forest, also the difference between the 0–250 and 250–500 mm classes was significant. For annual RC, however, no clear general trend with P over all land-use types can be noted (Figure 13b).

Weighted mean annual runoff (R), annual runoff coefficient (RC), annual soil loss (SL) and annual unit plot soil loss (SLu) for each land-use type in Europe and the Mediterranean as a function of the annual precipitation (P) interval.

Median trend, distribution and weighted mean of annual runoff (R) by annual precipitation (P) interval and land-use type in Europe and the Mediterranean. Significant differences (α=0.05) in R between subsequent P classes are indicated by full lines while insignificant differences are indicated by dashed lines. PL= number of plots, PY= number of plot-years.

Median trend, distribution and weighted mean annual of runoff coefficient (RC) by annual precipitation (P) interval and land-use type in Europe and the Mediterranean. Significant differences (α=0.05) in RC between subsequent P classes are indicated by full lines while insignificant differences are indicated by dashed lines. PL= number of plots, PY= number of plot-years.
With respect to annual SL (Figure 13c), there is a general trend towards higher annual SL with higher annual P in plots with crop cultivation (vineyard, tree crops, bare, cropland and fallow). Mean annual SL under bare conditions increases between 250–500 mm and 500–750 mm, but levels afterwards and even decreases between the 750–1000 mm and >1000 mm annual P classes. Similarly, annual SL in cropland generally increases with increasing annual P, but mean annual SL for the >1000 mm annual P class is higher than for the 750–1000 mm class. Annual SL increases gradually with increasing annual P for post-fire, rangeland and forest annual PL. For grassland, the highest annual SL are observed in the 500–750 mm class. For shrubland, annual SL first increases to a maximum in the 250–500 mm class, and then declines over the 500–750 mm and 750–1000 mm classes. Mean annual SL increases again between the 750–1000 and >1000 mm classes, but the latter class corresponds to only 3 PL so results are not conclusive. Significant differences in annual SL between subsequent P-classes were found between the 0-250mm and 250-500mm classes for shrubland and between the 250-500mm and 500-750mm classes for bare, cropland and shrubland plots. For cropland, also the 500–750 mm and 750–1000 mm classes were found to have significantly different distributions. Corrections for plot length and slope gradient resulted generally in annual SLu rates that are lower than the original annual SL rates. Nevertheless, the trends observed in the graph depicting the relation between annual P and annual SL (Figure 13c) were found to persist. Correlating the Modified Fournier Index with annual R, RC and SL did not yield clearer trends than using annual P measured on the plots (Figure 14). Median trend, distribution and weighted mean of annual soil loss (SL) by annual precipitation (P) interval and land-use type in Europe and the Mediterranean. Significant differences (α=0.05) in R between subsequent P classes are indicated by full lines while insignificant differences are indicated by dashed lines. PL= number of plots, PY= number of plot-years. Relations between Modified Fournier Index (MFI) and annual precipitation (P), annual runoff (R) and annual soil loss (SL) for cropland plots in Europe and the Mediterranean

IV Discussion
1 How representative are the available plot runoff and soil loss data for Europe and the Mediterranean?
The assessment of annual R and annual SL rates based on the review of data measured on runoff and soil loss plots is inherently biased since runoff plot experiments are generally set up to answer specific research questions and not to be representative for the entire range of actual hillslope conditions (Auerswald et al., 2009; Cerdan et al., 2006; Vanmaercke et al., 2012). Many plot-measuring sites are located in areas that experience interrill and rill erosion. Furthermore, no plots are located in Iceland or in Scandinavia above 61 degrees latitude (Figure 1), probably due to logistic problems in cold environments and the low population density in these areas. A comparatively small number of runoff plots was found for northeastern Europe (Table 3), since many of these data have not been published in international journals and are not easily accessible. While even more plot data are likely existing, the database compiled in this study is currently the largest compilation of field-measured annual R and SL data at plot scale and for Europe and the Mediterranean region.
The relatively short measuring periods for the plots (mean: 6 yrs, median: 4 yrs) can be attributed to the relatively short duration of research or PhD projects during which most of the plots are established. This also means that it is difficult to assess temporal variability in Ra and SLa on the longer term for a specific plot, which is nevertheless an important aspect for conservation planning (Bagarello et al., 2011). Figure 2 shows a clear decline of the number of plots in operation after 1996. This can be attributed to the fact that plot measurements are time- and labour-consuming and field-measured plot experiments are abandoned in favour of modelling studies, which have become more prevalent with increasing and cheaper computing power over the last 20–30 years (Merritt et al., 2003). A similar decline in catchment sediment yield studies was observed by Vanmaercke et al. (2011b), but was already starting between 1970 and 1980 when many large state-led monitoring schemes were abandoned. Most plot-scale studies are conducted by individual researchers however, and the number of publications on runoff and soil loss plots has not decreased in recent years (Figure 2). Moreover, recent publications review previously published data to assess erosion rates and variability under field conditions at national or regional scales (Table 1), to use these data in new spatial analyses (e.g. Cerdan et al., 2010) or to validate erosion models (e.g. Amore et al., 2004; Licciardello et al., 2009; Quinton, 1994; Tsara et al., 2005). Hence, there is an ongoing interest in field-measured annual R and SL data and a review and regular update of existing data sets offers opportunities to use these data in new analyses and thus to give added value to previously published studies (Baade and Rekolainen, 2006).
A large proportion of research has been carried out on bare and cropland plots (Table 5). Runoff and soil loss after wildfires have received considerable attention, despite annual R and SL rates still being low in comparison to those witnessed on some other land uses (Shakesby, 2011; Table 7). Permanent crops (tree crops and vineyards) have received more attention in plot studies in comparison to the areal percentage they represent in the CORINE land-cover map (Vanmaercke et al., 2012). Despite the relatively small area occupied by these land-use types, they are important contributors to total soil loss in Europe (Cerdan et al., 2010) and can be the dominant land-use type in certain regions. Furthermore, some of the highest annual RC (up to 40.2%) and SL rates (up to 151 Mg.ha-1.yr-1) recorded in the database occurred on permanent crops (i.e. vineyards and tree crops) (Figure 7, Table 7). The quantitative measurements currently available for these land-use types may still not be sufficient to make a comprehensive assessment of erosion risk in these land-use types (e.g. Gómez et al., 2008).
Also the distribution of plot lengths and slope gradients of the plots shows a research bias (Figure 4). Plots with a plot length between 20 and 30 m, i.e. close to the standard USLE plot (22.1 m; Renard et al., 1997), are by far the most frequent. Plot length is mainly determined by logistic limitations of the studies, and the relation to actual field slope lengths, for which very few data are available, is unknown. On the other hand, observed frequency distributions of the plot slope gradients for different land-use types (Figure 5) show both a difference between slope gradient distributions for different land-use types as they occur in the field as well as a research bias towards steeper slopes. This was demonstrated by Cerdan et al. (2006) who found that mean slope gradient of plots on grassland, forest and shrubland corresponded well to mean slope gradient for these land-use types on the reclassified CORINE land cover map, while mean slope gradients for arable, vineyards, orchard and post-fire plots were found to be steeper than the CORINE average. On the whole, the plots show a relatively steep slope gradient, with the majority of plots having a slope gradient above the 9% of the standard USLE plot (Figure 4). As was shown by Boardman (1998), careless extrapolation or generalization of these data can therefore lead to overestimations of annual SL rates. Therefore, the controlling environmental variables should always be accounted for when evaluating or extrapolating rates of annual R, RC and SL (Cerdan et al., 2010; de Vente et al., 2011). Hence, the rates of annual R and SL presented in this study may not be representative for the flat regions in Europe and the Mediterranean (Cerdan et al., 2010). Nevertheless, the extent of the database compiled in this study allows for the best representation of the relations between measured runoff, soil loss and annual precipitation for different land-use types currently available.
While the use of runoff and soil loss plots allows for a relatively easy assessment of annual runoff, runoff has received considerably less attention than soil loss in the literature as can be noted in the literature review (Table 1) and the number of plots and plot-years (Table 5). One of the problems underlying this discrepancy is the complex relationship between R and environmental variables (Wischmeier, 1966) and as a result also the relation between annual R and SL is less studied. Nevertheless, the assessment of annual R is an important part of many erosion models (Merritt et al., 2003) and hence a better understanding of runoff generation and runoff-soil loss relationships would contribute towards better erosion models (e.g. Kinnell and Risse, 1998). Runoff generation is also an important problem in itself, both on-site (e.g. loss of plant available water; Wallace, 2000) and off-site (e.g. flash floods). For instance, water is a key resource in the Mediterranean (Araus, 2004; Vanmaercke et al., 2011b) and annual RC are often higher than annual RC for comparable land uses in temperate regions, while annual SL rates tend to be lower in the pan-Mediterranean CZ than in the temperate CZ (Table 7, Figure 11), hence excessive runoff may be of more concern than soil loss.
2 Effects of plot length and slope gradient on annual runoff, runoff coefficient and soil loss
Annual SL on bare, cropland and tree crop plots are positively correlated to plot L-factor (Table 6). This finding concurs with results found by Cerdan et al. (2010). For these land-use types, also annual R and RC are positively correlated with plot length (Table 6). This can be attributed to the high connectivity along flow paths under these land-use types, allowing runoff to accumulate along the flow path over longer plot lengths and resulting in higher annual SL due to the increased detachment and transport capacity of the overland flow (Govers and Poesen, 1988; Prosser and Rustomji, 2000). A significant effect of plot S-factor on annual SL was only found for cropland and shrubland, although the correlations between the plot S-factor and annual SL are also relatively strong for bare plots and fallow plots (Table 6, p = 0.06).
With respect to climatic zone, Cerdan et al. (2010) found that for arable and bare plots the relation between LS-factor (equation 3) and annual SL was only significant for plot data collected outside the Mediterranean CZ and not for plots in the Mediterranean CZ, which was attributed to the higher surface rock fragment cover in the Mediterranean, especially on steeper slopes (Poesen et al., 1998). In this study, similar results are found for bare plots in the temperate CZ where the relation between the L-factor (rs = 0.62, p < 0.01), S-factor (rs = 0.25, p = 0.03) and LS-factor (rs = 0.44, p < 0.01) was significant, while for plots in the pan-Mediterranean CZ none of these relations were significant. However, for cropland in both the temperate and pan-Mediterranean CZ significant relations between L-factor and annual SL (rs = 0.36, p < 0.01 and rs = 0.20, p = 0.01, respectively) and between LS-factor and annual SL (rs = 0.33, p < 0.01 and rs = 0.17, p = 0.03, respectively) were found. This could be due to the inclusion of more pan-Mediterranean cropland plots in the database as compared to the database used by Cerdan et al. (2010). By including more plots, the probability of including low-frequency, high-intensity events for which there is a significant topographic effect increases. Nevertheless the correlation coefficient between LS-factor and annual SL remains smaller for cropland plots in the pan-Mediterranean (rs = 0.17, p = 0.03) than for cropland plots in the temperate zone (rs = 0.33, p < 0.01).
For grassland, no significant correlations were found between the L-, S- or LS-factor and annual SL. This can be explained by the high root density of the grass, which reduces soil erodibility (De Baets et al., 2006). The correlation between plot length and R was significantly positive (Table 6) however, which may be explained by the fact that runoff is more likely to converge on longer slopes, and the effect of grass cover on annual R retention is often less pronounced and more variable than the effect on annual SL, as is shown by studies on the effectiveness of grass buffer strips (e.g. Blanco-Canqui et al., 2004) and grassed waterways (e.g. Fiener and Auerswald, 2003). Contra-intuitively, a significant negative relation was noted between slope gradient and both R and RC for grassland. This is probably due to the fact that grasslands in the northern regions often lie on gentle slopes but are very wet due to a clayey subsoil. These soils often need drainage (e.g. Øygarden, 1996; Øygarden et al., 1997; Turtola and Paajanen, 1995; Turtola et al., 2007; Warsta et al., 2009). For plots under forest and post-fire a significant negative correlation between plot length and annual R and annual RC was found (Table 6). This is probably related to the heterogeneity of soil cover and macropore distribution in these land-use types, meaning that as plot length increases, runoff is more likely to flow through patches with increased roughness or infiltration capacity. For shrubland, a significant negative correlation between slope gradient and annual R and RC is observed, while the correlation between the S-factor and annual SL is positive (Table 6). An explanation for this can be offered by the reduced tendency in surface sealing of bare patches with steeper slopes (Poesen, 1984, 1986), reducing runoff but increasing splash erosion (Bradford et al., 1987). Nevertheless, plots under shrubland are characterized by a high spatial variability at plot scale and hence the establishment of relations between plot length and slope gradient and annual R, RC and SL is difficult (Cammeraat, 2002).
No significant effect of plot length on annual R or RC was found for vineyards, but the L-factor was significantly correlated with annual SL (Table 6). This is mainly due to the fact that high annual R and RC already occur on short plots, resulting in high runoff rates independent of plot length. For tree crops, which are also considered as a land-use type with permanent cultivation, a positive relation between plot length and annual R and RC was found. This difference with vineyards may be due to differences in soil types or the generally higher vegetation cover associated with tree crops, though no data were available to check these hypotheses. As was already indicated by Cerdan et al. (2010), our results illustrate that the relations between plot length and/or slope gradient and annual SL depend on the considered land use. In addition, our analyses show that also relations of plot length and slope gradient with R and RC depend on land-use type. Plot length and slope gradient may affect annual R and RC differently than they affect annual SL. Furthermore, a further subdivision of the major land-use types analysed by Cerdan et al. (2010) (bare, arable, permanent cultivation and permanent vegetation) reveals that also within these groups, relations between topographic factors and annual SL may vary. For instance, in the permanent vegetation group, no significant relation between LS-factor and annual SL was found for forest, grassland and rangeland plots, but there was a significant relationship between the LS-factor and annual SL for shrubland (Table 6).
3 Frequency distributions and relationships between annual runoff and soil loss for various land uses
Previous studies (e.g. Cerdan et al., 2010; Kosmas et al., 1997) showed the importance of vegetation cover as an important determinant of annual SL on a (sub)continental scale. This study further confirms this, as significant differences between the distributions of SL data were mostly found between land-use types with cultivation and land-use types with semi-natural vegetation. Construction sites have the highest observed annual SL rates (Figure 7, Table 7), although the limited amount of available data (Table 5) does not allow clear conclusions on the frequency distribution of the observed values. Bare plots and plots with crop cultivation (cropland, vineyards and tree crops) form a second cluster with high R, RC and SL rates (Figures 7 and 8). All these land-use types are characterized by severe anthropogenic disturbance (construction sites and bare PL) or intensive tillage of the soil for agriculture (cropland, tree crops and vineyards). On these plots, annual SL regularly exceeds tolerable soil loss rates (T-values) of 5–12 Mg.ha-1.yr-1 (Montgomery, 2007). The mean annual SL rate of 11.6 Mg.ha-1.yr-1 for tree crops in the pan-Mediterranean zone (Table 7) was found to be much higher than reported in previous studies such as Cerdan et al. (2010) who found an average annual SL rate for orchards in the Mediterranean of 1.67 Mg.ha-1.yr-1 and Fleskens and Stroosnijder who concluded that annual SL in olive orchards is unlikely to exceed 10 Mg.ha-1.yr-1. Annual SL in olive orchards is often the result of infrequent high-intensity rain events and depends strongly on spatial scale and land management (Gómez et al. 2008). Depending on which data are included in the analysis and the way the available data are analysed they may be deemed low (Fleskens and Stroosnijder, 2007) or frequently in excess of tolerable rates (Gómez et al., 2008). The use of a larger database (20 PL, corresponding to 80 PY) in this study indicates that annual SL in tree crops in the pan-Mediterranean zone can frequently exceed tolerable levels.
Vineyards show a high coefficient of variation for SLa in addition to a high average SLa value (Table 7). This is attributed to the high erosion susceptibility of this land-use type after vineyard establishment. Wicherek (1991) reports a sharp decline in soil losses over the first three years after vineyard establishment in the Aisne region, northern France (57.3, 28.7 and 1.4 Mg.ha-1.yr-1, respectively). Tropeano (1984) also reports the highest soil loss (8.3 Mg.ha-1.yr-1) to occur in the first year after vineyard establishment in the Piemonte region, northwestern Italy, which was reduced to 1.2 Mg.ha-1.yr-1 in the following year. This is also illustrated by Engels (personal communication 2009) who notes that soil loss from an old undisturbed vineyard in the Moselle region, western Germany, is considerably lower (0.5 Mg.ha-1.yr-1) than that from an adjacent vineyard where vines were removed and roots destroyed (4.4 Mg.ha-1.yr-1). Nevertheless, no clear trends in annual SL between 8 and 17 years after vineyard establishment were found for vineyards in the Douro region, northern Portugal (de Figueiredo, personal communication 2010; de Figueiredo and Gonçalves Ferreira, 1993; de Figueiredo and Poesen, 1998). Brenot et al. (2006) measured soil loss by vine stock unearthing (i.e. over the complete period since vineyard establishment) for vineyards of different ages (10–54 years) in Burgundy, France, but no trend in annual SL with respect to vineyard age is observed. In contrast to the high SLa observed for several vineyards, some vineyards established on very steep slopes (c. 45%) with very stony soils show comparatively very low mean SLa values in the Douro region, Portugal (0.2–0.5 Mg.ha-1.yr-1: de Figueiredo and Gonçalves Ferreira, 1993; de Figueiredo and Poesen, 1998) or the Moselle region (0.02–0.5 Mg.ha-1.yr-1: Richter, 1980). Although the number of available plots on which SL data were collected directly after vineyard (re)planting (n = 1) or removal of the vines (n = 1) is limited, these results indicate that vineyards are in general very vulnerable to soil loss after establishment or periodical replanting of the vines due to the often intense soil disturbances (Borselli et al., 2006). After a number of years, annual SL in vineyards often decreases to relatively low rates as is shown by the median annual SL of 0.9 Mg.ha-1.yr-1 for the whole of Europe and the Mediterranean (n = 39, Table 7). This is attributed to the stony soils often associated with vineyards which develop an erosion pavement with a high rock fragment cover which drastically reduces SL (Poesen et al., 1994).
Land-use types with a (semi-)natural vegetation generally have mean annual SL values well below the T-value, illustrating the dominant control of vegetation cover over annual SL (Montgomery, 2007). However, SL rates up to 7 Mg.ha-1.yr-1 can be observed on all land-use types and differences between land-use types are mainly situated in the frequency of occurrence for high annual SL values (i.e. above 10 Mg.ha-1.yr-1) (Figure 7). Scaling of annual SL values to unit plot soil loss rates results generally in somewhat lower annual SLu rates (Figure 6), but does not change the differences between different land-use types observed above (Figure 7, Table 7). Hence, while plot length and slope gradient may have an effect on annual SL values, they cannot fully account for the differences between different land-use types.
The distribution of annual R and RC values follows the same pattern as annual SL, with land-use types with a (semi-) natural vegetation cover having generally lower annual R and RC rates (Figure 7). However, fewer significant differences are observed between different land-use types combinations than for annual SL. There is a larger overlap between the annual R and RC data distributions (Figure 7) which indicates that the effect land-use types have on annual SL is more pronounced than they are on annual R and RC. No standard exists for tolerable annual R and RC levels, but excessive runoff also has several negative off-site effects such as gully formation (Poesen et al., 2006) and flooding (García-Ruiz et al., 2010).
With respect to the relations between annual R/RC and SL/SLu, the small differences between r2 values obtained for weighted and unweighted regressions (Table 8) can be explained by the distribution of the number of plot-years since a large part of the number of plots corresponds to a low number of plot-years and hence has similar weights (Figure 3). The generally better correlations between annual R and SL than between annual RC and SL (Table 8) can be explained by the fact that mainly the total runoff volume, rather than the runoff coefficient, determines the erosive power and transport capacity of the overland flow. Nevertheless, the inclusion of differences in annual precipitation makes it easier to compare RC rates between different plots. Hence, both annual R and RC are important variables to consider. The general trend in the relation between annual R and SL is that (semi-)natural land uses show annual SL rates of up to one order of magnitude less than SL values for land-use types with crop cultivation for the same annual R rate (Figure 9b). This is not the case for tree crops, however, where the increase in annual SL with increasing R is much stronger than for other land-use types (Figure 9a). However, this trend should be interpreted with caution as it is based on only 13 PL. For shrubland, clusters of data that correspond to individual plot-measuring sites show a good correlation between annual R and SL (Figure 9a), while the global regression is affected by more scatter in the data, which is most likely determined by local environmental factors that differ from plot site to site. Summer crops like maize, potatoes and sugar beet were found to result in R and SL rates which are among the highest for cropland (Figure 10).
4 Effects of land use on runoff coefficient and soil loss for different climatic zones
The higher annual RC observed in the cold zone may be attributed to the combination of snowmelt and a frozen soil in spring (e.g. Alström and Bergman, 1990; Wade and Kirkbride, 1998) and the generally lower annual evapotranspiration rate at high latitudes (Weiß and Menzel, 2008). Furthermore, annual R and RC rates in the pan-Mediterranean zone are generally higher than in the temperate CZ (Table 7, Figure 11). This is likely due to the combination of soil properties, the often smaller and more discontinuous natural vegetation cover in the pan-Mediterranean region, and the seasonality of the rainfall with a large fraction of the annual rainfall concentrated in a few important events during a short winter season (Altava-Ortiz et al., 2011; Mehta and Yang, 2008).
Nevertheless, these higher annual R and RC rates in the pan-Mediterranean do not result in high annual SL (Figure 11, Table 7). Smaller annual SL rates in the Mediterranean were also observed by Cerdan et al. (2010) and can be explained by the generally higher rock fragment cover in Mediterranean soils which is known to reduce soil loss rates (Poesen and Lavee, 1994; Poesen et al., 1994). Especially for vineyards, which are often located on stony soils, annual SL rates in the pan-Mediterranean CZ are much smaller than those in the temperate CZ (Figure 11). Nevertheless, Sanchis et al. (2008) observed that also for soils with a low rock fragment content (<10% content by mass) soil erodibility was lower in the Mediterranean. This is attributed to a dominance of clay rich soils in arable land in the Mediterranean and the associated low soil erodibility of these soils (Torri et al., 1997). Furthermore, very erodible soil types like loess-derived soils occur almost exclusively in the temperate CZ. Hence, it should be noted that differences between CZ include more than just a climatic effect as CZ is also a proxy for particular soil properties which affect annual R and SL like soil susceptibility to cracking. Apart from differences between precipitation and soil characteristics between the temperate and pan-Mediterranean CZ, differences in conventional land management like tillage frequency and depth in cropland may also account for part of the observed variability. Furthermore, for grassland and forest, there is no significant difference between median SL rates for plots in the pan-Mediterranean zone and plots in the temperate CZ, which again indicates the important effect of vegetation in controlling annual SL.
5 Effects of annual precipitation on annual runoff and soil loss for different land-use types
For all land uses, there is a consistent trend towards higher annual R with increasing annual P which is more pronounced for land-use types with crop cultivation than land-use types with (semi-) natural vegetation (Figure 13a). For all land-use types with crop cultivation, except for vineyards, there is a trend towards increasing annual RC with increasing precipitation, indicating that as rainfall increases there is a larger percentage of excess rainfall (Figure 13b). This may be related to distribution of rainfall patterns throughout the year, with areas with high annual P generally having a more uniform precipitation distribution throughout the year. This causes seasonal saturation of the soil and faster runoff formation (Ponce and Hawkins, 1996).
For annual SL, there is also an increase in annual SL rates with increasing annual P for most land-use types (Figure 13c) but, unlike annual R, the increase is more gradual. For bare and cropland plots, there is even a decrease in annual SL in the highest annual P classes (>750 mm for bare and >1000 mm for cropland). This is likely due to variations in seasonality of the rainfall. The relatively uniform rainfall distribution in regions with high annual P makes these regions less prone to unfrequent extreme rainfall events in periods of the year when the soil is vulnerable to erosion (Edwards and Owens, 1991; Larson et al., 1997).
Shrubland is generally limited to the pan-Mediterranean region and a maximum in SL is observed in the 250–500 mm annual P class (Figure 13c), and both the differences between the lower and higher precipitation class are significant. A similar trend was noted by Kosmas et al. (1997) who attributed this to an initial increase in annual SL with increasing annual P as the erosion potential of the rain also increases, combined with insufficient vegetation cover in drylands. However, as annual P futher increases the vegetation cover also increases, effectively reducing annual SL for higher annual P. This is similar to the mechanism and trend for sediment yield proposed by Langbein and Schumm (1958). However, this trend is not observed for other land-use types, as vegetation cover in shrubland can be expected to be the most sensitive to changes in annual P. Contrary to the study by Kosmas et al. (1997), this maximum is not noticeable for annual R, which may indicate that, at a Mediterranean-wide scale, runoff volume is more determined by other environmental characteristics than vegetation cover, such as surface sealing and rock fragment cover.
In general, significant differences in the frequency distribution of both annual R and annual SL were mostly found between the 250–500 mm and 500–750 mm P classes (Figures 13a and 13c). This may indicate that significant changes in the rainfall-runoff and rainfall-soil loss relations occur at around these annual P values which are likely related to changes in rainfall regime and distribution throughout the year. Hence, comparison of annual R and SL rates with measures that take rainfall distribution throughout the year into account could improve the analysis. Nevertheless, the use of MFI as an indicator of climatic effect on annual R, RC and SL did not yield better correlations than the use of annual P (Figure 14). The MFI values used in this analysis represent long-term average (i.e. climatological) values which may not be representative for the specific years of plot measurements. This explains why low annual R and SL rates occur regularly in zones with high MFI values. Furthermore, the relation between MFI and annual P (Figure 14) and MFI and EI30 values (Torri et al., 2006) is not straightforward.
VI Conclusions
While a decrease in the number of studies using runoff and soil loss plots under natural rainfall has occurred since the mid-1990s in Europe and the Mediterranean, there is still an interest in plot data. Recently, several studies have assessed the erosion problem by reviewing existing soil loss data from plots and by using these data in empirical studies or erosion model validations. Many of these studies investigated soil loss rates in relation to land use, topography and soil properties. However, plot-scale runoff has been largely neglected in these studies and the relation between runoff and soil loss has not been reviewed on a scale covering Europe and the Mediterranean. Nevertheless, runoff estimation and its relation to soil loss are important parts of many erosion models. On a (sub)continental scale also climatological differences are often not taken into account in overview studies. To address these issues, the largest data set of runoff and soil loss plot data for Europe and the Mediterranean region was compiled in this study, which includes for the first time both annual runoff and annual runoff coefficients at the continental scale. In general, soil loss studies using runoff plots are mainly focused on bare plots and cropland and fewer data are available for construction sites, tree crops and vineyards, in spite of the high soil losses that may be associated with these land-use types. Variation in annual runoff and soil loss rates observed on the plots ranges over several orders of magnitude. While land-use types with crop cultivation (cropland, fallow, tree crops and vineyards) have higher mean soil loss rates than land-use types under (semi-)natural vegetation (grassland, rangeland, shrubland, forest and post-fire), there are still large variations within each of these land-use types, which can only be partly accounted for by topographic (i.e. plot length and slope gradient) differences. Annual runoff rates follow the same pattern as annual soil loss rates, but differences between land uses are less clear. The generally good relations between annual runoff and annual soil loss illustrate the key importance of the relation between runoff and soil loss for a good assessment of soil loss rates. Further quantitative analysis of this relation may also contribute significantly to improve predictions in erosion models. Furthermore, the relation between annual runoff and soil loss also depends on climate, with comparatively high runoff coefficients in cold climates and lower soil losses in the pan-Mediterranean region. This indicates that runoff-soil loss relations may show important regional variations. Apart from the importance of runoff as a causal factor of soil loss, runoff in itself is associated with several problems such as flooding and plant-available water. Techniques specifically aimed at reducing runoff and runoff coefficients can contribute to a more efficient use of rainwater on-site to increase food production, especially in drier regions like the Mediterranean. As expected, annual rainfall was found to be related to annual runoff and soil loss, and the vegetation feedback effect for shrubland proposed by Kosmas et al. (1997) was also observed on a Mediterranean-wide scale for soil loss, but not for runoff. Nevertheless, a large part of the variation in runoff and soil loss rates remains unexplained and better relations between annual precipitation and runoff and soil loss can likely be obtained by accounting for rainfall erosivity, but use of a Modified Fournier Index did not yield better results than the use of the annual precipitation measured on the plots. Hence there are still possibilities to expand the research to better account for rainfall erosivity. Further research may also focus on the effects of several important soil characteristics (e.g. texture, organic matter content) and a more detailed analysis of several relations for which general trends were established in this study (e.g. the relation between annual runoff and plot length and plot slope gradient, the relation between annual runoff and soil loss). In conclusion, this meta-analysis of field plot data for Europe and the Mediterranean allows a quick assessment of the impact of land-use change scenarios on annual runoff, runoff coefficient and soil loss on a continental scale.
Footnotes
Appendix: List of abbreviations and symbols
Land-use types:
Climatic zones:
Funding
The research reported in this paper was conducted in the framework of the EC-DG RTD- 6th Framework Research Programme (sub-priority 1.1.6.3) - Research on Desertification- project DESIRE (037046): Desertification Mitigation and Remediation of land - a global approach for local solutions. M. Vanmaercke received grant-aided support from the Research Foundation Flanders (FWO), Belgium.
