Next Article in Journal
Ecological Engineering and Restoration of Eroded Muddy Coasts in South East Asia: Knowledge Gaps and Recommendations
Next Article in Special Issue
Settling Analyses and Investigations of Clayey Sediments to Storage in Geotubes
Previous Article in Journal
The Role of Financial Flexibility on Enterprise Sustainable Development during the COVID-19 Crisis—A Consideration of Tangible Assets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Investigation of Equilibrium Dynamics in Braided Gravel Beds for the Preservation of a Sustainable Fluvial Environment

by
Marilena Pannone
* and
Annamaria De Vincenzo
School of Engineering, University of Basilicata, 85100 Potenza, Italy
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(3), 1246; https://doi.org/10.3390/su13031246
Submission received: 17 December 2020 / Revised: 20 January 2021 / Accepted: 21 January 2021 / Published: 25 January 2021

Abstract

:
Gravel bars have an important role in the exchange between surface and subsurface waters, in preventing and mitigating riverbank erosion, in allowing the recreational use of rivers, and in preserving fluvial or riparian habitats for species of fishes, invertebrates, plants, and birds. In many cases, gravel bars constitute an important substrate for the establishment and development of ground flora and woody vegetation and guarantee higher plant diversity. A sustainable management of braided rivers should, therefore, ensure their ecological potential and biodiversity by preserving a suitable braiding structure over time. In the present study, we propose an analytical–numerical model for predicting the evolution of gravel bars in conditions of dynamical equilibrium. The model is based on the combination of sediment balance equation and a regression formula relating dimensionless unit bedload rate and stream power. The results highlight the dependence of the evolving sediment particles’ pattern on the ratio of initial macro-bedforms longitudinal dimension to river width, which determines the gradual transition from advective and highly braiding to diffusive transport regime. Specifically, the tendency to maintain braiding and flow bifurcation is associated with equilibrium average bed profiles and, therefore, equilibrium average stream power characterized by the maximum period that does not exceed transverse channel dimension.

1. Introduction

Braided rivers are typically organized in multithread patterns, with central or alternate bars dividing channels whose reach-scale statistical properties are controlled by discharge, bed slope, and grain size [1,2,3,4,5]. Braiding in gravel-bed streams may derive from one of several mechanisms, including chute cut-off of point or alternating bars, flow division around a transverse bar, dissection of multiple bars, and central bar formation [6]. Interactions among braided pattern, channel morphology, and bedload transport are significant [7]. Indeed, bedload transport in gravel-bed rivers is a very intricate phenomenon that is characterized by high variability in time and space [8,9,10,11]. Such variability derives from different mechanisms that act at different spatial scales [12,13,14,15], ranging from grain to reach: (1) the movement of single grains (e.g., [16,17,18,19]); (2) the creation of small bed forms, such as sediment waves [6], bed waves [20], and bedload sheets [21,22]; and (3) large-scale morphological processes, such as bar formation and migration, activation and abandonment of anabranches and side-channel shift [6,23,24]. The combination of all these processes results in a fluctuating nature of bedload transport, which was initially recognized by field surveys [25,26,27,28] and later confirmed by innovative measurement procedures [29,30,31,32]. As a matter of fact, several experimental studies performed under steady flow and sediment-supply conditions highlighted huge fluctuations of bedload transport rates both in single-thread [12,13,22,33] and braided channels [6,33]. At grain scale, channel morphological dynamics are the result of single-particle displacements from erosion sites and sediment diffusion in river bedload (e.g., [16,17,18,19]). Specifically, ref. [16] introduced a new conceptual model for the diffusion of moving-bed particles, which suggests that particle motion is diffusive in nature and includes at least three ranges of temporal and spatial scales characterized by different diffusion regimes. Later on, by means of particle tracers, ref. [34] confirmed the diffusive nature of particle displacement in real gravel-bed rivers and investigated gravel particle paths and their contribution to self-forming processes, such as gravel bar development. Particle paths were found to be mostly bank-line parallel, even in river bends. Additionally, transverse bedload was evaluated as 4% of the longitudinal counterpart. Results from physical models of braided rivers (e.g., [24]) highlighted the influence of bars in acting as deposition sites for moving particles, with over 80% of sediment tracers that were found at side, diagonal, and point bars. So far, very few examples of analytical or semi-analytical treatment of sediment continuity equation are available in the literature (e.g., [35]). In any case, none of these provides a solution in terms of actual bed elevation function. A two-dimensional fully numerical model was developed by [36] to simulate bar growth, channel widening, and scour holes at bar lee. In the present study, we make a system of two-dimensional sediment mass balance and a suitable constitutive relationship that relates dimensionless bedload rate and stream power in the case of braided gravel beds at dynamical equilibrium. Such a relationship was derived by the authors in a previous high-grain Reynolds numbers experimental study [3]. The crucial assumption enabling the analytical–numerical solution of the system in terms of actual bed elevation function is represented by the recognition of the critical state of flow along the braids and, as a consequence, by the equivalence of local section-average energy loss and local section-average bed slope (e.g., [3,37,38,39,40,41]). Ten experiments were carried out by [3] with fluvial gravel (D50 = 9 mm, D10 = 5 mm, and D90 = 13 mm) to analyze the relationship between flow stream power and bedload transport rate in braided rivers at high grain Reynolds numbers. The experimental runs were conducted in a rectangular, 1 m wide, and 16.7 m long flume with Plexiglas walls. The experimental device consisted of a water re-circulating/sediment feeding system. Gravel was added from a 2.3 m long sediment-supplying basin (located upstream of the flume and downstream of water inflow, with imposed average bottom slope) by means of a sediment feeder mounted on a travelling crane. The mixture of sediment and water coming from the supplying basin entered the flume and modified its slope until the achievement of equilibrium conditions. The runs differed from each other in terms of flow rate and initial slopes of flume and sediment-supplying basin bottom. See Table 1 for a summary of the experimental conditions, where S0 and Seq indicate initial and equilibrium flume bottom slopes, S0SB initial supplying basin bottom slope, and teq the running time needed to achieve equilibrium.
The equilibrium grain Reynolds number, which ranged from 601 to 710, was consistent with the condition of fully turbulent flow. Flow rates were measured by an electromagnetic current meter located on the recirculation pipe. For every run, water surface and bed elevations were measured by a manual point gauge mounted on a travelling carriage at 55 flume and 10 supplying basin cross-sections on average. Each of the flume cross-sections, which were spaced 30 cm apart, was monitored over time at 11 grid-points spaced 10 cm apart, with 605 points surveyed at the different sampling times. Sampling times number varied from 3 (R1) to 9 (R9), as a function of the total running time. As a consequence, during each single run, the total number of bed and free surface elevation measurements ranged from 1815 (for R1) to 5445 (for R9). Sediment flux at the flume inlet, which was proportional to the imposed sediment supplying basin bottom slope and water discharge, was estimated as the average volume of sediment eroded by water per unit time from the sediment-supplying basin. The average bedload output rate of the flume was measured by evaluating the cumulative weight of sediment conveyed into a flume downstream tank. The achievement of equilibrium conditions was determined by verifying that input and output bedload rates were equal and constant.
For each run, at equilibrium, section-averaged dimensionless unit bedload rate q ¯ b x * (where x indicates stream-wise direction, the over-bar indicates section-average and the asterisk dimensionless quantities) was evaluated at flume scale (which is representative of a fluvial reach) as a function of the section-averaged dimensionless unit stream power ω ¯ * [3]. Data were then analyzed and compared to those available from previous laboratory [2,42,43,44] and field studies [1,45,46,47]. The analysis highlighted how different combinations of grain Reynolds number and stream power range led to different types of equilibrium correlation laws (linear or power). In particular, for high grain Reynolds numbers (>601) and for ω ¯ *   < 1, which are the conditions investigated in the present work as typical field conditions, the experimental data could be well-interpreted by the following power law:
q ¯ b x *   =   0.1 ω ¯ * 2.02
Equation (1) was obtained by regression from our 10-run experimental laboratory data and field data collected along 66 braided gravel-bed rivers all over the world.
As reported by [3], in the case of braided gravel-bed rivers at high grain Reynolds numbers, Equation (1) estimates bedload transport rates more accurately than other formulas previously appeared in the literature [2,48,49,50]. Data of water surface elevations exhibited small oscillations about an average value during every run. Free surface oscillations confirmed the establishment of local critical conditions in the braids, as expected from the equilibrium Froude number of the experimental tests, which ranged between 1.15 and 1.33.

2. Methods

Let function η   =   η x , y , t represent bed elevation above a given reference level, with y and t indicating cross-sectional coordinate and time counted from equilibrium, respectively. We will proceed by subdividing η into section-average ( η ¯ ) and corresponding deviation ( η ):
η x , y , t   =   η ¯ x , t + η x , y , t   =   a i f e q x + N x , t + η x , y , t
In Equation (2) and in what follows, a indicates a constant, ifeq the reach-scale equilibrium bed slope, N the de-trended section-averaged bed profile (which reflects large-scale bedforms, such as bars) and η   accounts for small-scale bedforms, such as sediment waves or sediment sheets. Additionally,
η ¯ x , t   =   1 B 0 B η x , y , t d y
where B indicates total bed width. Unit-width sediment continuity reads:
1 p η t + · q b   =   s x , y , t
where p is sediment porosity, s indicates a source term (possibly attributable to slope or tributary sediment inputs) and qb the bedload rate vector per unit channel width:
q b   =   q b x , q b y
with:
q b x   =   q b x x , y , t   =   q ¯ b x x , t + q b x x , y , t
and
q b y   =   q b y x , y , t q b y x , y , t
Note that the hypothesis of negligible average transverse bedload and total transverse bedload practically coinciding with its deviation, expressed by Equation (7), is consistent with the experimental analysis carried out by [34], which highlighted how, in braided rivers, transverse bedload is only a small percentage of the longitudinal counterpart. The substitution of Equation (2) and Equations (5)–(7) into Equation (4) followed by section-averaging yields:
1 p N t + q ¯ b x x   =   S x , t
with:
S x , t   =   1 B 0 B s x , y , t d y
The difference of Equations (4) and (8) is:
1 p η t + q b x x + q b y y   =   s x , y , t
where:
s x , y , t   =   s x , y , t S x , t
As above summarized, from their previous laboratory experiment [3], the authors had found that in braided gravel beds at high grain Reynolds numbers and river reach-scale equilibrium:
q ¯ b x *   =   q ¯ b x g ϱ s ϱ ϱ d 3 0.1 ω ¯ * 2   =   0.1 ω ¯ 2 γ 2 g ϱ s ϱ ϱ d 3
where the average unit stream power is expressed by:
ω ¯   =   γ Q J ¯ B
In Equations (12) and (13) and in what follows, d indicates representative grain size, Q flow rate, g acceleration due to gravity, ρ water density (with γ = ρg), ρs sediment density, and j ¯ section-averaged energy loss. It is worth recalling here that the term “reach-scale equilibrium” refers to a situation in which the input bedload rate to the reach is equal to the output bedload rate from the reach, and the reach-scale bed trend slope is equal to the reach-scale free surface slope. Since in braided gravel beds water flows in critical conditions within each of the braids (e.g., [3,37,38,39,40,41], the corresponding section-averaged specific energy   E ¯ s , given by the sum of section-averaged kinetic energy U(x)2/2g (with U(x) indicating section-averaged velocity) and section-averaged flow depth H(x), is characterized by zero-mean longitudinal derivative (e.g., [51]):
d E ¯ s d x   =   d d x H x + U 2 x 2 g   =   d d H H x + Q 2 2 g b 2 H 2 x d H d x   =   1 F r 2 d H d x   =   0
where b indicates the equilibrium-constant total wetted width and Fr section-averaged Froude number ( F r   =   U / g H ). In these conditions, section-averaged energy loss is equal to section-averaged bed slope:
J ¯   =   E ¯ x   =   x η ¯ + E ¯ s   =   η ¯ x
Therefore:
J ¯   =   i f e q N x
and
q ¯ b x   =   Γ i f e q N x 2
with:
Γ   =   0.1 Q 2 B 2 g ϱ s ϱ ϱ d 3
Note that values of j ¯ < 0 do not constitute a physical contradiction in terms of energy valley-downstream variation because section-averaged bed elevation η ¯ does not have a physical but just a mathematical meaning. The same value of η ¯ can be representative of completely different transverse η-distributions. Water cannot be on top of bars and only flows within highly meandering braids. Equation (12) may also be thought as deriving from a formally analogous equation relating the actual (total) quantities:
q b x *   =   q ¯ b x + q b x g ϱ s ϱ ϱ d 3   =   0.1 ω ¯ + ω 2 γ 2 g ϱ s ϱ ϱ d 3
In this case,
q b x   =   q ¯ b x + q b x   =   Γ i f e q N x η x 2   = = Γ i f e q N x 2 + Γ η x 2 2 Γ i f e q N x η x
with q ¯ b x given by Equation (17) and
q b x   =   Γ η x 2 2 Γ i f e q N x η x
By analogy, we will assume (and we will justify later) that:
q b y   =   Γ η y 2 2 Γ i f e q N x η y
Substitution of Equation (17) into Equation (8) and Equations (21) and (22) into Equation (10), respectively, gives a variable-coefficient diffusion equation and a variable-coefficient advection-diffusion equation with source terms:
N t   =   2 Γ i f e q N x 1 p 2 N x 2 + S x , t 1 p
and
η t + 2 Γ 1 p 2 N x 2 η x   =   = 2 Γ 1 p i f e q N x η x 2 η x 2 + 2 Γ 1 p i f e q N x η y 2 η y 2 + s x , y , t 1 p
Note that we are consistently using the term diffusion equation and advection-diffusion equation for (23) and (24), respectively, because both N and η represent a length proportional to the section-averaged or local concentration of sediment particles. A linearization of Equations (23) and (24) will be achieved here by taking the reach-averaged values of N / x , η / x and η / y in the parentheses on their right-hand sides. The reason for doing this is that, given the experimentally detected wave-like behavior of functions N and η , the parentheses on the right-hand side of Equations (23) and (24) represent a sort of oscillatory diffusion coefficients, which can assume positive as well as negative values depending on the sign of local and section-average bed slope. Considering a generalized version of Fick’s law, “oscillatory diffusion coefficients” means that at small (in the case of N / x ) and very small (in the case of η / x and η / y ) scale, there is a pseudo-periodical alternation of zones characterized by a net flux of sediment towards higher concentrations and zones characterized by a net flux of sediment toward lower concentrations, which overall creates a mutual compensation. Note that the same type of compensation does not take place in the case of the oscillatory velocity 2 Γ 1 p 2 N x 2 in Equation (24), because advection moves the particles regardless of their local concentration. Based on the experimental data collected during the experimental runs performed by [3] and, as expected, due to the increasing frequency of decreasing-scale components of the bed elevation function, it results that:
i f e q N x ¯ R η x ¯ R , η y ¯ R
where the over-bar and subscript R jointly indicate reach-average. As an example, in the case of the reference run R4 (see Table 1), i f e q   =   0.015 , N / x ¯ R   =   0.002 , η / y ¯ R   =   1.538   × 10 4 , η / x ¯ R   =   1.446   × 10 19 . The conditions expressed by Relationship (25) enabled the following approximation of Equations (23) and (24), respectively:
N t   =   2 Γ i f e q 1 p 2 N x 2 + S x , t 1 p
and
η t + 2 Γ 1 p 2 N x 2 η x   =   2 Γ i f e q 1 p 2 η x 2 + 2 η y 2 + s x , y , t 1 p
It should be emphasized that the linear (and, in the above perspective, leading) parts of Equations (21) and (22) can respectively be re-written as:
q b x   =   2 Γ i f e q N x η x   =   0.2 F r 2 ε r 2 η x τ ¯ b * g ϱ s ϱ ϱ d 3
and
q b y   =   2 Γ i f e q N x η y   =   0.2 F r 2 ε r 2 η y τ ¯ b * g ϱ s ϱ ϱ d 3
where τ ¯ b * indicates dimensionless section-averaged bed shear-stress:
τ ¯ b *   =   ρ H i f e q N / x ρ s ρ d ,
and εr relative bed roughness ( ε r   =   d / H ). Note also that Equation (29) is fully compatible with that given in more general cases for mild secondary currents by [52] and later by [53], when put in the form:
q s y   =   α r η y τ s * g ϱ s ϱ ϱ d 3
In Equation (31), α is a dimensionless constant (equal to 4.93 in the study by [53]), and r is a coefficient given by the ratio of lift to drag forces. In Tealdi’s et al. model [53], qsy is the cross-sectional bedload rate along the sloping and counter-sloping side-wall (with the slope expressed by η / y ) and τ s * the dimensionless shear-stress at the middle bank of a trapezoidal cross-section. In the present model, τ ¯ b *   is obtained as a function of the section-averaged flow depth, which is practically measured starting from the average braided-bed level; additionally, q b y q b y represents the (zero-mean) cross-sectional bedload rate along the sloping and counter-sloping braid walls. Therefore, Equations (29) and (31) may be considered as coinciding with α = 0.2 and r, which represents the hydrodynamical engine of bedload transport, consistently given by the ratio of Froude number squared (a measure of lifting potential magnitude) to relative bed roughness squared (a measure of drag potential magnitude).
Equations (26) and (27) are two coupled partial differential equations. The first is a diffusion equation with constant diffusion coefficient D   =   2 Γ i f e q / 1 p   and space–time dependent source term S x , t / 1 p . The second is an advection–diffusion equation with constant diffusion coefficient D   =   2 Γ i f e q / 1 p , space–time dependent velocity v x , t   =   2 Γ / 1 p 2 N / x 2 and space–time dependent source term s x , y , t / 1 p . Variable-coefficients advection–diffusion equations were already solved by the authors in a stochastic framework in the case of river-flow transport of dissolved chemicals or fine suspended sediments in the presence of stationary bed morphology heterogeneity [54,55,56]. In the present case, due to the intrinsically non-stationary nature of bed morphology heterogeneity in braided rivers (e.g., [3]), no analytically treatable stochastic approach could apply. Note that (26) and (27) jointly reproduce the diffusive nature of particle displacement in gravel beds highlighted by [16] and later experimentally verified by [34] in the case of real gravel-bed rivers. From (26) and (27), and in the absence of external forcing ( S x , t   =   0 , s x , y , t   =   0 ) , one can write:
η t + 2 Γ 1 p 2 N x 2 η x   =   η t + 1 i f e q N t η x   =   2 Γ i f e q 1 p 2 η x 2 + 2 η y 2
The meaning of Equation (32) is that, besides the isotropic diffusion undergone by the small-scale bed undulations represented by η , the latter experience upstream migration where section-averaged bed elevation tends to decrease ( N / i f e q t   =   2 Γ / 1 p 2 N / x 2 < 0 ) and downstream migration where the section-averaged bed elevation tends to increase ( N / i f e q t   =   2 Γ / 1 p 2 N / x 2 > 0 ) . As an elementary example, let us start from a section-averaged bed elevation given by the superposition of the linear equilibrium trend and a single cosine, which may mathematically represent local changes in bed elevation induced by large-scale bedforms:
η ¯ x , 0   =   a i f e q x + N x , 0   =   a i f e q x + A cos π x L
where A and L indicate initial wave height and initial wave half-length, respectively. From Equation (26) with no source term, one obtains:
N x , t   =   Aexp π 2 D L 2 t cos π x L
where N = 0 (that is, at x = x0 = L/2), section-averaged bed elevation coincides with the elevation obtainable from the only reach-scale trend: η ¯ x 0 , t   =   a i f e q x 0 . For x < x0, N-profile is convex: 2 N / x 2 < 0   and the small-scale three-dimensional undulations represented by η migrate upstream (v(x, t) < 0) with an exponentially decaying velocity; conversely, for x > x0, N-profile is concave: 2 N / x 2 > 0   and the small-scale undulations migrate downstream (v(x, t) > 0) at the same gradually decreasing rate. Note that the convex N-profile is associated with N > 0 and η ¯ x , t > a i f e q x ; the vice versa holds for the concave counterpart. Considering that the critical state of flow implies i f e q i c , where i c indicates reach-scale critical bed slope, the quantity η ¯ c   =   a i f e q x may be assimilated to the reach-scale critical bed elevation. In this case, hinging on basic hydraulics (e.g., [51]), the condition η ¯ x , t > η ¯ c (i.e., the condition of local section-averaged bed elevation larger than the critical one) can be associated with a locally slow (sub-critical) sediment wave. The vice versa can be said in the case of concave N-profile, which would turn out to be associated with a locally fast (supercritical) sediment wave (Figure 1). Now, again from basic hydraulics, we know that small perturbations (here represented by η ) typically go upstream in the presence of sub-critical free-surface flows; conversely, they go downstream in the presence of supercritical flows. Therefore, that result makes our model perfectly compatible with a sediment-stream version of the classic gradually-varied flows theory. Specifically, we know that the approximate integration of de Saint-Venant equations for rectangular channel section [57]:
U H x + H t   =   0
H x + U g U x + 1 g U t   =   i f J ¯
for i f J ¯ leads to:
d U d H   =   ± g H
and
U ± g H H x + H t   =   0
Note that Equation (38) represents free-surface waves that propagate with celerity:
c   =   U ± g H
Therefore, for slow flows ( U / g H < 1 ) and d U / d H < 0 , c < 0 and the perturbations go upstream; on the other hand, for fast flows ( U / g H > 1 ) and d U / d H < 0 , c > 0, and the perturbations go downstream.

3. Results

To solve system (26)–(27) showing that it can explain the tendency of real rivers to maintain braiding under particular conditions and follow and monitor their evolution from an equilibrium configuration to the other, we resorted to a mixed analytical–numerical method. Specifically, we solved Equation (26) analytically for S(x,t) = 0 and an initial condition represented by an oscillatory pseudo-periodic function that is qualitatively consistent with our laboratory experimental observations [3]:
N x , 0   =   Δ c o s 2 π x Λ c o s 2 3 π x Λ
where Δ indicates the amplitude and Λ the characteristic length-scale, obtaining:
N x , t   =   Δ exp 4 π 2 D Λ 2 t cos 2 π x Λ exp 12 π 2 D Λ 2 t   cos 2 3 π x Λ
Then, we solved Equation (27) by particle-tracking (e.g., [58,59,60,61]) for s x , y , t = 0 and:
η x , y , 0   =   ξ h x m y
where ξ has the dimensions of a length, and h(x) and m(y) are two numerical coefficients defined as:
h x   =   1     0 x B 0   elsewhere
m y   =   1     B 2 β y B 2 + β 0   elsewhere
The reference macro-data set belongs to a river Basento reach (Basilicata region, southern Italy, Figure 2 and Figure 3) that can be considered as truly braided according to several classification criteria [62,63,64,65,66]: B = 150 m, Q = 305 m3/s, ifeq = 0.01, p = 0.24, d = 0.045 m, ρs = 2700 Kg/m3 and Fr = 1.09.
Additionally, we assumed Δ = 1 m, ξ = 2d = 0.09 m, and β = d/2 = 0.0225 m. Therefore, the initial bed elevation deviation turned out to be represented by a 9 cm-high double array of sediment particles instantaneously placed about the channel axis between x = 0 and x = B = 150 m (total number of particles: 6667), with a consistently almost zero average deviation:
1 B 0 B η x , y , 0 d y   =   2 β ξ h x B   =   2.7   x   10 5   0 x B 0   x < 0 ;   x > B
Provided that the linear structure of (26) and (27) enables the superposition principle, the aim of the present numerical experiment was to investigate how an elementary stream-wise perturbation evolves toward different types of braided or unbraided configurations as a function of the crucial ratio Λ/B. The application to real initial conditions will be the object of future research, when evolving post reach-equilibrium DEMs (Digital Elevation Models) will be available from the literature or our ongoing supplementary laboratory experiments.
For particle-tracking purposes, the barycenter of each sediment particle was virtually tracked based on the following discretized solution of the trajectory equation:
X t + d t   =   X t + v X t , t d t + 2 D d t G 0 , 1
Y t + d t   =   Y t + 2 D d t G 0 , 1
where X and Y, respectively, indicate the stream-wise and transverse component of the particle position vector, G (0,1) indicates a zero-mean and unit-variance Gaussian distribution, dt is the elementary time step and:
v x , t = 2 Γ 1 p 2 N x 2 = 0.8 Q 2 π 2 Δ B 2 Λ 2 1 p g ϱ s ϱ ϱ d 3 · 3 e x p 12 π 2 D Λ 2 t c o s 2 3 π x Λ e x p 4 π 2 D Λ 2 t c o s 2 π x Λ
D = 0.2 Q 2 i f e q B 2 1 p g ϱ s ϱ ϱ d 3
The elementary time step was chosen in such a way as to limit the length of the elementary purely diffusive displacement to a small percentage of channel total width:
d t   =   d n 2 D
with:
d n   =   B 2000
When a particle overcame a side boundary (representing the impervious bank), it was back-reflected by d Y   =   2 D d t 2 G 0 , 1 , where d t 2   =   d t d t 1 and dt1 was the time needed to reach the boundary itself. Note that, during dt2, the particle was not allowed to move along the longitudinal direction. It should be emphasized that, although our particle-tracking is 2-D in the horizontal projection of the channel plan view, the effect of lift and drag is inherent in the velocity and diffusion coefficient of the advection–diffusion equation that it solves. There are in the literature a few examples of 3D particle tracking applied to the study of saltating particles diffusion in turbulent water flow in more general conditions (e.g., [67]). The peculiarity of our study consists in the specific and more restricted field of application (braided gravel beds) and in the analytical nature of the model, allowed by the combination of the high grain Reynolds numbers quadratic relationship between dimensionless unit bedload rate and stream power and the typical critical state of flow, which determines the way the stream power is evaluated. Figure 4 shows the velocity distribution sampled for Λ = B/8 by the double sediment array at t = t0 and t = τd, with:
τ d   =   Λ 2 4 π 2 D
representing a characteristic diffusive time (i.e., the time needed by pure diffusion to cover a distance proportional to Λ). Figure 5 shows the solution of Equation (27) in terms of particles’ barycenter position at t = τd. Finally, Figure 6 shows the corresponding X-histogram. As one can see, the high-frequency oscillations of velocity induce a marked fragmentation of the initial sediment strip and a relatively reduced diffusion. The histogram displaying X-frequency (f) distribution exhibits an almost two-peaked shape with a central minimum. Figure 7, Figure 8 and Figure 9 show the same as Figure 4, Figure 5 and Figure 6 but referred to Λ = B/4. In this case, the initial velocity was less oscillating and overall characterized by a smaller magnitude. The number of sediment accumulation sites was reduced while their dimensions were increased due to the gradual transition from advective to diffusive regime. The X-histogram revealed a multi-peaked structure with higher concentrations of sediment in a lower number of longitudinal bands. Figure 10, Figure 11 and Figure 12 refer to Λ = B/2. The initial velocity was now rather low and only weakly oscillating. The corresponding particles’ distribution after the actual τ d   highlights the formation of three well-defined sediment accumulation sites in about 190 m, associated with three histogram peak zones. Note that the case-study in Figure 3, which represents a clearly braided branch of river Basento, and was characterized by Λ/B ≅ 1/2, include three large bars in about 220 m. The case Λ = B (Figure 13, Figure 14 and Figure 15) was the first characterized by a monotonic initially sampled velocity. From Figure 14, one can appreciate the dominant effect of diffusion (many sediment particles had already reached the boundaries at   t   =   τ d ) with considerable back-displacements (X < 0) and the slow transformation of the bed configuration into a two-dimensional sand–waves pattern. This tendency finds its maximum expression for Λ = 2B (Figure 16, Figure 17 and Figure 18), with channel width that is now almost uniformly occupied by sediment particles and the X-histogram that is rather close to a (skewed) Gaussian distribution, as in Taylor-like diffusion. Since the numerical experiments were conducted with the same B, Q, ifeq, ρs, p, and d and, therefore, with the same diffusion coefficient D, we can conclude that the tendency to maintain braiding and flow bifurcation is associated with equilibrium average bed profiles N and, therefore, equilibrium average stream power γ Q i f e q N / x characterized by maximum period Λ that does not exceed transverse channel dimension.
Finally, in Table 2, we summarize the main numerical parameters and results: the ratio Λ/B, the diffusion time τd, the Péclet number Pe = Λvmax/D—which is representative of the relative importance of advection and diffusion—and the braiding index BI, which was defined as the ratio of channel length occupied by sediment particles at t = τd to velocity characteristic length Λ. As one can see, the smaller the ratio Λ/B, the smaller the time needed for diffusion over Λ (τd), the larger Pe and the larger BI.
For the sake of completeness, and to analyze the free evolution of the equilibrium section-averaged bed profile, Figure 19, Figure 20, Figure 21, Figure 22 and Figure 23 show N(x) at t = t0 and t = τd for Λ = B/8, Λ = B/4, Λ = B/2, Λ = B and Λ = 2B, respectively. As one can see, due to the specific structure of (26), the pseudo-periodic behavior of the initial distribution tends to approach a perfectly periodic (other than considerably attenuated and characterized by lower frequency) oscillation, which is in turn affecting the local deviation evolution shown in the previous figures.

4. Discussion and Conclusions

In the present study, we propose an analytical–numerical methodology for the reconstruction of the equilibrium dynamics of gravel braided beds. The mathematical model is based on the combination of sediment continuity with or without forcing term and the practically quadratic relationship between dimensionless unit bedload rate and dimensionless unit stream power, which was derived by the authors in a previous work by regression, including ad hoc laboratory experimental measurements performed at high grain Reynolds numbers and real river data available in the literature. The crucial assumption enabling a semi-analytical solution in terms of bed elevation function was represented by the recognition of the critical state of flow, widely documented in the literature, which makes local section-averaged energy loss coinciding with local section-averaged bed slope. The first important result of the study is represented by the specific structure of the bed elevation deviation governing equation, coupled with a variable-coefficient diffusion equation holding for the section-average component. Note that the diffusive nature of sediment transport had already been highlighted by previous experimental investigations. A partial simplification of the two coupled differential equations was achieved in terms of diffusion coefficient, which was found to be proportional to the sum of reach-scale bed slope, section-averaged bed slope, and deviatory bed slope. As a matter of fact, due to the pseudo-periodical oscillation of section-averaged bed elevation and corresponding deviation, which was experimentally verified, the diffusion coefficients practically include terms that determine a mutual compensation of mass transfer from higher to lower concentrations of sediment and vice versa. For that reason, the dependence of these diffusion coefficients on the oscillatory part of the bed elevation function was neglected, and the only dependence on the constant reach-scale bed slope was maintained. A second important result, which basically proves the consistency of the present mathematical formulation, is represented by the direction of the small-scale bed forms migration according to the local section-averaged bed curvature. Specifically, where the section-averaged bed is convex, sub-critical bed waves would move upstream; conversely, where the section-averaged bed is concave, super-critical bed waves would migrate downstream. This conclusion is fully compatible with classic hydraulics and gradually-varied flows theory. The system of coupled equations that govern section-averaged bed elevation and corresponding deviation was analytically–numerically solved for a set of macro-data belonging to a real truly braided river reach of southern Italy, assuming an initially pseudo-periodical section-averaged bed. The advection–diffusion equation governing the bed elevation deviation was solved by particle-tracking, based on the constant diffusion coefficient and the velocity values locally sampled by each sediment particle, with an initial condition represented by a central double array of stones. The solution was plotted in terms of particles’ barycenter position and corresponding histogram after a characteristic diffusive time, for five different values of the ratio of (initial) longitudinal macro-bedforms maximum period to channel maximum width. The results indicate that the smaller the ratio, the more accentuated the tendency of the river to maintain the braided bed structure, with the formation of smaller and highly concentrated sediment heaps whose number is proportional to the frequency of the oscillatory section-averaged bed profile. Conversely, when the ratio is large (specifically, larger than 1), bed dynamics tends to approach a Taylor-like dispersive process (e.g., [68]) with asymmetric tails (e.g., [69]). This result means that the initial disturbance being the same, the longitudinally smoother the equilibrium braided bed (longer macro-bedforms or bars), the more marked the diffusive nature of transport, with the displacement of particles in all directions that are only bounded by river valley side boundaries. Conversely, the longitudinally more heterogeneous the equilibrium section-averaged bed (shorter macro-bedforms or bars), the more marked the advective nature of transport, with mainly longitudinal particle displacement and sediment mass fragmentation. The physical explanation for this behavior is that lower-frequency section-averaged bed waves are associated with lower-amplitude advective particle displacements and, therefore, relatively more intense diffusion. Finally, it should be noted that the system of coupled differential equations derived in the present study may easily be utilized (eventually in its complete version (23)–(24)) for a finite-differences/elements solution, starting from the real DEM of the river bed at t = 0. As a matter of fact, it is widely recognized that human-induced factors, such as embankments, roads, bridges, and dams construction, can affect sediment transport and bed morphology of the regulated river [70,71,72], and can critically impact on location and extent of gravel bars [73], which are important natural structures for fluvial or riparian habitats and, therefore, for the protection and conservation of species of fishes, invertebrates, plants, and birds [74,75,76]. Several fish species live in braided rivers, and their connected wetlands and gravel bed grounds constitute suitable spawning areas. Additionally, in many cases, gravel bars represent an important substrate for the establishment and development of ground flora and woody vegetation and guarantee higher plant diversity [77]. Sustainable management of braided rivers should, therefore, ensure their ecological potential and biodiversity by preserving braiding structure over time [1]. The analytical–numerical solution proposed by the present study, which models gravel braided bed evolution from an equilibrium configuration to the other for imposed flow and bedload rate, may, therefore, be applied for the analysis of the impact of fluvial engineering works on river habitat and its biotic components, as well as for the planning and the monitoring of sustainable river restoration.

Author Contributions

Conceptualization: M.P. and A.D.V.; Data Curation: A.D.V.; Methodology: M.P.; Software: M.P.; Validation: A.D.V. Both authors have read and agreed to the published version of the manuscript.

Funding

This research has been supported by the MIUR PON R&I 2014–2020 Program (project MITIGO, ARS01-00964); and by projects L.IDRO.AM.BIO (Laboratory of environmental and biological hydrodynamics) U.P.B. Pannone17 and S.PO.T.S.A.B. (Stream power and solid transport in braided beds) U.P.B. De Vincenzo16—University of Basilicata.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the available data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mosley, M.P. Analysis of the effect of changing discharge on channel morphology and instream uses in a braided river, Ohau River, New Zealand. Water Resour. Res. 1982, 18, 800–812. [Google Scholar] [CrossRef]
  2. Bertoldi, W.; Ashmore, P.; Tubino, M. A method for estimating the mean bed load flux in braided rivers. Geomorphology 2009, 103, 330–340. [Google Scholar] [CrossRef]
  3. De Vincenzo, A.; Brancati, F.; Pannone, M. An experimental analysis of bed load transport in gravel-bed braided rivers with high grain Reynolds numbers. Adv. Water Resour. 2016, 94, 160–173. [Google Scholar] [CrossRef]
  4. Redolfi, M.; Tubino, M.; Bertoldi, W.; Brasington, J. Analysis of reachscale elevation distribution in braided rivers: Definition of a new morphologic indicator and estimation of mean quantities. Water Resour. Res. 2016, 52, 5951–5970. [Google Scholar] [CrossRef] [Green Version]
  5. Mirauda, D.; De Vincenzo, A.; Pannone, M. Statistical characterization of flow field structure in evolving braided gravel beds. Spat. Stat. 2019, 34, 100268–doi10. [Google Scholar] [CrossRef]
  6. Ashmore, P. Channel morphology and bed load pulses in braided, gravel-bed streams. Geogr. Ann. 1991, 73, 37–52. [Google Scholar] [CrossRef]
  7. Germanoski, D.; Schumm, S.A. Changes in Braided River Morphology Resulting from Aggradation and Degradation. J. Geol. 1993, 101, 451–466. [Google Scholar] [CrossRef]
  8. Pitlick, J.; Mueller, E.R.; Segura, C.; Cress, R.; Torizzo, M. Relation between flow, surface-layer armoring and sediment transport in gravel-bed rivers. Earth Surf. Process. Landf. 2008, 33, 1192–1209. [Google Scholar] [CrossRef]
  9. Garcia, C.; Laronne, J.B.; Sala, M. Continuous monitoring of bedload flux in a mountain gravel-bed river. Geomorphology 2000, 34, 23–31. [Google Scholar] [CrossRef]
  10. Habersack, H.; Hauer, C.; Liedermann, M.; Tritthart, M. Modelling and monitoring and management of the Austrian Danube. Proc. Water 2008, 21, 29–31. [Google Scholar]
  11. Downs, P.W.; Soar, P.J.; Taylor, A. The anatomy of effective discharge: The dynamics of coarse sediment transport revealed using continuous bedload monitoring in a gravel-bed river during a very wet year. Earth Surf. Process. Landf. 2016, 41, 147–161. [Google Scholar] [CrossRef] [Green Version]
  12. Gomez, B.; Naff, R.L.; Hubbell, D.W. Temporal variations in bedload transport rates associated with the migration of bedforms. Earth Surf. Process. Landf. 1989, 14, 135–156. [Google Scholar] [CrossRef]
  13. Singh, A.; Fienberg, K.; Jerolmack, D.J.; Marr, J.; Foufoula-Georgiou, E. Experimental evidence for statistical scaling and intermittency in sediment transport rates. J. Geophys. Res. 2009, 114, F01025. [Google Scholar] [CrossRef] [Green Version]
  14. Heyman, J.; Mettra, F.; Ma, H.B.; Ancey, C. Statistics of bedload transport over steep slopes: Separation of time scales and collective motion. Geophys. Res. Lett. 2013, 40, 128–133. [Google Scholar] [CrossRef] [Green Version]
  15. Ma, H.; Heyman, J.; Fu, X.; Mettra, F.; Ancey, C.; Parker, G. Bed load transport over a broad range of timescales: Determination of three regimes of fluctuations. J. Geophys. Res. Earth Surf. 2014, 119, 2653–2673. [Google Scholar] [CrossRef]
  16. Nikora, V.; Habersack, H.M.; Huber, T.; McEwan, I.K. On bed particle diffusion in gravel-bed flows. Water Resour. Res 2002, 38, 17/1–17/9. [Google Scholar] [CrossRef]
  17. Pyrce, R.S.; Ashmore, P.E. Particle path length distributions in meandering gravel-bed streams: Results from physical models. Earth Surf. Process. Landf. 2003, 28, 951–966. [Google Scholar] [CrossRef]
  18. Ganti, V.; Meerschaert, M.M.; Foufoula-Georgiou, E.; Viparelli, E.; Parker, G. Normal and anomalous diffusion of gravel tracer particles in rivers. J. Geophys. Res. 2010, 115. [Google Scholar] [CrossRef]
  19. Ancey, C.; Heyman, J. A microstructural approach to bed load transport: Mean behaviour and fluctuations of particle transport rates. J. Fluid Mech. 2014, 744, 129–168. [Google Scholar] [CrossRef] [Green Version]
  20. Hoey, T.B. Temporal variations in bedload transport rates and sediment storage in gravel-bed rivers. Prog. Phys. Geogr. 1992, 16, 319–338. [Google Scholar] [CrossRef]
  21. Kuhnle, R.A.; Southard, J.B. Bed load transport fluctuations in a gravel bed laboratory channel. Water Resour. Res. 1988, 24, 247–260. [Google Scholar] [CrossRef]
  22. Recking, A.; Frey, P.; Paquier, A.; Belleudy, P. An experimental investigation of mechanisms involved in bed load sheet production and migration. J. Geophys. Res. Earth Surf. 2009, 114, F03010. [Google Scholar] [CrossRef]
  23. Hoey, T.B.; Cudden, J.R.; Shvidchenko, A.B. The consequences of unsteady sediment transport in braided rivers. In Gravel-Bed Rivers V.; Mosley, M.P., Ed.; The New Zealand Hydrological Society: Thorndon, Wellington, New Zealand, 2001; pp. 121–142. [Google Scholar]
  24. Kasprak, A.; Wheaton, J.M.; Ashmore, P.E.; Hensleigh, J.W.; Peirce, S. The relationship between particle travel distance and channel morphology: Results from physical models of braided rivers. J. Geophys. Res. Earth Surf. 2015, 120, 55–74. [Google Scholar] [CrossRef]
  25. Ehrenberger, R. Direct bedload mesurements on the Danube at Vienna and their results to date. Die Wasserwirtsch. 1931, 34, 1–9. [Google Scholar]
  26. Mühlofer, L. Investigation into suspended load and bedload of the river Inn, near Kirchbichl, Tirol. Die Wasserwirtsch. 1933, 26, 1–6. [Google Scholar]
  27. Einstein, H.A. The calibration of bedload trap used in the Rhine. Schweuzerische Bauztg. 1937, 110, 29–32. [Google Scholar]
  28. Nesper, F. Results of bedload and silt movement observations on the Rhine at the Brugg Bridge. Schweuzerische Bauztg. 1937, 110, 143–148, 161–164. [Google Scholar]
  29. Emmett, W.W. The Channels and Waters of the Upper Salmon River Area, Idaho; USGS Professional Paper 870A; United States Government Printing Office: Washington, DC, USA, 1975. [CrossRef]
  30. Reid, I.; Frostick, L.E.; Layman, J.T. The incidence and nature of bedload transport during flood flows in coarse-grained alluvial channels. Earth Surf. Process. Landf. 1985, 10, 33–44. [Google Scholar] [CrossRef]
  31. Whiting, P.J.; Dietrich, W.E.; Leopold, L.B.; Drake, T.G.; Shreve, R.L. Bedload sheets in heterogeneous sediment. Geology 1988, 16, 105. [Google Scholar] [CrossRef]
  32. Cudden, J.R.; Hoey, T.B. The causes of bedload pulses in a gravel channel: The implications of bedload grain-size distributions. Earth Surf. Process. Landf. 2003, 28, 1411–1428. [Google Scholar] [CrossRef]
  33. Hoey, T.B.; Sutherland, A.J. Channel morphology and bedload pulses in braided rivers: A laboratory study. Earth Surf. Process. Landf. 1991, 16, 447–462. [Google Scholar] [CrossRef]
  34. Liedermann, M.; Tritthart, M.; Habersack, H. Particle path characteristics at the large gravel-bed river Danube: Results from a tracer study and numerical modelling. Earth Surf. Process. Landf. 2013, 38, 512–522. [Google Scholar] [CrossRef]
  35. Bolla Pittaluga, M.; Luchi, R.; Seminara, G. On the equilibrium profile of river beds. J. Geophys. Res. Earth Surf. 2014, 119, 317–332. [Google Scholar] [CrossRef]
  36. Jang, C.L.; Shimizu, Y. Numerical simulation of relatively wide, shallow channels with erodible banks. J. Hydraulic Eng. 2005, 131, 565–575. [Google Scholar] [CrossRef] [Green Version]
  37. Fahnestock, R.K. Morphology and hydrology of a glacial stream—White River, Mount Rainier, Washington. U.S. Geol. Surv. 1963, 422, 70. [Google Scholar]
  38. Schumm, S.A.; Khan, H.R. Experimental study of channel patterns. Geol. Soc. Am. Bull. 1972, 83, 1755–1770. [Google Scholar] [CrossRef]
  39. Boothroyd, J.C.; Ashley, G.M. Processes, bar morphology, and sedimentary structures on braided outwash fans, Northeastern Gulf of Alaska. In Glaciofluvial and Glaciolacustrine Sedimentation; Jopling, A.V., McDonald, C., Eds.; Special publication—Society of Economic Paleontologists and Mineralogists: Tulsa, OK, USA, 1975; pp. 193–222. [Google Scholar]
  40. Ergenzinger, P. Chaos and order: The channel geometry of gravel bed braided rivers. In Geomorphological Models; Ahnert, F., Ed.; Catena Verlag: Cremlingen-Destedt, Germany, 1987; Catena Suppl. 10; pp. 85–98. [Google Scholar]
  41. Grant, G.E. Critical flow constrains flow hydraulics in mobile-bed streams. Water Resour. Res. 1997, 33, 349–358. [Google Scholar] [CrossRef]
  42. Viparelli, M.; Pica, M. Streams in Large Alluvial Beds of High Slope, xii iahr. International Congress, Fort Collins. 1967. Available online: http://www.diia.unina.it/collana02.html (accessed on 15 December 2020).
  43. Young, W.; Davies, T. Bedload transport in braided gravel-bed river model. Earth Surf. Process. Landf. 1991, 16, 499–511. [Google Scholar] [CrossRef]
  44. Zarn, B. Einfluss der flussbettbreite auf die wechselwirkung zwischen abfluss, morphologie und geschiebetransportkapazität. In VAW Mit-Teilung 154; Vischer, D., Ed.; Laboratory of Hydraulics, Hydrology and Glaciology (VAW): ETH Zurich, Switzerland, 1997. (In German) [Google Scholar]
  45. Church, M.; Rood, K. Catalogue of alluvial river channel regime data, Natural Sciences and Engineering Research Council of Canada. Rep. Dept. Geol. Univ. Br. Columbia 1983. Available online: http://www.nced.umn.edu/ (accessed on 15 December 2020).
  46. Hey, R.D.; Thorne, C.R. Stable channels with mobile gravel beds. J. Hydraulic Eng. 1986, 112, 671–689. [Google Scholar] [CrossRef]
  47. Viparelli, C. Corsi d’acqua naturali e leggi che ne regolano il modellamento. In Quaderno n.4; Istituto di Idraulica e Costruzioni Idrauliche: Napoli, Italy, 1972. (In Italian) [Google Scholar]
  48. Meyer-Peter, E.; Müller, R. Formulas for bed-load transport. In Meeting of the International Association for Hydraulic Structures Research, 2nd ed.; International Association Hydraulic Research: Stockholm, Sweden, 1948; pp. 39–64. [Google Scholar]
  49. Hunziker, R.P. Fraktionsweiser Geschiebetransport. Mitteilungen der Versuch-sanstalt fur Wasserbau. In Hydrologie und Glaziologie, 138; Eidgenossische Technische Hochschule (ETH): Zurich, Switzerland, 1995; p. 209. (In German) [Google Scholar]
  50. Wong, M.; Parker, G. Re-analysis and correction of bed-load relation of Meyer-Peter and Muller using their own database. J. Hydraulic Eng. 2006, 132, 1161–1168. [Google Scholar] [CrossRef] [Green Version]
  51. Chow, V.T. Open-Channel Hydraulics; Mc Graw-Hill: New York, NY, USA, 1959; p. 680. [Google Scholar]
  52. Parker, G.; Andrews, E.D. Sorting of bed load sediment by flow in meandering bends. Water Resour. Res. 1985, 21. [Google Scholar] [CrossRef]
  53. Tealdi, S.; Camporeale, C.; Ridolfi, L. Long-term morphological river response to hydrological changes. Adv. Water Resour. 2011, 34, 1643–1655. [Google Scholar] [CrossRef]
  54. Pannone, M. Transient Hydrodynamic Dispersion in Rough Open Channels: Theoretical Analysis of Bed-Form Effects. J. Hydr. Eng. 2010, 136, 155–164. [Google Scholar] [CrossRef]
  55. Pannone, M. Longitudinal Dispersion in River Flows Characterized by Random Large-Scale Bed Irregularities: First-Order Analytical Solution. J. Hydr. Eng. 2012, 138, 400–411. [Google Scholar] [CrossRef]
  56. Pannone, M.; De Vincenzo, A.; Brancati, F. A mathematical model for the flow resistance and the related hydrodynamic dispersion induced by river dunes. J. Appl. Math. 2013, 432610. [Google Scholar] [CrossRef] [Green Version]
  57. De Saint-Venant, B. Comptes Rendus De L’académie Des Sciences; Sciences Academy of France Institute: Paris, France, 1871; pp. 3–4. [Google Scholar]
  58. Pannone, M. Effect of nonlocal transverse mixing on river flows dispersion: A numerical study. Water Resour. Res. 2010, 46, W08534. [Google Scholar] [CrossRef]
  59. Pannone, M. Predictability of tracer dilution in large open channel flows: Analytical solution for the coefficient of variation of the depth-averaged concentration. Water Resour. Res. 2014, 50, 2617–2635. [Google Scholar] [CrossRef]
  60. Pannone, M.; De Vincenzo, A. Stochastic numerical analysis of anomalous longitudinal dispersion and dilution in shallow decelerating stream flows. Stoch. Env. Res. Risk Assess. 2015, 29, 2087–2100. [Google Scholar] [CrossRef]
  61. Pannone, M.; Mirauda, D.; De Vincenzo, A.; Molino, B. Longitudinal Dispersion in Straight Open Channels: Anomalous Breakthrough Curves and First-Order Analytical Solution for the Depth-Averaged Concentration. Water 2018, 10, 478. [Google Scholar] [CrossRef] [Green Version]
  62. Leopold, L.B.; Wolman, M.G. River channel patterns: Braided, meandering, and straight. U.S. Geol. Surv. Prof. Pap. 1957, 282-B, 39–84. [Google Scholar]
  63. Schumm, S.A. Patterns of alluvial rivers. Annu. Rev. Earth Planet. Sci. 1985, 13, 5–27. [Google Scholar] [CrossRef]
  64. Rosgen, D.L. A Classification of Natural Rivers. Catena 1994, 22, 169–199. Available online: http://www.science.earthjay.com/instruction/HSU/2016_fall/GEOL_332/labs/lab_04/Rosgen_1994_Classification_Rivers.pdf (accessed on 15 December 2020). [CrossRef] [Green Version]
  65. Montgomery, D.R.; Buffington, J.M. Channel-reach morphology in mountain drainage basins. Geol. Soc. Am. Bull. 1997, 109, 596–611. [Google Scholar] [CrossRef]
  66. Church, M. Bed material transport and the morphology of alluvial rivers. Annu. Rev. Earth Planet. Sci. 2006, 34, 325–354. [Google Scholar] [CrossRef]
  67. Bialik, R.J.; Nikora, V.I.; Rowinski, P.M. 3-D Lagrangian modeling of saltating particles diffusion in turbulent water flow. Acta Geophys. 2012, 60, 1639. [Google Scholar] [CrossRef]
  68. Taylor, G.I. The dispersion of matter in turbulent flow through a pipe. Proc. R. Soc. Lond. 1954, 223, 446–468. [Google Scholar] [CrossRef]
  69. Pannone, M. On the exact analytical solution for the spatial moments of the cross-sectional average concentration in open channel flows. Water Resour. Res. 2012, 48, W08511. [Google Scholar] [CrossRef]
  70. Mirauda, D.; De Vincenzo, A.; Pannone, M. Simplified entropic model for the evaluation of suspended load concentration. Water 2018, 10, 378. [Google Scholar] [CrossRef] [Green Version]
  71. De Vincenzo, A.; Molino, A.J.; Molino, B.; Scorpio, V. Reservoir rehabilitation: The new methodological approach of Economic Environmental Defence. Int. J. Sediment Res. 2017, 32, 288–294. [Google Scholar] [CrossRef]
  72. De Vincenzo, A.; Covelli, C.; Molino, A.J.; Pannone, M.; Ciccaglione, M.; Molino, B. Long-Term Management Policies of Reservoirs: Possible Re-Use of Dredged Sediments for Coastal Nourishment. Water 2019, 11, 15. [Google Scholar] [CrossRef] [Green Version]
  73. Assani, A.A.; Petit, F. Impact of hydroelectric power releases on the morphology and sedimentology of the bed of the Warche River (Belgium). Earth Surf. Process. Landf. 2004, 29, 133–143. [Google Scholar] [CrossRef]
  74. Ock, G.; Gaeuman, D.; McSloy, J.; Kondolf, G.M. Ecological functions of restored gravel bars, the Trinity River, California. Ecol. Eng. 2015, 83, 49–60. [Google Scholar] [CrossRef]
  75. Zeng, Q.; Shi, L.; Wen, L.; Chen, J.; Duo, H.; Lei, G. Gravel Bars Can Be Critical for Biodiversity Conservation: A Case Study on Scaly-Sided Merganser in South China. PLoS ONE 2015, 10, e0127387. [Google Scholar] [CrossRef] [PubMed]
  76. Rempel, L.L.; Healey, M.; Richardson, J.R. Macroinvertebrate community structure along gradients of hydraulic and sedimentary conditions in a large gravel-bed river. Freshwat. Biol. 2000, 45, 57–73. [Google Scholar] [CrossRef]
  77. Gilvear, D.; Francis, R.; Will, N.; Gurnell, A. 26 Gravel bars: A key habitat of gravel-bed rivers for vegetation. Dev. Earth Surf. Process. 2007, 11, 677–700. [Google Scholar] [CrossRef]
Figure 1. Schematic of small-scale sand-waves propagation direction domains.
Figure 1. Schematic of small-scale sand-waves propagation direction domains.
Sustainability 13 01246 g001
Figure 2. Basento river basin.
Figure 2. Basento river basin.
Sustainability 13 01246 g002
Figure 3. Study-case reach of Basento river—from Google Earth (modified). The two transverse red lines circumscribe the river reach under investigation. Their length is respectively proportional to Bmin and Bmax. Wave-length Λ is roughly represented by bars’ average longitudinal dimension. Flow direction is from left to right.
Figure 3. Study-case reach of Basento river—from Google Earth (modified). The two transverse red lines circumscribe the river reach under investigation. Their length is respectively proportional to Bmin and Bmax. Wave-length Λ is roughly represented by bars’ average longitudinal dimension. Flow direction is from left to right.
Sustainability 13 01246 g003
Figure 4. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/8.
Figure 4. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/8.
Sustainability 13 01246 g004
Figure 5. Particles distribution at t = τd for Λ = B/8.
Figure 5. Particles distribution at t = τd for Λ = B/8.
Sustainability 13 01246 g005
Figure 6. Histogram of particles longitudinal position at t = τd for Λ = B/8.
Figure 6. Histogram of particles longitudinal position at t = τd for Λ = B/8.
Sustainability 13 01246 g006
Figure 7. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/4.
Figure 7. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/4.
Sustainability 13 01246 g007
Figure 8. Particles distribution at t = τd for Λ = B/4.
Figure 8. Particles distribution at t = τd for Λ = B/4.
Sustainability 13 01246 g008
Figure 9. Histogram of particles longitudinal position at t = τd for Λ = B/4.
Figure 9. Histogram of particles longitudinal position at t = τd for Λ = B/4.
Sustainability 13 01246 g009
Figure 10. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/2.
Figure 10. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B/2.
Sustainability 13 01246 g010
Figure 11. Particles distribution at t = τd for Λ = B/2.
Figure 11. Particles distribution at t = τd for Λ = B/2.
Sustainability 13 01246 g011
Figure 12. Histogram of particles longitudinal position at t = τd for Λ = B/2.
Figure 12. Histogram of particles longitudinal position at t = τd for Λ = B/2.
Sustainability 13 01246 g012
Figure 13. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B.
Figure 13. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = B.
Sustainability 13 01246 g013
Figure 14. Particles distribution at t = τd for Λ = B.
Figure 14. Particles distribution at t = τd for Λ = B.
Sustainability 13 01246 g014
Figure 15. Histogram of particles longitudinal position at t = τd for Λ = B.
Figure 15. Histogram of particles longitudinal position at t = τd for Λ = B.
Sustainability 13 01246 g015
Figure 16. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = 2B.
Figure 16. Velocity sampled by the sediment strip at t = t0 (full line) and t = τd (dotted line) for Λ = 2B.
Sustainability 13 01246 g016
Figure 17. Particles distribution at t = τd for Λ = 2B.
Figure 17. Particles distribution at t = τd for Λ = 2B.
Sustainability 13 01246 g017
Figure 18. Histogram of particles longitudinal position at t = τd for Λ = 2B.
Figure 18. Histogram of particles longitudinal position at t = τd for Λ = 2B.
Sustainability 13 01246 g018
Figure 19. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/8.
Figure 19. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/8.
Sustainability 13 01246 g019
Figure 20. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/4.
Figure 20. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/4.
Sustainability 13 01246 g020
Figure 21. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/2.
Figure 21. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B/2.
Sustainability 13 01246 g021
Figure 22. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B.
Figure 22. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = B.
Sustainability 13 01246 g022
Figure 23. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = 2B.
Figure 23. De-trended section-averaged bed profile at t = t0 (full line) and t = τd (dotted line) for Λ = 2B.
Sustainability 13 01246 g023
Table 1. Summary of the experimental laboratory conditions for the derivation of the equilibrium relationship between dimensionless bedload rate and stream power [3].
Table 1. Summary of the experimental laboratory conditions for the derivation of the equilibrium relationship between dimensionless bedload rate and stream power [3].
RunQ (m3/s)S0S0SBSeqteq (h)
R10.04440.0140.0840.00746
R20.0310.0130.0670.012167
R30.0360.01350.0670.011168
R40.02540.0270.05660.015261
R50.03040.0140.0410.015208
R60.0400.010.03310.01296
R70.0450.01920.04910.009144
R80.0480.010.03310.011197
R90.0530.01040.0320.01261
R100.0350.0220.01450.01572
Table 2. Summary of numerical simulations characteristic parameters. In red the output (Braiding Index).
Table 2. Summary of numerical simulations characteristic parameters. In red the output (Braiding Index).
RunΛ/Bτd (s)PeBI
R11/832.74210.528.8
R21/4130.96105.334.5
R31/2523.8352.672.8
R412095.3426.472.3
R528381.3513.241.8
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pannone, M.; Vincenzo, A.D. Theoretical Investigation of Equilibrium Dynamics in Braided Gravel Beds for the Preservation of a Sustainable Fluvial Environment. Sustainability 2021, 13, 1246. https://doi.org/10.3390/su13031246

AMA Style

Pannone M, Vincenzo AD. Theoretical Investigation of Equilibrium Dynamics in Braided Gravel Beds for the Preservation of a Sustainable Fluvial Environment. Sustainability. 2021; 13(3):1246. https://doi.org/10.3390/su13031246

Chicago/Turabian Style

Pannone, Marilena, and Annamaria De Vincenzo. 2021. "Theoretical Investigation of Equilibrium Dynamics in Braided Gravel Beds for the Preservation of a Sustainable Fluvial Environment" Sustainability 13, no. 3: 1246. https://doi.org/10.3390/su13031246

APA Style

Pannone, M., & Vincenzo, A. D. (2021). Theoretical Investigation of Equilibrium Dynamics in Braided Gravel Beds for the Preservation of a Sustainable Fluvial Environment. Sustainability, 13(3), 1246. https://doi.org/10.3390/su13031246

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