Applying the RUSLE and ISUM in the Tierra de Barros Vineyards (Extremadura, Spain) to Estimate Soil Mobilisation Rates

: Spain is one of the largest wine producers in the world, with Extremadura (south-west Spain) being its second-largest producing region after Castilla La Mancha. Within Extremadura, the most traditional and productive viticulture region is the Tierra de Barros, which boasts an annual production of 3 × 10 6 litres. However, no soil erosion assessment has been undertaken in any vineyard in the region to ascertain environmental sustainability. Therefore, the Improved Stock Unearthing Method (ISUM) and the Revised Universal Soil Loss Equation (RUSLE) were applied to assess the long-term soil erosion rates. Both methods were applied using an experimental plot (2.8 m × 148.5 m) encompassing 99 paired vines in a 20-year-old vineyard under a tillage management system and on bare soils throughout the year. The ISUM and RUSLE found total soil mobilization values of 45.7 Mg ha − 1 yr − 1 and 17.4 Mg ha − 1 yr − 1 , respectively, a di ﬀ erence of about 5 times. Mapping techniques showed that soil surface declined to an average of − 6.2 cm, with maximum values of − 28 cm. The highest values of soil depletion were mainly observed in the upper part and the form of linear features following the hillslope direction. On the other hand, under the vines, the soil surface level showed accumulations of up to + 2.37 cm due to tillage practices. Our study demonstrated the potential of high soil erosion rates occurring in conventional vineyards managed with tillage in the inter-row areas and herbicides under the vines within the Tierra de Barros. Also, we demonstrated the elevated di ﬀ erences in soil mobilisation rates using the ISUM and RUSLE. Therefore, further research must be conducted in other vineyards to determine the suitability of the models for assessing soil erosion rates. Undoubtedly, soil conservation measures must be designed and applied immediately due to high erosion rates.


Introduction
Soil erosion directly affects natural and anthropogenic resources, biogeochemical cycles, and ecosystem functions [1]. This environmental issue must be mitigated to achieve neutrality of land degradation [2,3] according to the Sustainable Development Goals of the United Nations [4]. For the ongoing issues in agriculture and land management, soil erosion is one of the most relevant concerns [5]. Lands have shown erosion rates of up to 10 Mg ha −1 yr −1 [6]. For example, high soil erosion rates were were applied to achieve this goal. These methods share a common feature of estimating the average soil mobilisation for a certain timeframe. The RUSLE is traditionally used with pre-established parameters such as rainfall erosivity, length of the slope, or soil management [37]. On the other hand, ISUM was recently updated from the Stock Unearthing Method (SUM) considering in situ measurements of the graft union as a passive indicator of the soil surface level. Therefore, the application of both methods will shed light on the soil erosion rates in the Tierra de Barros wine production region to help assess the environmental sustainability of the vineyards.

Study Area
A representative vineyard located at 38 • 50 26" N / 06 • 32 49" W in the so-called Tierra de Barros in the province of Badajoz (Extremadura, south-west Spain) was selected. The study area is situated on the farm called La Agraria, a few kilometres away from the municipality of Arroyo de San Serván ( Figure 1A). The experimental area is a conventional vineyard with vines mounted on a trellis with a support irrigation system and a total extension of about 20 ha. The espalier system allowed grape production three years after exposing an American rootstock graft union, which was planted after the Phylloxera crisis. The studied vineyard seems to have a complex topography, and the investigated portion has a double slope direction, with a concave profile ( Figure 1B,C). The parent materials are gravel and cobles with a fluvial origin. The soil texture varies among several horizons: (1) loamy (Ap), (2) loamy-clay (Bt and 2Bt, argic), and (3) loam-silty (Ck, calcic). The soil is classified as a Calcic-Luvisol [38] ( Figure 1D). The observed soil profile has a useful soil depth of 142 cm. On the surface, there is moderate stoniness (5%-15%), being more abundant in the Ap horizon, where it can reach up to 50% with gravel from 2 to 6 cm. The average inclination on this farm is 7.1%. The climate is Mediterranean with an average annual temperature of 16.5 • C (Arroyo de San Serván meteorological station; REDAREX, 38 • 51 29" N/06 • 28 22" W). January (8.0 • C) is the coolest month and July (26.9 • C) is the warmest. November is the wettest month (55.9 mm) and July the driest (4.4 mm); the average annual total precipitation is about 425.2 mm. The research site is located within a 20-year-old vineyard plantation with the Cayetana grape variety (local variety), which follows a planting frame of 2.8 m × 1.5 m ( Figure 1E,F). Each row (210 m long) contains approximately 140 vines. The whole farm receives the same tillage treatment (three times a year) depending on the humidity conditions, trying to preserve the soil bare all year round with herbicides. To obtain a correct representation of the amount of soil mobilisation in the study area, a total of 2871 measurements were taken along the 99 opposite paired-vines ( Figure 2).

The Improved Stock Unearthing Method (ISUM)
The ISUM is a low-cost and easy-to-apply method for estimating the current soil surface and the soil mobilisation rate considering the graft union of the vine as biomarkers at the time of planting and involves taking measures in the inter-row areas [36]. The graft of the vine indicates the initial surface level that relates to the soil at the time of planting. The grafts of the strains are only a few centimetres from the ground to avoid plant losses by freezing or by contact with herbicides ( Figure 2A). Also worth noting is that the vertical growth of the graft is negligible since it is the new vine that undergoes the evolution. In our study site, after interviewing the farmers and owners that planted the field, it was confirmed that the graft unions were located 2 cm above the surface. Possible errors from 1 to 2 cm due to extreme rainfall events or modifications by the tractors after the plantation took place are assumed in this research but are relatively small. For instance, an error of up to 3.3 cm was calculated with a similar method by Biddoccu et al. [39]. The farmers explained that the error is very small because of the gentle slopes, flat terrain due to tillage before planting the vines, and previous estimations using GPS and total stations. Moreover, we also visited recent plantations close to the study area in the same viticultural region to confirm the 2 cm above the surface. For this reason, we are "estimating" and not "quantifying" soil mobilisation rates. Similar errors were also researched in vineyards using the ISUM in other Spanish vineyards, for example in the La Mancha or Valencia regions. The vines were planted from the nursery while housed within a plastic cylinder between 30 and 40 cm in length and 10 cm in diameter that prevented damage by rabbits and protected the plants from herbicides. This protection makes it easier to find the graft unions.
To perform the ISUM measurements, a meter tape was stretched between each opposite paired-vine and measurements were taken every 10 cm along with tape ( Figure 2B). To avoid problems with buried graft unions, another meter tape was used to take the measurements along the inter-row area transect at 20 cm above the graft. Some concerns could still exist concerning this methodology because we did not provide any information on the correction of the slope of the tape used for measuring the soil surface height in the inter-row areas. However, it was confirmed by the farmers that the soil was flattened and that the paired vine rows are on the same slope inclination, so it was not necessary to apply a slope correction. In this work, 29 measurements were taken between each opposite paired-vine, resulting in a total of 2871 measurements along the 99 opposite paired vine rows ( Figure 2C-F). In this way, greater accuracy can be obtained [40] than in measurements recorded in previous works [36] where parts of the surface between the rows were not measured.
We applied the erosion-deposition (ER) equation proposed by Paroissien et al. [41] to estimate the total soil mobilisation rate (Equation (1)): where the volume (Vol, m 3 ), the total area of the measured field (St; ha), the years of the vines (Av, 20 years), and the apparent bulk density (BD; gr cm −3 ) were used as input elements. The soil bulk density is 1.39 g cm −3 , which represents the average of measurements collected in a 100 cm 3 cylinder at 12 different points along the row and the inter-row area. The method allows for an estimation of soil mobilisation in the two areas of "rows" and "inter-rows" with a considered width of 150 cm and 10 cm, respectively.

The Revised Universal Soil Loss Equation (RUSLE) Model
A wide range of parametric models that allow estimation of soil loss can be incorporated into a geographic information system (GIS). The most widely used method is the Revised Universal Soil Loss Equation [35] which was derived from the Universal Soil Loss Equation (USLE) previously proposed by Wischmeier and Smith [37]. The RUSLE model is composed of five factors expressed as (Equation (2)): A = R * K * LS * C * P where A is the average rate of soil loss (Mg ha −1 y −1 ), R is the rainfall erosivity factor (MJ mm ha −1 h −1 yr −1 ), K is the factor of soil susceptibility to erosion, LS is the factor of length (m) and magnitude of the slope (%), C is the cover and management factor, and P is the factor of soil conservation practices.
To create the RUSLE map, we used the RUSLE tool implemented in ArcGIS 10.5 (ESRI, Colorado, USA). It is necessary to add a Digital Elevation Model (DEM) (http://centrodedescargas.cnig.es/ CentroDescargas/index.jsp) and the values of each of the above-mentioned factors. In Spain, a DEM with a resolution greater than 5 meters is not available. Therefore, the available 5-m DEM was converted into points and then the 3-m DEM (3 m × 3 m pixels) was derived using the Topo to Raster technique implemented in the ANUDEM software (ANU Fenner School of Environment and Society and Geoscience, Australia) developed by Michael Hutchinson.
The R factor has been obtained based on the adapted universal equation of soil loss defined by ICONA [42] as: where e represents the basis of the Napierian algorithms, PMEX is the average value of the annual series of maximum monthly rainfall (mm) using the meteorological station of Extremadura Irrigation Advisory Network (REDAREX, 38º51% 29" N/06º28 22" W), MR is the average precipitation of the October-May period (mm), and F24 is the ratio between the annual maximum 24-h rainfall squared and the sum of the maximum 24-h rainfall of all months of the same year.
The K factor has been calculated using the equation proposed by Wischmeier and Smith [30] as: where OM represents the percentage of the organic matter determined through the method proposed by Walkley and Black [43]. M is the product of the percentage fractions of sand, clay, and silt. The granulometric analysis of the samples was done with the laser diffraction particle analyser (Beckman Coulter LS 13 320), using the liquid aqueous module with an optical model of Fraunhofer.rf780d. s and p are, respectively, the soil structure and permeability classes. We used the nomogram established by Wischmeier and Smith [37] to determine the s and p parameters. We used a RUSLE tool implemented in ArcGIS, which requires only one single DEM as input, and automatically calculates the LS factor based on the USPED (Unit Stream Power-based Erosion Deposition) equation: where U is the upslope area per unit width (a measurement of water flow) in meters (m 2 /m), β is the slope angle in degree, 22.1 is the length of the standard USLE plot in meters, and 0.09 = 9% = 5.15 • is the slope of the standard USLE plot. The range of values of the exponents is for m = 0.2−0.6 and n = 1.0−1.3, where the lower values are used for prevailing sheet flow and higher values for prevailing rill flow. The value of the C factor suggested by ARPAV (Agenzia Regionale per la Prevenzione e Protezione Ambientale del Veneto) for vineyard crops was used in this study. In the case of the P factor, the value proposed by Panagos et al. [6] for Spain was adopted. Table 1 shows the RUSLE parameter values used.

Mapping the Current Soil Surface Level for ISUM
The same person took the measurements between the initial distance of the graft unions and the current soil surface in May 2018 to avoid any bias ( Figure 2G,H). The data obtained from the measurements were incorporated into ArcMap version 10.5 software (ESRI, Colorado, USA) for its representation. First, a grid point was created with the "fishnet" tool using a shapefile format. Subsequently, the Digital Elevation Model (DEM) was generated through interpolation methods using the Geostatistical Wizard extension incorporated in ArcGIS. Ordinary Kriging (OK) emerged as the best interpolation method in terms of the lowest Root Mean Square Error (RMSE) and the highest coefficient of determination (R 2 ), after applying geostatistical tools to obtain the most accurate scale for the map. Table 2 compares the performance statistics of the mean, RMSE and R 2 . A final map was generated with a resolution of 0.03 m pixel size that allows an accurate visualisation. Since the possible error can range from 1 to 2 cm, two maps were also included with ±1−2 cm to show this second possible error detected after performing the OK. These calculations also included the total soil mobilisation rates. The main goal was to obtain the potential flow paths and soil mobilisation patterns along the studied area after extreme rainfall events ( Figure 2I,J).

Current Soil Surface using ISUM
The average, median, and the 5th and 95th percentile values obtained between each opposite paired-vine (n = 99) are depicted in Figure 3. The average soil surface is −6.24 cm, which represents a soil mobilisation of 3.12 mm yr −1 . The maximum soil depletion reaches up to −28 cm. On the other hand, the maximum soil accumulation is +19 cm. Considering the soil surface level between each paired-vine, the highest soil depletion is observed from the intervals of 60 cm to 250 cm with an average value of −7.98 cm and a maximum value of −28 cm. In the intervals 0-50 and 260-290, the average values are close to 2.40 cm, reaching values up to −22 cm. As expected under the influence of tractor wheels and tillage, the highest soil depletion values are found in the inter-row areas and the highest accumulations occur close to the vines.
Paying attention to different intervals, the highest accumulations and the lowest soil depletions happen within the intervals 0-50 and 260-290 cm. The average accumulation values reach +14.76 cm, with maximum values higher than +19 cm. Meanwhile, the maximum average values reaching +8.3 cm are recorded from the intervals of 60 to 250 cm. The highest absolute maximum value was higher than +18.6 cm and located between the intervals 210 and 250 cm. A map of the current soil surface level of the study area generated using the Ordinary Kriging procedure is presented in Figure 4A. Figure 4B depicts the longitudinal profile curve of the topographical changes. It is observed in Figure 4A that soil accumulation occurs in the lower areas. The highest values of soil depletion are observed in the upper reaches and the form of linear features following the slope direction. Under the vines, the soil surface level accounts for an average accumulation value of +2.37 cm. However, soil accumulations between 20 cm and 25 cm are observed in the middle and the foot reaches the plot. The maps considering the possible error of 1 and 2 cm are also included in Figure 5. In these maps, we can observe that the mean average values of soil depletion increased, showing even higher soil mobilisation processes along with the studied plot.   Figure 6 and Table 3 present the erosion map generated using the RUSLE tool implemented in ArcGIS and the total area per mobilised soil rates, respectively. The results show values of soil losses ranging from 0 to 45.8 Mg ha −1 yr −1 in the studied row and inter-row areas, and close to them. The mean values were 17.4 Mg ha −1 yr −1 . These values are not homogeneous in the study area as we observed that there are differences among hillslope positions. In some areas, coinciding with the ISUM map, we observed similar hotspots with the highest soil mobilisation reaching up to 45 Mg ha −1 yr −1 . If we regard the study area in a more general context, we can observe that the paired-rows and inter-row areas are crossed by a linear feature similar to a gully or ephemeral river channel. This situation could explain the general soil mobilisation along with the plot, which is noted against the tractor pass direction. The spatial erosional pattern that was identified by RUSLE is comparable with the results obtained by the ISUM: high erosion in the location corresponding to the gully/ephemeral gully channel along the middle of the studied area.    Table 4 shows the soil mobilisation estimates using the ISUM and RUSLE methods. The differences between the ISUM and RUSLE models are considerable, obtaining soil mobilisation rates of 45.7 Mg ha −1 yr −1 and 17.4 Mg ha −1 yr −1 , respectively. Moreover, through the use of ISUM, the values in various units of erosion for total soil loss and differences between the row (0-50 cm) and inter-row (50-240 cm) areas are also summarised in the table. In the inter-row areas, a soil mobilisation of 45.1 m 3 is found, which corresponds to 62.7 Mg ha− 1 for 20 years. Conversely, under the vines, soil mobilisation reaches 5.9 m 3 , which corresponds to a total of 8.2 Mg ha −1 during the same period.

Discussion
This study is based on a single vine row with no spatial replication, thereby making it necessary to remark that this is a preliminary study for an initial comparison of ISUM and RUSLE. We cannot pretend to establish strong robustness in terms of representativeness of the whole vineyard. The main aim is to show whether both methods can or cannot be used to model soil mobilisation rates and make the first measurements in this viticultural area. It is necessary to increase the number of measuring points in the rows and inter-rows in the future to check the robustness in terms of representativeness.
Our research applied the traditional and well-established RUSLE methodology for estimation of soil losses. The ISUM was also used as a more recent approach based on direct measurements in the field. The ISUM and RUSLE show considerably different soil erosion rates of 45.69 Mg ha −1 y −1 and 17.4 Mg ha −1 y −1 , respectively, for the study area of Arroyo de San Serván in the Tierra de Barros, Extremadura. In the literature, soil erosion estimates using the RUSLE are generally significantly different, sometimes more than one order of magnitude from those of other models. A clear example is a study carried out by Busacca et al. [44] who found an average erosion rate of 31.4 Mg ha −1 yr −1 using the RUSLE while the Cs137 technique gave an estimated reduction of 11.6 Mg ha −1 y −1 . Also, the RUSLE gave erosion rates higher than other laboratory techniques or in situ experiments with rainfall simulation or USLE plots. It was observed that the RUSLE performed poorly in estimating the long-term erosion rates in the Tierra de Barros. Our data showed a soil lowering during the last 20 years, and this fits better with the observed soil depletion during the fieldwork and recent soil mobilisation that occurred during extreme rainfall events that happened during the last year of this study ( Figure 2I,J). Other studies apply rainfall simulation on experimental plots, which do not show the whole picture of the spatial and temporal variability of soil erosion [45].
In Spain, the use of the RUSLE and direct measurements have shown contrasting responses. López Bermúdez [46] found that the USLE, as applied in the Segura river in Southeastern Spain, registered 105.5 Mg ha −1 yr −1 , while the reservoir bathymetries obtained a value of 7.5 Mg ha −1 yr −1 , on average, subsequently reaching an overestimation of 14 times. Soto et al. [47] applied the USLE to all the watersheds in Spain and estimated a soil loss of 23.44 Mg ha −1 yr −1 , whereas the reservoir bathymetries yielded 14 times lower soil loss.
The overestimation of the USLE method is more apparent when comparing the whole dataset developed by the researchers. López Bermúdez [46] found that the USLE measured values from 7 to 302 Mg ha −1 yr −1 at the Spanish basins, while bathymetry measured values from 2 to 4 Mg ha −1 yr −1 , and the soil erosion plots indicated values between 0.08 and 3.1 Mg ha −1 yr −1 . Thus, the USLE overestimated soil erosion rates in both cases. Bagarello et al. [48] attributed these discrepancies to scale issues. Among other authors, Jardí et al. [49] estimated soil erosion rates in road embankments of 122-864 Mg ha −1 yr −1 with USLE, while plot measurements yielded a rate of 0.9 Mg ha −1 yr −1 . This can be attributed to the fact that USLE shows an average rate for a long period but the plots only show the period of measurement, and this can cause the soil erosion rates to differ by two orders of magnitude over a decade.
Another example of overestimation of soil losses by USLE is given by Avellanas et al. [50] for measurements in pine and oak forest (burnt and control) in Zaragoza, Spain. Using the USLE they estimated soil losses of 0, 1.61, 11.52, and 64. We found that the USLE and ISUM estimated different soil erosion rates at the "Tierra de Barros" study site. Previous research in the Basque Country estimated soil loss of 16-70 Mg ha −1 yr −1 with the USLE and 18 to 90 Mg ha −1 yr −1 by erosion pins. However, once the forest was cut, the soil erosion rates increased to 36-116 Mg ha −1 yr −1 with USLE and 55-105 with erosion pins [53], which can lead to an overestimate of soil erosion rates. Due to its ease of use and low cost, the ISUM has been applied for many different goals such as to assess differences among parent materials, vineyards with different ages or slope positions. This method allows for determining in great detail the micro-topographical variations of the vine plantations, which can be helpful to design correct soil management systems and where to apply them.
It is confirmed that the key factor for obtaining the best resolution in the generation of the final model depends on the number of measurements [40]. Therefore, in this study, measurements were carried out every 10 cm by the same person. Moreover, we have to keep in mind some limitations and potential sources of errors during the field campaign, including considering the distance of the graft to the surface, since at the time of planting not all strains were at the same distance and in the same direction. On the other hand, the RUSLE model has been used for decades in numerous studies, consolidating itself as one of the most widely used methods in estimating soil losses. The difference concerning the use of ISUM lies mainly in the scale of measurement and the parameters that are required for the generation of the models.
The RUSLE model is a quick and easy application method in many agricultural areas, but the obtention of some factors such as rainfall erosivity or erodibility can be complicated without a long-term database. The results obtained with the RUSLE model in this study do not show great differences regarding other research works developed in vineyards in the Mediterranean basin or elsewhere. Soil erosion rates are determined by a multitude of factors, and range from 6.4 Mg ha −1 yr −1 to more than 22 Mg ha −1 yr −1 in sloping vineyards, decreasing in time from 19.46 Mg ha −1 yr −1 to 11.28 Mg ha −1 yr −1 in old vineyards or between 53.9 and 69.5 Mg ha −1 yr −1 depending on the mechanisation practices [54,55]. As already mentioned, the resolution (and also the scope) of the two methods are very different, especially as applied in this work. Thus, we include both maps and estimations to compare the average values (if considered representative of the whole selected vineyard) with the average RUSLE value, as a validation of the RUSLE. The RUSLE is very sensitive to some factors (e.g. LS or C factors), so using general factors to estimate soil erosion on a small plot could lead to very misleading results. In this research, we intend to show the readers the obtained results and how we can improve this research in the future.
The soil erosion rates recorded with the ISUM in this study differ from those found in the literature. The most important difference between the volumes of soil loss could be related to the number of measurements between the opposite pair-vines, since the years of the vines in many cases are similar. In vineyards of 25 years of age, and taking five measurements, the erosion rates were estimated as 2.5 Mg ha −1 yr −1 [36]. In another work, the rates of erosion appear conditioned by the source of material which reached an erosion rate of −87.7 Mg ha −1 yr −1 in lands dominated by loams in an 8-year old vineyard [54]. However, the differences are reduced if we compare our results with those of Rodrigo-Comino et al. [40], in which different numbers of measurements were compared, and the erosion rates for 18 points were −40.1 Mg ha −1 yr −1 . The high soil erosion rates observed in this study requires that policymakers, farmers, and companies involved in wine production develop erosion control strategies to achieve environmental sustainability. Increasing the use of vegetation covers and mulches such as straw or cover crops, and reducing the use of herbicides, should be urgently considered.
It is pending to apply the ISUM and RUSLE for vineyards plantations with organic management systems to confirm that the use of soil erosion control measures is effective to achieve land degradation neutrality [56]. Previous research carried out with other methodologies demonstrates that organic farming results in lower erosion rates (due to soil quality improvement) as compared to chemical farming, mainly due to biota activation, and the reduction in the connectivity of the flows [57,58].

Conclusions
The ISUM and RUSLE were applied to assess soil erosion rates in one of the most important European viticultural regions: Tierra de Barros, Extremadura (south-west Spain). We concluded that this vineyard in the Tierra de Barros shows non-sustainable soil erosion rates under its conventional soil management system characterised by the use of machinery to till the soil and the application of herbicides to keep the soil bare. The ISUM and RUSLE estimated soil mobilisation rates of 45.7 Mg ha −1 yr −1 and 17.4 Mg ha −1 yr −1 , respectively. Generally, RUSLE erosion estimates are more than one order of magnitude higher or lower in comparison with other soil erosion measurement methods. In our case, the ISUM and RUSLE did not obtain similar rates but indicated similar patterns, although the RUSLE is a very sensitive model to some topographical and cover factors, which greatly affect a direct comparison. Both mapping techniques showed that the soil surface reached on average Funding: The study was financed by Junta de Extremadura and FEDER (Project IB16052). The landscapes of Extremadura inspired this research.