Assessing and Modeling Soil Detachment Capacity by Overland Flow in Forest and Woodland of Northern Iran

Land use has significant effects on the erosion process, since it influences the soil detachment capacity by causing an overland flow (Dc). The effects of different land uses on the rill detachment capacity have not been explained in depth, and the hydraulic parameters providing accurate estimates of this soil property have not been completely identified. This study quantifies Dc at low flow rates in woodland and forestland, compared to two other land uses (cropland and grassland), in the Saravan watershed (Northern Iran), and develops prediction models of Dc and rill erodibility (Kr). Dc was measured on undisturbed soil samples, collected in the four land uses, and characterized in terms of the main physico-chemical properties in a flume experiment, simulating five slopes and five shallow water flows. The results showed that Dc was significantly lower in woodland and forestland compared to cropland and grassland, as the consequence of the changes in the main soil properties and the more developed vegetation cover and structure. Dc was positively correlated to clay and silt contents of soils, and negatively correlated to sand content, aggregate stability, root density, and organic matter. The stream power and unit stream power were found to be very accurate predictors of Dc in woodland and forestland, respectively. Kr values, which assumed the lowest values in woodland and forestland, were provided by interpolating Dc and the shear stress of water flow. Overall, this study has confirmed that vegetation cover and improved soil properties in forestland and woodland may help to reduce erosion in delicate environment ecosystems, such as the forests of Northern Iran.


Introduction
Soil erosion consists of detachment, transportation, and deposition of sediments due to rainfall and surface runoff [1,2]. Soil detachment is the separation of soil particles from the matrix over the watershed (Northern Iran). We hypothesize that the more developed vegetation cover and structure of woodland and forestland are able to reduce the rill detachment capacity, compared to the other land uses. Moreover, regression equations are proposed to predict D c and K r from hydraulic parameters under controlled conditions of flow rate and terrain slope.
The study is a contribution to understand the relationships between the soil erodibility and both the soil properties and the hydraulic parameters under different land uses.

Study Area
In Iran, the anthropogenic activities (e.g., deforestation, intensive agriculture, farm abandonment), which play important influences on the natural soil pedogenesis, have been very intense in the past 50 years. In the northern part of this country, the reduction of forestland cover, due to the expansion of arable land and intensive cultivation, population growth, and technological development [41], has led to soil degradation in large areas [42]. Therefore, this region seems to be suitable for the study aims, since degradation of soil quality and increased erosion rate have been the main results of land use changes.
The study area is the Saravan Forestland Park, which is one of the oldest forestlands in Guilan province. The experiments were carried out in the Saravan watershed, which covers an area of 14.87 km 2 (outlet coordinates 37 • 08 04 N, 49 • 39 44" E) between 50 and 250 m above the mean sea level ( Figure 1).
Forests 2020, 11, x FOR PEER REVIEW 3 of 21 structure of woodland and forestland are able to reduce the rill detachment capacity, compared to the other land uses. Moreover, regression equations are proposed to predict Dc and Kr from hydraulic parameters under controlled conditions of flow rate and terrain slope. The study is a contribution to understand the relationships between the soil erodibility and both the soil properties and the hydraulic parameters under different land uses.

Study Area
In Iran, the anthropogenic activities (e.g., deforestation, intensive agriculture, farm abandonment), which play important influences on the natural soil pedogenesis, have been very intense in the past 50 years. In the northern part of this country, the reduction of forestland cover, due to the expansion of arable land and intensive cultivation, population growth, and technological development [41], has led to soil degradation in large areas [42]. Therefore, this region seems to be suitable for the study aims, since degradation of soil quality and increased erosion rate have been the main results of land use changes.
The study area is the Saravan Forestland Park, which is one of the oldest forestlands in Guilan province. The experiments were carried out in the Saravan watershed, which covers an area of 14.87 km 2 (outlet coordinates 37°08′04′′ N, 49°39′44″ E) between 50 and 250 m above the mean sea level ( Figure 1). The climate of the area is typically Mediterranean, Csa type, according to the Köppen-Geiger classification [43]. The mean annual temperature and precipitation are 16.3 °C and 1360 mm, respectively [44].  The climate of the area is typically Mediterranean, Csa type, according to the Köppen-Geiger classification [43]. The mean annual temperature and precipitation are 16.3 • C and 1360 mm, respectively [44].
Forestland, cropland, grassland, and woodland in the Saravan watershed were the land uses under investigation in this study. The representative tree and plant species of the forestland are Carpinus betulus, Parrotia, willow, and pussy willow. Woodland is covered by smaller trees with lower density than forestlands. Grassland is covered by herbaceous vegetation and shrubs. Cropland consists of planted black berry, greengage, and sloe, subjected only to planting and harvesting as tillage operations.

Soil Sampling and Analyses
D c was measured between August and September 2018 after collecting undisturbed soil samples of the four land uses in areas with similar slope gradients and elevation. In cropland, soil was sampled from flat patches.
The sampling procedure was adopted in accordance with the work of Zhang et al. [4,10]. A steel ring (0.1 m in diameter and 0.05 m in height) was used to extract the samples from soil. Before sampling, rocks, weeds, and litter were removed over the soil surface. In other to ensure minimum disturbance of the soil sample, roots and other debris surrounding the steel ring were clipped, while gently pressing the ring into the soil. After that, the steel ring was inserted into the soil and the sample bottom was carefully trimmed to remove the excess soil. Two cotton cloths were used to cushion the top and bottom of the steel ring, in order to avoid disturbance while the sample was transported to the laboratory.
Bulk density, aggregate stability in water, organic matter content, and root weight density were determined on additional 100 soil samples (25 samples × 4 land uses), randomly collected in the same location of each land use. Soil bulk density and aggregate stability were measured by the oven-drying and the wet-sieving methods, respectively. Soil organic matter was estimated by using the potassium dichromate colorimetric method. Root weight density was measured using a washing method over a sieve with a 1-mm mesh and subsequently roots of each soil sample were weighted after oven-drying for 24 h at 65 • C.

The Laboratory Flume
To generate rill flows, a laboratory scale flume was used ( Figure 2). Combinations of soil samples of each land use and different values of flow rate and terrain slope were simulated (hereinafter each combination will be indicated as "experiment").
More detailed information about the experimental flume is reported in Asadi et al. [45,46]. To summarize, this flume, with a rectangular cross section, was 3.5-m long and 0.2-m wide. In order to reproduce the natural roughness of the soil, a 5-mm layer was put on the flume bed [14].

The Experimental Design
For each experiment, a sample of soil was packed in a steel ring (0.1 m in diameter, 0.05 m in height). The surface of each soil sample, sieved through a 2-mm mesh, was wetted by light spraying. The wetted sample in the steel ring was then placed in a hole (0.1 m in diameter) of the flume bed, located far from the flume outlet [10], making sure that the sample upper surface was flushed by water flowing over the flume bed surface [14,18]. Then, the flume slope and flow rate were adjusted to the desired values. The slope of the flume, set manually, could range between 0 and 38%. Water was poured into the flume from an upstream tap and the flow rate measured at the outlet, as explained in Section 2.6.
Before starting the experiment, a cover panel was used to prevent soil samples from scouring. Once the water flow started and stabilized, the panel was removed and D c was measured together with the hydraulic parameters of the water flow. Each experiment was carried out for no more than five minutes. In order to reduce the influence of uneven detachment within the sample, the time of the experimental test was controlled by the scouring depth (0.02 m) for each soil sample [14,18]. The experiment ended when the depth of the eroded soil in the steel ring was 0.015 m. Finally, the wet soil of the sample was oven-dried at 105 • C for 24 h to measure its dry weight.

Determination of the Rill Detachment Capacity and Erodibility
During each experiment, Dc [kg s −1 m −2 ] was calculated as the average value of four replicates, using the following equation: Soil detachment in rills occurs when shear stress of the water flow acting on the soil (τ) exceeds the critical shear stress (τc) and the sediment load is lower than sediment transport capacity [17].  For each land use, five water flow rates (0.22, 0.33, 0.44, 0.56, and 0.67 L m −1 s −1 ) and as many slope gradients (3.4%, 8.5%, 17.6%, 26.8%, and 36.4%) were simulated in the flume. Overall, 400 soil samples (4 land uses × 5 water flow rates × 5 slope gradients × 4 replications) were subjected to the experiments.

Determination of the Rill Detachment Capacity and Erodibility
During each experiment, D c [kg s −1 m −2 ] was calculated as the average value of four replicates, using the following equation: where: • Soil detachment in rills occurs when shear stress of the water flow acting on the soil (τ) exceeds the critical shear stress (τc) and the sediment load is lower than sediment transport capacity [17]. According (2)

Determination of the Hydraulic Parameters
When the flow stabilized into the flume, the water depth was measured in two cross sections (placed 0.4 m and 1 m upstream of the flume outlet, respectively), using a level-probe with an accuracy of 1 mm. For each cross section, three points were considered for these measurements, of which one in the center of the flume and two others at 0.01 m from each side of the flume. Based on these six measurements, the average value of the flow depth was calculated for each experiment.
In accordance with the equations developed from several authors [47] for streams, rivers, and open channels, a rectangular cross section was assumed for the rills and its hydraulic radius (R, [m]) was calculated as follows: where h (equal to 0.2 m in this study) and p [m] are the flow width and depth, respectively. The flow rate per unit width (q, [L m −1 s −1 ]) was measured five times per experiment, collecting the water into a graduated plastic cylinder. The flow velocity (v, [m s −1 ]) m was estimated using a fluorescent dye technique in ten replicated measurements. The water temperature (T, [ • C]) was also recorded in order to calculate its viscosity. Based on these data, the Reynolds number (Re, [-]) was calculated. Furthermore, the mean water flow velocity (V, [m s −1 ]) was determined reducing the water velocity by 0.6, 0.7, or 0.8, when the flow was laminar, transitional or turbulent, respectively [48].

Statistical Analysis
Hypothesizing that the samples were independent and the distribution was normal, the t-test was used to assess the statistical significance of the soil properties among the land uses (at p level lower than 0.05). The normal distribution hypothesis was checked using QQ-plots.
Then, a Principal Component Analysis (PCA) was applied to the soil parameters, in order to find correlations (using Pearson's method), [50] among the soil properties, to identify the existence of meaningful derivative variables (the Principal Components, PCs) and to group soil samples according to the investigated land uses.
The relationships between D c (considered as the dependent variable) and the hydraulic parameters (i.e., V, τ, ω, φ and E, assumed as independent variables) were analyzed by a non-linear regression method based on power functions, according to the literature e.g., reference [25]. The accuracy of the regression equations was evaluated using quantitative indexes, as the relative error (RE), the coefficient of determination (r 2 ) and the coefficient of efficiency (NSE) [51]. The optimal values of these indexes are one for r 2 and E as well as zero for RE, while the acceptance limits of r 2 and NSE are 0.50 and 0.35, respectively. The prediction capacity of a model is considered to be good when r 2 and E are higher than 0.75 [52][53][54].

Variability of Soil Properties and Detachment Capacity among Land Uses
The characterization of soil texture shows that all samples were clay loamy (according to SSDS, 2017, classification) [55], with about 33%-37% of clay, 47%-49% of silt, and 14%-20% of sand ( Table 2). The clay and sand contents were slightly different in cropland compared to the other land uses, although not to a significant degree (p < 0.05).
In general, the other soil properties varied among the investigated land uses. The cultivated soil had the highest bulk density and the lowest aggregate stability, root weight density and organic matter content. All these properties were significantly different from those detected in the other land uses ( Table 2).
The aggregate stability was the maximum possible level in woodland (about 3-fold and 2-fold the values measured in cropland and grassland, respectively) and slightly lower (−25%) in forestland (although not significantly). The soils sampled in woodland and forestland showed the highest root weight density (with a maximum of 0.63 ± 0.03 kg m −3 in forestland) and organic matter content (with a maximum of 1.87 ± 0.02% in forestland). The differences in these soil properties were statistically significant compared to the two other land uses ( Table 2).
As shown in Figure 3, the rill detachment capacity was significantly influenced by the land use. On average, this soil property was the highest in cropland (0.049 kg m −2 s −1 ) and the lowest in forestland (0.008 kg m −2 s −1 ). D c of grassland and woodland was closer to the latter land use rather to cropland. The rill detachment capacity of cropland was 2-fold, 4-fold, and 6-fold greater than in grassland, woodland and forestland, respectively. Moreover, the variability of this soil properties was larger than in the other land uses (standard deviation of 0.045 kg m −2 s −1 against 0.024, 0.01, and 0.005 kg m −2 s −1 of grassland, woodland, and forestland) (Figure 3).  The analysis of Pearson's matrix shows that rill detachment capacity was positively correlated with silt and clay contents of soil and negatively correlated with bulk density (although not significantly), root weight density, sand content, aggregate stability, and organic matter. Many other significant correlations were also found among the soil properties, except for bulk density, The analysis of Pearson's matrix shows that rill detachment capacity was positively correlated with silt and clay contents of soil and negatively correlated with bulk density (although not significantly), root weight density, sand content, aggregate stability, and organic matter. Many other significant correlations were also found among the soil properties, except for bulk density, which was not correlated with D c and organic matter (Table 3). Using PCA, two principal components (PCs) were identified, explaining together about 80% of the total variance of the soil properties. In more detail, PC1 explained 64% of this variability, whereas PC2 explained another 15%. D c and all soil properties (except for the bulk density) had a significant loading on the first PC, while the bulk density had a significant weight on the second PC (Table 4).
In other words, D c was associated with low values of organic matter, root weight density, sand content, and aggregate stability of soil, as well as with high values of silt and clay contents (see PC 1 ), while no direct association was found between D c and the bulk density (Table 4 and Figure 4a).   Table 4 reports the power equations correlating the rill detachment capacity with individual hydraulic parameters, while the values of the indexes evaluating the prediction capacity of these equations are shown in Table 5. In general, all the equations simulated Dc quite well, as shown by the satisfactory values of r 2 (over 0.50, except for the equation Dc-E in woodland) and NSE (always over 0. 35). NSE was only negative in the regression Dc-τ for cropland, (which means that the model Plotting the sample scores on the first two PCs, evident differences in soil properties emerged among land uses. Three well-differentiated groups, one for each land use, were evidenced besides a clear overlapping of points corresponding to FL and WL (Figure 4b). Table 4 reports the power equations correlating the rill detachment capacity with individual hydraulic parameters, while the values of the indexes evaluating the prediction capacity of these equations are shown in Table 5. In general, all the equations simulated D c quite well, as shown by the satisfactory values of r 2 (over 0.50, except for the equation D c -E in woodland) and NSE (always over 0. 35). NSE was only negative in the regression D c -τ for cropland, (which means that the model performance is very poor) while the RE was very high (>85%) in this case. The best predictor of D c was the stream power for cropland, grassland and woodland, since, by using these equations, the best combination of RE, r 2 and NSE (always good) was achieved, while the rill detachment capacity was better simulated by the unit stream power for forestland (RE = −0.1, r 2 = 0.92 and NSE = 0.92) ( Table 6). Conversely, the D c values were more scattered around the regression line using flow velocity, shear stress, and unit energy compared to stream power and unit stream power ( Figure 5).  Table 6. Values of the criteria adopted for evaluating the accuracy of equations in Table 5  performance is very poor) while the RE was very high (>85%) in this case. The best predictor of Dc was the stream power for cropland, grassland and woodland, since, by using these equations, the best combination of RE, r 2 and NSE (always good) was achieved, while the rill detachment capacity was better simulated by the unit stream power for forestland (RE = −0.1, r 2 = 0.92 and NSE = 0.92) ( Table 6). Conversely, the Dc values were more scattered around the regression line using flow velocity, shear stress, and unit energy compared to stream power and unit stream power ( Figure 5). Different linear regression equations with expression [2] were established between D c and τ for the different land uses, always with high coefficients of determination (r 2 > 0.90) between D c and τ ( Figure 6). Cropland showed the maximum rill erodibility (0.0092 s m −1 ) and critical shear stress (1.71 Pa). The K r value for cropland was 1.8, 4.6, and 9.2 times greater than in grassland, woodland, and forestland soil, respectively. Thus, the rill erodibility was the lowest in forestland (0.001 s m −1 ), while woodland showed a slightly higher value (0.002 s m −1 ). The critical shear stress was higher in cropland (1.71 Pa) and grassland (1.65 Pa), and lower in forestland (1.0 Pa) and woodland (0.70 Pa) ( Table 7). Notes: RE = relative error (%); r 2 = coefficient of determination; NSE = coefficient of efficiency of Nash and Sutcliffe.

Relationships between the Rill Detachment Capacity, Rill Erodibility and Hydraulic Parameters
Different linear regression equations with expression [2] were established between Dc and τ for the different land uses, always with high coefficients of determination (r 2 > 0.90) between Dc and τ ( Figure 6). Cropland showed the maximum rill erodibility (0.0092 s m −1 ) and critical shear stress (1.71 Pa). The Kr value for cropland was 1.8, 4.6, and 9.2 times greater than in grassland, woodland, and forestland soil, respectively. Thus, the rill erodibility was the lowest in forestland (0.001 s m −1 ), while woodland showed a slightly higher value (0.002 s m -1 ). The critical shear stress was higher in cropland (1.71 Pa) and grassland (1.65 Pa), and lower in forestland (1.0 Pa) and woodland (0.70 Pa) ( Table 7).

Impacts of Land Use on Soil Properties and Rill Detachment Capacity
Literature studies have demonstrated that land uses influence the physico-chemical and hydrological characteristics of a soil e.g., reference [33,35,[56][57][58][59]. It is also known that soil detachment capacity affected by overland flow varies widely under different land uses [26,60]. The characterization of the soil properties under the different land uses of this study confirms this variability [61][62][63]. The clay loam texture of the investigated soils was the same in all studied land uses, indicating the homogeneity of soil forming processes and the similarity of parent materials [61]. Therefore, the plant survival and growth in woodland, forestland and grassland, combined with the agricultural practices used in cropland, may explain the variability of the soil properties and particularly the variations detected for the rill detachment capacity.
Agricultural cultivation and its associated biomass extraction, and the use of machinery are responsible for the decrease in organic matter and the compaction of the soil. The higher bulk density values and the lower aggregate stability of the cultivated soil could be due to these actions. Organic matter content is significantly lower compared to the values measured in forestland and woodland and is in close accordance with other recent studies [26,27,64]. As a matter of fact, it is well known that soils with a noticeable vegetation cover and structure-such as woodland and forestland of this study-show a higher content of organic matter [65][66][67], which, in turn, increases root weight density, micro-porosity, and macro-porosity and thus increases the aggregate stability of soils. Plant roots bind soil particles at the soil surface [29,68], increasing the soil infiltration rate [69] and providing additional surface roughness [70].
The changes in soil properties also reflect its hydrological parameters, such as the higher water infiltration capacity as well as the reduced runoff generation ability and low soil erodibility in vegetated soil [71][72][73]. The latter property can be directly linked to the rill detachment capacity, which, as has been shown in this study, was significantly influenced by the land use. The lower rill detachment capacity among the four land uses was detected in forestland and woodland soil, corroborating outcomes of some recent studies e.g., references [4,7,27,40]. As a matter of fact, in cropland the surface soil was disturbed by farming operations (e.g., weed mechanical removal, tillage), which often left the soil bare and disrupted the aggregates, thus increasing its potential for erosion. Conversely, the higher vegetal cover and structure of forestland and woodland determined a lower rill detachment capacity. This is in accordance with many authors e.g., references [74][75][76][77][78], who concluded that vegetation cover is an important factor in the reduction of soil erosion, although some studies suggested that vegetation cover might not be so important or is even insignificant for rill erosion [7,69].
The correlation analysis and PCA basically confirm the relationships between some important soil properties and its erodibility. As a matter of fact, these statistical techniques have demonstrated that, for the investigated land uses, the rill detachment capacity is directly associated with the organic matter, root density, and aggregate stability, as well as the textural properties of sampled soils. In other words, when the organic matter (and thus the aggregate stability and root density) as well as the sand content decrease, the rill detachment capacity increases. Conversely, the influence of bulk density on rill detachment capacity is not evident under the experimental conditions. Moreover, in soils with higher silt and clay contents, particle detachability increases (as shown by the positive correlations between D c , silt and clay contents). As found by Li (2015), an increase in silt-in loess soils with already very high silt content-implies an increase in D c while cohesion works in the opposite direction, which is consistent with numerous results (e.g., references [4,60,64]). The susceptibility of soils with finer particles is an important outcome of this study, since it confirms that soils with prevalent clay or silt contents are more erodible compared to soils with a proper vegetal cover, such as woodland or forestland. Therefore, land management activities favoring an evolution from agriculture to forestland land use may be suggested and, at the same time, deforestation of woody areas should be avoided in order to avoid intolerable soil erosion rates. This is a worldwide priority [79] due to the control vegetation exerts on soil particle detachment.
As a response to the associations between soil properties and its erodibility found in this study, two land uses (cropland and grassland) were discriminated from woodland and forestland, as shown by the evident clusters arranged along a clear gradient on the score plot. This is in accordance with a previous investigation carried out in the same environment, which showed that soils with different land uses can be grouped in separate clusters, depending on the changes induced on their physico-chemical properties caused by land use [35]. There are close relationships between soil detachment capacity on one hand and soil bulk density, size, silt content, cohesion, and root density [27], and bulk density on the other [27,29].

Relationships between the Rill Detachment Capacity, Erodibility and Hydraulic Parameters
Beside the soil properties, soil detachment in rills depends on the hydraulic characteristic of the overland flow [18]. The regression analysis carried out in this study has demonstrated that the rill detachment capacity can be estimated from common hydraulic parameters using power equations with different accuracy levels, which are generally appreciable. Also, many other authors [18][19][20][21][22][23][24]80] showed that flow discharge, shear stress, stream power, unit stream power, and unit energy had significant positive correlations with the rill detachment capacity and thus they can be assumed as being D c predictors.
Moreover, the simple regression models found in this study were more accurate in forestland and woodland soils. The response curves of D c versus hydraulic parameters also showed greater proximity of woodland in relation to forestland than the other land uses. The study has shown that, when the best predictor should be chosen among the analyzed parameters, the stream power is the most representative hydraulic variable for cropland, grassland, and woodland, while the unit value of this parameter is suggested for forestland. Conversely, the use of the flow velocity, shear stress, and unit energy should be avoided to estimate rill detachment capacity, since in some cases this may lead to large overestimation or underestimation of the soil erodibility. The higher accuracy of stream power and unit stream power may be due to the fact that this hydraulic parameter simultaneously takes into account the flow velocity (such as V) and slope (such as E and τ), which both have a clear influence on soil detachability due to the overland flow.
The evaluation of the effects of land use on rill erodibility and critical shear stress is important, because K r and τ c are two of the most important parameters reflecting soil resistance to rill erosion [17]. It has been shown above that the lower rill detachment capacity detected in forestland and woodland, compared to the other land uses, may be due to the interactive effects of vegetation and soil properties, and this also influences the values of rill erodibility and critical shear stress, which are mathematically linked to the D c (as shown by the accuracy of the regression models developed in this study). This accuracy is consistent with findings of Zhang et al. [4], who stated that the linear regression functions between soil detachment and shear stress can be used to fit the observed data satisfactorily, with coefficients of determination of up to 0.9 for each land use. The lower rill erodibility and critical shear stress is in accordance with previous studies, showing that cropland is subject to the highest erosion rates among the investigated land uses [4,7,27,40].

Conclusions
Few studies have explained the quantitative relations among the soil detachment capacity due to overland flow and the properties of forest soils. Moreover, the methods for its estimation from hydraulic parameters are quite scarce in forestland and in woodland. This study has quantified the rill detachment capacity and proposed prediction equations of this parameter for woodland, forestland, grassland, and cropland. It has been demonstrated how and to what extent the rill detachment capacity is affected by the different land uses in the case study of Saravan watershed in Northern Iran. Although the textural characteristics of the investigated soils were quite similar among the four studied land uses, the hydraulic parameters and the physico-chemical properties were significantly different. The mean rill detachment capacity of soils was lower in woodland and forestland than in cropland and grassland. This outcome confirms the working hypothesis that the more developed vegetation cover and structure of woodland and forestland are able to reduce the rill detachment capacity, compared to other land uses. Rill detachment was positively correlated with clay and silt contents and negatively correlated with sand content, aggregate stability, organic matter, and root density, while no statistically significant correlations were found with bulk density.
For all land uses, a power function was used to estimate the rill detachment capacity from the hydraulic parameters and a generally good accuracy (in particular for woodland and forestland) was achieved for the proposed equations. The best predictor of rill detachment capacity was the stream power for cropland, grassland, and woodland, and the unit stream power for forestland. The use of flow velocity and unit energy should be avoided when estimating the rill detachment capacity. The rill erodibility and critical shear stress, which were the lowest in forestland and woodland out of all of the investigated land uses, were predicted with good accuracy along with the slope and intercept of linear regressions between the rill detachment capacity and shear stress.
Overall, the case study of the Saravan watershed has provided a contribution in understanding the linkages between the rill detachment capacity on one side, and the soil properties and hydraulic parameters on the other side under different land uses. It has been finally demonstrated that forestland and woodland can be suitable land uses (mainly in soils with prevalent clayey and silty content) when the soil conservation purposes are achieved by land planners and watershed managers. Therefore, afforestation in the croplands should be prioritized when possible, in order to reverse or mitigate possible erosion-degradation.