Comparing Hydric Erosion Soil Loss Models in Rainy Mountainous and Dry Flat Regions in Portugal

: Soil erosion is a severe and complex issue in the agriculture area. The main objective of this study was to assess the soil loss in two regions, testing different methodologies and combining different factors of the Revised Universal Soil Loss Equation (RUSLE) based on Geographical Information Systems (GIS). To provide the methodologies to other users, a GIS open-source application was developed. The RUSLE equation was applied with the variation of some factors that compose it, namely the slope length and slope steepness (LS) factor and practices factor (P), but also with the use of different sources of information. Eight different erosion models (M1 to M8) were applied to the two regions with different ecological conditions: Montalegre (rainy-mountainous) and Alentejo (dry-ﬂat), both in Portugal, to compare them and to evaluate the soil loss for 3 potential erosion levels: 0–25, 25–50 and >50 ton/ha · year. Regarding the methodologies, in both regions the behavior is similar, indicating that the M5 and M6 methodologies can be more conservative than the others (M1, M2, M3, M4 and M8), which present very consistent values in all classes of soil loss and for both regions. All methodologies were implemented in a GIS application, which is free and available under QGIS software.


Introduction
Soil erosion is a severe and complex issue in the agricultural area, and can be defined as a process of detachment of particles that are transported or deposited in other areas by water or wind actions. The Food and Agriculture Organization of the United Nations (FAO) in its main report starts with "Soils are fundamental to life on Earth" [1], while other research indicates that erosion is endangering several areas such as natural resources, flooding and affecting water reservoirs [2,3]. Given this, it is fundamental to assess and control the dimension of erosion and its socio-economic and environmental impacts.
Several environmental factors affect soil erosion rates, such as the characteristics of the soil, the slope, the precipitation, the vegetation coverage and type, etc. The slope is one of the most influential factors in soil erosion estimation [4].
In Portugal, as in many other countries of the European Mediterranean, soil erosion is one of the main environmental problems affecting the sustainability of agro-forestry systems [5]. The soil loss problem is concentrated in arable crop areas due to the adopted exploitation of agricultural management practices and intense farming that increases soil loss [6]. Several studies have been performed in the Mediterranean watersheds, since this region is particularly prone to erosion, due to the long dry periods, heavy erosive rainfall events, steep slopes, low organic matter, low nutrient contents, and rapid landuse changes [7][8][9]. For instance, Rellini et al. [9] applied the RUSLE model in Portofino In this work, we apply the RUSLE equation with the variation of some factors that compose it, namely the factors LS and P, but also with the use of different sources of information (for R factor). GIS and Remote Sensing (RS) data/techniques have been used to retrieve the different factors. Several studies used GIS software to produce erosion risk maps to assess soil loss [15,[25][26][27][28][29][30].
The different erosion models were applied to two study cases, Montalegre and Alentejo (in Portugal), to compare them, to evaluate the soil loss and answer this question: which methodology is more appropriate to infer soil erosion risk in a rainy mountainous zone (Montalegre) and in a dry and flat zone (Alentejo)?
The open-source GIS application was developed using an open-source programming language, and it can be adapted to other countries by adjusting the rules. The application is free and open-source and available to any user to use and improve it.

Methodology
The developed application relies upon the implementation of eight methodologies combining different factors. Each methodology allowed the generation of a different final erosion risk map, considering the different combination of the factors involved in the RUSLE equation. The RUSLE equation allows estimation of the mean annual soil loss based on six erosion factors (Equation (1)) [12]: where A is the computed estimation of the average annual soil loss (ton/km 2 ·year), R is the erosivity factor, calculated over one year (MJ mmk/m 2 ·h·year), K is the soil erodibility factor (ton·km 2 ·h/km 2 ·MJ·mm), LS (unitless) is the combination of slope length (L) and slope steepness (S), C (unitless) is the land cover factor, and P is the practices used for erosion control (unitless, ranging between 0 and 1). In this work, the areas with a high risk of soil water erosion were delimited based on the identification of potential soil erosion, through the application of Equation (2), which considers the factors of soil erosivity (R), soil erodibility (K), slope length and steepness (LS), and practices used for erosion control (P) translated into the expression: Equation (2) is widely used, with different methodologies used to obtain the R, K and LS values. In this work, we tested and compared four different methodologies to obtain the LS factor and two different sources to obtain R values. The K values were expressed based on Pimenta [31] tables and the P value was also considered in the potential erosion risk, which reflects the agricultural practices that support the terraces and/or earth-retaining walls, with a crucial role in the reduction of the soil erosion. In all the methodologies applied, the annual soil loss was estimated from Equation (2), multiplied by 0.6, 0.8 or 1, corresponding to the P-factor value. The following sub-sections present the different sources and methodologies employed to estimate the factors.

R Factor
The R values were estimated based on two different sources: based on the precipitation erosivity map obtained from the SNIAmb of the Portuguese Environment Agency, I.P and based on JRC work [32]. In the first source, the information used to determine the R factor was taken from the map available in SNIAmb of the Portuguese Environment Agency, I.P (Precipitation erosivity-estimated based on 449 50.8 mm udometric stations-R. Wischmeier [4]), which is found in units of the English System (American foot/acre·tons). The R map is provided in raster format with 250 m of spatial resolution and in the Lisbon Datum Hayford Gauss with false origin-Militar Coordinates (EPSG: 20790). As the units of the map are in the English System, they should be converted to the International System (IS). This conversion was performed based on Foster et al. [23]. In JRC work, the R factor was estimated based on Panos Panagos et al. [33]. The values were computed from the precipitation data with different temporal resolutions and then normalized to the R factor with 30 min of temporal resolution using linear regression. The Gaussian Process Regression (GPR) model was used to interpolate the precipitation erosivity values of the stations.

K Factor
The K factor can be obtained using soil characteristics, referring to sand content, silt, clay, organic matter, permeability and structure. Pimenta [31] estimated the soil erodibility values based on the Wischmeier and Smith [4] abacus. In this abacus is presented the clay percentage, the sand percentage, the organic matter percentage, the structure and the permeability. To obtain the soil erodibility, Pimenta [31] established Equation (3): where K is the K factor (ton·km 2 ·h/km 2 ·MJ·mm); M is the parameter that defines the particle size correspondent to the product (% of lime + fine sand) × (100 − %clay); OM is the organic matter (%); s is the soil structure; and p is the permeability.

LS Factor
The LS factor was tested and compared using four different methodologies, based on Wischmeier and Smith [4], Mitasova [13], Moore and Wilson [17] and McCool et al. [24], respectively.
The LS factor calculation based on Wischmeier and Smith [4] is estimated from Equation (4) where LS is the slope length and steepness parameter (unitless), A is the slope length (in meters), and θ is the slope inclination (in degrees). The m value varies from 0.2 to 0.5, depending on the slope percent: slope ≥ 5%, m = 0.5; slopes of 3.5 to 4.5%, m = 0.4, slopes of 1 to 3, m = 0.3, slope < 1% or uniform gradients, m = 0.2. The LS factor can also be obtained from Mitasova [13], according to Equation ( where A is the slope length, θ is the slope inclination (in degree), and m and n are physical coefficients. The annual soil loss can be estimated according to Equation (2), but using the LS estimated from the Mitasova [13] model. The third approach to estimate the LS factor was based on Moore and Wilson [17] and translated by Equation (6): where flow_accumulation is the raster corresponding to the flow accumulation, and θ is the slope inclination (in degrees). In the last methodology tested, the LS factor was based on McCool et al. [24]. In this estimation, for slope values less than 9%, the LS was estimated according to Equation In slope values higher than 9%, the LS factor can be estimated as Equation (8): Land 2021, 10, 554

of 18
In both equations, the m constant value is estimated from Equations (9) and (10): where A is the slope length (in meters), θ is the slope inclination (in degrees), m is a constant value that depends on slope value, and β is the ratio between runoff erosion and erosion caused by the impact of raindrops. Table 1 presents the four methodologies used to estimate the LS factor and the different sources used to estimate the R and K values. In Table 1, each equation has one factor R, K and LS. The objective of this study is to combine all the factors (obtained from different methodologies) in the two regions of study, resulting in eight combinations (2(R × K) × 4LS). Then, these combinations will be analyzed in terms of three levels of P values, 1, 0.8, and 0.6, to infer the impact of conservative practices in each erosive model tested ( Table 2). Table 2. Eight methodologies tested considering the different factors.

GIS Application
An open-source GIS application was already developed to estimate the annual soil loss [30]. However, this application implemented the RUSLE equation according to the methodology proposed by Wischemeier and Smith [4]. In this work, a new open-source GIS application, Erosion Risk, was developed under QGIS 3.18 Zurich software [34], implementing different methods for each factor involved to estimate soil loss, giving the user the possibility of choosing the best method to use (see Supplementary Materials). The application was developed using the Python programming language and several algorithms from Processing Toolbox packages, such as GRASS, SAGA, and GDAL, among others [35][36][37]. In addition, modules from Qt API, QGIS API and PyQGIS API were used to develop the code [38,39]. Due to the authors' experience in developing QGIS plugins and due to the integration of several external algorithms in this software [40], QGIS was chosen to implement the application. The Plugin Builder was used to start the creation of the plugin. The Erosion Risk application is composed of a single button that opens the main window, and the GUI was created through Qt Designer framework. The main window is composed of a toolbox widget, which provides a series of pages or compartments in a toolbox [39]. Figure 1 presents the Erosion Risk main window (with the LS factor page opened).
The application was developed using the Python programming language and se algorithms from Processing Toolbox packages, such as GRASS, SAGA, and GDAL, a others [35][36][37]. In addition, modules from Qt API, QGIS API and PyQGIS API were to develop the code [38,39]. Due to the authors' experience in developing QGIS pl and due to the integration of several external algorithms in this software [40], QGIS chosen to implement the application. The Plugin Builder was used to start the creati the plugin. The Erosion Risk application is composed of a single button that open main window, and the GUI was created through Qt Designer framework. The window is composed of a toolbox widget, which provides a series of pag compartments in a toolbox [39]. Figure 1 presents the Erosion Risk main window the LS factor page opened). Each one of the factors that composes the RUSLE equation was distributed b pages/compartments of the main window. In this regard, the toolbox is composed o pages, including a final one to create the final erosion risk map: R, K, LS, P and RU On each page, the methods previously described for the LS factor were implemented the user can choose which one they intend to use. In R factor, two sources are availa define the input: the R values based on the SNIAmb precipitation erosivity map or b on JRC; in K factor, if the soil characteristics are available, they should be defined i characteristics Table; on the another hand, if these characteristics are not availabl user can confer the Pimenta [31] tables available in the application (through Open K button) and assign the rating values; in LS factor, the four methodologies previ defined are available to choose in a combo box; in the P factor, a table is availab assign the P values to the input land use land cover map. The open-source application is free and the code is available to users, and the developers can chan according to other methodologies and adapt it to other regions. Each one of the factors that composes the RUSLE equation was distributed by the pages/compartments of the main window. In this regard, the toolbox is composed of five pages, including a final one to create the final erosion risk map: R, K, LS, P and RUSLE. On each page, the methods previously described for the LS factor were implemented and the user can choose which one they intend to use. In R factor, two sources are available to define the input: the R values based on the SNIAmb precipitation erosivity map or based on JRC; in K factor, if the soil characteristics are available, they should be defined in Soil characteristics Table; on the another hand, if these characteristics are not available, the user can confer the Pimenta [31] tables available in the application (through Open K tables button) and assign the rating values; in LS factor, the four methodologies previously defined are available to choose in a combo box; in the P factor, a table is available to assign the P values to the input land use land cover map. The open-source GIS application is free and the code is available to users, and the developers can change it according to other methodologies and adapt it to other regions.

Case Study
The hydric erosion methodologies involved in the application were tested in two different regions in Portugal: (i) mountainous and rainy region (Montalegre, North of Portugal); and (ii) dry and flat region (Alentejo, South of Portugal; Figure 2). The hydric erosion methodologies involved in the application were tested in two different regions in Portugal: (i) mountainous and rainy region (Montalegre, North of Portugal); and (ii) dry and flat region (Alentejo, South of Portugal; Figure 2). In Montalegre, the grassland system includes permanent semi-natural meadows (called lameiros) and common lands, which represent around 70% of the agricultural area [41], and bovine livestock is the main agricultural activity. In this region the Atlantic climate favors high precipitation occurrence (1531 mm/year), mainly occurring from autumn to spring, autumn-winter temperatures lower than 12 °C and mean monthly temperatures ranging from 3.5 °C to 17.2 °C. The high precipitation values occurred from autumn to spring, combining with high slope values and promoting soil loss. According to Koppen's classification [42] for the period (1980-2010), the region belongs to group Csb (temperate, with dry summer, which is not very hot but extensive), while Thornthwaite's rational climate classification [43] describes it as A B′1 s b′4 super-humid (hydric index 137; A), primer-mesothermic (thermal efficiency index: 622 mm evapotranspiration; B'1), with a moderate shortage of water in the summer (aridity index: 16.8; s) and moderate thermal efficiency summer concentration index: 48%; b′4 (above 20% ≥ typically continental).
The Alentejo study zone was chosen as having opposite characteristics to Montalegre. The study zone is located in the South of Portugal, and it is a protected zone, a National Hunt zone and Forestry Perimeter, that belongs to the western part of the Serra Morena, in the municipality of Moura, Beja District, Portugal (Figure 2). It is named Herdade da Contenda (HC), and it extends along 19 km of the Portuguese-Spanish border with an area of 52.7 km 2 . In this region, the elevation values range from 246 m to 543 m, with the South and Southeast part presenting the highest elevation values (Figure 2). The  Figure 2). In Montalegre, the grassland system includes permanent semi-natural meadows (called lameiros) and common lands, which represent around 70% of the agricultural area [41], and bovine livestock is the main agricultural activity. In this region the Atlantic climate favors high precipitation occurrence (1531 mm/year), mainly occurring from autumn to spring, autumn-winter temperatures lower than 12 • C and mean monthly temperatures ranging from 3.5 • C to 17.2 • C. The high precipitation values occurred from autumn to spring, combining with high slope values and promoting soil loss. According to Koppen's classification [42] for the period (1980-2010), the region belongs to group Csb (temperate, with dry summer, which is not very hot but extensive), while Thornthwaite's rational climate classification [43] describes it as A B 1 s b 4 super-humid (hydric index 137; A), primer-mesothermic (thermal efficiency index: 622 mm evapotranspiration; B'1), with a moderate shortage of water in the summer (aridity index: 16.8; s) and moderate thermal efficiency summer concentration index: 48%; b 4 (above 20% ≥ typically continental).
The Alentejo study zone was chosen as having opposite characteristics to Montalegre. The study zone is located in the South of Portugal, and it is a protected zone, a National Hunt zone and Forestry Perimeter, that belongs to the western part of the Serra Morena, in the municipality of Moura, Beja District, Portugal (Figure 2). It is named Herdade da Contenda (HC), and it extends along 19 km of the Portuguese-Spanish border with an area of 52.7 km 2 . In this region, the elevation values range from 246 m to 543 m, with the South and Southeast part presenting the highest elevation values (Figure 2). The annual average temperature is around 17.0 • C (maximum 24.3 • C and minimum 9.3 • C). The mean annual precipitation value is 729.5 mm with January and February presenting the highest values. During summer, the precipitation presents very low values.
The two regions have climates that clearly distinguish themselves, since Montalegre is a humid and rainy region, while Alentejo is a dry and flat region (Figure 2).

Data Acquisition
The data were obtained from different sources with different coordinate systems assigned. All the data were projected to the European Terrestrial System 1989 Portugal Transverse Mercator 2006 (ETRS89-PTTM06) coordinate system. The obtained information was also limited to the study zone, using algorithms from QGIS, such as Clip layer by mask. This process was implemented in the application.
In order to obtain the R factor, the raster file with the precipitation erosivity values available in the SNIAmb da Agência Portuguesa do Ambiente, I.P. [44] was used. The precipitation erosivity was estimated based on 449 udometric stations with rainy events of 50.8 mm. A rainy event is preceded by an interval of at least 6 h without the occurrence of precipitation [4]. The data is in American tons/acre units. In this case, the units should be converted to the IS corresponding to MJ·mm/ha·h.ano. Table 3 presents the data sources for R, K and LS factors used in this work. The K values, corresponding to soil erodibility, can be obtained from the soil map and land suitability map for agriculture, when available. In the cases that this information is not available, the user can assign the P values according to Pimenta [31] tables, which are available in the GIS application. A Digital Elevation Model (DEM), obtained from Shuttle Radar Topography Mission (SRTM) [45] with a spatial resolution of 25 m, was used to derive the LS factor. The data was resampled to 25 m in both regions, Montalegre and Alentejo, following the mosaic created from DEM.

Data Processing
The maps were created using the GIS application and considering the different methodologies adopted for the LS factor, the different sources for R values and with different P values. The final erosion risk maps were classified in three categories: 0-25; 25-55; and >55 (ton/ha·year). In order to evaluate the results obtained, the percentages of each area in each of the three categories were obtained, using the QGIS r.mapcalc.simple tool. Some statistical variables such as average, standard deviation and confidence interval were also estimated. Additionally, to compare the results obtained between the two regions, and between the different categories, several graphics were created and analyzed.

Results and Discussion
The different methodologies were applied to each study area to compare them, to evaluate and answer the question "which methodology is more appropriate to study the soil erosion risk in a mountainous and rainy zone such as Montalegre municipality and, on the other hand, which methodology is more appropriate for a dry and flat zone such as Alentejo"? The final maps were created for the eight methodologies proposed in this study. Figures 3 and 4 present the maps obtained.

Results and Discussion
The different methodologies were applied to each study area to compare them, to evaluate and answer the question "which methodology is more appropriate to study the soil erosion risk in a mountainous and rainy zone such as Montalegre municipality and, on the other hand, which methodology is more appropriate for a dry and flat zone such as Alentejo?" The final maps were created for the eight methodologies proposed in this study. Figures 3 and 4 .  Figures 3 and 4 illustrate the spatial distribution of the potential soil erosion risk map of Montalegre municipality and Alentejo region, respectively. The annual average soil loss was grouped in three categories: 0-25 (green), 25-55 (yellow) and >55 (red) (ton/ha·year). The maps were created using different combinations between eight methodologies. Additionally, different P values were assigned (1, 0.8, and 0.6). Overall, in the Montalegre municipality, which is characterized as a rainy and sloped region, the category >55 (ton/ha·year) was mainly observed in the Southwest zone of the municipality in all maps of soil erosion rate for all the methodologies, emphasized by the M5 and M6 methodology. This result is expected since these zones are influenced by higher slopes, and areas with high slope values are related to a high level of erosion [30]. The category with less soil erosion, 0-25 (ton/ha·year), mostly covered the Northeast zone of Montalegre. This was also verified in all the methodologies.
In the Alentejo region, a region with smother slopes when compared to Montalegre, the soil erosion rates were mostly identified in the South/center zone of the area. In methodologies M5, M6 and M7, the area with moderate/high erosion risk is not as evident.
In the remaining methodologies, the spatial distribution of the soil erosion rates is similar. As can be observed from Figures 4 and 5, the spatial variations in length-slope factors by the different methodologies applied have significant effects on the total erosion estimation. Regarding Figures 3-5 we considered for evaluation two approaches: (i) for each region, we evaluated the average value (area) in each category, in all methodologies and the methodologies that are above or below the average value, to understand the most appropriate methodology to be used in each region; and (ii) considering both regions, we evaluated if the results present the same tendency/pattern.
Considering the overall tested methodology, in Montalegre (82,330 ha), 33.8% of the land surface was classified as having low soil erosion (0-25 ton/ha·year), 15.6% was classified as having a moderate soil erosion risk (25-55 ton/ha·year) and 50.6% of the area was classified as having high soil erosion risk (>55 ton/ha·year). In the first and in the second categories (Figure 5a,b), is possible to verify that all the methodologies are above the average value, except for M5 and M6, indicating that these methodologies are more  Table 4 presents the statistical values estimated from the percentage of area, such as average, standard deviation, minimum value, maximum value, and confidence interval, regarding the two different regions. Table 4. Statistical variables for the three levels (0-25, 25-55 and >55 ton/ha·year) including: the total area value (in ha) followed by the percentage area ± confidence interval, the standard deviation, minimum value, and maximum value, for each region, Montalegre and Alentejo. Considering the average of all tested methodologies, the soil loss in the Montalegre region is 27,830 (33.8 ± 4.8) ton/ha·year, where 33.8% of soil losses occur in the category of 25 to 55 ton/ha·year. In this region and according to all methodologies, between 42.9% (35,320 ha) and 58.3% (47,998 ha) of the total area are at high risk of erosion (>55 ton/ha·year). In Montalegre, the areas with the lowest risk of erosion (class 0 to 25 ton/ha·year) represent 11.3% (9303 ha) to 19.9% (16,384 ha).

Soil
Considering the Alentejo region, between 34.4% (1662 ha) and 49.8% (2405 ha) are considered to be in the low erosion risk areas. In the category of moderate soil erosion (25-55 ton/ha·year), 800 (16.6 ± 7.8) ton/ha·year is the amount of soil loss and 27.8% (1343 ha) to 54.8% (2647 ha) of the areas are at high risk of erosion. Figure 5 presents the graphics generated with the percentage of area in each category in the two different regions (bars) and the average value for each category in the two regions, Montalegre (orange line) and Alentejo (grey line).
Regarding Figures 3-5 we considered for evaluation two approaches: (i) for each region, we evaluated the average value (area) in each category, in all methodologies and the methodologies that are above or below the average value, to understand the most appropriate methodology to be used in each region; and (ii) considering both regions, we evaluated if the results present the same tendency/pattern.
Considering the overall tested methodology, in Montalegre (82,330 ha), 33.8% of the land surface was classified as having low soil erosion (0-25 ton/ha·year), 15.6% was classified as having a moderate soil erosion risk (25-55 ton/ha·year) and 50.6% of the area was classified as having high soil erosion risk (>55 ton/ha·year). In the first and in the second categories (Figure 5a,b), is possible to verify that all the methodologies are above the average value, except for M5 and M6, indicating that these methodologies are more conservative, compared with M1, M2, M3, M4, M7 and M8. Considering the third category (Figure 5c), where the erosion risk is higher, the opposite of the previous conclusion is observed; M5 and M6 are lower than the average value.
Considering the Alentejo region (4830 ha) and all methodologies, 42.1% of the areas have low soil erosion risk (0-25 ton/ha·year), 16.6% have moderate soil erosion risk (25-55 ton/ha·year) and 43.3% of the area was identified as having high soil erosion risk (>55 ton/ha·year). In the first category (Figure 5a), M5 is above the average value, indicating that this is a more conservative methodology. Considering the second category ( Figure  5b), M5, M6 and M7 are above the average value, and in the third category ( Figure 5c) the opposite happens, where M5 is below the average, and M6 and M7 are in the line of the average value; the remaining methodologies are above.
Considering the relation between both regions, Montalegre and Alentejo, it is expected that the Alentejo region presents lower erosion values compared to Montalegre, a rainy region, especially in the third category. Considering this category, in M1, M3, M5 and M7, it is observed that its percentage of area in Alentejo is lower than in Montalegre. In the remaining methodologies, this does not happen. In these methodologies, the R values were obtained from the precipitation erosivity values available in the SNIAmb da Agência Portuguesa do Ambiente, I.P [44] and estimated based on Wischemeier and Smith [4]. The opposite is expected to be observed in the first category, where the percentage of areas are related to zones with lower soil erosion risk, meaning that Alentejo areas should be higher than Montalegre. Observing Figure 5a, this hypothesis is confirmed again in the M1, M3, M5 and M7 methodologies. The same situation as before is observed in this category, meaning that the R value source can influence the results in the methodologies.
Regarding all methodologies, in both regions the behavior is similar, indicating that M5 and M6 methodologies are more conservative than the others. The referred methodologies are based on Moore and Wilson [17] LS factor estimation.
The K factor, which expresses the susceptibility of a soil to erosion, is related to soil properties such as organic matter content, soil texture, soil structure and permeability and stone covers. The mean value of soil erodibility USLE-K factor in the water erosion area was calculated as 0.0321 with a standard deviation of 0.0107 [46]. However, Mediterranean countries (Italy, Spain, Greece, and Portugal) have mean K factor values between 0.039 and 0.042 (ton·ha·h/ha·MJ·mm), which is higher than the mean erodibility in Europe [47].
The obtained mean K value of water erosion areas in Montalegre was 0.032 and was concentrated in the range of 0.028-0.036 (from K values obtained from method seen in [4]) and 0.020 in the range of 0.013-0.035 (from K values obtained from the JRC method). In the Alentejo region, the obtained mean K value was 0.039, which ranged between 0.029-0.039 according to K values obtained with Wishchmeier and Smith [4], and the mean value 0.040 ranged between 0.034-0.047, according to the JRC method. In the Alentejo region, the K value is higher than in the Montalegre region.
Regarding the two methodologies for calculating the K factor, these present different behaviors of potential erosion according to the regions. The K factor based on JCR methodology, when compared to the Wischemeier and Smith [4] methodology, presents higher values of potential erosion in the Alentejo region (lower slopes and less rainy), but more conservative erosion values in sloping and rainy regions (Montalegre).
This difference of K values among regions could be related to the incorporation of the protective effect of surface stone cover, which is usually not considered for soil erodibility calculations. The effect of stone covers in K-factor calculations is likely to result in an overestimation of soil erosion, particularly for Montalegre, where the highest percentages of surface stone cover were observed. Soil erodibility, together with management practices (P factor) and vegetation cover (C factor) can be influenced by agricultural practices. Therefore, the K factor estimation can provide important information for mitigation of soil erosion through the implementation of improved conservation practices (e.g., increase or preserve soil organic carbon in areas prone to high levels of soil erosion risk or adaption of soil management at areas of high risk) [46].
Considering the question proposed in the objectives: which methodology is more appropriated to infer on soil erosion risk in a rainy mountainous (Montalegre) and in a dry and flat zones (Alentejo), within this work it is possible to conclude that the M1, M2, M3, M4 and M8 methodologies present very consistent values in all classes and in both regions. The maps resulting from M5 and M6 were very different from the previous group, with M7 in an intermediate situation.

Conclusions
The presented study provides an open-source GIS application, Erosion risk, developed to produce soil erosion maps implementing different methodologies and combinations of factors. The application was developed for QGIS software, using open-source programming language and libraries. The GIS application is free and respects open-source freedoms, providing the possibility to be improved and/or adapted to other methodologies. Four methodologies for the creation of the LS factor were implemented in the application and were tested using two regions with different characteristics: a rainy and mountainous region in the North of Portugal and a dry-flat region in the South of Portugal. The combination with different sources of R and P values allows taking in consideration some scenarios providing options to policy makers for managing soil erosion issues most efficiently for prioritization of regions with different characteristics. The erosion soil maps created allowed us to assess soil loss through the implementation of different methodologies that combined different factors of the RUSLE equation.
In the future, new functionalities will be incorporated in the open-source GIS application, such as: (i) tools that allow the user to dynamically create the maps involved and perform comparisons at the moment of the creation; (ii) new methods for R, K, and LS factors will be added to each page; and (iii) statistical functionalities will also be added to compare, for instance, areas using statistical measures such as average, standard deviation and confidence interval, as performed in this work. After the validation and testing of these methodologies regarding other situations, a web GIS application was developed by the authors and these methodologies will be incorporated in the web GIS application, in order to reach more users.
To conclude, further research in the area of soil erosion estimation (with particular focus on the existent LS estimation methods) using GIS is required.