Evaluating the Effects of the Rill Longitudinal Proﬁle on Flow Resistance Law

: In this paper, for the ﬁrst time, the effect of the longitudinal proﬁle shape of the rill (uniform, concave, and convex) on ﬂow resistance law was studied. The ﬁrst part of the paper is based on a theoretical equation to estimate the Darcy–Weisbach friction factor f , deduced from the power velocity distribution and rill measurements performed on a plot. At ﬁrst, the equation to estimate the Γ parameter of the velocity proﬁle was calibrated using all available measurements. Then an analysis of the hydraulic characteristics at reach scale, for comparable values of discharge, was carried out, comparing the different proﬁle shapes. To assess the inﬂuence of the rill proﬁle shape on ﬂow resistance law, this calibration was also carried out using the data categorized by proﬁle. In the second part of the paper, an analysis of the scour depth and eroded rill volume was developed for four rills of each conﬁguration. The results showed that an accurate estimate of f can be obtained by calibrating the ﬂow resistance equation for each proﬁle shape. The component of the Darcy–Weisbach friction factor that is due to the proﬁle shape varied from 0.68 to 14.6% of the overall friction factor for the concave proﬁle, and from 3.4 to 26.9% for the convex proﬁle. The analysis also showed that, for the convex proﬁle, the scour was concentrated downstream of the slope change, while for the uniform and concave proﬁles, it was uniformly distributed. Furthermore, the scour depth measured in correspondence of the rill thalweg had a generally increasing trend, with discharge for all the investigated proﬁle shapes. Total eroded rill volume of the concave proﬁle was lower than those detected for the uniform and convex proﬁles and was characterized by a reduction of 57.9%, as compared to the uniform proﬁle.


Introduction
Many researchers have studied the effects of length and slope gradient on runoff and soil loss using uniform plots, i.e., having a constant slope, while a small number of studies have focused on the effect of profile shape (concave, convex) [1,2].
Young and Mutchler [3,4] investigated experimentally the effect of slope shape (convex, concave, uniform) on runoff and soil loss at plot scale and showed that concave hillslopes tends to have reduced total sediment loss compared to a uniform profile slope.
Rieke-Zapp and Nearing [2] carried out laboratory experiments to determine the relationship between slope shapes and soil erosion using five slope shape treatments (convex-linear, concave-linear, nose slope, head slope and uniform). These experiments demonstrated that rill patterns, sediment yield and runoff, all changed with slope shape. In particular, these authors reported a decrease of 75% of soil loss for concave slopes, as compared to uniform ones with similar surface area. Sensoy and Kara [5] carried out an experimental investigation using nine experimental field plots with different profile shapes, established on a 30% hillslope, subjected to natural rainfall. This field experiment indicated that slope shape (uniform, concave, convex) metal cylinder of known volume, and determining its weight after drying. A uniform soil moisture condition for all the experiments was achieved by wetting the soil until infiltration stopped. Then, rills were manually incised along the maximum slope direction of the plot and were shaped using a clear flow discharge of 0.1 L s −1 , applied for 3 min. The rills had the same longitudinal profile as the plot. For each rill, nine longitudinal segments, bounded by two cross-sections spaced 0.624 m apart (Figure 3), were determined.   The plot was prepared with three different profiles (concave, convex and uniform) ( Figure 1) using a mean plot slope s p of 18%. In particular, for the concave and convex profiles, the plot was divided into two 3.27 m long parts, and two slope values (12 and 24%) were set up to obtain a mean value of 18% ( Figure 2). The soil was a clay loam characterized by 32.7%, 30.9% and 36.4% of clay, silt, and sand, respectively. metal cylinder of known volume, and determining its weight after drying. A uniform soil moisture condition for all the experiments was achieved by wetting the soil until infiltration stopped. Then, rills were manually incised along the maximum slope direction of the plot and were shaped using a clear flow discharge of 0.1 L s −1 , applied for 3 min. The rills had the same longitudinal profile as the plot. For each rill, nine longitudinal segments, bounded by two cross-sections spaced 0.624 m apart (Figure 3), were determined.   For all the profile shapes, soil was manually hoed and raked and then compacted by a woody bar. The bulk density of the soil was determined by sampling the soil, using a metal cylinder of known volume, and determining its weight after drying. A uniform soil moisture condition for all the experiments was achieved by wetting the soil until infiltration stopped. Then, rills were manually incised along the maximum slope direction of the plot and were shaped using a clear flow discharge of 0.1 L s −1 , applied for 3 min. The rills had the same longitudinal profile as the plot. For each rill, nine longitudinal segments, bounded by two cross-sections spaced 0.624 m apart (Figure 3), were determined. A constant clear inflow discharge, ranging from 0.18 to 0.63 L s −1 , was used for the experiments.
Measurements of water depth h, cross-section area σ, wetted perimeter C, bed slope s, and flow velocity V were carried out for rill reaches (i.e., the channel between a given measured cross-section and the rill end) on uniform (99 reaches), concave (71 reaches), and convex (72 reaches) rills. The indirect measurements of Reynolds number Re = V h/vk, in which νk is the kinematic viscosity, Froude number F = V/√(g h), in which g is the acceleration due to gravity, and Darcy-Weisbach friction factor f = 8 gRs /V 2 , where R is the hydraulic radius, were also obtained.
For each experimental run, a 3D-DTM of the plot area was built from a set of 70 photographs using a 3D-photo reconstruction technique with the image-processing software Agisoft Photoscan Professional (Version 1.1.6, Agisoft, St. Petersburg, Russia). The slope gradient of each rill segment was determined using the rill thalweg imaged by the 3D-DTM. The slope of each rill reach was calculated by averaging the values measured in the considered reach.
The mean water depth h and the mean hydraulic radius R in the reach were determined using the method developed by Di Stefano et al. [36]. This method couples an accurate ground survey of the rill channel, obtained by close-range photogrammetry, with the survey of the water tracks inside the channel marked by a dye solution. The measured water depths and geometric cross-sections obtained by the 3D-DTM, at a distance interval of 6.2 cm within the reach, were used to calculate the values of the hydraulic cross-section area σ and wetted perimeter C in the reach and the values of h and R using the following equations [36]: A constant clear inflow discharge, ranging from 0.18 to 0.63 L s −1 , was used for the experiments.
Measurements of water depth h, cross-section area σ, wetted perimeter C, bed slope s, and flow velocity V were carried out for rill reaches (i.e., the channel between a given measured cross-section and the rill end) on uniform (99 reaches), concave (71 reaches), and convex (72 reaches) rills. The indirect measurements of Reynolds number Re = V h/ν k , in which ν k is the kinematic viscosity, Froude number F = V/ √ (g h), in which g is the acceleration due to gravity, and Darcy-Weisbach friction factor f = 8 gRs /V 2 , where R is the hydraulic radius, were also obtained.
For each experimental run, a 3D-DTM of the plot area was built from a set of 70 photographs using a 3D-photo reconstruction technique with the image-processing software Agisoft Photoscan Professional (Version 1.1.6, Agisoft, St. Petersburg, Russia). The slope gradient of each rill segment was determined using the rill thalweg imaged by the 3D-DTM. The slope of each rill reach was calculated by averaging the values measured in the considered reach.
The mean water depth h and the mean hydraulic radius R in the reach were determined using the method developed by Di Stefano et al. [36]. This method couples an accurate ground survey of the rill channel, obtained by close-range photogrammetry, with the survey of the water tracks inside the channel marked by a dye solution. The measured water depths and geometric cross-sections obtained by the 3D-DTM, at a distance interval of 6.2 cm within the reach, were used to calculate the values of the hydraulic cross-section area σ and wetted perimeter C in the reach and the values of h and R using the following equations [36]: in which w is the surface width and N is the number of the cross-sections in the reach. A methylene blue solution was applied as a dye tracer to measure the flow velocity [28,[37][38][39][40]. The measured surface velocity was converted to the mean flow velocity V by a correction factor of 0.8 [41][42][43].
For four rills of each profile shape, Digital Elevation Models (DEMs) at the ends of the shaping phase (D1) and the experimental runs (D2) were obtained from the 3D-DTM. The DEM of the difference (DoD) between D2 and D1 was used to measure the scour depth values inside the whole rill channel and calculate the total eroded rill volume RV. Considering that the minimum level of detection of the used DEMs was ±0.005 m, variations of bed elevation in the range ±0.005 m were not detectable.
Finally, for the same rills, the scour depth values SD were measured, at 0.005 m intervals along the rill thalweg, as the difference in height between the longitudinal profiles at the end of the shaping phase (D1) and at the end of the experimental runs (D2) (an example for one rill with a convex profile shape is plotted in Figure 4). For each rill, the mean scour depth SD m was also calculated as the mean of the SD values.
in which w is the surface width and N is the number of the cross-sections in the reach. A methylene blue solution was applied as a dye tracer to measure the flow velocity [28,[37][38][39][40]. The measured surface velocity was converted to the mean flow velocity V by a correction factor of 0.8 [41][42][43].
For four rills of each profile shape, Digital Elevation Models (DEMs) at the ends of the shaping phase (D1) and the experimental runs (D2) were obtained from the 3D-DTM. The DEM of the difference (DoD) between D2 and D1 was used to measure the scour depth values inside the whole rill channel and calculate the total eroded rill volume RV. Considering that the minimum level of detection of the used DEMs was ±0.005 m, variations of bed elevation in the range ±0.005 m were not detectable.
Finally, for the same rills, the scour depth values SD were measured, at 0.005 m intervals along the rill thalweg, as the difference in height between the longitudinal profiles at the end of the shaping phase (D1) and at the end of the experimental runs (D2) (an example for one rill with a convex profile shape is plotted in Figure 4). For each rill, the mean scour depth SDm was also calculated as the mean of the SD values.

The Rill Flow Resistance Equation
For an open channel flow, the local flow velocity profile v(y) along a given vertical is represented by the following functional relationship [44][45][46]: , , , , * , , , , = where ϕ is a functional symbol, v is the local velocity, y is the distance from the bottom, d is a characteristic bed particle diameter, u*= √(g R s) is the shear velocity, ρ is the water density and μ is the dynamic water viscosity.

The Rill Flow Resistance Equation
For an open channel flow, the local flow velocity profile v(y) along a given vertical is represented by the following functional relationship [44][45][46]: where ϕ is a functional symbol, v is the local velocity, y is the distance from the bottom, d is a characteristic bed particle diameter, u * = √ (g R s) is the shear velocity, ρ is the water density and µ is the dynamic water viscosity. Applying the Π-Theorem of the dimensional analysis and using as dimensional independent variables y, u * and µ, Equation (3) can be expressed by the following dimensionless functional relationship [30]: in which ϕ 1 is a functional symbol.
Hypothesizing the Incomplete Self-Similarity in u * y/ν k [25,[47][48][49][50], neglecting the influence of Re (turbulent flow regime) and considering that the flow Froude number accounts also for the ratio h/d [26], the integration of Equation (4) yields to the following velocity profile: where Γ(s, F) is a function to be defined by velocity measurements and the exponent δ is calculated by the following theoretical equation [51,52] The following expression of the Darcy-Weisbach friction factor f is deduced [25,48,53] by integrating Equation (5): From the velocity profile (Equation (5)) the following estimate Γ v of Γ function [25,53] is obtained by setting equal to α h the distance y at which the local velocity is equal to the cross-section average velocity V: The coefficient α, which is less than 1, takes into account that both V is located below the water surface and the mean velocity profile in the cross-section is considered. The coefficient α has to be calculated by the following theoretical relationship deduced by Ferro [25]: For this investigation α varied in a narrow range and the mean value of 0.124 was assumed. Considering that, according to Equation (5), Γ theoretically depends only on s and F [26,54], Γ v can be estimated using the following power equation: where a, b and c are coefficients to be determined from experimental measurements.

Rill Flow Resistance for Different Profile Shapes
The experimental runs were carried out for Reynolds numbers corresponding to turbulent flow (2764 ≤ Re ≤ 11015 for the uniform profile, 2282 ≤ Re ≤ 14346 for the concave one, and 2486 ≤ Re ≤ 9487 for the convex one) and for Froude number values corresponding to subcritical and supercritical flow conditions (0.54 ≤ F ≤ 1.70 for the uniform profile, 0.64 ≤ F ≤ 1.63 for the concave one, and 0.44 ≤ F ≤ 1.71 for the convex one) ( Figure 5). In other words, the investigated ranges of Re and F are very similar for the three different profiles.

profiles.
Using all the 242 measurements, i.e., neglecting the influence of the profile shape, Equation (10) was calibrated obtaining the following relationship: characterized by a coefficient of determination equal to 0.9894. The comparison between the 242 measured Γv values, obtained by Equation (8) and those calculated applying Equation (11) is plotted in Figure 6a. Introducing Equation (11) into Equation (7), the following flow resistance law was obtained: A good agreement, characterized by a root mean square error RMSE equal to 0.1743, between the measured Darcy-Weisbach friction factor values and those calculated by Equation (12) was detected (Figure 6b), even if Equation (12) tended to slightly overestimate the f values for the uniform profile. The friction factor values calculated by Equation (12) were characterized by errors less than or equal to ±20% for 96.3% of cases and less than or equal to ±10% for 78.5% of cases. Figure 6 highlights that the f measurements obtained for the convex profile covered a different range than those obtained from the uniform and concave profiles. Different characteristics of the convex profile are also revealed by the frequency distributions of the hydraulic variables F, f, s and V plotted in Figure 7. Using all the 242 measurements, i.e., neglecting the influence of the profile shape, Equation (10) was calibrated obtaining the following relationship: characterized by a coefficient of determination equal to 0.9894. The comparison between the 242 measured Γ v values, obtained by Equation (8) and those calculated applying Equation (11) is plotted in Figure 6a. Introducing Equation (11) into Equation (7), the following flow resistance law was obtained: A good agreement, characterized by a root mean square error RMSE equal to 0.1743, between the measured Darcy-Weisbach friction factor values and those calculated by Equation (12) was detected (Figure 6b), even if Equation (12) tended to slightly overestimate the f values for the uniform profile. The friction factor values calculated by Equation (12) were characterized by errors less than or equal to ±20% for 96.3% of cases and less than or equal to ±10% for 78.5% of cases. Figure 6 highlights that the f measurements obtained for the convex profile covered a different range than those obtained from the uniform and concave profiles. Different characteristics of the convex profile are also revealed by the frequency distributions of the hydraulic variables F, f, s and V plotted in Figure 7.  As this analysis showed that profile shape influences flow characteristics, the theoretical flow resistance law was calibrated for each investigated condition. Table 1 lists the a, b and c coefficients of Equation (10) obtained for the uniform, concave, and convex profiles.  As this analysis showed that profile shape influences flow characteristics, the theoretical flow resistance law was calibrated for each investigated condition. Table 1 lists the a, b and c coefficients of Equation (10) obtained for the uniform, concave, and convex profiles.  Figure 8a shows the comparison between the measured Γv values obtained by Equation (8) and those calculated applying Equation (10), with the specific a, b and c coefficients (Table 1) for the different profiles. The coefficient of determination of the obtained relationships is equal to 0.9812, 0.9979 and 0.9882 for the uniform, concave, and convex profiles, respectively. Figure 8b shows good agreement, characterized by a RMSE equal to 0.1234, between the measured values of the Darcy-Weisbach friction factor and those calculated by Equations (7) and (10), with a, b and c coefficients obtained for each profile shape. In this case, the calculated friction factor values are characterized by errors always less than or equal to ±20% and less than or equal to ±10% for 86.8% of cases.  Figure 8a shows the comparison between the measured Γ v values obtained by Equation (8) and those calculated applying Equation (10), with the specific a, b and c coefficients (Table 1) for the different profiles. The coefficient of determination of the obtained relationships is equal to 0.9812, 0.9979 and 0.9882 for the uniform, concave, and convex profiles, respectively. Figure 8b shows good agreement, characterized by a RMSE equal to 0.1234, between the measured values of the Darcy-Weisbach friction factor and those calculated by Equations (7) and (10), with a, b and c coefficients obtained for each profile shape. In this case, the calculated friction factor values are characterized by errors always less than or equal to ±20% and less than or equal to ±10% for 86.8% of cases. file on flow resistance, was calculated as the difference f-fu between the measured f value of the concave or convex profile and that of, fu, calculated for the uniform profile in the same hydraulic conditions (Equations (7) and (10), with a, b and c coefficients obtained for the uniform profile). This fmor value is almost equal to zero for 9.9% (concave) and 4.2% (convex) of the measured values. The ratio fmor/f varies from 0.68 to 14.6% for concave profile and from 3.4 to 26.9% for convex profile (Figure 9). (8) and those calculated by Equation (10) with the shape profile-specific a, b and c coefficients (a) and measured Darcy-Weisbach friction factor values and those calculated by Equations (7) and (10) with the shape profilespecific a, b and c coefficients (b).  (8) and those calculated by Equation (10) with the shape profile-specific a, b and c coefficients (a) and measured Darcy-Weisbach friction factor values and those calculated by Equations (7) and (10) with the shape profilespecific a, b and c coefficients (b).

Figure 8. Comparison between the measured Γv values obtained by Equation
Moreover, to study the influence of profile shape on rill flow resistance, the additive method of shear stresses [55], already used to investigate the effects of sediment transport for rills [29,56], and of riparian vegetation for a vegetated channel [57], was applied. In particular, the Darcy-Weisbach friction factor, f mor , corresponding to the effect of the profile on flow resistance, was calculated as the difference f -f u between the measured f value of the concave or convex profile and that of, f u , calculated for the uniform profile in the same hydraulic conditions (Equations (7) and (10), with a, b and c coefficients obtained for the uniform profile). This f mor value is almost equal to zero for 9.9% (concave) and 4.2% (convex) of the measured values. The ratio f mor /f varies from 0.68 to 14.6% for concave profile and from 3.4 to 26.9% for convex profile (Figure 9). Water 2022, 14, x FOR PEER REVIEW 11 of 19 Figure 9. Frequency distribution of the ratio fmor/f. For all three investigated profile shapes, Figure 10 also shows that rills are generally scoured by flow. The blue areas of this figure identify rill zones where the variation of bed elevation is lower than the minimum level of detection (±0.005 m), or lack variation. Moreover, for the convex profile, the scour is concentrated downstream of the slope change, while for uniform and concave profiles, it is uniformly distributed.  For all three investigated profile shapes, Figure 10 also shows that rills are generally scoured by flow. The blue areas of this figure identify rill zones where the variation of bed elevation is lower than the minimum level of detection (±0.005 m), or lack variation. Moreover, for the convex profile, the scour is concentrated downstream of the slope change, while for uniform and concave profiles, it is uniformly distributed.

Analysis of Rill Scour Depth and Eroded Volume for Different Profile Shapes
The knowledge of the scour depth values obtained by the DoDs allows calculation of the total eroded volume RV. Figure 11 shows the comparison between the applied discharge Q and RV using an example of four rills analysed for each profile. The evident trend, as expected, is that RV increases with Q for all the investigated profiles. Finally, for the uniform, concave, and convex profiles, the sum of the four RV values resulted in 0.0644, 0.0271 and 0.0654 m 3 , respectively.
For the forementioned four experimental runs for each profile shape, the frequency distribution of the scour depth values SD was plotted in Figure 12. This figure highlights that scour depth, measured in corresponding rill thalwegs, has a generally increasing trend with discharge for all investigated profile shapes. Water 2022, 14, x FOR PEER REVIEW 12 of 19 The knowledge of the scour depth values obtained by the DoDs allows calculation of the total eroded volume RV. Figure 11 shows the comparison between the applied discharge Q and RV using an example of four rills analysed for each profile. The evident trend, as expected, is that RV increases with Q for all the investigated profiles. Finally, for the uniform, concave, and convex profiles, the sum of the four RV values resulted in 0.0644, 0.0271 and 0.0654 m 3 , respectively.
For the forementioned four experimental runs for each profile shape, the frequency distribution of the scour depth values SD was plotted in Figure 12. This figure highlights that scour depth, measured in corresponding rill thalwegs, has a generally increasing trend with discharge for all investigated profile shapes.  Figure 13 shows the frequency distribution of the SD/SD m ratio for the investigated discharges. This figure demonstrates that the variable SD/SD m can be considered independent of discharge. Figure 14 shows the empirical frequency distribution of SD/SD m for each profile shape, using an example for a single discharge value (≈0.35 L s −1 ). The overlapping of the three distributions suggested that SD/SD m does not depend on profile shape. Water 2022, 14, x FOR PEER REVIEW 13 of 19 Figure 11. An example of the comparison between the applied discharge Q and the total eroded rill volume RV, for four rill channels for each profile shape.  Figure 13 shows the frequency distribution of the SD/SDm ratio for the investigated discharges. This figure demonstrates that the variable SD/SDm can be considered independent of discharge. Figure 11. An example of the comparison between the applied discharge Q and the total eroded rill volume RV, for four rill channels for each profile shape.
Water 2022, 14, x FOR PEER REVIEW 13 of 19 Figure 11. An example of the comparison between the applied discharge Q and the total eroded rill volume RV, for four rill channels for each profile shape.  Figure 13 shows the frequency distribution of the SD/SDm ratio for the investigated discharges. This figure demonstrates that the variable SD/SDm can be considered independent of discharge.   Figure 14 shows the empirical frequency distribution of SD/SDm for each profile shape, using an example for a single discharge value (≈0.35 L s −1 ). The overlapping of the three distributions suggested that SD/SDm does not depend on profile shape.    Figure 14 shows the empirical frequency distribution of SD/SDm for each profile shape, using an example for a single discharge value (≈0.35 L s −1 ). The overlapping of the three distributions suggested that SD/SDm does not depend on profile shape.  The relationship between SD m and the discharge was determined for each profile shape ( Figure 15). The mean scour depth increased with discharge according to a power relationship depending on the profile. The relationship between SDm and the discharge was determined for each profile shape ( Figure 15). The mean scour depth increased with discharge according to a power relationship depending on the profile.  Figure 5 shows that the rill flows for the three examined profiles were characterized by comparable Reynolds numbers, while for the convex profile, they were characterized, on average, by Froude numbers lower than those related to the uniform and concave profiles, and less than 1 for 93% of cases. For the latter, the pairs (Re, F) were generally overlapped. The detected trend for Froude number confirmed that the profile shape influences the kinematic flow characteristics (Figure 7a), while the comparable Reynolds number values ( Figure 5) suggest that the profile does not affect the flow regime. Figure 7 confirmed this overlapping trend for the measured flow velocity (Figure 7d) and Froude number (Figure 7a) on uniform and concave profiles. The slope values of the concave profile were less than those of the uniform profile, while the s values of the convex profile were greater than those measured on the uniform profile (Figure 7c). For the convex profile, the occurrence of the highest slope values and the lowest velocity values can be explained by the relevant erosion phenomena. Indeed, the highest RV values detected for the convex profile at the single rill and plot scale, demonstrate that for this profile, the erosion phenomena determine a major increase of grain roughness and sediment transport that lower flow velocity and increase the Darcy-Weisbach friction factor ( Figure  7b). The overlapping values of flow velocity and friction factor f for uniform and concave profiles can be explained by the similar distribution of erosion phenomena occurring in the rill channels (Figure 10a,b).

Rill Flow Resistance for Different Profile Shapes
The comparison between Figure 6b and 8b and the error statistics highlight that a single flow resistance equation (Equation (12)) is applicable to all investigated profiles, despite the differences highlighted by the analysis of hydraulic characteristics at reach scale. However, the f estimate is more accurate if specific calibration coefficients (a, b and c) are applied for each profile. Figure 9 demonstrates that the ratio fmor/f was greater than 10% for 39.1% and 69.6% of the investigated cases for the concave and convex profiles, respectively. Thus, for the examined conditions, the effect of profile shape on flow resistance can be considered relevant.  Figure 5 shows that the rill flows for the three examined profiles were characterized by comparable Reynolds numbers, while for the convex profile, they were characterized, on average, by Froude numbers lower than those related to the uniform and concave profiles, and less than 1 for 93% of cases. For the latter, the pairs (Re, F) were generally overlapped. The detected trend for Froude number confirmed that the profile shape influences the kinematic flow characteristics (Figure 7a), while the comparable Reynolds number values ( Figure 5) suggest that the profile does not affect the flow regime. Figure 7 confirmed this overlapping trend for the measured flow velocity (Figure 7d) and Froude number (Figure 7a) on uniform and concave profiles. The slope values of the concave profile were less than those of the uniform profile, while the s values of the convex profile were greater than those measured on the uniform profile (Figure 7c). For the convex profile, the occurrence of the highest slope values and the lowest velocity values can be explained by the relevant erosion phenomena. Indeed, the highest RV values detected for the convex profile at the single rill and plot scale, demonstrate that for this profile, the erosion phenomena determine a major increase of grain roughness and sediment transport that lower flow velocity and increase the Darcy-Weisbach friction factor (Figure 7b). The overlapping values of flow velocity and friction factor f for uniform and concave profiles can be explained by the similar distribution of erosion phenomena occurring in the rill channels (Figure 10a,b).

Rill Flow Resistance for Different Profile Shapes
The comparison between Figures 6b and 8b and the error statistics highlight that a single flow resistance equation (Equation (12)) is applicable to all investigated profiles, despite the differences highlighted by the analysis of hydraulic characteristics at reach scale. However, the f estimate is more accurate if specific calibration coefficients (a, b and c) are applied for each profile. Figure 9 demonstrates that the ratio f mor /f was greater than 10% for 39.1% and 69.6% of the investigated cases for the concave and convex profiles, respectively. Thus, for the examined conditions, the effect of profile shape on flow resistance can be considered relevant.
The convex profile was characterized by the highest values of the ratio f mor /f (3.4-26.9%) (Figure 9), which demonstrates that the contribution of the component of total flow resistance due to the profile was more pronounced than for the concave profile. This result can be explained as higher eroded volumes were detected for the convex profile and f mor accounts for sediment transport and roughness effects on the flow resistance.

Analysis of Rill Scour Depth and Eroded Volume for Different Profile Shapes
Previous studies [23,58] considered a slope value of 18% as a threshold to discriminate between the gentle (GS) and the steep slope (SS) condition. This distinction is explained as differences in hydraulic (flow depth, velocity, Reynolds number, Froude number) and sediment transport (flow transport capacity, actual sediment load) variables occur. For the SS condition, results by Nicosia et al. [58] supported the conclusion that the sediment transport capacity T c was so high that it was not a limiting factor for the actual sediment transport, while for the GS case, the actual sediment transport could be limited by both T c and the soil particle detachability and transportability.
For the convex profile, Figure 10c shows that the scour was concentrated, while it was evenly distributed for uniform ( Figure 10a) and concave profiles (Figure 10b). In particular, the noticeable scour of the rills in convex profiles was located downstream of the slope change from a gentle (12%) to a steep (24%) slope condition, because in the upstream half of the rill, the actual sediment transport was limited by T c , while in the downstream half, it was not. Notwithstanding, the spatial distribution of the scour was similar for uniform and concave profiles, the scour depth values of concave profiles were lower than those of uniform and convex profiles, which instead were similar ( Figure 10). Consequently, also the eroded volume of the rill for the concave profile was lower than the eroded volumes for the uniform and convex profiles ( Figure 11). In particular, the concave profile resulted in 57.9% less soil loss than the uniform profile. This result confirms the literature findings [2,[6][7][8], which have reported 50 to 80% less soil loss from concave hillslopes than from uniform slopes. Figure 11 also shows that increasing values of flow discharge determine a higher RV. This result agrees with that obtained by Di Stefano et al. [56] who found that when discharge increased, the contribution of sediment transport to flow resistance and the eroded rill volume also increased.
The results shown in Figures 13 and 14 suggest that SD/SD m was independent of discharge and profile shape, respectively. In other words, the mean value SD m accounted for the effects of both discharge and profile shape on scour depth.
Finally, Figure 15 demonstrates that profile shape determined a different increasing relationship between mean scour depth along the rill thalweg and discharge. The mean scour depth for the uniform profile was always higher than that obtained for the convex profile. Instead, the total eroded volume obtained for the uniform profile was slightly lower than that obtained for the convex profile. This result can be explained as the DoD information accounts for the spatial distribution of scour depth, while SD m represents the mean value of measurements performed along a single linear flow path.

Conclusions
In this paper, the effect of the rill profile on flow resistance law was evaluated using measurements obtained from uniform, concave, and convex profiles.
At first, the equation to estimate the velocity profile parameter Γ was calibrated by all the available measurements and an accurate estimate of the Darcy-Weisbach friction factor (Equation (12)) was obtained. Then, to assess the influence of each profile on flow resistance law, the equation to estimate Γ was also calibrated when the data for each profile shape was distinguished. This choice determined an appreciable improvement of the friction factor estimate.
The Darcy-Weisbach friction factor corresponding to the effect of the profile on flow resistance, varied from 0.68 to 14.6% of the overall f factor for the concave profile, and from 3.4 to 26.9% of f for the convex profile.
The developed analysis also showed that for convex profiles, scour was concentrated downstream of the slope change, while for uniform and concave profiles, it was uniformly distributed. Furthermore, the scour depth measured in the corresponding rill thalwegs, had a generally increasing trend with discharge for all investigated profiles.
The total eroded volume for the concave profile was lower than those obtained for uniform and convex profiles and was characterized by being 57.9% less than for the uniform profile.
Finally, the analysis indicated that the mean value of scour depth in the rill thalweg accounted for the effects of both discharge and profile shape, and increases with increased discharge according to a relationship varying with the profile shape.
The main limitations of the study are the use of a single soil type and one mean plot slope. These can be overcome by other experiments using different soils and values of mean plot slope.
The present investigation highlights a markedly less soil loss from concave slopes than from either uniform or convex slopes. Moreover, the proposed approach allows estimating rill flow velocity, which is an essential component of process-based soil erosion models.