Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Biodegradation of Cyanide Using Soda Lake-Derived Alkaliphilic Microbial Consortia
Previous Article in Journal
The Impact of Cobalt Species on the Hazardous Characteristics of Cobalt-Leaching Residue: A Case Study from Guangdong Province, China
Previous Article in Special Issue
Characteristics of Soil Physical Properties and Spatial Distribution of Soil Erosion on Ridge-Slope Farmland in the Black Soil Areas of Northeast China
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 16
Issue 20
10.3390/w16202955
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

One-Dimensional Numerical Cascade Model of Runoff and Soil Loss on Convergent and Divergent Plane Soil Surfaces: Laboratory Assessment and Numerical Simulations

by
Babar Mujtaba
1,
João L. M. P. de Lima
2,3,* and
M. Isabel P. de Lima
2,3
1
Research Institute on Mines and Environment (RIME), Université du Québec en Abitibi-Témiscamingue (UQAT), 445 Blvd. De l’Université, Rouyn-Noranda, QC J9X 5E4, Canada
2
Department of Civil Engineering, Faculty of Sciences and Technology, University of Coimbra, Rua Luís Reis Santos, 3030-788 Coimbra, Portugal
3
MARE—Marine and Environmental Sciences Centre/ARNET—Aquatic Research Network, University of Coimbra, 3030-790 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Water 2024, 16(20), 2955; https://doi.org/10.3390/w16202955
Submission received: 1 September 2024 / Revised: 11 October 2024 / Accepted: 14 October 2024 / Published: 17 October 2024
(This article belongs to the Special Issue Mechanism and Simulation of Water Erosion and Nutrient Loss on Hillslopes)
Browse Figures
Figure 1
<p>Sketch and notation used for the cascade of <span class="html-italic">n</span>-planes representing converging (sloping to the left) or diverging (sloping to the right) cascades.</p> "> Figure 2
<p>(<b>a</b>) Laboratory setup used in the experiments, consisting of a rainfall simulator and a soil flume (top and bottom left); the flume has outlets on both ends for the converging plane (bottom middle) and diverging plane surface (bottom right) experiments. (<b>b</b>) Dashed lines represent contour lines relative to an arbitrary datum (ground level) for a 20% flume slope, while solid lines indicate the border wall of the flume’s geometry.</p> "> Figure 3
<p>On the <b>right</b>, a schematic sketch illustrates the approximate representation of the soil flume surface planar geometry using converging and diverging cascade planes, each with an equal length <math display="inline"><semantics> <mrow> <mi>L</mi> </mrow> </semantics></math> of 0.5 m, measured along the direction of the plane slope. The width <math display="inline"><semantics> <mrow> <mi>W</mi> </mrow> </semantics></math> of each plane is also shown, which is measured in the perpendicular direction. On the <b>left</b>, the variation in mean rainfall intensity across the planes of the converging and diverging cascades is shown.</p> "> Figure 4
<p>Seven plane soil surfaces: one rectangular, three converging with different convergence angles θ, and three diverging with different divergence angles θ. The schematic sketch also illustrates the approximate representation of the soil surface planar geometry using converging and diverging cascade planes of equal length <span class="html-italic">L</span> (1 m). The width <span class="html-italic">W</span> of each plane is also shown. The dimensions of the rectangular plane soil surface (4 m × 2 m) are provided as well.</p> "> Figure 5
<p>Observed and numerically simulated hydrographs for converging and diverging surfaces with different combinations of rainfall intensity (<span class="html-italic">I</span>) and slope (<span class="html-italic">S</span>). The Nash–Sutcliffe (<math display="inline"><semantics> <mrow> <mi>N</mi> <mi>S</mi> </mrow> </semantics></math>) coefficients are also shown.</p> "> Figure 6
<p>Observed and numerically simulated sediment graphs for converging and diverging surfaces with different combinations of rainfall intensity (<span class="html-italic">I</span>) and slope (<span class="html-italic">S</span>). Note that the Y-axis scale is not the same on the top-left figure. The Nash–Sutcliffe (<math display="inline"><semantics> <mrow> <mi>N</mi> <mi>S</mi> </mrow> </semantics></math>) coefficients are also shown.</p> "> Figure 7
<p>Simulated hydrographs for plane soil surfaces with different angles (θ) of convergence (<b>left</b>) and divergence (<b>right</b>). Section A represents a rectangular surface where θ = 0°.</p> "> Figure 8
<p>Runoff peaks for plane soil surfaces with different angles (θ) of convergence (<b>left</b>) and divergence (<b>right</b>). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference.</p> "> Figure 9
<p>Simulated sediment graphs for plane soil surfaces with different angles (θ) of convergence (<b>left</b>) and divergence (<b>right</b>). Section A represents a rectangular surface where θ = 0°. Note that the Y-axis scales are not the same.</p> "> Figure 10
<p>Peak soil loss for plane soil surfaces with different angles (θ) of convergence (<b>left</b>) and divergence (<b>right</b>). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference. Note that the Y-axis scales differ.</p> ">
Versions Notes

Abstract

:
A one-dimensional numerical overland flow model based on the cascade plane theory was developed to estimate rainfall-induced runoff and soil erosion on converging and diverging plane surfaces. The model includes three components: (i) soil infiltration using Horton’s infiltration equation, (ii) overland flow using the kinematic wave approximation of the one-dimensional Saint-Venant shallow water equations for a cascade of planes, and (iii) soil erosion based on the sediment transport continuity equation. The model’s performance was evaluated by comparing numerical results with laboratory data from experiments using a rainfall simulator and a soil flume. Four independent experiments were conducted on converging and diverging surfaces under varying slope and rainfall conditions. Overall, the numerically simulated hydrographs and sediment graphs closely matched the laboratory results, showing the efficiency of the model for the tested controlled laboratory conditions. The model was then used to numerically explore the impact of different plane soil surface geometries on runoff and soil loss. Seven geometries were studied: one rectangular, three diverging, and three converging. A constant soil surface area, the rainfall intensity, and the slope gradient were maintained in all simulations. Results showed that increasing convergence angles led to a higher peak and total soil loss, while decreasing divergence angles reduced them.

1. Introduction

Rainfall-induced soil erosion poses a serious environmental threat, leading to the reduced fertility of arable lands, lower crop yields, intensified surface runoff, and the transportation of harmful pollutants to freshwater bodies. Consequently, soil conservation is crucial for sustainable agriculture, necessitating an accurate estimation of soil loss and the identification of critical areas for implementing soil conservation measures, especially on complex soil surface geometries.
Mathematical models offer cost-effective tools to understand rainfall-induced soil erosion processes, evaluate soil and water conservation practices under different environmental conditions, and assist in policy development for future climate scenarios. However, no universal model can accurately predict soil erosion across all temporal and spatial scales.
Characterizing the rainfall-induced runoff response is essential for estimating soil erosion. The kinematic wave (KW) theory is widely used for modeling overland flow, with one-dimensional (1D) KW approximations proving effective on small slopes [1,2].
Sediment transport data, which are crucial for testing soil erosion models, can be obtained from field plots and laboratory soil flumes, taking into account factors such as rainfall characteristics, slope geometry, soil properties, and land cover. While numerous studies have evaluated, under controlled laboratory conditions, physically based soil erosion and transport models on simple geometries such as square or rectangular plane soil surface geometries [3,4,5,6], complex soil surfaces, which are common in nature, are less studied.
Only a few researchers [7,8,9,10,11] have examined erosion on curved rectangular surfaces (convex, plane, and concave). These types of geometries are similar to those often found in rocky catchments where landslides are common [12,13]. To the best of our knowledge, to date, no experimental or numerical studies have assessed the impact of varying converging and diverging angles of plane soil surfaces on rainfall-induced runoff and soil loss. Additionally, no physically based 1D model exists for estimating erosion on these surfaces.
This paper has two objectives: first, to develop a 1D cascade model for estimating rainfall-induced runoff and soil erosion on converging and diverging plane surfaces, and to test it using controlled laboratory experiments. The physically based mathematical model, grounded in mass conservation principles, simulates sediment detachment, transport, and deposition processes [14,15]. Laboratory experiments were conducted using a rainfall simulator and a soil flume with a planar circular sector geometry to create converging and diverging flow conditions. Second, this paper evaluates the influence of different plane surfaces’ converging and diverging angles on runoff and soil loss through numerical simulations on seven plane surfaces, including one rectangular/linear configuration (angle along the flow direction θ = 0°), and surfaces with three unique converge and divergence angles (θ = 8.5°, 16.7°, and 24.2°). The soil surface area, rainfall intensity, and slope gradient were kept constant across all numerical simulations.

2. Proposed Cascade Model

Soil erosion is a complex phenomenon influenced by the geometry and geomorphology of soil surfaces. The 1D cascade model, used in this work to estimate rainfall-induced runoff and soil erosion, represents the soil surface geometry as a series of cascading rectangular planes. The model consists of three components: (i) Horton’s equation for soil infiltration; (ii) the KW approximation of the 1D Saint-Venant shallow water equations for overland flow generation; and (iii) the 1D sediment transport equation for soil loss calculation.

2.1. Infiltration Model

Horton [16,17,18] introduced an exponential decay equation to describe how the infiltration rate varies over time after a soil surface is exposed to constant rainfall:
f c t = f 1 + f 0 − f 1 e − k t
where f c ( t ) is the infiltration rate at any time t, f 0   is the initial infiltration rate when rainfall starts (at t = 0), f 1 is the final or equilibrium infiltration rate, and k is the exponential time decay coefficient. Cumulative infiltration F , at any time, t, is:
F = ∫ 0 t f c t d t = f 1 t + f 0 − f 1 k   1 − e − k t
Using Equation (1) to eliminate t in Equation (2), the resultant equation becomes the following:
F = f 0 − f c k − f 1 k ln ⁡ f c − f 1 f 0 − f 1
Before ponding occurs, cumulative infiltration is given by F = I t , where I is the mean rainfall intensity. Ponding occurs when the infiltration rate equals the rainfall intensity I . Setting f c = I in Equation (3), one obtains the following:
F p = f 0 − I k − f 1 k ln ⁡ I − f 1 f 0 − f 1
where F p is the cumulative infiltration at ponding time t p [19,20], which is the following:
t p = F p I = f 0 − I k I − f 1 k I ln ⁡ I − f 1 f 0 − f 1

2.2. Overland Flow Model

Overland flow on any given plane is described by the KW approximation of the 1D Saint-Venant equation [1,2,21]:
∂ h ∂ t + 1 W ∂ Q ∂ x = r
where h is the water depth, W is the runoff width, Q is the runoff discharge, r is the rainfall excess intensity, t is time, and x is the distance in the flow direction.
The KW equation assumes that the friction slope equals the bed slope. By using existing open channel flow equations, a relationship between discharge Q and water depth h at any point on a plane can be expressed as follows:
Q = α W h β
where α is the KW resistance parameter and β is the Bakhmeteff dimensionless coefficient, assumed to be equal to 5/3. The parameter α is related to Manning’s roughness coefficient n , and it can be calculated as α = n   S 1 / 2 , where S is the bed slope of the plane.
Substituting Equation (7) in Equation (6), the KW equation can be rewritten as follows:
∂ h ∂ t + β α h β − 1 ∂ h ∂ x = r
A cascade of rectangular n-planes representing converging or diverging surfaces exposed to rainfall is shown in Figure 1. Conceptually, these n-planes could vary in their lengths, widths, and slopes.
Here, the numerical scheme for overland flow computation is described for the converging cascade planes as an example. The governing equation for overland flow (Equation (8)) is derived using the second-order single-step Lax–Wendroff numerical scheme, expressed in finite difference form as follows [21]:
h j i + 1 = h j i + ∆ t r j i − β α h j + 1 i β − 1 + h j − 1 i β − 1 2   h j + 1 i − h j − 1 i 2 ∆ x + ∆ t 2 2   r j i + 1 − r j i ∆ t − β α ∆ t 2 2 ∆ x   [ h j + 1 i β − 1 + h j i β − 1 2   r j + 1 i + r j i 2 − β α h j + 1 i β − 1 + h j i β − 1 2 h j + 1 i − h j i ∆ x − h j i β − 1 + h j − 1 i β − 1 2 r j i + r j − 1 i 2 − β α h j i β − 1 + h j − 1 i β − 1 2   h j i − h j − 1 i ∆ x ]
where i and j represent the grid points in time and space, respectively. For the downstream boundary of each plane, the following first order scheme is used [21]:
h j L i + 1 = h j L i + ∆ t r j L i − β α   h j L i β − 1 + h j L − 1 i β − 1 2   h j L i − h j L − 1 i ∆ x , for   x = L
where h j L refers to the last point in space x and L is the length of each plane along the flow direction. The initial condition for every plane is given as follows:
h x , 0 = 0 ,       for   0 ≤ x ≤ L    and   t = 0
The upstream boundary condition for the first plane k is given as follows:
h k 0 , t = 0 ,       for   x = 0   and   0 ≤ t ≤ ∞  
For the remaining cascade planes, the discharge leaving the downstream boundary of each plane enters the upstream boundary of the next plane, serving to establish the boundary condition required for the computation of the water depth of that plane.
As an example, the computation of the upstream boundary condition for plane n is given as follows, for the upstream boundary entering the discharge and water depth, respectively:
Q n − l x , t = α n − l W n − l h n − l β x , t ,     for   x = L n − l   ( plane   n l )   and   0 ≤ t ≤ ∞  
h n 0 , t = Q n − l α n W n 1 β ,     for   x = 0   ( plane   n )   and   0 ≤ t ≤ ∞
The numerical scheme needs to maintain the Courant number for stability, i.e.:
Δ t Δ x < 1 α β h β − 1 ( x , t )

2.3. Soil Erosion Model

Rainfall-induced soil erosion involves the processes of sediment detachment, transportation, and deposition [22]. Sediment detachment and transportation occur due to the impact of raindrops and the shear stress of overland flow, known as interrill and rill erosion. Deposition refers to the settling of eroded sediments. These processes are described by the 1D continuity equation of sediment transport [22,23]:
∂ ( c h ) ∂ t + ∂ ( c q ) ∂ x = e R − d + e I
where h is water depth, q is local discharge per unit width, c is sediment concentration in overland flow, e R is the rill erosion rate, d is the sediment deposition rate, e I is the interrill erosion rate, t is time, and x is the distance in the flow direction. The interrill erosion rate e I is given by the following equation [24]:
e I = K I   I   r
where K I   is the interrill coefficient, I is rainfall intensity, and r is rainfall excess intensity. Rill erosion occurs when sediments cannot resist the shear stress of overland flow, and its rate can be expressed as follows [22,25]:
e R = K R   τ e b
where K R is a soil detachability factor for shear stress, τ e is the average “effective” shear stress assuming broad shallow flow, and b is an exponent typically ranging from 1.0 to 2.0.
The sediment deposition rate d is expressed as follows [26]:
d = ∈ c   V s
where ∈ is a dimensionless coefficient depending on soil and fluid properties, c is the sediment concentration in overland flow, and V s is the settling velocity of the sediment particle, which can be expressed as follows:
V s = W o γ s − γ γ g d s
with
W o = 2 3 + 36 ν 2 g   d s 3 γ s γ − 1 − 36 ν 2 g   d s 3 γ s γ − 1
where γ s and γ are the specific weights of sediment and water, respectively, ν is the kinematic viscosity of water, g is the acceleration due to gravity, and d s is the size of the settling sediment.
The numerical procedure proposed to solve the sediment continuity equation (Equation (16)) uses the four-point implicit scheme [14,27]:
ϕ ∆ t c h j + 1 i + 1 − c h j + 1 i + 1 − ϕ ∆ t   c h j i + 1 − c h j i + ω ∆ x c q j + 1 i + 1 − c q j i + 1 − ω ∆ x c h j + 1 i − c q j i = ω ϕ e R j + 1 i + 1 + 1 − ϕ e R j i + 1 + 1 − ω   ϕ e R j + 1 i + 1 − ϕ e R j i − ω ϕ d j + 1 i + 1 + 1 − ϕ   d j i + 1 − 1 − ω   ϕ d j + 1 i + 1 − ϕ   d j i + ω e I i + 1 + 1 − ω   e I i
where ϕ is the weighting factor for space and ω is the weighting factor for time. If ϕ = 1 and ω = 0, the numerical scheme becomes explicit and needs a Courant number (Equation (15)) to be maintained for stability. If ϕ = 0.5 and ω > 0.5, the scheme becomes implicit and is unconditionally stable; however, for accuracy, the Courant number must still be maintained. For our simulations, we used ϕ = 0.5 and ω = 0.6, which made the numerical scheme stable for the computation of small volumes of runoff and sediment, typically associated with laboratory conditions, which were used in this study to test the performance of the model.
The initial and upstream boundary conditions for each cascade plane are as follows [14]:
c x , t p = K I   I   ( t p )   r ( t p ) ∈ V s + r ( t p ) ,   for   x ≥ 0   and   t = t p
c 0 , t = K I   I   ( t )   r ( t ) ∈ V s + r ( t ) ,   f o r   x = 0   a n d   t ≥ t p
These boundary conditions imply that the numerical scheme will start computing the sediment concentration from the ponding time t p . The downstream boundary condition was considered open for each cascade plane.

2.4. Calibration

To achieve a more accurate representation of the real-world system being studied and ensure that the model outputs align closely with the observed data, several model parameters can be adjusted, in particular infiltration parameters ( f 0 , f 1 , and k , in Equation (1)) and two soil erosion parameters ( K I   and K R in Equations (17) and (18), respectively). This process enhances the reliability and predictive capability of the model for the specific conditions under which it was calibrated.

3. Materials and Methods

The performance of the model proposed in Section 2 was evaluated by comparing numerical results with laboratory data from experiments using a rainfall simulator and a soil flume. This chapter presents the model testing procedure (Section 3.1), including a description of the laboratory experiments and the corresponding scheme used in the numerical model. Additionally, it describes the numerical simulations used to explore the effect of several parameters on the model’s response, along with the scenarios adopted for this purpose (Section 3.2).

3.1. Model Testing Procedure

3.1.1. Experimental Data

Data for testing the model were obtained from laboratory experiments conducted on a soil flume under controlled conditions, using simulated rainfall (Figure 2). The rainfall simulator’s sprinkler was equipped with a 3.58 mm orifice diameter HH-22 FullJet nozzle, installed 2.2 m above the soil flume surface, producing a full cone spray. This simulator has been used in several other studies [9,28]. Carvalho et al. [28] calculated the raindrop mean diameter and mean fall velocity as 0.74 mm and 1.92 m s−1, respectively, at a constant nozzle operational pressure of 50 kPa for the same nozzle used in this study. To maintain this constant pressure, a hydraulic system was employed, pumping water from a constant-head reservoir with a submerged pump. The rainfall spatial distribution was measured by placing rain gauges 0.15 m apart on the soil flume and measuring the collected rainfall volume during a 2-minute rainfall event. The uneven rainfall spatial distribution on the flume was a direct consequence of using a single full-cone nozzle, positioned in the middle of the flume surface approximately 1 m from the outlet, along the slope.
The soil flume had a planar circular sector geometry with a length of 2.02 m, limited by arc lengths of 0.10 m and 2.10 m, leading to the formation of converging and diverging flow conditions. The flume’s depth was 0.12 m, and it was constructed from zinc-coated iron metal sheets. A geotextile rug was placed over the iron mesh at the bottom of the soil flume to allow for free drainage. However, drainage flux was not measured, as the focus was solely on surface runoff measurements needed to validate the numerical model’s overland flow generation component. The soil flume had runoff collection outlets at both ends, and by tilting the flume in opposite directions, experiments could be performed on both converging and diverging plane surfaces.
After sieving through a 5 mm mesh, sandy loam soil was spread manually in the flume in layers and gently compacted with a steel plate. The soil surface was levelled using a wooden blade, with a final soil depth of 0.1 m occurring. Before each experiment, the soil was gently saturated with water (using a hose connected to tap water) until ponding occurred, followed by 1.5 h of drainage. After each experiment, the soil was removed, replaced with air-dried soil, and the saturation procedure was repeated to ensure consistent soil moisture across all experiments.
A total of four experiments were conducted: two on converging and two on diverging plane surfaces. Each experiment was performed once without repetition. The data on applied mean rainfall intensity and the slope of the converging and diverging plane surfaces are given in Table 1. The rainfall duration for each experiment was 7 min. High rainfall intensity (≈66 mm/h) and a steep gradient (20%) created conditions prone to high soil erosion, while lower rainfall intensity (≈37 mm/h) with a gentle gradient (5%) ensured low soil loss conditions on both surface geometries. This approach was used to test the model’s ability to simulate runoff and erosion for converging and diverging plane surfaces under these contrasting scenarios.
Runoff samples were collected at the flume outlet every minute for 10 s, from the initiation of runoff until the end of rainfall. Two additional samples were collected at 10-s intervals after the rainfall stopped. The runoff samples were weighted, oven-dried at 80 °C for 24 h, and weighted again to estimate runoff and soil loss.

3.1.2. Numerical Simulations Data

For the numerical simulations aiming at calibrating and testing the model (Section 2) based on the laboratory data, the surface was represented by a cascade of four planes (Figure 3), with a total area of 2.2 m2, which was approximately equal to the surface area of the soil flume (2.24 m2). Each plane had a constant length ( L ) of 0.5 m along the flow direction, while the widths ( W ) of the planes varied. For a given cascade (either converging or diverging), all planes had a similar slope. Since the rainfall spatial distribution over the flume was not uniform, the simulated rainfall distribution on the cascade planes was also uneven (Figure 3).
The infiltration parameters ( f 0 , f 1 , and k ) and two soil erosion parameters ( K I   and K R ) were optimized by trial and error to fit the observed hydrographs and sediment graphs (Table 1). The remaining soil erosion parameters (b, ∈, γs, γ, ds, and ν) were kept constant for all simulations (Table 2), which is justified as some are physical parameters, while others are expected to have a minimal impact on the model results.

3.1.3. Model Evaluation

The model’s skill in simulating runoff peak, total runoff volume, soil loss peak, and total sediment mass was assessed using the Percentage of Deviation, D e v ( % ) :
D e v ( % ) = S i m − O b s O b s × 100
where O b s and S i m represent the observed and simulated data, respectively.
The goodness of fit for the shape of the hydrographs and sediment graphs was evaluated using the Nash–Sutcliffe coefficient of efficiency ( N S ):
N S = 1 − ∑ t = 1 n O b s t − S i m t 2 ∑ t = 1 n O b s t − O b s m 2
where O b s t   and S i m t represent the observed and simulated data at time t, O b s m is the mean of the observed data, and n is the total number of data points (n = 9 in the experiments). The N S values range from −∞ to 1. As the N S value approaches 1, the model’s predictive efficiency improves. A value of 1 indicates a perfect fit between the model’s prediction and the observed data.

3.2. Numerical Simulation Scenarios for Varying Surface Convergence and Divergence Angles

A total of seven soil geometries with plane surface profile curvatures were numerically investigated, as shown in Figure 4. These included one rectangular/linear geometry, three diverging geometries, and three converging geometries. Each soil surface geometry had a constant area of 8 m2. Three distinct angles of convergence or divergence (θ along the flow direction), specifically 8.5°, 16.7°, and 24.2°, resulted in three different geometries for both converging and diverging configurations. For the numerical analysis, the converging and diverging soil surface geometries were modeled using a cascade arrangement composed of four planes. Each individual plane maintained a consistent length of 1 m along the flow direction, with variations introduced in the widths of these planes. The rectangular surface was modeled as a single plane measuring 4 m in length and 2 m in width. For any given cascade arrangement, all planes had an identical slope of 20%. The total simulation time for each soil surface geometry was 9 min, with a uniform rainfall intensity of 45 mm/h simulated for 5 min on all cascade planes. The infiltration ( f 0 , f 1 , and k ) and soil erosion ( K I   and K R ) parameters used were kept constant for all numerical simulations and are shown in Table 2. The remaining soil erosion parameters ( b , ∈ , γ s , γ , d s and ν ) used are also listed in Table 2, and coincide with the values used to test and calibrate the model (Section 3.1.2). In all numerical simulations, identical conditions were maintained to ensure that only the variation in soil surface geometry influenced runoff generation and soil loss estimation. In the converging surface simulations, the flow generated by rainfall moved from larger to smaller cascade planes, while in the diverging surface simulations, it moved in the opposite direction. The initial, upstream boundary and downstream boundary conditions for all planes in the model cascade are described in Section 3.2, particularly in Equations (10) to (14).

4. Results

4.1. Calibration and Testing of the Numerical Model

The infiltration parameters f 0 , f 1 , and k , as well as the two soil erosion parameters K I and K R , determined during the model calibration (Section 2.4) are listed in Table 1. The remaining soil erosion parameters ( b , ∈ , γ s , γ , d s , and ν ; see Section 2.3 and Section 3.1.2) used are listed in Table 2.
Simulated and observed hydrographs and sediment graphs for different rainfall intensity and slope combinations on converging and diverging plane surfaces are presented in Figure 5 and Figure 6, respectively. The numerical model’s simulated hydrographs aligned well with the observed ones, as indicated by the N S values close to 1 (Figure 5). The rising limbs of the simulated hydrographs matched the observed values closely, although there was a delay in the recession limbs.
Overall, the model effectively reproduced the main characteristics of sediment transport. All simulated sediment graphs showed good agreement with the observed data, supported by N S values consistently higher than 0.7 (Figure 6). However, the cascade model revealed a better capability to simulate runoff rates compared to soil loss rates.
Table 3 and Table 4 present the simulated and observed runoff peaks, runoff volumes, soil loss peaks, and total soil losses, along with the corresponding D e v ( % ) values, indicating that the model performed reasonably well. There were two exceptions: the total runoff in experiment 4 and the total soil loss in experiment 2, where D e v ( % ) values were higher. These cases involved lower rainfall intensities (≈37 mm/h) and gentle slopes (5%), which resulted in lower runoff volumes and sediment transport yields. Further adjustment of the model’s parameters could reduce the D e v ( % ) values, but this might decrease the N S values for the corresponding hydrographs and sediment graphs. Although the N S value, used as a model performance indicator, suggests that the model may face challenges in accurately estimating runoff volumes and total soil loss for smaller quantities, further testing of the model should be conducted under a broader range of conditions and parameters.
The observed soil loss on converging surfaces was substantially greater than on diverging surfaces due to the increased stream power generated under converging flow conditions [9]. For the higher intensity rainfall simulation (≈66 mm/h) and the steeper slope (20%), the total soil loss increased by 788%, whereas for the lower intensity (≈37 mm/h) and the milder slope (5%), the increase was 269%. In the converging flow experiment that combined both high rainfall intensity and a steep slope (experiment 1, in Table 3 and Table 4), small but distinct rills developed on the soil surface. Note that the total runoff volume for the converging flow conditions was also greater than that for the corresponding diverging flow conditions, with an increase of 39% for the higher rainfall intensity and steeper slope, and 178% for the other scenario. Overall, the numerical simulations align with this behavioral trend (Section 4.2).

4.2. Numerical Simulations for Plane Soil Surfaces with Varying Converging and Diverging Angles

The infiltration and soil erosion parameters used for each numerical simulation were constant and are listed in Table 2. The simulated values for peak runoff and peak and total soil loss across all simulated plane surface geometries are listed in Table 5.
Among the three converging geometries, the soil surface with a 24.2° angle exhibited the lowest runoff peak, while displaying the highest peak and total soil loss. Conversely, within the diverging geometries, the soil surface with an 8.5° angle recorded the lowest runoff peak and the highest peak and total soil loss.
Figure 7 displays simulated hydrographs for seven different plane soil surface geometries. There are differences between the rising and recession limbs of the hydrograph for the rectangular surface compared to those of the converging and diverging surfaces. The discharge increases faster under converging conditions in comparison to the rectangular and diverging surfaces.
The trend in the runoff peaks for all plane soil surface geometries are better depicted in Figure 8. Runoff peaks for converging surfaces were consistently lower than those for the rectangular surface; for the experimental data tested, the results were 0.5–1.4% lower. Notably, an increase in the converging angle of the plane soil surface geometry led to a linear decrease in the runoff peak (R2 = 0.99). Conversely, for diverging soil surface geometries, the runoff peaks were higher than those for the rectangular one; for the data tested, the results were 0.5–1.5% higher. An increase in the diverging angle of the soil surface geometry resulted in a linear increase in the runoff peak (R2 = 0.99).
Figure 9 and Figure 10 present the simulated sediment graphs and peak soil loss for the various plane soil surface geometries, respectively. These outputs also indicate a behavioral trend. For the varied tested scenarios, which led to a wide range of results, peak soil losses for converging surfaces were remarkably higher by 36–657% compared to the rectangular surface. The peak soil loss exhibited an exponential increase with the increasing converging angle of the soil surface geometry (R2 = 0.91). In contrast, diverging surfaces had a peak soil loss 20–41% lower than the rectangular surface, and an increase in the diverging angle of the soil surface geometry led to an exponential decrease in the peak soil loss (R2 = 0.99). The soil loss peaks of converging surfaces exceeded those of the diverging surfaces by a substantial 168–1250% for the converging/diverging angles studied. Compared to the rectangular surface, the total soil loss on converging surfaces was between 36% and 672% higher, while on diverging surfaces, it was 20% to 40% lower (Table 5), with angles of convergence or divergence along the flow direction varying from 8.5° to 24.2°.

5. Discussion

In this study, the concept of representing soil surfaces with convergent and divergent geometries using a cascade of n-planes accommodates spatial changes in these surfaces. Although the model is used here to simulate processes on plane surfaces, it can be adjusted for curved surfaces by varying the slope of the model planes, which broadens its applicability.
For the model version explored, the results reveal that when the flow leaves the downstream boundary of one plane and enters the upstream boundary of the next plane, it experiences an abrupt change (shock wave phenomenon [29]) due to the change in cascade geometry. This issue, observed with a model consisting of four cascading planes, can also affect the shape of the hydrographs. Therefore, the number of cascade planes used for any specific study area needs to be increased until such abrupt oscillations are minimized.
Furthermore, the numerical model was tested with observed laboratory data from single events (without replicates). Therefore, experimental errors (such as those from soil preparation, rainfall simulator handling, and runoff sampling) cannot be ruled out. Future work should consider varying the conditions (such as rainfall intensities and duration, slopes, surface shapes and dimensions, and soil types and covers) to further evaluate the model’s performance and its applicability at the hillslope scale in heterogeneous landscapes. Additionally, a larger number of observed simulations should be tested.

6. Conclusions

A 1D cascade model for estimating rainfall-induced runoff and soil erosion on various soil surface geometries (rectangular/plane, converging planes, and diverging planes) was proposed and tested using a soil flume and a rainfall simulator. Overall, for the tested slopes and rainfall intensities, the simulated hydrographs and sediment graphs compared well with the observed ones, as suggested by the Nash–Sutcliffe coefficients of efficiency ( N S ). However, the model’s ability to simulate the observed hydrographs ( N S > 0.87) was better than that to simulate the sediment graphs ( N S > 0.70).
To examine how varying the converging or diverging angle of plane soil surfaces affects runoff and soil loss, a series of simulations using the proposed model were conducted with seven different geometric configurations. These configurations included a rectangular/linear surface (angle along the flow direction θ = 0°), three surfaces each with distinct converging angles, and three surfaces each with unique diverging angles. The studied converging and diverging angles were 8.5°, 16.7°, and 24.2°. The converging and diverging plane surfaces were modeled with a four-plane cascade. Throughout all the simulation runs, the soil surface area, rainfall intensity, slope gradient, and infiltration and soil parameters remained constant. As the convergence angle of the plane soil surface increased, there was a reduction in the runoff peak but an increase in both the peak and total soil loss. Conversely, increasing the divergence angle of the plane soil surface led to an increase in the runoff peak, while causing a decline in both the peak and total soil loss. Compared to the rectangular soil surface, all converging geometries showed higher peak and total soil loss, while the diverging geometries displayed lower values. This indicates that converging plane soil surfaces are the most hazardous in terms of soil loss generation.
The presented model is thought to be a useful step towards simulating overland flow generation and sediment transport on convergent and divergent plane surfaces, avoiding the complexity of more advanced 2D models. Future studies should evaluate the ability of this model to simulate runoff and sediment loss on various curvatures (convex and concave) of different soil surface shapes (converging, linear, and diverging) under varying rainfall intensities, soil types, and slopes. The application of the model to larger areas, while also considering a greater number of planes in the cascade model, should be explored, as this will allow for the modeling of runoff and soil loss at the hillslope scale and in more heterogeneous landscapes. Additionally, the impact of different soil covers, which is prevalent at the hillslope scale [30,31], also needs to be taken into consideration.

Author Contributions

Conceptualization, J.L.M.P.d.L. and B.M.; methodology, J.L.M.P.d.L. and B.M.; modeling, J.L.M.P.d.L. and B.M.; validation, J.L.M.P.d.L. and B.M.; formal analysis, B.M., J.L.M.P.d.L., and M.I.P.d.L.; writing—original draft preparation, B.M. and J.L.M.P.d.L.; writing—review and editing, J.L.M.P.d.L. and M.I.P.d.L.; and funding management, J.L.M.P.d.L. and M.I.P.d.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly financed by National Funds through the FCT—Fundação para a Ciência e a Tecnologia (Portuguese Foundation for Science and Technology), credited to MARE https://doi.org/10.54499/UIDB/04292/2020 and https://doi.org/10.54499/UIDP/04292/2020 and to ARNET https://doi.org/10.54499/LA/P/0069/2020.

Data Availability Statement

Data that support the findings are available from the first author upon reasonable request.

Acknowledgments

Experiments were conducted at the Laboratory of Hydraulics, Water Resources and Environment of the Department of Civil Engineering of the Faculty of Sciences and Technology of the University of Coimbra (Portugal).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. de Lima, J.L.M.P.; Singh, V.P. The influence of the pattern of moving storms on overland flow. Adv. Water Resour. 2002, 25, 817–828. [Google Scholar] [CrossRef]
  2. Isidoro, J.M.G.P.; de Lima, J.L.M.P. An analytical closed form solution for 1D kinematic overland flow under moving rainstorms. J. Hydrol. Eng. 2013, 18, 1148–1156. [Google Scholar] [CrossRef]
  3. Liu, Q.Q.; Xiang, H.; Singh, V.P. A simulation model for unified interrill erosion and rill erosion on hillslopes. Hydrol. Process. 2006, 20, 469–486. [Google Scholar] [CrossRef]
  4. Deng, Z.Q.; de Lima, J.L.M.P.; Jung, H.S. Sediment transport rate-based model for rainfall-induced soil erosion. Catena 2008, 76, 54–62. [Google Scholar] [CrossRef]
  5. An, Y.; Liu, Q.Q. Two-dimensional hillslope scale soil erosion model. J. Hydrol. Eng. 2009, 14, 690–697. [Google Scholar] [CrossRef]
  6. Arguelles, A.C.C.; Jung, M.; Mallari, K.J.B.; Pak, G.; Aksoy, H.; Kavvas, L.M.; Eris, E.; Yoon, J.; Lee, Y.; Hong, S. Evaluation of an erosion-sediment transport model for a hillslope using laboratory flume data. J. Arid Land 2014, 6, 647–655. [Google Scholar] [CrossRef]
  7. Rieke-Zapp, D.H.; Nearing, M.A. Slope Shape Effects on Erosion. Soil Sci. Soc. Am. J. 2005, 69, 1463–1471. [Google Scholar] [CrossRef]
  8. de Lima, J.L.M.P.; Isidoro, J.M.G.P.; de Lima, M.I.P.; Singh, V.P. Longitudinal Hillslope Shape Effects on Runoff and Sediment Loss: Laboratory Flume Experiments. J. Environ. Eng. 2018, 144, 04017097. [Google Scholar] [CrossRef]
  9. Mujtaba, B.; de Lima, J.L.M.P.; de Lima, M.I.P.; Vargas, M.M. Importance of the hydraulics of converging and diverging flows on soil erosion override effects of rock fragment cover patterns: Experimental assessment. Z. Geomorphol. 2019, 62, 183–197. [Google Scholar] [CrossRef]
  10. Mombini, A.; Amanian, N.; Talebi, A.; Kiani-Harchegani, M.; Rodrigo-Comino, J. Surface roughness effects on soil loss rate in complex hillslopes under laboratory conditions. Catena 2021, 206, 105503. [Google Scholar] [CrossRef]
  11. Chen, T.; Shu, J.; Han, L.; Tian, G.; Yang, G.; Lv, J. Modeling the effects of topography and slope gradient of an artificially formed slope on runoff, sediment yield, water and soil loss of sandy soil. Catena 2022, 212, 106060. [Google Scholar] [CrossRef]
  12. Goodman, R.E.; Shi, G. Block Theory and Its Application to Rock Engineering; Prentice Hall: Englewood Cliffs, NJ, USA, 1985. [Google Scholar]
  13. Raghuvanshi, T.K. Plane failure in rock slopes—A review on stability analysis techniques. J. King Saud Univ. 2019, 31, 101–109. [Google Scholar] [CrossRef]
  14. Lopes, V.L. A Numerical Model of Watershed Erosion and Sediment Yield. Ph.D. Thesis, University of Arizona in Tucson, Tucson, AZ, USA, 1987. [Google Scholar]
  15. Smith, R.E.; Goodrich, D.C.; Woolhiser, D.A.; Unkrich, C.L. KINEROS—A kinematic runoff and erosion model. In Computer Models of Watershed Hydrology; Singh, V.P., Ed.; Water Resources Publications: Littleton, CO, USA, 1995; pp. 697–732. [Google Scholar]
  16. Horton, R.E. The role of infiltration in the hydrologic cycle. Eos Trans. Am. Geophys. Union 1933, 14, 446–460. [Google Scholar]
  17. Horton, R.E. The interpretation and application of runoff plane plot experiments with reference to soil erosion problems. Soil Sci. Soc. Am. J. 1938, 1, 401–437. [Google Scholar] [CrossRef]
  18. Horton, R.E. An Approach Toward a Physical Interpretation of Infiltration-Capacity. Soil Sci. Soc. Am. J. 1941, 5, 399–417. [Google Scholar] [CrossRef]
  19. Chow, V.T.; Maidment, D.R.; Mays, L.W. Applied Hydrology; McGraw Hill: New York, NY, USA, 1988. [Google Scholar]
  20. Brutsaert, W. Hydrology: An Introduction; Cambridge University Press: New York, NY, USA, 2005. [Google Scholar]
  21. Singh, V.P. Kinematic Wave Modelling in Water Resources: Surface-Water Hydrology; John Wiley and Sons Ltd.: Chichester, UK, 1996. [Google Scholar]
  22. Foster, G.R. Modeling the erosion process. In Hydrologic Modeling of Small Watersheds; Haan, C.T., Johnson, H.P., Brakensiek, D.C., Eds.; American Society of Agricultural Engineers: St. Joseph, MI, USA, 1982; Volume 5, pp. 259–380. [Google Scholar]
  23. Bennett, J.P. Concepts of mathematical modeling of sediment yield. Water Resour. Res. 1974, 10, 485–492. [Google Scholar] [CrossRef]
  24. Lane, L.J.; Shirley, E.D. Erosion and sediment yield equations: Solutions for overland flow. In Soil Erosion and Conservation, Proceedings of the Workshop on USLE Replacement, West Lafayette, IN, USA, 7–9 January 1985; Trimble, S.W., Ed.; USDA, National Soil Erosion Laboratory, Purdue University: West Lafayette, IN, USA, 1985; 22p. [Google Scholar]
  25. Croley, T.E. Unsteady overland sedimentation. J. Hydrol. 1982, 56, 325–346. [Google Scholar] [CrossRef]
  26. Mehta, A.J. Characterization tests for cohesive sediments. In Proceedings of the Conference on Frontiers in Hydraulic Engineering, ASCE/MIT, Cambridge, MA, USA, 9–12 August 1983; pp. 79–84. [Google Scholar]
  27. Lopes, V.L.; Lane, L.J. Modeling sedimentation processes in small watersheds. IAHS Publ. 1988, 174, 497–508. [Google Scholar]
  28. Carvalho, S.C.P.; de Lima, J.L.M.P.; de Lima, M.I.P. Using meshes to change the characteristics of simulated rainfall produced by spray nozzles. Int. Soil Water Conserv. Res. 2014, 2, 67–78. [Google Scholar] [CrossRef]
  29. Kibler, D.F.; Woolhiser, D.A. The Kinematic Cascade as a Hydrologic Model; Hydrology Papers 39; Colorado State University: Fort Collins, CO, USA, 1970. [Google Scholar]
  30. Bueno-Hurtado, P.; Seidou, O. Empirical and physical modelling of soil erosion in agricultural hillslopes. J. Hydrol. Hydromech. 2024, 72, 279–291. [Google Scholar] [CrossRef]
  31. Wang, C.; Fu, X.; Zhang, X.; Wang, X.; Zhang, G.; Gong, Z. Modeling soil erosion dynamic processes along hillslopes with vegetation impact across different land uses on the Loess Plateau of China. Catena 2024, 243, 108202. [Google Scholar] [CrossRef]
Water 16 02955 g001
Figure 1. Sketch and notation used for the cascade of n-planes representing converging (sloping to the left) or diverging (sloping to the right) cascades.
Figure 1. Sketch and notation used for the cascade of n-planes representing converging (sloping to the left) or diverging (sloping to the right) cascades.
Water 16 02955 g001
Water 16 02955 g002
Figure 2. (a) Laboratory setup used in the experiments, consisting of a rainfall simulator and a soil flume (top and bottom left); the flume has outlets on both ends for the converging plane (bottom middle) and diverging plane surface (bottom right) experiments. (b) Dashed lines represent contour lines relative to an arbitrary datum (ground level) for a 20% flume slope, while solid lines indicate the border wall of the flume’s geometry.
Figure 2. (a) Laboratory setup used in the experiments, consisting of a rainfall simulator and a soil flume (top and bottom left); the flume has outlets on both ends for the converging plane (bottom middle) and diverging plane surface (bottom right) experiments. (b) Dashed lines represent contour lines relative to an arbitrary datum (ground level) for a 20% flume slope, while solid lines indicate the border wall of the flume’s geometry.
Water 16 02955 g002
Water 16 02955 g003
Figure 3. On the right, a schematic sketch illustrates the approximate representation of the soil flume surface planar geometry using converging and diverging cascade planes, each with an equal length L of 0.5 m, measured along the direction of the plane slope. The width W of each plane is also shown, which is measured in the perpendicular direction. On the left, the variation in mean rainfall intensity across the planes of the converging and diverging cascades is shown.
Figure 3. On the right, a schematic sketch illustrates the approximate representation of the soil flume surface planar geometry using converging and diverging cascade planes, each with an equal length L of 0.5 m, measured along the direction of the plane slope. The width W of each plane is also shown, which is measured in the perpendicular direction. On the left, the variation in mean rainfall intensity across the planes of the converging and diverging cascades is shown.
Water 16 02955 g003
Water 16 02955 g004
Figure 4. Seven plane soil surfaces: one rectangular, three converging with different convergence angles θ, and three diverging with different divergence angles θ. The schematic sketch also illustrates the approximate representation of the soil surface planar geometry using converging and diverging cascade planes of equal length L (1 m). The width W of each plane is also shown. The dimensions of the rectangular plane soil surface (4 m × 2 m) are provided as well.
Figure 4. Seven plane soil surfaces: one rectangular, three converging with different convergence angles θ, and three diverging with different divergence angles θ. The schematic sketch also illustrates the approximate representation of the soil surface planar geometry using converging and diverging cascade planes of equal length L (1 m). The width W of each plane is also shown. The dimensions of the rectangular plane soil surface (4 m × 2 m) are provided as well.
Water 16 02955 g004
Water 16 02955 g005
Figure 5. Observed and numerically simulated hydrographs for converging and diverging surfaces with different combinations of rainfall intensity (I) and slope (S). The Nash–Sutcliffe ( N S ) coefficients are also shown.
Figure 5. Observed and numerically simulated hydrographs for converging and diverging surfaces with different combinations of rainfall intensity (I) and slope (S). The Nash–Sutcliffe ( N S ) coefficients are also shown.
Water 16 02955 g005
Water 16 02955 g006
Figure 6. Observed and numerically simulated sediment graphs for converging and diverging surfaces with different combinations of rainfall intensity (I) and slope (S). Note that the Y-axis scale is not the same on the top-left figure. The Nash–Sutcliffe ( N S ) coefficients are also shown.
Figure 6. Observed and numerically simulated sediment graphs for converging and diverging surfaces with different combinations of rainfall intensity (I) and slope (S). Note that the Y-axis scale is not the same on the top-left figure. The Nash–Sutcliffe ( N S ) coefficients are also shown.
Water 16 02955 g006
Water 16 02955 g007
Figure 7. Simulated hydrographs for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). Section A represents a rectangular surface where θ = 0°.
Figure 7. Simulated hydrographs for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). Section A represents a rectangular surface where θ = 0°.
Water 16 02955 g007
Water 16 02955 g008
Figure 8. Runoff peaks for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference.
Figure 8. Runoff peaks for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference.
Water 16 02955 g008
Water 16 02955 g009
Figure 9. Simulated sediment graphs for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). Section A represents a rectangular surface where θ = 0°. Note that the Y-axis scales are not the same.
Figure 9. Simulated sediment graphs for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). Section A represents a rectangular surface where θ = 0°. Note that the Y-axis scales are not the same.
Water 16 02955 g009
Water 16 02955 g010
Figure 10. Peak soil loss for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference. Note that the Y-axis scales differ.
Figure 10. Peak soil loss for plane soil surfaces with different angles (θ) of convergence (left) and divergence (right). A rectangular surface is represented by θ = 0°. Trend lines fitted to the data are included for reference. Note that the Y-axis scales differ.
Water 16 02955 g010
Table 1. Rainfall intensity, slope, and manually optimized infiltration and soil erosion parameters used in numerical simulations, determined during model calibration to fit the hydrographs and sediment graphs observed in the laboratory experiments, for convergent and divergent plane soil surfaces.
Table 1. Rainfall intensity, slope, and manually optimized infiltration and soil erosion parameters used in numerical simulations, determined during model calibration to fit the hydrographs and sediment graphs observed in the laboratory experiments, for convergent and divergent plane soil surfaces.
Exp.Soil
Surface		Rainfall
Intensity (mm h−1)		Slope
(%)		Infiltration ParametersSoil Erosion Parameters
f 0
(m s−1)		 f 1
(m s−1)		 k
(-)		 K I
(kg s m−4)		 K R
(kg m N−1.5 s−1)
1Conv. plane66.1202.50 × 10−53.0 × 10−60.0601.2 × 1070.52
237.451.25 × 10−53.5 × 10−60.0091.5 × 1070.80
3Div. plane65.8205.00 × 10−55.5 × 10−60.03015.0 × 1070.45
437.351.34 × 10−54.0 × 10−60.0071.5 × 1073.12
Table 2. Infiltration and soil erosion parameters used for all numerical simulation scenarios involving different surfaces’ geometries. This set of parameters partially differs from those used for calibrating and testing the model based on the laboratory data.
Table 2. Infiltration and soil erosion parameters used for all numerical simulation scenarios involving different surfaces’ geometries. This set of parameters partially differs from those used for calibrating and testing the model based on the laboratory data.
ModelParameterValueUnits
Infiltration f 0 2.5 × 10−5m s−1
f 1 3.0 × 10−6m s−1
k 0.02-
Soil Erosion K I 1.7 × 107kg s m−4
K R 0.1kg m N−1.5 s−1
b 1.5-
∈ 0.5-
γ s 15,650N m−3
γ 9980N m−3
d s 0.0004m
ν 1.3 × 10−6m2 s−1
Table 3. Observed and numerically simulated runoff peak and volumes for all experiments on converging and diverging surfaces. The rainfall duration was 7 min. The Percentage of Deviation D e v ( % ) of the numerical simulations relative to the laboratory experiments is also presented.
Table 3. Observed and numerically simulated runoff peak and volumes for all experiments on converging and diverging surfaces. The rainfall duration was 7 min. The Percentage of Deviation D e v ( % ) of the numerical simulations relative to the laboratory experiments is also presented.
Exp.Soil
Surface		Rain Intensity (mm h−1)Slope
(%)		DataRunoff Peak (mL s−1) D e v (%)Total Runoff Volume (L) D e v (%)
1Conv. plane66.120Obs.32.97−0.610.420.6
Sim.32.7610.48
237.45Obs.15.281.02.48−12.2
Sim.15.432.17
3Div. plane65.820Obs.29.66−5.17.508.2
Sim.28.168.12
437.35Obs.12.49−3.20.8949.6
Sim.12.091.32
Table 4. Observed and numerically simulated soil loss peak and total soil loss for all experiments on converging and diverging surfaces. The rainfall duration was 7 min. The Percentage of Deviation D e v ( % ) of the numerical simulations relative to the laboratory experiments is also presented.
Table 4. Observed and numerically simulated soil loss peak and total soil loss for all experiments on converging and diverging surfaces. The rainfall duration was 7 min. The Percentage of Deviation D e v ( % ) of the numerical simulations relative to the laboratory experiments is also presented.
Exp.Soil
Surface		Rain
Intensity (mm h−1)		Slope
(%)		DataSoil Loss Peak
(g m−2 s−1)		 D e v (%)Total Soil Loss
(g)		 D e v (%)
1Conv. plane66.120Obs.2.30−12.31598.140.8
Sim.2.021611.48
237.45Obs.0.16−3.870.59−27.7
Sim.0.1651.06
3Div. plane65.820Obs.0.27−5.7179.93−3.4
Sim.0.25173.87
437.35Obs.0.066.919.12−3.7
Sim.0.0718.41
Table 5. Numerically simulated runoff peak and peak and total soil loss for seven plane soil surfaces: one rectangular, and the others each with distinct converging or diverging angles. All planes had an identical slope of 20% and were subjected to a uniform rainfall intensity of 45 mm/h for 5 min. The Percentage of Deviation D e v ( % ) relative to the rectangular surface is also presented.
Table 5. Numerically simulated runoff peak and peak and total soil loss for seven plane soil surfaces: one rectangular, and the others each with distinct converging or diverging angles. All planes had an identical slope of 20% and were subjected to a uniform rainfall intensity of 45 mm/h for 5 min. The Percentage of Deviation D e v ( % ) relative to the rectangular surface is also presented.
SimulationSoil Surface
Geometry		Conv./Div. Angle
(Degrees)		Runoff Peak
(mL s−1)		 D e v (%)Soil loss Peak
(g m−2 s−1)		 D e v (%)Total Soil Loss
(g)		 D e v (%)
1Rectangular plane074.25 0.20 280.09
2Converging plane8.573.90−0.50.2735.7381.5336.2
316.773.56−0.90.44121.0625.85123.4
424.273.24−1.41.50657.12162.25672.0
5Diverging plane8.574.600.50.16−19.7224.83−19.7
616.774.961.00.13−32.2190.40−32.0
724.275.331.50.12−40.9167.26−40.3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mujtaba, B.; de Lima, J.L.M.P.; de Lima, M.I.P. One-Dimensional Numerical Cascade Model of Runoff and Soil Loss on Convergent and Divergent Plane Soil Surfaces: Laboratory Assessment and Numerical Simulations. Water 2024, 16, 2955. https://doi.org/10.3390/w16202955

AMA Style

Mujtaba B, de Lima JLMP, de Lima MIP. One-Dimensional Numerical Cascade Model of Runoff and Soil Loss on Convergent and Divergent Plane Soil Surfaces: Laboratory Assessment and Numerical Simulations. Water. 2024; 16(20):2955. https://doi.org/10.3390/w16202955

Chicago/Turabian Style

Mujtaba, Babar, João L. M. P. de Lima, and M. Isabel P. de Lima. 2024. "One-Dimensional Numerical Cascade Model of Runoff and Soil Loss on Convergent and Divergent Plane Soil Surfaces: Laboratory Assessment and Numerical Simulations" Water 16, no. 20: 2955. https://doi.org/10.3390/w16202955

APA Style

Mujtaba, B., de Lima, J. L. M. P., & de Lima, M. I. P. (2024). One-Dimensional Numerical Cascade Model of Runoff and Soil Loss on Convergent and Divergent Plane Soil Surfaces: Laboratory Assessment and Numerical Simulations. Water, 16(20), 2955. https://doi.org/10.3390/w16202955

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop