Calculation Model to Predict the Static Armor Layer Size Distribution After the Reconstruction of a Gravel River Bed

The size distribution of armor layer in mountain rivers is an important factor that affects the stability of the river bed. However, there are relatively few studies on the prediction of armor layer size distribution in the reconstruction process after the previous static armor layer becomes unstable. In response to the above challenges, this study considers the incipient probability of sediment particles as the starting point, and comprehensively considers the coupling relationship between the initial bed materials, bed structure, armor ratio, and flow intensity, using a simple calculation model for predicting the static armor layer size distribution after reconstruction of a gravel-river bed. This paper introduces the concept of critical incipient particle size D c , and considers that sediment particles smaller than D c will incipient easily, resulting in a higher probability of being washed out, whereas the incipient probability of sediment particles larger than D c (this part of the large-particle sediment includes not only the original particles on the bed surface, but also the large sediment particles exposed by the erosion of the bed subsurface) is relatively small. At the same time, this model also uses the armor ratio to reflect the impact of the bed surface structure. This study cites data from five sets of laboratory flume experiments to verify the calculation model, and the experimental results show that the model calculation results are in good agreement with the experimentally measured data, especially in predicting the median diameter D 50 of the static armor layer. Our calculation model provides theoretical guidance for the study of mountain riverbed stability, earthquake prevention and disaster reduction.


INTRODUCTION
Natural disasters in mountain rivers profoundly affect the healthy development of rivers and, consequently, human survival. A scientific understanding of the impact of natural disasters on the evolution and development of rivers is particularly important for dealing with natural disasters in mountain rivers (Xie et al., 2018;Wang et al., 2019;Wang et al., 2021). Globally, mountainous area accounts for 30% of the total land area. Mountain water disasters have occurred frequently in the past two decades. For example, secondary disasters (such as landslides, mudslides, and collapses) caused by the 2008 Wenchuan earthquake in China have had a huge impact and profoundly changed the evolution and development of rivers Dai et al. (2011);Li et al. (2014); Fan et al. (2016), and these disasters have resulted in a huge loss of life and property in the local area. Mountain rivers are characterized by high terrain, large drops, rapid floods, and wide bed composition, and these characteristics make related scientific research challenging. For gravel-bed rivers, secondary disasters caused by the earthquake greatly affected the water and sediment supply, and the bed armor layer Hassan et al. (2006); Mao et al. (2011); Wang et al. (2020a) has been in the dynamic cyclic process of formation-destruction-reformation for a long time (Vericat et al., 2006;Orrú et al., 2016). The instability and reconstruction process of the armor layer will have an important impact on river embankment projects, water diversion and sediment prevention projects, and navigation projects.
The surface layer of a mountain river bed has a wide composition, and the particle size distribution ranges from fine sand to boulders. In addition, a gravel river bed is a morphology often present in mountain rivers, and a static armor layer is typically formed on the riverbed surface (Hassan et al., 2006;Wang et al., 2020b). Under the conditions of insufficient upstream sediment supply, the sorting effect of flow often causes fine-grained sediment to scour and move out, leaving larger particles of sediment on the river bed, which then form a static armor layer (Chin et al., 1994;Church et al., 1998;Parker and Sutherland, 1990;Wang et al., 2020a;Wang et al., 2020b). In the past, various scholars have studied the bed size distribution of the static armor layer formation process through flume experiments Little and Mayer, (1972); Shen and Lu, (1983) or field experiments Rovira and Núñez-Gonz); Gessler (1971) first developed a mixed sand and gravel size distribution that considered the coarsening process. His method only considered the pulsating effect of flow, and ignored the random distribution of the incipient drag force of sediment particles on the bed surface. He believed that the pulsation of the drag force follows the law of normal error; therefore, the probability that the sediment stays on the bed surface will not be incipient, and then the size distribution of the static armor layer can be obtained. On this basis, other scholars have begun to study the of static armor size distribution prediction model. Shen and Lu (1983) used three regression equations to calculate the surface size distribution of a river bed based on Gessler's sediment size distribution prediction model and the research results of Einstein's hidden coefficient Einstein, (1950) ;Gessler, (1971) from the perspective of its influence on the incipient heterogeneous sediment. Of course, some scholars believe that the scope of application of Shen and Lu's model has certain limitations. For example, some scholars assert that this method has poor accuracy for beds that include fine sand Kellerhals and Church, (1977) ;Ettema, (1984); Odgaard (1984) discussed and analyzed the above-mentioned problems. He chose a normal curve to describe the size distribution of the static armor layer, and the results were supported by both theory and experiment. Although the method is simple in form, his conclusion is that the size distribution after the formation of the armor layer is unrelated to the initial bed materials, which does not seem to fit the actual situation very well. Garde (2006) believes that a common shortcoming of previous studies is that they are all based on the model of (Little and Mayer, 1972). Moreover, he believes that the applicability of this model depends on whether bed sampling technology is sufficiently accurate. He introduced the parameter Kramer's M, and coupled it with the initial bed material and the flow strength to obtain a model for calculating the median particle size of the armor layer, used a wider range of laboratory flume data and channel data (median particle size range is 0.8-18.12 mm, and geometric standard deviation range is 1.5-15.76) for verification, and was able to obtain good results. He et al. (2002) believed that when the incipient probability was determined, the bed load had a maximum incipient particle size. When calculating the stable armor layer size distribution, the maximum incipient particle size of the bed surface must first be determined. Under the constraint of a certain incipient probability, sediment particles larger than this size cannot be moved. On this basis, a calculation model is proposed for the static armor layer under the condition of no upstream sediment supply. Although there are already many results that can be applied to the prediction of the bed surface size distribution of the static armor layer in the formation stage, the influencing factors that affect the static armor layer size distribution are complicated, such as exposure He et al. (2002); Bai et al. (2013); Zhang et al. (2016), riverbed structure Powell et al. (2016); Bertin and Friedrich, (2018); Wang et al. (2020b), cyclical changes in current intensity (Hassan et al., 2006;Mao, 2012). Although many achievements have been made at present, the scope of application of these results has certain limitations. For example, the results of Gessler, Shen and Liu, etc. are not suitable for a static armor layer of fine sediment, and Garde's model (Grade et al., 2006) pays more attention to the median diameter of the armor layer. The results reported by He et al.(2002) can be used to accurately calculate the bed surface size distribution when the armor layer is formed under the condition of clear water scour, but this method is not suitable for predicting the bed surface size distribution when it forms after the previous static armor layer is broken.
To improve the relatively immature calculation model for the surface layer size distribution of the new static armor layer that forms after the previous static armor layer becomes unstable, this study proposes a new model that is based on the idea of incipient probability and the critical incipient particle size, introduces the armor ratio and other parameters, reconstructs the calculation model of the bed surface bed size distribution when the new static armor layer is formed again after the previous armor layer becomes unstable, and verifies it using laboratory flume experimental data.

DATA USED FOR ANALYSIS
Experiments were conducted in a 37 m long and 1 m wide flume ( Figure 1A) at the State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, China. The initial bed materials comprised a unimodal gravel with a median grain size D 50 4 mm, and distribution truncated at 1 and 16 mm ( Figures 1B,C). Figure 1B shows a partial picture of the initial bed materials, and Figure 1C shows the size distribution curve.
During the experiment, the water level was measured using four automatic water level meters set at distances of 13, 17, 21, and 25 m from the starting point of the flume, which can measure the change in water level real time. An automatic bed load and sediment transport rate measurement system was set at a position 34 m away from the starting point; this system can continuously measure the bed load during the experiment and measure the cumulative weight of the bed load every 1 s. During the test, the topography of riverbed was measured without water during the necessary period. The topography was measured using a total station (Nikon-DTM-352c). The riverbeds at 17, 21, and 25 m were manually sampled and screened (a 20 cm × 20 cm square was sprayed with red spray paint on the target section, and colored particles were separated from the bed surface. Specifically, tweezers were used to remove the coarse-grained sediment, and then a fine-bristle brush was used to remove fine-grained sand to obtain the surface gradation information of the river bed. Herein, six experiments were conducted, which contained the formation of a static armor layer and the reformation of a new static armor layer after the previous one was broken (Wang et al., 2020a). In each case, the initial bed materials and bed slopes were the same. The experiments differed in terms of discharge (flow intensity). Here, we used data from five experiments to verify our calculation model, and the hydraulic conditions of each set of laboratory flume experiments used are shown in Table 1. The entire suite of experiments has been discussed by (Wang et al., 2020a). The grain size distributions obtained through this method are weight-by-area samples, and they require conversion to weight-by-volume samples to directly compare the grain sizes using the Bunte and Abt (2001) conversion. Wang et al. (2020a) reported a study on the critical breakup condition of a static armor layer, and concluded that when the bed shear stress increases to 1.2 times that of the bed shear stress of the prearmored bed, the previous armor layer broke. In our experiments, we used the bed shear stress τ (τ ρgRJ, where R is the hydraulic   The experimental data for this study were selected from Wang et al. (2020a).
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 660216 radius and J is the energy slope) to express the flow intensity, in which τ form indicates the flow intensity during the static armor layer formation phase, and φ τ/τ form indicates the relative flow intensity. For natural rivers, the bed surface size distribution that determines the boundary conditions can be predicted using related models. As early as 1971, Gessler. (1971) conducted related studies on the surface size distribution of beds that determined flow conditions based on probability statistics. The calculation formula is as follows: where, p 0i is the sediment weight percentage of i-the sediment class in the bed materials, D i is the particle size of i-the sediment class, and D min and D max are the minimum and maximum particle sizes of sediment in the bed materials,respectively; F i is the sediment weight percentage of the bed materials with a particle size smaller than D i ; q Di is the probability that a particle with a size of remains immobile ; τ and τ c represent the current shear stress and critical flow shear stress of i-the sediment class respectively, σ is the mean square deviation of the instantaneous flow velocity; and x is the integral variable. He et al. (2002) believed that when the incipient probability was determined, the bed load had the largest particle size D cmax . When calculating the stable static armor layer size distribution, the maximum incipient particle size of the bed surface must first be determined. Under a certain incipient probability constraint, a sediment larger than this size cannot be moved. The bed load motion includes different forms such as sliding, rolling, and jumping, but rolling is the main form of motion, so the incipient critical shear stress Θ c obtained from rolling balance is expressed as follows: where, B is the comprehensive coefficient; ρ s and ρ are the sediment and fluid densities, respectively; g is the gravity acceleration; ξ is the exposure factor, which can be calculated by an empirical formula Misri et al. (1984); Samaga et al. (1986); D m is the mean particle size; β is the incipient standard; α is the relationship coefficient between the flow friction velocity and the average flow velocity; L W is the force arm of the gravity of the sediment particles under water; C d and L d are the drag force andarm coefficient, respectively; and C L and L L are the upward lift force and moment arm coefficient, respectively. From Eq. 3, it can be concluded that the relative exposure coefficient ξ → 0 at the time when the critical incipient shear stress and the exposure of the bed surface are maximum. According to the weak motion standard, when the bed surface sediment particles are maximally exposed, the incipient critical shear stress is 0.024. The maximum incipient particle size of the sediment particles is calculated as follows: Based on the above theory, He et al. (2002) obtained a calculation model for the surface size distribution of a static armor layer under no upstream sediment supply: where, p i is the percentage of particles D i in the static armor layer, m is the maximum incipient particle size in the initial bed material size distribution with the number of groups corresponding to the smallest particle size being one and the number of groups corresponding to D cmax , being 1 ≤ m < n, where n is the total number of groups of the sediment size distribution. The two models described above can provide better results for predicting the surface size distribution of the static armor layer formed by the initial bed materials. In particular, the model proposed by He et al. (2002)has the advantages of simple calculation and strong practicability. Figure 2 shows a comparison between the predicted and measured results of the bed surface size distribution after the initial static armor layer is formed in Experiments 3-0-6-0. As shown in the figure , the predicted and measured curves are in good agreement. Our data show that the relative error in predicting the median diameter D 50 of the riverbed was less than 14%.
As mentioned above, although Gessler (1971) pioneered a model for calculating the bed surface size distribution of the static armor layer, and the method has simple calculations, clear parameters, and good generalization, there are some obvious shortcomings. First, regardless of the strength of the flow intensity , different sediment sizes will remain on the river bed, which is not in line with the actual situation. Second, the calculation of incipient shear stress uses uniform sediment, which is inconsistent with actual rivers. Finally, the calculations of incipient shear stress and of incipient probability are nested rather than independent of each other, and thus, there is great irrationality. Although Gessler's model has the abovementioned shortcomings, it provides a very good theoretical basis for calculation of the coarsening of the bed surface size distribution. He et al. (2002) also obtained a set of methods for calculating the bed surface size distribution based on incipient probability, but these methods also have obvious defects, especially in calculating the bed surface size distribution after the previous static armor layer becomes unstable. Their theory includes the following three hypotheses: 1) The probability of particles smaller than D cmax on the bed changes depending on the particle size, and the probability of coarser particles on the bed is higher. Therefore, the probability of not being washed out can be expressed as (D i /D cmax ) θ , where is the constant coefficient. 2) The Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 660216 composition of the bed surface also affects the retention of sediment particles on the riverbed. The larger the particle size, the greater the probability of staying on the bed, which is represented by (D m /D cmax ) θ . 3) Particles larger than D cmax on the bed surface include two parts: one is particles larger than D cmax on the initial bed material itself, and the other is particles exposed on the subsurface layer. The above three hypotheses are very important for the prediction model of the static armor layer formation process, but for the new process of forming a new static armor layer after the previous one is destroyed, it is not sufficient to rely solely on these three hypotheses, and other influencing factors need to be considered.

CALCULATION MODEL OF THE STATIC ARMOR LAYER SIZE DISTRIBUTION AFTER THE RECONSTRUCTION
The current calculation method Gessler, (1971); Shen and Lu.  Bertin and Friedrich, (2018); therefore, the factors considered are not comprehensive enough. In particular, when a new static armor layer is formed after the previous armor layer is destroyed, the bed structure has a significant influence on the movement of the bed load (Vericat et al., 2006;Wang and Liu, 2009;Orrú et al., 2016;Wang et al., 2020a, Wang et al., 2020b. Therefore, it is necessary to fully consider the influence of the riverbed structure when calculating the riverbed surface size distribution after the static armor layer becomes unstable. The bed surface coarsening process is not only a simple sediment particle sorting process, but also is affected by many influencing factors such as flow intensity, bed structure and texture, sediment particle shape and density, and upstream sediment supply. The bed stability of gravel-bed rivers is primarily due to bed coarsening, which is comprised two parts: 1) bed surface coarsening and 2) bed surface structure (Wilcock, 1987;Wang et al., 2020a). To study the relationship between bed structure and bed composition more intuitively, Wolcott (1989) defined the bed structure as comprising "textural structures" (including armoring and particle interlocking) and "geometric structures" (including clusters and imbrication). At present, the technology for obtaining gravel riverbed compositions is relatively mature and includes both direct sampling methods Bunte and Abt, (2001) and image processing methods (Nie et al., 2015). According to the stacking patterns between different particles, Wittenberg. (2002) described and defined the current riverbed structure and particle structure. The arrangement of sediment particles on the bed surface can reflect the performance of the riverbed structure to a certain extent. Figure 3A shows the bed morphology after the bed was roughened. After the coarsening is completed, the sediment particles on the riverbed transition from their initial free-spreading state to a large-particle sediment stacking structure. There are some significant structural characteristics in the bed, such as cluster structures (red dashed frame in Figure 3B) and linear structures (blue dashed frame in Figure 3B). These bed structures will significantly increase the stability of the bed, have a huge impact on the bed load transport rate, and then affect changes in the surface composition of the bed. In natural rivers, when the shear stress τ of the flow exceeds the critical shear stress τ c , sediment transport occurs. Eq. 8 can be used to pre-calculate the bed load transport rate (Dietrich et al., 1989): where, q b is the bed load transport rate of single-width by weight, and k and c are empirical coefficients (where c is approximately 1.5). Dietrich et al. (1989) obtained a dimensionless sediment transport ratio q p , which is the transport rate for the coarse surface normalized by the transport rate for a surface as fine as the subsurface or load: where, τ cs and τ ct are the critical boundary shear stresses of the surface and the subsurface, respectively; a is gravel with a uniform specific gravity of unity; and D 50 and D 50l are the median grain size of the bed surface and load, respectively. To estimate the influence of bed structure on the process of particle entrainment and sediment transport, Hassan and Church. (2000) established a set of methods for estimating the corresponding shear stress of the bed structure based on the total shear stress, which has been widely used. Their method quantifies the bed structure from a unique perspective, which also provides us with a comprehensive perspective for understanding the influence of the bed structure on the stability of the static armor layer. In this method, the parameter D 50 /D 50t , which is the armor ratio, is defined as the ratio of the surface median diameter to the subsurface median diameter (Hassan et al., 2006). Eq. 9 can be rewritten in the following form: where, τ cs 0.045g(ρ s − ρ)D 50 , and g is the acceleration due to gravity. The Shield number is close to 0.045 when dealing with a widely graded sediment mixture, so it can be approximated as τ cs /τ ct D 50 /D 50s . When q p 0, there is no sediment supply upstream, which is part of the process of static armor layer formation. Figure 4 shows the loci for various q p and experimental data form the experiments (Church et al., 1998;Hassan and Church, 2000;Wang et al., 2020b). After the formation of the static armor layer, the bed load movement intensity was low, the mobilization of the river bed was small, and the river bed was in a stable state. In this analysis, our experimental data were concentrated around q p 0.1. The degree of riverbed coarsening was positively correlated with τ/τ ct . The higher the degree of riverbed coarsening, the stronger the stability of the riverbed, and the greater the influence of the bed structure. Especially in the pebble river channel, which had no upstream sediment supply, the "stone cells" formed on the bed surface will greatly promote the stability of the river bed (Church et al., 1998). As the current mathematical quantification method of bed structure is still difficult to apply, in the previous static armor layer size distribution calculation models, few scholars have fully considered the influence of the bed structure. Based on the above discussion, we believe that the riverbed armor ratio can be used as a reflection of the riverbed structure and, consequently, it was introduced into our model.
This study is also based on the incipient probability, by comprehensively considering the immobility probability of the large sediment particles in the riverbed and combining the influence of the bed structure on the coarsening of the gravel riverbed to establish a new calculation model. To reasonably calculate the bed surface size distribution of the static armor layer that formed afterinstability of the previous armor layer, it is necessary to determine the criterion for the large-grained sediment on the bed surface, which is like the study of He et al. (2002). However, our understanding of the critical particle size D c is different from that of He et al. (2002). They believe that sediment particles larger than the maximum incipient particle size on the river bed cannot be moved, and we believe that the sediment on the bed surface should be divided into particles that are larger and smaller than the critical size D c . Out of the two particle size groups, the sediment particles larger than the critical size have a small incipient probability, while sediment particles smaller than the critical size have a higher incipient probability. Therefore, the choice of critical particle size D c has become an important factor affecting the accuracy of the model's prediction. The calculation of the critical incipient shear stress obtained with bed load rolling as the main form of motion is shown in Eq. 3 and Eq. 4. It can be shown that this parameter is related to the flow intensity (using shear stress τ to characterize flow strength) and the particle size on the bed surface. 0.045 or 0.047 can be used when dealing with wide-graded nonuniform Shields numbers (Andrews, 1984;Wilcock and McArdell, 1993;Hassan and Church, 2000). In this study, we continue to use this standard, and the critical dimensionless Shields number is 0.047. From this, we determined that the critical incipient particle size in the model is D c τ/[0.047g(ρ s − ρ)]. Figure 5 shows the calculation results of the bed surface size distribution of the new static armor layer formed after instability of the previous armor layer with different critical particle sizes in Experiment 2. It can be seen from the figure that when the critical dimensionless Shields number Θ c 0.047, the bed size distribution obtained is closer to the measured result than Θ c 0.024, FIGURE 4 | Plot of Eq. 10, showing the correlation between D 50 /D 50s and τ b /τ ct for various q p . The solid circles represent the data of Exp. HM1-HM6 Hassan and Church, (2000), the open five-point stars represent HM7-HM9 Hassan and Church, (2000), and the crosses represent data from (Church et al., 1998). The open circles represent values of q p computed from Eq. 10 for our data.
FIGURE 5 | The calculation results of the bed surface size distribution of the new static armor layer formed after instability of the previous armor layer with different critical particle sizes in Experiment two.
Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 660216 which also indicates that the value of 0.047 for this parameter is relatively accurate. Moreover, the value of θ in the prediction model of He et al. (2002) does not significantly consider the influence of the riverbed structure; in this model, if the linear relationship is set to 1, the non-linear relationship is set to 0.5. However, we believe that the value of this parameter will affect the prediction accuracy of the model, and we believe that the bed structure has a greater impact on the bed surface layer after the formation of the initial static armor layer. To address this, we have added a parameter known as "armor ratio" that reflects the influence of the bed structure in the model, so when considering the parameter θ, its value can be increased appropriately. This is more in line with the actual situation. After increasing the coefficient θ from 0.5 to 0.6, the predicted configuration was found to be closer to the measured size distribution.
We believe that the main factors affecting the movement of sediment on a river bed include flow conditions, characteristics of the bed materials (composition of the bed surface and subsurface Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 660216 8 layers), and bed structure. The flow intensity is the main driving force affecting the sorting of the riverbed. For the same bed materials, the greater the flow intensity, the surface particles and the greater the corresponding D c . Therefore, the influence of flow intensity can be determined by the value of D c .The armor ratio can reflect the influence of subsurface bed materials. The bed surface sediment also affects the size distribution calculation, which can be reflected by the ratio of D i to D c and is related to thearmor ratio. After the bed materials were scoured by flow, the sediment particles formed stacked structures of different sizes, such as linear structures, triangular structures, and cluster structures (see Figure 3). The influence of these structures greatly affects the calculation of bed surface size distribution. Moreover, the influence of bed structure cannot be characterized by the characteristic particle diameter of the riverbed. However, the final result of the cluster structure formed on the surface of the river bed is to make the bed surface more stable, so we can change the value of D c and add the armor ratio to reflect the influence of the bed structure. Based on the above consideration of the factors affecting the surface sediment particles of the river bed, we propose a method for calculating the percentage of the static armor surface layer as follows: The accuracy of Eqs 11, 12 were ascertained by comparison with the measured data form laboratory flume experiments. Figure 6 shows the comparison between the predictions of our model and the data measured by laboratory flume for the size distribution curve of the new static armor layer that formed after the previous static armor layer was broken. The results show that the values calculated by this model are consistent with the measured data. Figure 7 shows a comparison between our model calculations and the actual measured static armor layer surface cumulative weight percentage. It can be seen from the figure that our predictive model still has significant deficiencies in the prediction of fine-grained sediment (the black solid squares in Figure 7, D i 2 mm), but it has a good ability to predict the proportion of larger particles (open squares and circles in Figure 7, D i 4 and 8 mm). The median diameter of the bed used in the static armor layer formation and breakup experiments was close to 8 mm, which means that our model is effective in calculating and grading the composition of the bed surface, especially in predicting the median diameter D 50 ; however, there are still some errors in the prediction of small sediment particles. We believe that this is due to the difficulty in sampling fine particles during the sampling process. Large errors will occur in the removal of fine particles from the river bed by fine brushes, which will lead to large deviations between the measured data and the calculated results.

CONCLUSION
The influence of the bed structure on the bed surface size distribution during the formation and destruction of the static armor layer was explored, and on this basis, a model that calculated the bed surface size distribution after the previous static armor layer became unstable was established and verified. The development degree of the bed structure has an impact on the formation and reconstruction of the static armor layer. After the previous static armor layer formed, the bed surface formed a structure with a different degree of particle Frontiers in Earth Science | www.frontiersin.org May 2021 | Volume 9 | Article 660216 aggregation, which increased the bed stability. The degree of influence of the bed structure on the bed stability can be initially obtained by stripping off the two factors that affect the stability of theriverbed (bed structure and the bed coarsening). The ratio of the median particle size between the bed surface and the subsurface D 50 /D 50t (armor ratio) was introduced to reflect the effect of the bed structure. The calculation method for the critical incipient particle size was obtained by reasonable analysis of the value of the dimensionless incipient shear stress Θ c . The index correction of the ratio of the fractional particle size to the critical particle size determines the influence of the initial bed materials on the formation process of the static armor layer. Combining the above theories and considering the incipient probability of sediment particles on the bed surface as the starting point, a simple prediction model for calculating the static armor layer size distribution after reconstruction of a gravel riverbed was established. The results of the calculation model are in good agreement with the results of the laboratory flume experiment data, especially for predicting the median particle size. Although our model still has certain errors in predicting small sediment particles, itstill provides information that can be used for the study of gravel riverbed stability, and provides theoretical guidance for the development and protection of mountain rivers.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.