Soil Erosion Processes in European Vineyards: a Qualitative Comparison of Rainfall Simulation Measurements in Germany, Spain and France

Small portable rainfall simulators are considered a useful tool to analyze soil erosion processes in cultivated lands. European research groups in Spain (Valencia, Málaga, Lleida, Madrid and La Rioja), France (Reims) and Germany (Trier) have used different rainfall simulators (varying in drop size distribution and fall velocities, kinetic energy, plot forms and sizes, and field of application) to study soil loss, surface flow, runoff and infiltration coefficients in different experimental plots and Mosel-Ruwer valley in Germany). The measurements and experiments developed by these research teams give an overview of the variety of methodologies used in rainfall simulations to study the problem of soil erosion and describe the erosion features in different climatic environments, management practices and soil types. The aims of this study are: (i) to investigate where, how and why researchers from different wine-growing regions applied rainfall simulations with successful results as a tool to measure soil erosion processes; (ii) to make a qualitative comparison about the general soil erosion processes in European terroirs; (iii) to demonstrate the importance of the development of standard method for measurement of soil erosion processes in vineyards, using rainfall simulators; and (iv) and to analyze the key factors that should be taken into account to carry out rainfall simulations. The rainfall simulations in all cases allowed infiltration capacity, susceptibility of the soil to detachment and generation of sediment loads to runoff to be determined. Despite using small plots, the experiments were useful to analyze the influence of soil cover to reduce soil erosion, to make comparisons between different locations, and to evaluate the influence of different soil characteristics. The comparative analysis of the studies performed in different study areas points out the need to define an operational methodology to carry out rainfall simulations, which allows us to obtain representative and comparable results and to avoid errors in the interpretation in order to achieve comparable information about runoff and soil loss.


Introduction
The concept of terroir defines a vineyard with particular regional vitivinicultural practices (land management, soil tillage, crop management, etc.) and identifiable bio-physical environmental conditions (soil, climate, landscape and topography) with direct influences on grape composition [1][2][3][4].
Several authors evidenced the importance of vineyard degradation during the last few decades, induced by applying chemical weeding, seasonal intensive tillage, green pruning and the use of heavy machinery [5][6][7][8][9].All these activities notably affect one of the most important components of the terroir: the soil.Pedological process research in terroirs has become one of the most relevant topics nowadays in geosciences [10].Vineyards are notably one of the land uses, in particular in Mediterranean environments, in which high erosion rates are being recorded [11].There is a growing awareness of the need to avoid high soil erosion rates and to reduce the transfer of pollutants downstream [12,13].To meet this goal, new agricultural measures are being applied through the development of new agri-environments with conceptual models of public spending efficiency [14] in the recovery of soil functions and services [15][16][17][18].The occurrence of extreme rainfall events, soil characteristics and in some cases steep slopes contribute to soil degradation.Several studies note that these factors cause soil degradation in vineyards and show how they affect soil surface characteristics [19,20] and can be relevant indicators to classify and analyze actual soil degradation processes in vineyards [13,[21][22][23][24][25][26].According to this, the most important question to investigate has become the spatial and temporal variability of soil patterns caused by geomorphological dynamics over only a few hectares [27][28][29][30].
To enable quantification of these spatially variable parameters, small portable rainfall simulators are useful tools, because they provide information about the detachment, transport and deposition of soil particles.They are the key to understanding soil erosion and hydrological process, and to design strategies to control the non-sustainable soil losses found in vineyards.Several authors have applied these experiments to measure rainfall-runoff processes, splash effect, soil erosion, infiltration, permeability, irrigation or nutrient movement [22,[31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48].However, despite the effort and time invested by researchers in the understanding of soil degradation using rainfall simulations, one of the most important problems is the lack of a standard rainfall simulator for an effective transfer of knowledge and generation of comparable data.
The characteristics of the different rainfall simulators, such as the spatial rainfall distribution, drop sizes and velocities, drop kinetic energy and plot forms and sizes are different among the rainfall simulators applied in the compared studies.In addition, different land covers and different experimental simulations characteristics (rainfall intensities, experimental time or plot sizes) make a successful comparison difficult [49,50].Moreover, climate is different at the different locations, which means that the rainfall simulators should produce rain of different characteristics.Therefore, it is important to have a fleet of rainfall simulators with diverse rainfall properties as they can all provide required types of experiments that are needed to match the environmental conditions of different regions.
However, a number of publications on rainfall simulations in vineyards asks for one standard methodology to be able to compare the results in different climatic settings and under different management.The data obtained from all the studies of rainfall simulations in vineyards could be of great significance to compare the simulated processes, and could also serve as a source of inputand validation-data for soil erosion modeling and extrapolations to larger scales [51,52].This could be an important and useful step for a reliable assessment of vulnerability, risk levels and control of soil degradation in European terroirs.These insights can in turn serve as a basis to develop more sustainable management of vineyards around Europe.
Therefore, the aims of this research are: (i) to investigate where, how and why researchers from different wine-growing regions applied rainfall simulations with successful results as a tool to measure soil erosion processes; (ii) to make a qualitative comparison about the general soil erosion processes in European terroirs; (iii) to demonstrate the importance of the development of a standard method for studying soil erosion processes in vineyards, using rainfall simulators and; and (iv) to analyze the key factors that should be taken into account to carry out rainfall simulations.

Study Area
The selected study areas are located in traditional and conventional European wine-growing regions (Figure 1 and Table 1).The area of study in Germany is located close to the village of Waldrach (49.7418N; 6.7524E), which is in the Ruwer-Mosel valley (Rhineland-Palatinate, Germany), a tributary of the Mosel River.These vineyards (with the Riesling grape variety) are characterized by conventional management, including soil tillage with machinery and grass/pruning cover.The site is located on steep slopes (25-50%) and the soil is characterized by a high rock fragments percentage on the surface, a silt proportion of 64.7% and a soil organic matter content of 6.1%.Total average annual rainfall is about 780 mm and the annual average temperature is 9.3 ˝C.For this area, Trier and Málaga universities studied the impacts of rainfall and human influences on soil erosion and determined with botanic marks of the graft unions soil loss rates of 3.4 Mg¨ha ´1¨yr ´1, high geometrical variations of the rills and elevated infiltration [53,54].
In France, the study site is located close to Paris, in the Montagne de Reims (49.16N; 4.12E).Total annual rainfall is similar to the German study area with 757 mm and maximum monthly events of about 72.9 mm.The most used grape varieties are Pinot noir, Pinot meunier and Chardonnay.In this case, several different soil treatments for the vineyards were applied: grass/pruning cover, plowing, and use of herbicides or without herbicide.Soils are characterized by similar percentage of silts and sands and relativity high soil organic matter (3.9-8.3%).Morvan et al. [20] analyzed the impact of different land management on soil in this area using rainfall simulations.They noted that ecological treatments are more effective against soil erosion than conventional management.
For Spain, five study sites were selected.The first one is the Montes de Málaga in Málaga (36.76N; ´4.39E).Martínez-Murillo and Ruiz Sinoga [55] analyzed the influence of the soil surface components on these hillslope vineyards with Muscat of Alexandria grape variety.This conventional plantation (with handmade and animal works) registers an annual average rainfall of 586.1 mm on soils with a high stoniness (about 52.8%) and low porosity (33.3%).
The second selected study area in Spain is Celler del Roure (38.7833N; 0.87E), in the Les Alcusses valley (Moixent municipality), which is located in Valencia region at 550 m.a.s.l. with a production of Monastrell grape variety.Sandy soils (>60% sand) with low soil organic matter content (1.01%) induced soil erosion dynamics in the form of rills and sheet wash.Annual average rainfall is about 420 mm with minimal monthly averages of 5 mm in summer.The rainfall simulations data were produced by the Valencia and Wageningen Universities, following their published methods [26, [31][32][33][34][35]56].
Els Hostalest de Pierola (41.59N, 1.77E) in Barcelona province, whose vineyards belong to the Penedès wine-growing region, is the third selected study area.This area (with mainly white varieties such as Parellada, Macabeo, Xarello and Chardonnay grape variety production) is characterized by soils with high contents of silt (25.8%) and sands (63.6%) and low soil organic matter (1.2%).Total average annual rainfall is 520 mm and the annual average temperature is 15 ˝C.Lleida University carried out several studies in that wine region with different topics, mostly related with soil, climatic and phenologic dynamics [57][58][59][60][61][62][63].
The fourth selected study area is Villamediana de Iregua (La Rioja) (42.26N; 2.23E).Soils have high clay (19.9%) and sand (39.9%) contents on terraces about 1-2 m height.This traditional vitivicultural area, with Tempranillo grape variety production is managed with tillage and herbicides using machinery.Annual average rainfall is 419 mm and annual average temperature is 14.6 ˝C, with minimum temperatures of 1.2 ˝C in winter.The Instituto Pirenaico de Ecología (CSIC) and La Rioja University performed several studies on soil erosion, land degradation and wheel tracks signals on these soils [39,64].
Finally, Campo Real in Madrid (40.35N; ´3.37E) with organic vineyards of Tempranillo grape variety production is the last experimental area.This is a semi-arid area having less than 400 mm of annual rainfall [65][66][67][68].As with the previous study area, soils have high sand content (58%), high clay content (24%) and low soil organic matter (1.2%).These soil characteristics, together with soil tillage with machinery, favor soil degradation processes.The regional extension services (IMIDRA) and researchers from the Autonomous Madrid University studied the effects of different vegetation cover (Brachypodium distachyon and annual barley) between the vine rows as protection measures against the land degradation processes in these vineyards.

Qualitative Comparison with Different Data of Rainfall Simulations in Vineyards
All the vineyard data from the different research groups were organized in a database in order to homogenize the information.An analysis was done to assess the differences among the experiments due to rainfall simulator design and characteristics, plot size and time intervals for recording runoff during the experiments (Table 2).The results were interpreted taking all parameters into account, which can play an important role in the final observed soil dynamics processes.
The first technical characteristic is the dropper system.In Moixent (Valencia, Spain), Ruwer-Mosel valley (Trier, Germany) and Villamediana de Iregua (La Rioja, Spain), the same nozzle (Lechler 460.608) was used.In the other locations, spraying systems (Campo Real, Madrid), Hardy nozzle (Montes de Málaga, Spain), a dropping system controlled by a water column (Els Hostalest de Pierola, Spain) and a Teejet nozzle (Montagne de Reims, France) were the principal applied brand.
The second characteristic was the size and form of the plot.In four rainfall simulations, ring plots with sizes between 0.25 and 0.45 m 2 were used (Moixent, Ruwer-Mosel valley, Montes de Málaga and Villamediana de Iregua).In the other locations, the plots were rectangular with different dimensions: 2 m 2 in Campo Real, 0.6 m 2 in Els Hostalest de Pierola, 0.25 m 2 in Montagne de Reims.With respect to the third characteristic, the droppers or nozzle height, there was more uniformity, as they ranged from 2 and 2.5 m above the soil surface.
The last characteristic was the operational method, where several differences were observed.For example, the total duration of the experiments ranged between 15 minutes in the vineyards of Campo Real and 90 minutes in Montagne de Reims, and the time intervals for collecting runoff samples ranged between 1 and 10 minutes.
For more concrete information about the different methods (rainfall intensities, study areas, duration of the experiments, etc.), the articles for each research group were also added to Table 2.
In order to homogenize the information, all data of the rainfall simulations experiments (data sheets) were organized into five-minute time steps.However, because some authors used different intervals in different experiments (between 5 and 15 minutes), results must be homogenized, calculating the suspended sediment load (SSL) per hour (g¨m 2 ¨h´1 ), surface flow (L¨m 2 ¨h´1 ) and suspended sediment load concentration (SSC in g¨L¨m 2 ¨h´1 ), using: where A is the area of the rainfall simulator plot and RFt is the duration of the rainfall experiment.After transforming the data, the obtained results were not exactly equals to the values of the published papers, but these calculations were necessary to make the graphics and commentaries similar.
For each site, soil erosion results are presented separately (Figures 2-8).First, bar and point graphs show the averages per intervals of 5 minutes of the suspended sediment load (SSL in g¨m ´2¨h ´1) and the suspended sediment load concentration (SSC in g¨L¨m ´2¨h ´1).The second graphic presents the surface flow (L¨m ´2¨h ´1), runoff (%) and infiltration (%) coefficients.Furthermore, tables for each study site were added to show: (i) the previous environmental plot characteristics during the experiments; and (ii) the averages pxq, totals, maximum (max) and minimum (min) values of the SSL, surface flow, SSC, runoff coefficient (%) and infiltration coefficient (%).To qualitatively compare the initial soil erosion processes with the different rainfall simulators characteristics and plot sizes, a final table with the results of the first 15 minutes of each rainfall simulations was presented.This time was chosen, because it represents the minimum duration of all the experiments.
Finally, a Spearman correlation test was performed to detect which factor shows the best trend (with the increasing or decreasing) with the SSL, surface flow and SSC (averages and maximum).

Ruwer-Mosel Valley (Trier, Germany)
Rainfall simulations at the Ruwer-Mosel valley were applied to study soil erosion processes before, during and after the vintage in conventional vineyards.Four rainfall simulations in August (2008) and five between September and December (2013) were performed.The experiments were done between August and September and coincided with tillage practices before the harvest.The rest of the simulations were done after the harvest, coinciding with the period of decreasing rainfall amounts and temperatures.Results of these rainfall simulation experiments were obtained with a rainfall intensity of 40 mm¨h ´1.
The plots were characterized by steep slopes (27.5 ˘5.6 ˝), with high vegetation cover (44.7 ˘33.8%) and high stone cover (57.8 ˘33.5%).The aim to this study was to measure the impact of rainfall and tillage practices on the soil before, during and after the harvest.
The average SSL was 1.25 ˘0.05 g¨m ´2¨h ´1 with a maximum value of 2.09 g¨m ´2¨h ´1 during the interval of 15-20 minutes.Surface flow was 0.21 ˘0.1 L¨m ´2¨h ´1 with a maximum value of 0.32 L¨m ´2¨h ´1.SSC showed higher values at the beginning of the experiment than at the end of the simulation, with an average of 6.8 ˘3.1 g¨L¨m ´2¨h ´1 and maximum of 11.7 g¨L¨m ´2¨h ´1.

Montagne de Reims (Champagne, France)
Rainfall simulations in Montage de Reims were applied in vineyards with three commonly used cultivation practices in Champagne: bare soil (1), bark and vine pruning (3), and grass cover (8).The goal was to quantify the influence of cultivation practices in the inter-row of vines and to determine the influence of the density of grass cover on soil loss and surface runoff.Twelve rainfall simulations were carried out during about 90 minutes with different rainfall intensities ranging between 20 and 76 mm¨h ´1 (Figure 3).
Environmental characteristics included a slope of 5.2 ˘0.8 ˝, a vegetation cover of 44.9 ˘43.2% and soil moisture of 19.4 ˘0.6%.Results were highly heterogeneous depending on the cultivation practice.Average SSL was 0.13 ˘0.1 g¨m ´2¨h ´1 with maximal values of 0.15 g¨m ´2¨h ´1.Average surface flow showed values of 0.15 L¨m ´2¨h ´1.0.94 ˘0.003 g¨L¨m ´2¨h ´1 was the average of the suspended sediment load concentration.Finally, in this case, we also observed high average infiltration (92.5 ˘2%) and low average runoff coefficient (7.5 ˘2%).

Els Hostalest de Pierola (Penedès)
The aim of the simulations was to evaluate the variability of infiltration along the slope due to disturbances created by land leveling operations.The simulation carried out in this experiment was quite similar to that performed in Málaga: six rainfall simulations (Figure 5) on three different points along the slope (top, middle and foot), with two replications in each of them.The simulations were done in June (2010) with a rainfall intensity ranging between 48 and 70 mm¨h ´1, with a duration of 60 minutes, although after 40 minutes a constant runoff volume was reached.slope of the study area was 9.8 ˘4˝w ith scarce soil cover (vegetation cover of 8.7 ˘21.2% and stone cover of 11 ˘1.6%).Soil moisture was 23.1 ˘11.2% and roughness was 1.7 ˘0.8%.
The average soil loss was 2.73 ˘0.86 g¨m ´2¨h ´1 with maximum values of 3.54 g¨m ´2¨h ´1.Surface flow was 0.21 ˘0.06 L¨m ´2¨h ´1.These results meant suspended sediment load concentrations of 13.27˘0.6g¨L¨m ´2¨h ´1, with the highest values being14.77g¨L¨m ´2¨h ´1.
Finally, runoff coefficient showed an average value of 43.8 ˘16% with a maximum value of 62.7% and infiltration averaged 56.3 ˘16% with values of up to 80%.

Moixent (Valencia, Spain)
The aim of this study was to measure the immediate effect of barley straw mulch on soil erosion and to detect the runoff processes.Sixteen rainfall experiments during the first days of July 2015 were conducted (Figure 6).Similar to the vineyard in Ruwer-Mosel valley, the experiment applied in this Mediterranean vineyard consisted of thirty minutes duration and intervals of five minutes to collect samples.This period corresponded with the driest period of the summer and before the vintage.Rainfall intensity was 40 mm¨h ´1.All experiments were carried out on lands with low slopes (2.5 ˘1.03 ˝), low vegetation cover (1.5 ˘2%) and low stone cover (21.1 ˘3.8%).
At the beginning, during the first 10 minutes, any runoff was noted.After that time, the average SSL was 6.86 ˘7.5 g¨m ´2¨h ´1 with a maximum of 16.36 g¨m ´2¨h ´1.The average surface flow was 1.17 ˘1.4 L¨m ´2¨h ´1 with a maximum value of 3.35 L¨m ´2¨h ´1.This resulted in an average SSC of 6.2 ˘0.9 g¨L¨m ´2¨h ´1 with a maximum of 6.91 g¨L¨m ´2¨h ´1.
In this study, in most cases, high infiltration rates (81.7 ˘21.3%) were observed.However, on a few occasions, runoff coefficients higher than 50% were also registered.

Villamediana de Iregua (La Rioja)
The of the rainfall simulations carried out in La Rioja was to measure the response of the soil to the effects of wheel tracks.Twenty-nine rainfall simulations were carried out in September, coinciding with the lowest soil moisture values (5-6%).The simulations had a rainfall intensity between <50 and >70 mm¨h ´1.
There was no vegetation cover, while a stone cover of 12.7 ˘16.2% was noted.The slope of the plots was 4.6 ˘3˝a nd with a SSW aspect.
Results did not show any surface flow and soil loss during the first five minutes of the experiment (Figure 7).However, after that time period, the average soil loss was 10.4 ˘8.3 g¨m ´2¨h ´1 with a maximum of 25.4 g¨m ´2¨h ´1.With respect to surface flow, 1.9 ˘1 L¨m ´2¨h ´1 (max. of 2.4 L¨m ´2¨h ´1) was observed.As a result, a suspended sediment load concentration of 4.1 ˘2.1 g¨L¨m ´2¨h ´1 was obtained.

Campo Real (Madrid)
Nine rainfall simulations were run on test plots with different land managements: (i) three with conventional tillage with tractors; (ii) three with secondary vegetation cover of Brachypodium distachyon; and (iii) three with annual barley (Hordeum vulgare L.).The experiments were conducted in summer after the spring soil tillage and in autumn before the fall soil tillage.The rainfall simulation durations were 15 minutes with intervals of one minute.Furthermore, measurements were also performed between 15 and 20 minutes to record the possible remaining runoff and soil loss after rainfall ended (inertia).Rainfall intensity was 130 mm¨h ´1.
Other envirommental plot characteristics during the experiments were 14 ˝slope, 40.5 ˘33.4% average vegetation cover and 17% stone cover.
Results (Figure 8) showed an average SSL of 2.13 ˘1.18 g¨m ´2¨h ´1 with a maximum of 3.15 g¨m ´2¨h ´1.Surface flow average was 2.54 ˘1 g¨L¨m ´2¨h ´1 with high values of 3.76 L¨m ´2¨h ´1.An average sediment concentration of 0.93 ˘0.6 g¨L¨m ´2¨h ´1 was noted.High infiltration (86.4 ˘7.6%) and a low runoff coefficient (13.6 ˘7.6%) were also observed.

Final Comparison
Different results were observed during the first 15 minutes of each rainfall simulation experiment (Table 3).High soil losses were observed in Montes de Málaga and Villamediana de Iregua (2.44 and 2.58 g¨m ´2¨h ´1, respectively).However, Waldrach, Moixent and Montagne de Reims obtained low rates due to the high infiltration coefficients (near 100%).Finally, high values of runoff were observed for Villamediana de Iregua and Els Hostalest de Pierola (between 22 and 54%).
Finally, a Spearman's rank correlation coefficient was calculated for these data to analyze the possible relationship between independent (plot, rainfall intensity and environmental conditions) and dependent (SSL, SF, SSC, RC, and IC) variables (Table 4).
First, it can observe that the plot size obtained a correlation between 0.519 and 0.593 with the surface flow (SF) and the runoff coefficient.Furthermore, correlation between an increase of the SSL and an increase in runoff and a decrease in infiltration were noted (0.607 and ´0.607, respectively).SF and SSC showed relationships with a decrease of the soil moisture inside the plot, ´0.522 and ´0.725, respectively.Finally, with a decrease in the percentage of slope and vegetation cover (´0.613 and ´0.667), an increase of SSC was noted.

Discussions
Different questions arose referring to the representativeness of rainfall intensity used in the experiments, the influence of the plot size on soil losses and runoff, or the duration of the experiments.It was also observed that there is a high variability in soil detached in each experiment, which ranged from an average value of 0.13 g¨m ´2¨h ´1 in Montagne de Reims (Champagne, France) to 10.4 g¨m ´2¨h ´1 in Villamediana de Iregua (La Rioja, Spain), but with maximum soil loss up to 25 g¨m ´2¨h ´1 (Villamediana de Iregua, La Rioja).The average surface flow was between 0.21 L¨m ´2¨h ´1 at Ruwer-Mosel valley (Trier, Germany) and 1.93 L¨m ´2¨h ´1 at Villamediana de Iregua (La Rioja, Spain).Differences were observed for the values of SSC.Els Hostalest de Pierola (Penedès, Spain) showed high values (13.3 g¨L¨m ´2¨h ´1).The lowest results were found in Campo Real (Madrid, Spain) and in Montagne de Reims with (Champagne, France) with 0.9 g¨L¨m ´2¨h ´1, although they applied high rainfall intensities.
After understanding the context of the rainfall simulation experiments of each investigation group (rainfall simulators, methods and study areas), homogenizing their data sheets in intervals of five minutes, and transforming all values into m 2 ¨h´1 , the problematic of soil erosion processes in European vineyards has been observed.
The experiments carried out under different soil tillage conditions and those carried out with different soil covers allowed us to confirm not only the effect of cover to reduce soil erosion [20], but also the utility of using rainfall simulations to evaluate these effects [39,57,60,61,64,66].Apart from these findings, one important question for the general discussion has not yet been clearly answered: how representative are the data for a general diagnostic about the most sustainable practices in European vineyards?Three causes for this problem could be proposed.
First, how many replicates must be done and where should the experiments be performed in an experimental area?High standard deviations in all experiments and an elevated variability of the trends have been shown.In many cases, runoff coefficients did not show a clear trend to delimit this parameter.Therefore, it may be important to carry out different repetitions or to increase the duration until a constant steady rate of runoff and suspended sediment load is achieved.
Second, the variability of soil surface components and soil tillage practices were determinant too.When and where should we conduct rainfall simulations: (i) before, during or after the vintage?(ii) On the top, middle or foot slope?(iii) In the rows or in the inter-rows?(iv) Which environmental characteristics must be considered on the plots?These responses would be possible: (i) to characterize at least all the different surface components (and their spatial and temporal variability during the year); (ii) to study with soil analysis if there are differences between the different parts of the slope and in the rows or inter-rows; and (iii) to control exactly for the work schedule of the vine-growers.Third, we have observed that after applying different rainfall simulators sizes and velocities, drop kinetic energy, plot forms and sizes, field of application, methodologies, etc.), homogenizing the intervals of the rainfall and understanding the concrete specifications from each study area, it would be difficult directly to conduct any quantitative diagnostics and comparisons.With the Spearman coefficient, for example, correlations between rainfall intensity and plot size were noted.Therefore, as Iserloh et al. observed [49,50], there is influence between the results and the technical characteristics of the rainfall simulators.It is clear that the most important reason for this study is to demonstrate the necessity of a standard rainfall simulator and methodology.
In the future, following this investigation, it is necessary to: (i) make quantitative analyses and comparisons between different study areas; (ii) calculate how many simulations and under which environmental and human conditions must be appropriated to develop rainfall simulation campaigns; (iii) spend time recollecting information, sampling methods and data treatments; and (iv) avoid interpretation errors and information losses after calculating averages and grouping the different intervals.

Conclusions
This collection of rainfall simulations in different contexts offers a wide range of aspects that can be analyzed concerning the magnitude of soil erosion processes in European wine-growing areas.The results shown in these rainfall simulation experiments allowed us to confirm that due to soil characteristics in which the experiments were done: (i) high infiltration and low runoff coefficients may be recorded, as well as the dynamic within the rainfall; and (ii) surface and sub-surface flow dynamics controlled the time at which runoff started and total runoff recorded; however, (iii) higher runoff rates have been indicated in other studies in the same vitivinicultural areas, when soil have higher silt contents.In addition to soil characteristics, the rainfall intensity and duration were additional control factors.

Figure 1 .
Figure 1.Localization of the study areas.
sediment load concentration (SSC) " SSL Surface flow

1 .
Montes de MálagaThe aim of this study was to compare rainfall impact on soil without vegetation cover and high stoniness.The rainfall simulations in Montes de Málaga were carried out in June (2003) in three different locations along a hillslope of 23 ˝in a conventional vineyard (Figure4): at the top (2), middle (2) and foot (2) slope.Two replications were done at each location.All experiments were carried out with 63-66 mm¨h ´1 rainfall intensity with a duration of 60 minutes.No vegetation cover and high stoniness (52.8%) characterized the plot.The average SSL was 3.13 ˘2.69 g¨m ´2¨h ´1 with maximum values of 6.85 g¨m ´2¨h ´1.Surface flow was 0.35 ˘0.15 L¨m ´2¨h´1(with maximum values of 0.62 L¨m ´2¨h ´1).The results gave rise to high SSC, (8.71 ˘0.56 g¨L¨m ´2¨h ´1), which reached 21.88 g¨L¨m ´2¨h ´1.Finally, runoff coefficients of 5.4 ˘2.5% and 94.6 ˘2.5% of infiltration were measured.

Figure
Figure Rainfall simulations in Campo Real (Madrid, Spain).

Table 1 .
Study areas with their characteristics.

Table 2 .
Rainfall simulators and methods.

Table 3 .
Rainfall simulation results during the first 15 minutes.

Table 4 .
Spearman correlation coefficient of whole dataset.