A Shear Reynolds Number-Based Classiﬁcation Method of the Nonuniform Bed Load Transport

: The aim of this study is to introduce a novel method which can separate sand- or gravel-dominated bed load transport in rivers with mixed-size bed material. When dealing with large rivers with complex hydrodynamics and morphodynamics, the bed load transport modes can indicate strong variation even locally, which requires a suitable approach to estimate the locally unique behavior of the sediment transport. However, the literature offers only few studies regarding this issue, and they are concerned with uniform bed load. In order to partly ﬁll this gap, we suggest here a decision criteria which utilizes the shear Reynolds number. The method was veriﬁed with data from ﬁeld and laboratory measurements, both performed at nonuniform bed material compositions. The comparative assessment of the results show that the shear Reynolds number-based method operates more reliably than the Shields–Parker diagram and it is expected to predict the sand or gravel transport domination with a <5% uncertainty. The results contribute to the improvement of numerical sediment transport modeling as well as to the ﬁeld implementation of bed load transport measurements.


Introduction
In terms of typical geomorphological features, rivers can be divided into three main section types: upper, middle and lower course [1]. One of the decisive differences among these river features is the erosion capacity of the flow, which constantly decreases along the river. In fact, it can be stated that the erosion prevails along the upper course rivers, yielding coarser gravel bed material and significant bed load transport. As to the middle course, all the three characteristic sediment transport processes, such as the local erosion of the bed, the transport of coarse and fine particles together with river bed aggradation take place. Finally, the lower course type rivers can be generally characterized with the deposition of fine sediments transported from upstream. Due to the spatially and temporally varying bed load transport capacity of the flow and the changing river planforms (straight, meandering, or braided river channel patterns [2]), the morphological features can also show significant variability.
Many classification methods can be found in the literature to describe the morphological properties of river sections (e.g., [3][4][5][6][7]). The main goal of these tools is to typify and predict the dominant morphological properties and processes, such as defining the predominant sediment transport mode (bed load, or suspended load) [8], the channel pattern type [1,9], the specific bed material grain size [7], the erosion capacity [10], or the bed armor measure [6]. According to these methods, the morphological properties and processes can be estimated as the function of a few hydro-morphological variables, e.g., the water depth, the channel width, the mean flow velocity, the mean grain size of the bed material, or the longitudinal bed slope.
For instance, there are classification methods, suggesting that the dominant grain size of the bed load (sand or gravel) can be obtained based on the d 50 grain size of the bed material [8]. In most cases, such a classification is straightforward due to the fact that a dominant bed material fraction can be chosen, considering, e.g., sand-bed (d 50 < 2 mm) or gravel-bed (d 50 > 2 mm) streams [8]. However, there are several situations, when such a clear distinction cannot be made as both sand-and gravel-dominated zones appear even in shorter river reaches. Moreover, the locally dominant bed load fraction can alter according to the current flow conditions. So, only the motion of finer particles of a nonuniform bed material can occur during low-water regime. In turn, the higher flood waves increase the sediment transport capacity of the flow, resulting in probably gravel-dominated bed load. That is, because the bed load is mainly fed by the bed itself [11], very different grain sizes can characterize the bed load transport within river sections with mixed-size bed material. These sections are typical to middle-course rivers, at the transition sections in terms of bed material.

Case Study
Such a river reach is the upper Hungarian Danube reach between rkm 1798 and rkm 1795 ( Figure 1). At this section, the middle course resulted in a meandering pattern type, with nonuniform bed material (d 50 ranges between 0.32 and 70.5 mm) [9,12]. Furthermore, conventional river regulation measures (e.g., groin fields) were installed, which further enhanced the diversity of the morphological properties and processes [13][14][15][16][17]. The following parameters characterize the river: the main channel width at mean water regime ranges between 150 m and 350 m; the average water surface gradient is around 0.0002-0.00025 [17], whereas the characteristic flow discharges are Q m = 2000 m 3 /s (mean flow), Q bf = 4300-4500 m 3 /s (range of bankfull discharge), Q 100 = 10,400 m 3 /s (100 year flood event) [18] and Q 200 =~10,800 m 3 /s (200 year flood event) [19]. Here, the dominant fraction in the bed material, and hence in the bed load, shows a strongly varying spatial distribution. That is, bed armoring and gravel bed load takes place in the main channel, while clear sand transport can be observed in the near-bank zones and in the groin fields. Finally, the accumulation and erosion of the gravel-sand mixture is detectable at some places, resulting in the formation of side bars ( Figure 1). The herein-presented study focuses on this issue: how the dominant fraction of the bed load transport can be determined, when both the river hydrodynamics and the morphological features show strong heterogeneity within the studied reach. When local or reach-scale morphodynamic analysis is carried out, it is far from evident what sort of empirical sediment transport formula to use. There can be found formulas for finer fractions (e.g., [20]), as well as for coarser fractions (e.g., [21]), both having limitations when considering mixed-size bed material. However, no such straightforward method can be found in the literature, which suggests the locally dominant bed load transport mode and thus the more suitable formula for nonuniform bed material. An attempt is made in this study to fill this gap, providing a methodology that works more accurately compared to earlier published methods.

Problem Statement
Significant variability of the dominant grain size in the bed load can take place in nonuniform bed material, e.g., in the case of the investigated Danube reach. Figure 2 shows the density function of a bimodal bed material sample, collected from the lower side bar (see Figure 1, green columns, left vertical axis). The black line (and the right vertical axis) represents the critical dimensionless bed shear stress values (or critical Shields number, Equation (1)), as the function of the grain size, calculated by the Shields curve [22]. If the dimensionless bed shear stress exceeds the critical value, the given particle is in motion, otherwise it stays calm.
where τ* is the Shields number, τ is the bed shear stress, g is the acceleration due to gravity, ρ s is the sediment density, ρ w is the water density and d is the grain diameter. As the fraction rates indicate, sand and gravel also appear well separately, both with significant amount (>30%). The grain size distribution indicates a d 50 of 7 mm, which would refer to gravel transport. However, examining the critical Shields numbers, it turns out that the dominant bed load grain size seems to be inconstant, varying depending on the flow conditions. If the dimensionless bed shear stress is in a lower range (τ * <~0.05, or equivalently τ <~2 N/m 2 ), only sand (<2 mm) movement can be expected. In turn, increasing dimensionless bed shear stress (τ * > 0.05, or equivalently τ >~2 N/m 2 ) results in the motion of even the coarser gravels (>2 mm), while the motion of the finer sand fraction tends to be rather suspended, like rolling [8], resulting in even coarser bed load. Concluding the above, the d 50 grain size of the bed material itself is not able to indicate the current dominant grain size in the bed load, especially not in the case of nonuniform bed material under varying flow conditions (for further similar study cases see, e.g., [23,24]).
When implementing computational morphodynamic models in cases like the above-presented mixed-size bed material, the variable dominant fraction in the bed load can lead to inappropriate application of sediment transport models. Nevertheless, the numerical simulation of morphodynamic processes can contribute to the better performance of river engineering activities, e.g., when planning restoration measures, and it is considered a widely applied investigation tool with many benefits. However, the reliable numerical description of the interaction between river hydrodynamics and the morphological features is a challenging task in rivers with uniform bed material composition, but shows even higher difficulties in the case of mixed-size bed material conditions (e.g., [25][26][27][28][29]). When the grain sizes cover a smaller range in the bed material, the potential erosion results also in a quite uniform grain size range of the bed load transport. For such conditions, many bed load sediment transport models have already been developed and successfully validated [22,[30][31][32][33]. On the other hand, the choice of the bed load sediment transport formula is still a challenge when the characteristic grain sizes can significantly change both in time and space [23,25,[34][35][36]. Accordingly, the most suitable numerical formula can change in time and space. For instance, as many studies present [34,[37][38][39][40][41][42][43][44], the van Rijn model [20] could provide better estimation of sediment transport in sand-dominated zones, whereas the Wilcock and Crowe model [21] suits gravel-dominated mixed-size cases better. Török et al. [34] also proved that the combined and parallel application of these two models provides better results when simulating sediment transport in mixed-size bed material. In the referred study, a physical indicator was defined as a parameter to decide which formula is to be applied. It is, however, not evident what this parameter should be. In such cases, the prediction precision of the dominant fraction in the bed load transport plays an emphatic role in the accuracy of the sediment transport calculation possibilities.

Existing Method for Bed Load Classification
One applicable method, the Shields-Parker river sedimentation diagram, can be used as a comprehensive bed-material and sediment transport nature classification tool for a particular river reach, with uniform bed content [7]. The diagram indicates the Shields number (τ*, or equivalently τ* bf , Equation (1)) based on the d 50 specific mean grain size versus the so-called explicit particle Reynolds number (Re p , Equation (2), where d = d 50 in the case of the Shields-Parker diagram). The diagram includes the well-known Shields curve (continuous black curve), which shows if the given particle is in motion (points above the curve), or remains still on the bed surface (points below the curve) [22].
where R = (ρ s − ρ W )/ρ W , ρ W is the water, ρ s is the sediment density, and ν is the kinematic viscosity of the fluid. Parker stated [8] that alluvial rivers can usually be divided into two types based on the function of particle Reynolds number (Re p , Equation (2)): sand-bed and gravel-bed rivers. According to Parker, a distinction can be defined based on Re p , which actually depends only on d 50 (therefore it can, in fact, be considered as the dimensionless substitution of the grain size). That is, considering that d = 2 mm is the threshold value between the sand and gravel fractions, the critical value for Re p also can be calculated according to Equation (1), where d = 2 mm results in Re p = 360 (Re p > 360 → gravel bed, Re p < 360 → sand-bed), which critical value is depicted by a vertical narrow in Figure 3, within the Re p range of 100-1000. Consequently, the method basically yields a grain size-based classification method for near-uniform bed materials. This simplified approach is not capable of classifying characteristic local morphodynamic processes for complex situations in terms of river bed morphology and inhomogeneous bed material composition [7] (e.g., in the case of bed material in Figure 2).
The basic idea here is to keep the simplicity of the classification method shown by Parker, but at the same time to involve a new indicator, which accounts for the local hydrodynamic effects that better characterize the interaction between the water flow and the river bed, i.e., the bed shear stress. Thus, based on literature review, field and laboratory measurements, a conjecture is introduced in this paper. For a better understanding of the herein-proposed methodology, the most relevant findings of the related publications are summarized in the following.

Introducing the Novel Method for Bed Load Classification
Shields published a well-known study regarding sediment transport investigation in 1936 [22]. The research work investigated the critical state of the bed-forming particles, applying a mathematical model of forces acting on a solitary particle. Shields pointed out that the rate of the acting forces is not constant, but the relation of the effective force of the flow to the resistance of a grain can be expressed with a universal function of the shear Reynolds number (Re*, Equation (3)). The interpretation of the shear Reynolds number is expressed in Equation (3): the ratio of grain size to the thickness of the laminar boundary layer.
where u* is the shear velocity (u * = ( τ ρ ) 0.5 , where τ is the local bed shear stress). The following table demonstrates the difference between the shear Reynolds number and the explicit particle Reynolds number. The constant and the variable parameters in Equations (2) and (3) are depicted in Table 1. The last column indicates the ratio of Re * to Re p . The main outcome of the table is that Re * implicates not only the grain size (d), but also a term regarding the acting forces of the flow (u*). This remark confirms the Re * -based investigation of Shields, instead of the Re p -based method. Table 1. The constant and the variable parameters in the equations of the shear Reynolds number (Re * ) and the explicit particle Reynolds number (Re p ) (Equations (2) and (3)).

Constants 1
Variables The Shields curve, actually, shows the relation between the shear Reynolds number and the acting forces [22,45]. The investigation proves that the acting forces and the stability of a grain depend not on the grain sizes, but on the hydraulic regime formed directly above the grains [7,46,47] (Considering the possible development of different kinds of turbulent flow pattern, e.g., in sand-bed rivers bed-form-induced turbulent flow, while in gravel-bed rivers, friction-induced turbulent pattern can occur [20,32,[48][49][50][51][52]). It leads to two important consequences: (i) as the effective acting force depends not only on the grain size, but also on the acting forces of the flow, it is suggested that the sediment transport natures need to be classified based on the shear Reynolds number, too; (ii) due to this, the validation limits of the sediment transport-related mathematical models (that is, the sediment transport models) should be defined as the function of the shear Reynolds number.
The investigation of the applicability of the van Rijn bed load transport formula [20] also enhances the benefits of the shear Reynolds number-based investigation, against the grain-size-based (or characteristic grain-size-based). As one of the most widely tested and applied sediment transport formulas for river engineering applications, the van Rijn bed load transport equation [20] considers the application of the Shields curve. The transport model was developed and validated for sand transport measurements, and, therefore, the application limit was defined accordingly (d < 2 mm) [7]. However, the validity range of the base mathematical model can be defined by the shear Reynolds number, and the grain size is not a suitable variable for this reason. Namely, the drag coefficient (C D ), the lift coefficient (α L ), or the effect of the spinning motion depends on the shear Reynolds number [53][54][55][56], and the recommended values of these variables in the van Rijn model were calibrated and validated for lower shear Reynolds numbers (Re* < 400) [20,57,58]. Therefore, the determination of the applicability of the van Rijn model is reasonable as the function of the shear Reynolds number, which is the hydraulic smoother (Re* < 400) regime.
At the same time, the lower Re* range implies lower bed shear velocity (u*) and smaller particles (d) (see Equation (3)), assuming the domination of the less-resistant finer fractions in the bed load and stable coarser grains. This is proved by the about 600 measurement data [20,57,58]. Also, Parker mentions [8] that, at rougher hydraulic regime (Re* > 1000), the motion of the sand fractions is rather suspended like rolling. Consequently, finer fractions disappear from the bed load and the coarser gravels remain and become dominant. Accordingly, the van Rijn bed load model is not applicable in rougher hydraulic (Re* > 1000) regime.
The laboratory experiments, where coarser (d 50 > 5.3 mm) and nonuniform bed material conditions were investigated by Wilcock et al. [26], are also taken into consideration here. Figure 4 shows the rate of the gravel load to the sand load as the function of the shear Reynolds number. Within any mixtures, the higher the shear Reynolds number, the coarser the composition of the bed load. It is visible that the gravel load is equal to the sand load when the Re* is between 300 and 400, around 350 ( Figure 4). Above this, the gravel transport is dominant. In turn, below 350 the sand transport prevails, which underlines the earlier assumption, i.e., the fine sediment bed load transport takes place and dominates in the smoother hydraulic range. Consequently, despite that the d 50 characteristic grain size indicates only gravel transport, sand-dominated sediment transport was also found. Based on the above-described summary, the following findings are emphasized

•
Shields has already presented [22] that the shear Reynolds number (Re*) is a suitable variable (instead of the grain sizes) to define the relation of the acting forces and thus to classify the sediment transport in the near-surface zone. In other words, the effective force of the flow (together with the resistance of the grain) must be taken into account for the proper investigation of the stability and mobility of bed-forming particles.

•
The applicability of both mathematical models of Shields [22] and van Rijn [20] can be described by the shear Reynolds number, but not by grain sizes. This also supports the fact that the shear Reynolds number-based classification of the sediment transport is expected to be more reliable than a grain-size-based method.

•
The investigation of the measurements, for which the van Rijn bed load formula was validated and calibrated [20,57,58], suggests that the sand bed load can occur at lower Re* range, Re* < 400. Complementing this by the statement of Parker [8], it is assumed that in rougher hydraulic regime (Re* >~1000), the sand transport appears only in suspended form, while the bed load consists mainly of gravel.

•
The Re*-based assessment of the laboratory experiments of Wilcock et al. [26] shows that the Re* range~3-400 seems to be a critical zone, where the yield of the sand and gravel are roughly the same. A characteristic grain size cannot refer to such phenomena.
Accordingly, the following conjecture was stated: The Re* parameter is a more suitable parameter for the prediction and distinction of the sand-or gravel-dominated bed load transport, than either simply a characteristic grain size (d 90 or d 50 ), or the Re p . Moreover, Re* between~300 and 400 seems to be a critical zone: above-gravel, below-sand domination can be expected in the bed load transport.

Description of the Validation Measurements
In order to investigate accuracy of the conjecture, an assessment considering characteristic local bed grain sizes, local bed shear stress values as well as calculated Reynolds numbers was performed using datasets from recent laboratory and field experiments of the authors. First, the results from the laboratory experiments of Török et al. [59] were used. In this experiment, morphological processes (scouring, bed armoring, and aggradation) were investigated around a single groin, using mixed-size bed material with an initial d 50 of 5.16 mm. The local grain size distributions were determined by an automated image-based grain detection software tool BASEGRAIN [60]. Local bed shear stress values were estimated according to the TKE-Method [61] using the near-bed point-velocity measurements, carried out with a 3D Acoustic Doppler Velocimeter (ADV) [62]. The measurements were carried out at quasi-equilibrium bed geometry. Then, negligible bed changes were measured in the flume (∆z < ± 1 mm/h); however, weak bed load transport (Q b < 0.3 kg/h) of mainly the sand particles (d 50 = 1.3 mm) could still be observed at the outlet of the flume, where the d 50 of the bed material ranged between 4 and 6 mm [59].
Second, field data from a section of a large river was also assessed. The field experiments were carried out in the Hungarian section of the Danube River at mean flow regime, where the main characteristic morphological parameters were already introduced in the Introduction point. Forty-seven bed material samples were taken by a drag-bucket sampler, for which the grain-size distributions of the samples were determined by sieving analysis [63]. Parallel to the local bed material samplings, time-averaged flow velocity profiles were measured by an acoustic Doppler current profiler (ADCP) (WorkHorse Rio Grande 1200 kHz) at the sampling locations. Based on these long-term, fixed-boat measurements, the local bed shear stress values could be estimated using the turbulent wall law [64]. No bed load sampling could be performed during the two measurement campaigns. However, the bed load rating curve at the closest monitoring station (6 km upstream) [17] represents that low, but not negligible, bed load transport took place during both measurements (Q b < 0.1 kg/s).
The data regarding the field and laboratory measurements can be found in Appendix A (Table A1).

Measurement Data Processing for Validation Purpose of the Stated Conjecture
Field and laboratory data from the above-introduced reach of the Danube River in Hungary and laboratory experiments were conducted to confirm the stated conjecture. In order to investigate the shear Reynolds number-based classification potential, the samples were sorted into three classes according to their origin place (Figures 5 and 6): the black squares refer to the sandy bed material samples, which were collected from the near-bank zones. Here, the bed material consists predominantly of silt and sand, gravel transport can be expected only during flood waves. The red stars originate from the side bars (e.g., bed material illustrated in Figure 2). As it was introduced in Chapter 2 and based on Figure 2, the dominant grain size in the bed load is variable (both sand and gravel can dominate) and it is strongly flow-dependent. Finally, the samples from the main channel are depicted by the blue diamonds. Based on the bed material analysis, only gravel transport can occur at this region of the channel, because sand can be vanishingly found. However, as the measurements indicate stable bed surface, pure gravel transport can be expected at mean flow conditions. The investigation examines whether the shear Reynolds number is a proper variable to separate the different classes of the points.  Furthermore, the comparison of the shear and particle Reynolds number-based investigations was also examined. First, the scattering of the three point classes as the function of Re* and Re p were compared. Based on this comparison, a more accurate separation method becomes outlined visually.
Finally, in order to perform a quantitative assessment on the performance of the Re*and Re p -based approaches, a statistical method was applied, which provides information on the representativity of the Reynolds number ranges for the characteristic bed load transport processes. Here, the log-normal distributions of the Re* and Re p values were calculated, separately for the three above-distinguished groups. The reason for using the log-normal distribution was the fact that asymmetric distributions were expected (based on the scattering of the points), where Re numbers cannot take negative values [65]. The probability density and the cumulative distribution functions are described in the Appendix A. Based on the cumulative distribution functions, the Re*and Re p -based separatenesses of the different point groups become qualitatively characterizable.

Results
First, the field measurement data were analyzed in Figure 5, separately. The scattering of the points reflects a clear tendency: the bed shear stress increases with higher grain sizes. Nevertheless, when indicating a separation line, which represents a Re* of 350, the points from the three well-distinguishable regions indeed show different features. Points representing the shallower, sandy regions fall below the line. As most of these sampling locations can be characterized with mean grain sizes in the gravel range (d 50 > 2 mm), the Shields-Parker diagram would consequently suggest clear gravel transport. However, considering the relatively low bed shear stress values (<2 N/m 2 ), the Shield curve predicts only sand motion (see the findings based on Figure 2, where τ * = 0.05 → τ ≈ 2 N/m 2 ). Indeed, the d 50 of the bed material samples at the outlet of the flume (4-6 mm) would refer to gravel domination in the bed load, but the corresponding d 50 of the bed load samples show clear sand domination (1.3 mm) [59].
Second, the points distributed along the line belong to the bars. At this part of the points, the domination of sand or gravel cannot be well distinguished, but most probably a similar amount of them is transported. Accordingly, two phenomena can occur: (i) the bed shear stress values are low and the bed materials are quite uniform and coarse, resulting in resting particles (red stars, τ <~1 N/m 2 ); (ii) higher flow transport capacity can cause almost the same order of magnitude in the sand and gravel discharge, if the bed material contains both fractions (red stars, where τ >~2 N/m 2 and see also Figure 4).
Finally, the blue diamonds located above the line represent the stable [17,59,64] main stream with coarse grains. Because of the armored bed surface, weak sediment transport occurs during mean flow discharge. Only higher flood waves can result in significant bed load transport [64], containing of local bed material, i.e., only gravel grains.
The constant particle Reynolds number of 360 values are indicated with a vertical dotted line. Unlike the separation line of Re* = 350, the particle Reynolds number does not distinguish the different point classes.
Assessing the relationship between local Shields number and the shear Reynolds number for both the laboratory and field experiments ( Figure 6A), the statements in the previous paragraphs are corroborated regarding the Re*-based classification of the landforms. Again, the Re* ≈ 300-400 range outlines as a suitable indicator for the identification whether the gravel or the sand particles dominate.
Contrarily, no clear distinction of the dominating morphodynamic processes can be made if the particle Reynolds number is considered ( Figure 6B). Most of the points, representing all three different transport modes, range between 1000 < Re p < 10,000, clearly below the critical condition provided by Shields. Accordingly, a stable bed surface is expected both in the laboratory and in the field study. However, this is only the case for the armored sections of the study sites (blue diamonds) [17,59,64], whereas the bed load measurements of the laboratory experiments and the bed load rating curve at the study site in the Danube River demonstrate that weak bed load transport took place in both cases. Moreover, the bed load material investigation of the laboratory experiment suggested that the bed contained mainly sand fractions. Overall, despite the fact that the Shields curve would suggest stable bed, weak and sand-dominated bed load transport was monitored during the field and laboratory measurements.
Finally, the results of the log-normal-distribution-based probability examination are presented. As a summary of the calculations, the probability data was summarized for the six log-normal distributions, considering Re* = 300, Re* = 400, and Re p = 360 as critical values. Table 2 shows that 5% of the near-bank points are in the Re* > 300 range, 60% of the gravel bar points take Re* between 300 and 400, while only 3% of the points from the main channel are expected below 400. It can also be seen that, below the critical particle Reynolds number value (d = 2 mm; Re p = 360), only the fine sediment transport nature can be expected. However, according to the log-normal distribution, 85% of the near-bank points are over 360. Thus, the Re p = 360 indicates well that below this value mainly finer sediment transport is expected. In turn, the sediment transport nature cannot be determined obviously above this border. Table 2. Probability values of the shear Reynolds number and the particle Reynolds number, calculated based on the fitted log-normal distribution for grouped data.
The probability values indicated for the three main morphodynamic processes emphasize the suitability of Re* for the determination of the locally dominant grain size range in the bed load transport. On the other hand, the application of Re p suggests no clear classification for the same categorization, when spatially varied bed material is present, such as at the investigated river reach.

Discussion
Previous studies report that the dominant sediment transport nature in rivers with uniform bed material can be reliably determined as the function of the so-called explicit particle Reynolds number (Re p ), but yields less accurate estimations in cases when the bed material composition shows strong variability even along shorter reaches.
Based on the available and accessible data, namely 70 bed material samples and related local bed shear stress values gathered both from recent own laboratory and field experiments, it could be confirmed that, instead of the utilization of the explicit particle Reynolds number (Re p ), or the application of characteristic grain sizes, e.g., d 50 or d 90 , the shear Reynolds number (Re*) is a more adequate parameter to assess the locally dominant sediment transport nature. The probability calculations indicate that the sand-dominant bed load is estimable by 95% accuracy, while the gravel majority transport by 97%. In contrast, the Re p -based former method shows significantly greater inaccuracy. Although the gravel-dominated locations are estimated precisely, the method calculates the prediction of fine sediment majority by 85% inaccuracy.
Besides the authors' own experimental data, it could also be shown that the Re*-based approach works reliably for the laboratory and field data based on which the widely used van Rijn [20,57,58] and the Wilcock and Crowe [26] bed load formulas were validated. The method is able to indicate the fine particle domination for all of the validation and calibration data of the van Rijn model correctly. Moreover, the transitional range between the sand and gravel majority becomes detectable in the case of sand-gravel mixture bed material measurements, with which the Wilcock and Crowe model was validated.
As another example, the authors assessed the datasets, consisting of 45 coupled bed shear stress -d 50 characteristic grain size value pairs from field measurements, published by Mueller et al. [66]. In that study, the median surface grain sizes of the gravel-bed streams and rivers were reported to vary between 0.027 and 0.21 m, which indicates quite coarse bed material and armored bed surface. Accordingly, the calculated Re* are consequently high, all estimated Re* values exceed 400.
Besides the promising validation results, the theory must be further proven by direct bed load field measurements. Also noted, that the local bed load content depends not only on the local bed material, but also on the arriving bed load material [67]. This obvious effect is not taken into account by the introduced method. Additionally, the benchmark measurements were carried out at weak bed load discharges. Thus, the validation needs to be extended to a wider range of the flow regime. However, the authors believe that the presented classification indicates meaningful progress in the investigation methods of the nonuniform bed load.

Conclusions
The herein-presented classification method of the locally dominant fraction in bed load transport can contribute to better implementation of different sediment transport investigation methods applied in large rivers. For instance, a well-known issue related to field sediment transport monitoring is the high uncertainty of direct bed load sampling methods. Having information on the local bed material and on the flow field, the offered Re*-based approach can suggest where and what sort of sampling techniques would be the most suitable to collect reliable sediment information.
The results can also contribute to the development of improved computational modeling tools. Indeed, instead of applying one specific empirical sediment transport formula in a simulation to calculate local bed load transport, the Re*-based approach can be utilized to distinguish between several formulas, each of them having a certain application range. Thus, the combined use of several formulas becomes possible.
Author Contributions: G.T.T. conceived the study, collected the data, conducted the data analysis, preformed the statistical analysis and drafted the manuscript. G.T.T. and S.B. discussed the results. S.B. worked on subsequent drafts of the manuscript. S.B. and J.J. participated in the design and coordination of the study.
Funding: This research was supported by the ÚNKP-17-3 and 17-4 New National Excellence Program of the Ministry of Human Capacities, which is highly appreciated, and was partly supported by MTA TKI of the Hungarian Academy of Sciences. Support of grant BME FIKP-VÍZ by EMMI is also kindly acknowledged. The authors acknowledge the funding of the OTKA FK 128429 grant.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
The probability density function for the log-normal distribution is expressed as: where the parameters are −∞ < µ < ∞ (scale parameter) and σ > 0 (shape parameter), can be obtained by the fitting of the function to the known points. The integration of the log-normal distribution curves results in the cumulative distribution function, which is: where Φ is the Laplace integral [68,69]. Table A1. Measured morphodynamic data of field [17,64] (at the introduced reach of the Danube River in Hungary) and laboratory [59] measurements and of the referenced laboratory experiments [17,59,64].