Hydrodynamic Controls of Particulate Metals Partitioning along the Lower Selenga River—Main Tributary of The Lake Baikal

: In this study, the downstream e ﬀ ects of pollutants spreading due to hydromorphological gradients and associated changes in sediment transport conditions along the braided-meandering and deltaic distributary reach of a large river downstream section are discussed. We demonstrate the signiﬁcance of hydrodynamic control for sediment-associated metal partitioning along the river. Typically, the downward decline of the sediment and metals spreading towards Lake Baikal is observed due to bu ﬀ er e ﬀ ects in the delta. During peak ﬂow, the longitudinal gradients in heavy metal concentration along the distributary delta reach are neglected due to higher concentrations delivered from the upper parts of the river. In particular, signiﬁcant variations of heavy metal concentrations associated with the river depth are related to sediment concentration and ﬂow velocity proﬁles. Various particulate metal behavior in silt-sand delta channels and the sand–gravel Selenga main stem emphasize the importance of near-bottom exchange for particles spreading with the river ﬂow. Using empirically derived Rouse numbers, we found quantitative relationships between the ratio of particulate metals sorting throughout depth in a single river channel and the hydrodynamic conditions of sediment transport.


Introduction
River sediments carry tremendous quantities of chemicals and thus act as the main driver of pollution to the recipient lakes and seas [1,2]. Most of the existing studies of suspended sediment geochemistry rely on sediment samples taken at the surface of the river channel, at one single sampling time [3,4]. However, particulate element sorting may exist due to hydrodynamic effects (e.g., velocity, transport capacity) or due to composition (chemistry, mineralogy, size) of grains-both impacts are unknown. Chemical effects of hydrodynamic sorting within the river water column may significantly influence estimates of riverine fluxes and thus require applications of hydrodynamic methods [5,6].
Hydrodynamic sorting of sediments abruptly changes along the lower reaches of large rivers where a decrease of sediment and metal spreading is observed due to channel gradient decline towards the receiving water body related to buffer effects in the delta. For this study, we took the case of the main tributary of the Lake Baikal-the Selenga River, its delta, and adjacent upstream 100 km reach. Here, the shift from single-channel reach located upstream from the delta to the huge distributary system occurs, which provides an opportunity as a unique field laboratory to explore varying hydrodynamic effects on sediment transport and composition. Physical controls on spatial transitions in bed material from gravel to sand have been discussed before [7,8], demonstrating relations of distribution of sediment grain size in deltaic channels and hydrodynamics and sediment transport conditions. At the same time, the transition of sediment fluxes, including partitioning between suspended load and bedload, as well as the hydrodynamic drivers of geochemical fluxes, remains unclear.
In this paper, we investigate hydromorphological gradients and associated changes in sediment transport conditions in the braided-meandering and distributary reach of the large Selenga river downstream section. We argue that the well-known Rouse model [2] of sediment transport [2], both with the suspended particulate matter concentration profiles, predicts suspended load to bedload partitioning and can be used as a proxy to describe particulate metal transport as a function of height above the bed from a limited number of measured parameters: suspended sediment concentration at a given depth, particle size distribution, and bottom shear velocity. The particular focuses of the paper are (1) downstream variations of sediment transport conditions and partitioning between suspended load and bedload, (2) the effects of these processes into particulate concentrations of metals during a particular flood event of 2018, and (3) the behavior of chemical elements with respect to hydrodynamic sorting within the river water column along the river course.

Data and Methods
The case study is the Selenga River, originating in Mongolia, which contributes about 50% of the total inflow into Lake Baikal and is considered as a key sustainability factor of the largest freshwater lake in the world [9]. Due to the lithochemical features of the lower Selenga basin, its waters are enriched with chemical elements (Fe, Zn, Mo, Cu), with concentrations corresponding to world-average values [10]. The main sources of river pollution are associated with the poorly treated industrial and municipal wastewaters from Ulan-Ude, Selenginsk, and Kabansk cities [11], as well as from cities of Mongolia [12]. The large inland delta (over 500 km 2 ) is considered as an important geochemical buffer to protect Lake Baikal [13,14]. It is a fluvially dominated fresh-water system that is characterized by up to eight orders of distributary channels [7,[15][16][17]; as such, the delta region has developed large lakes and widespread wetland regions that are adjacent to the channels [8,13,[18][19][20]. Here, a significant part of the suspended and dissolved forms of heavy metals and metalloids, coming from various sources, are stored in the delta [13,14,18,21,22]. Previously, a decrease in concentrations of Cu, Pb, Zn in delta reaches, and Mn and Fe in bottom sediments was observed [13,23].

Sampling and Samples Processing
A hydrochemical screening campaign was carried out between 27 July and 1 August 2018. The measurement transects were distributed 150 km along the lowermost reach of the Selenga River ( Figure 1). This reach is specific for varying conditions with respect to flow hydraulics. The main course of the Selenga River in its lower reach (100-150 km from the Lake Baikal) crosses the Khamar-Daban ridge, flowing in a narrow V-shaped valley in a north direction, and rarely braiding. The top of this reach is where the main source of water pollution from domestic effluent and industrial wastewater-Ulan-Ude city-is located. Next, the river course rotates to the west and enters a vast flatland. Active braiding starts about 40 km from Lake Baikal, where the river forms a delta. The reach upstream from the delta will hereafter be referred to as "main channel" (transects from S18-17 to S18-28 in Figure 1) in contrast to "delta channels" (transects from S18-4 to S18-16 in Figure 1) and "delta edge" (transects S18-1, S18-2, and S18-3 in Figure 1). The latter transects correspond to the lowermost locations within delta channels in the interface between river and lake.
To assess the chemical concentrations and fluxes associated with the bed and suspended sediments, water sampling was carried out on 28 transects along with the simultaneous discharge measurements with an acoustic doppler current profiler (ADCP) to account for local hydraulic conditions. We used a Teledyne RDInstruments RioGrande WorkHorse 600 kHz ADCP (Teledyne RD Instruments USA, 14020 Stowe Drive, Poway, CA 92064, USA) unit mounted on a moving boat. Water samples were pumped out with a filterless submersible 12V pump from three layers (top, midsection, and near-bottom) to account for the vertical distribution of the suspended sediment. For each sample in a depth profile, the boat was repositioned at its original location, and sampling was performed while drifting at the river water velocity.
Bottom sediment samples were taken using a tethered flow-through sampler: an O-shaped metal pipe with a textile filter attached at the rear end. Three samples were taken at each transect-two under each bank and one at the deepest section.
The sampling was carried out during the peak of a flood event ( Figure 2). The maximum streamflow discharge calculated from the stage-discharge curve at the Kabansk gauging station (5 km upstream of the delta apex, see Figure 1

Hydrodynamic Modelling of Sediment Transport
The hydrodynamic conditions of sediment transport and particulate geochemistry evolution along the Lower Selenga River were estimated using based on one-dimensional momentum balance approach, which enables us to calculate bed shear stress condition [24]: where U is measured depth-averaged flow velocity in the flow direction (x), τ b f is bankfull shear stress, ρ is the density of water, g is gravitational acceleration, and ∂η ∂x is the change in bed elevation in the x-direction. Assuming incompressible, quasi-steady, one-dimensional flow ( where S f is the friction slope and Fr is the Froude number, equal to U/ gH. Shear velocity (V * ) is then equal to The precise estimates of bed shear velocity u * and shear stress τ b were done based on ADCP [25]: where u is the integrated velocity of the flow at the sampling site taken from the ADCP measurements, k is von Karman's constant assumed equal to 0.40, H the water depth at the sampling site, e is the base of natural logarithms, z 0 is bed roughness calculated as z 0 = 0.1·d 84 . Local boundary shear stress values τ b at the sampling verticals was based on the method by [26]: where ρ is the water density assumed unity, k is von Karman's constant assumed equal to 0.40, u is the near-bed velocity, z is the depth at the sampling point, and z 0 is the bottom roughness scale calculated as z 0 = 0.1·d 84 . A dimensionless Rouse number (Ro) is used to determine conditions of sediment transport and was used as a proxy to explain sediment geochemistry: where k is von Karman's constant (0.41), ω s is settling velocity (m/s) of the bottom sediments, which is a function of grain size, shape, and density) [27]. Settling velocity (ω s ) of natural sediment particles is based on methods presented by Dietrich and is estimated with the following equation: where ρ s and ρ are densities (kg m −3 ) of grain and fluid respectively, g is the gravitational acceleration (m s −2 ), v is the kinematic viscosity of the fluid and ω * is dimensionless particle settling velocity (−) which is estimated with an empirical equation: where empirical coefficients R 1 , R 2 , and R 3 are calculated as where the values of Corey shape factor CSF (−) and the Powers roundness P (−) were set constants equal to 0.7 and 3.5 (typical coarse sand), and D * is dimensionless nominal particle diameter calculated with the following formula: where D n is the nominal diameter of a particle (m). Suspended load flux was estimated for each ADCP transects as an average of suspended sediment concentration (SSC). For this estimate, we did not consider any and assumed that the water samples (near-surface, middle and near-bottom layers) capture the representative suspended load flux: Bedload fluxes (q G ) were estimated using a simplified formula for bedload transport, proposed by [28,29] and recalibrated on natural sand-bed river bedload transport data sets [30] The bedload in each cross-section was counted based on Q G as Further partitioning between bed load and sediment load was assessed as Q R Q R +Q G for each transect.

Grain Size and Hydrogeochemical Analyses
The water samples were then filtered for suspended material through a 0.45-µm membrane filter. The concentrations of elements were derived by inductively coupled plasma mass-spectrometry (ICP-MS) and atomic-emission spectroscopy (ICP-AES) methods. For the present study, we consider insoluble major elements (Fe, Al) and metals (Cd, Cr, Cu, Mo, Ni, Pb, As, Mn, Sr), which concentrations increase with decreasing grain size and hence are subject to hydrodynamic sorting. In addition, these elements are important regional pollutants [31].
Bottom sediments coarser than 1 mm were sieved manually to determine grain size. The finer sediment grain size was measured with a Fritsch Analysette 22 NanoTec Laser particle sizer (FRITSCH GmbH, Industriestrasse 8, 55743 Idar-Oberstein, Germany). All grain sizes were classified into 3 categories: clay (grain sizes d < 5 µm), silt (d = 5-50 µm), and sand (>50 µm). The average 50% (d50) and 84% (d84) sizes were also calculated for each sample. To calculate the d50 and d84 sediment particle diameters, the measured grain size distribution curves were linearized in double logarithmic coordinates, and the corresponding values were interpolated.
Particulate concentrations of metals c (µg/L) were determined based on the relative concentrations of elements within a sample C 0 (µg/g) accounting measured sediment concentrations SSC (g/m 3 ): Further suspended load to bedload S-B analyses was used to characterize metals distribution over channel depth from the bottom (B) to suspended (S) sediments. Analyses of hydrodynamic sorting of a chemical within a vertical were done based on a ratio between concentrations of element x in the subsurface layer C 0surf (µg/g) and near-bottom layer C 0bot (µg/g): Hence, Kx > 1 refers to hydrodynamic sorting of element x, whereas Kx < 1 is related to element x enrichment in the subsurface layers. We considered the sorting as homogenous when Kssc varies as 1 ± 0.1, and heterogeneous is >1 ± 0.3. The Kx was also applied to physical parameters such as velocity V, d50, SSC, and sand fractions share. This quantitative description of metals partitioning between the bed and suspended load is referred further as the suspended load to bedload analyses (SB-analyses). This approach develops DS-analyses related to partitioning between dissolved (D) and suspended (S) modes [31]. The exact quantitative descriptions of partitioning coefficients related to DS partitioning and SB partitioning are published initially in [32].

Downstream Variations of Sediment Transport
Values of V * corresponding to flow conditions during the time period of measurements presented herein (Figure 3b) demonstrate the shifted conditions of sediment transport. Both V * , τ, and Ro demonstrate a decreasing trend along the river course, with an abrupt decline along the distributary river system (Figure 3).
The maximum streamflow discharge of the Selenga River that was measured by the ADCP unit at the delta apex on 28 July reached 2200 m 3 /s. A comparison with the corresponding value at the gauging station showed only an 8% discrepancy between the gauge and the ADCP measurements. Figure 3 shows a significant alteration of hydraulic conditions between the main channel and the delta. All stream parameters rapidly decrease towards the delta edge. The smallest measured channels conveyed only 2%-4% of the main channel discharge.
The main channel (upstream of the delta), on the contrary, shows little streamflow variability-only 15% of the mean streamflow discharge (1900 m 3 /s). The mean water velocity also showed high stability-a 6% variance of the mean 1.8 m/s. The described contrast between the hydraulic conditions in the delta and the main Selenga River channel allows for a comparison of the sediment partitioning between the two reaches. Rouse number varies between maximum values within the transects S18-09 Ro = 5.09 in the upper part of the delta and Ro = 0.01 at one of the outer transects at the downstream part of the delta (transect S18-02). The values are ranged at certain parts of the river channel: they were between 0.01 and 1.04 at the vertical at delta edge transects, 0.06 and 5.21 (average Ro = 1.55) at delta channels transects, and from 0.11 to 3.84 (average Ro = 2.81) at the main channel upstream from the delta (transects from S18-17 to S18-28).
With changes in hydraulic conditions, the sediment transport mode is shifted along the three considered reaches. Assessment of the measured Q G +Q R for each transect showed the downstream decrease of the sediment flux, therefore implying that some of the sediment must be depositing. In the main channel upstream from the delta, the Q R /(Q G +Q R ) ratio variation lies between 10% and 30%, implying the dominant bedload transport conditions. In the most downstream part of the delta, the Q R /(Q G +Q R ) ratio reaches 90%-100%. This increase of suspended sediment contribution to total transport along with the decline of the Rouse number value (see Figure 3b) is associated with the general increase of sediment sorting downstream. These variable conditions of sediment transport are clearly seen in the various grain-size sorting conditions observed in the delta (Figure 4) compared to upstream from the delta. In most profiles, particle distribution with depth displays a relatively homogenous fashion. Particle coarsening with increasing depth is observed only at the delta edge where clay and silt particles are in the range of the critical sizes to be transported by the flow [8] and hence are settled. Clay particles are rather poorly sorted, especially along the upper part of the river, due to low settling velocities.

Sediment Grain Sizes and Geochemistry
The suspended sediments are poorly sorted in a range from 0.5 to 500 µm, with patterns specific to the delta channels ( Figure 5a) and in the main channel (Figure 5b), with an average of 15 µm. Even though the average during the survey remained constant along a reach, some changes in the particular classes can be seen. The grain size mode is increasing from the surface to bottom sediments. In both patterns, grain sizes of suspended sediments are characterized by two dominant classes (bimodal pattern), both in near-surface, middle, and near-bottom samples ( Figure 5). The first peak is associated with the 1-10-µm class, and the second belongs to 50-100 µm. Both in the delta and main channel, an increase of the coarser fractions (10-50 µm) is observed in the near-bottom sediments. Bed material samples generally have large mean diameters (several hundreds of µm) with a single-mode and evenly distributed along the spectrum from 1 to 100.000 µm. Finer bed material exists at the most downstream delta branches at the outlet, where due to sand and gravel termination upstream, only clay and silt particles are presented in the bottom layer (Figure 5b). Interestingly, the latter have the same grain size patterns as suspended sediments. These grain size effects are clearly seen within the relative concentrations of elements C 0 (µg/g). The abundance of coarse particles, such as quartz or albite grains in the bottom sediments, decreases other element content. This dilution is maximal in the bottom where coarse-grained quartz concentration is maximum (Table 1). This implies that bottom sediments are enriched in elements such as Si, while suspended sediments are enriched in Al, Fe, and other elements carried in small or platy minerals such as micas or clays and associated with Al concentrations. Among studied elements, only Sr is not influenced by sediment sorting.  The majority of the elements' concentrations mostly demonstrate an increasing trend along the studied reach. This might be due to higher discharge during sampling dates in the delta reach (Table 2). Under streamflow conditions of 2000-2200 m 3 /s, the concentrations of most elements in the delta branches are 2-3 times higher compared to the adjacent upstream reach surveyed at 1800-2000 m 3 /s. The most pronounced increase is related to near-bottom concentrations, which increased downstream, mostly with a factor of 3 ( Table 2). A statistically significant increase in the concentration of metals in suspended sediments from the surface to the river bottom was revealed as an average for the whole considered river reach ( Figure 6). The vertical distribution ratio (coefficient Kx, see Section 2.3) indicates that the physical characteristics of the flow are the most uniformly distributed over the depth parameters among all-flow velocity V, m/s (Kv = 1.04). Suspended sediment concentration SSC, mg/l, is characterized as heterogeneous and is lying within the range of metals concentrations value (K SSC = 1.12). The metals can be ranked by the increase of the Kx value within the vertical profile in the following sequence: Pb-As-Mn-Mo-Cr-Ni-Cu-Cd. The value of Kx varies from K Pb = 1.04 to K Cd = 1.42. These effects are different for the particular river reaches. Rather homogenous metals distribution or even decrease of concentrations (Kx < 1) over depth are observed in the upper channel reach upstream from the delta consistently with an observed increase of the bedload transport and coarsening of the near-bottom suspended sediments. In the delta where no sand or gravel transport is observed, the concentrations of metals increase in the near-bottom layer, both with increases of metal-bearing clay and silt fractions. The role of hydrodynamic sorting of grain sizes due to changing runoff conditions is clearly seen by the relationship between particulate concentrations and grain sizes ( Table 3). The distribution of the metals within the depth column is mostly due to hydrodynamic sorting of grain sizes. Most chemical elements (both major elements-Na, Mg, K, Ca, Al, Fe, and Mn-and metals) show strong relations between grain-size and chemical composition [33][34][35][36]. Under the conditions of the hydrodynamic sorting, we observed significant correlations (Rcor > 0.5, p-values < 0.001) between size fraction 0.5-1 µm for most of the elements (Table 3). Most elements show strong relations between grain-size and chemical composition (Figure 7). This can be explained both by typical enrichment of clay fractions by the metals (as have been shown by various regional studies demonstrating the crucial role of <16-µm sediments in transporting metals [37]. Also elevated correlations (Rcor > 0.60, p < 0.001, n = 33) were observed with the 50-100-µm size fraction. The latter emphasize the possible role of coarser sediments (silt-sand) due to other speciation of metals related to minerals. Examples of bimodal distribution with sediment grain sizes can be found in literature in relation to Fe-Mn oxides, which can drive variation of relative distribution percentage for each metal speciation (Cu, Ni, Pb, Cr, and Zn) at a primary binding phase in sediment particles [36].

Bed and Suspended Sediment Fluxes
The calculations indicate that sediment flux decreases downstream with an increasing bifurcation in the delta [8], which generally corresponds with the previous studies and also the hypothesis of gravel termination [20]. The latter can be responsible for the abrupt decline in bedload flux, which almost stops in the delta outlets. Shear stress and Rouse number are also reduced downstream due to water partitioning among the bifurcating channels. As one can see from the obtained datasets, the shear stress, which is not sufficient to produce gravel transport for the delta channels, also induces depositional patterns in the suspended sediment flux: the increase of near-bottom concentrations is the evidence of downward fluxes of the most suspended particles, also including most fine clay particles.
The hydrodynamic gradient lays within the threshold for bedload transport conditions. Generally, it is assumed that for the rivers with a similar size as Selenga, values lower than the critical Rouse number Ro* = 2.5 indicate that a sediment particle begins to contribute to the suspended load [38,39], whereas higher values indicate that a sediment particle is most likely transported as part of the bedload [40]. To test the relationship between the observed Ro numbers, we plotted our results (Figure 8) of the relationship between the ratio of suspended load to total load and Rouse number based on the laboratory data for the mixed load from Guy, Simons, and Richardson [41]. Lines shown are those from the Einstein integrals of Guo and Julien, as obtained by Shah-Fairbank [42] for h/d 100 and 100.000. It can be clearly seen that the measured conditions in the Selenga River are generally close to fit line for the rivers at h/d50 = 100.000. At the same time, this shows that one can obtain an extremely large variability in sediment concentration and thus the Q R /(Q G + Q R ) ratio in deep sand-bed rivers when the Rouse number is fairly large (Ro > 0.5). The case study on Selenga River is quite instructive in this regard, which reflects conditions of the abrupt shift from a mixed load and bedload dominated channel (Ro < 2) to a suspended sediment dominated channel (Ro > 2).

Particulate Metals Behavior and SB-Analyses
Both hydrodynamic (stream velocities) and hydrological (water discharge rates) impacts can explain the patterns of particulate metals along the Lower Selenga River. Absolute concentrations along the river course are mostly driven by temporal changes of water discharge during flood propagation. The latter significantly influences the increase of particulate concentrations in the delta compared to the upper adjacent main channel reach. This fact, which contradicts previous observations on suspended sediment retention within a delta [8,19] and storage of particulate and dissolved chemicals [18], can be explained by the extensive role of storm events in geochemical flows. As it was shown in the previous regional studies [44], relatively fine particles (62-250 µm; particle size 1) are eroded preferentially over coarser particles during the high flow events, implying their prevalence in suspension during these high flow conditions. Vice versa, these fine fractions are most important in the transport of metals [6,45], and, thus, overall trends of concentrations increase are consistent among the different heavy metals, although the magnitude of change differs. This implies an important consideration that the delta function as a geochemical buffer is time-and discharge-dependent. This statement is in general agreement with the contrasting sedimentation patterns observed via satellite images [19], which reveal sediment storage (a downward decline of sediment transport) in the Selenga delta during high discharges (>1500 m 3 /s). Sediment concentrations increase can be seen under lower streamflow conditions (<1500 m 3 /s). This might be particularly important due to observed influence of climate change on the seasonal discharges of the Selenga River [46][47][48], which cause a significant impact on channel processes in the lower reaches of the river [32] and chemical composition due to the drought period observed for the last 20 years [49,50].
We performed a quantitative description of metals partitioning between the bed and suspended fluxes (SB-analyses) using the ratio between concentrations of element x in the subsurface layer and near-bottom layer Kx (see Section 2.3). This parameter significantly varied along the transition zone from the sand-gravel of the Selenga River to the sandy-silt channels of the delta area, where the depositional pattern led to higher sorting of fine metal-bearing fractions. For element concentrations, Kx increased significantly in the delta channels ( Table 4). The hydrodynamic structure of the flow remains largely unchanged (a decrease in velocities from the surface to the bottom is characterized by a similar ratio Kv = 0.94 upstream of the delta and Kv = 0.96 in the delta), while the average velocity decreases steadily downstream the delta area (see Figure 3). The leading factor in the formation of the heavy metal transfer pattern as a suspended load should be considered in the peculiarities of the mass transfer between the channel flow and bottom sediments, which differs depending on the type of sediment that the channel consists of. On average, in the delta sandy-silt channels, the increase in heavy metals concentration in suspended sediment is 30%-40% higher than in the Selenga River main channel. These gradients can be associated with a larger concentration of fine silt at the bottom of the floodplain channels, and the deposition of heavier sand mineral fractions containing metals in their crystal lattice. At the same time, the SSC gradient in the delta channels increases even more significantly (Kssc = 1.75 in the delta versus Kssc = 0.99 in the main channel). The tremendous role of larger mineral fractions sedimentation in the formation of vertical heterogeneity is confirmed by the sand fraction (>50 µm) content increase in the composition of the bottom sediments of the delta channels (Ksand (%) = 1.75). Finally, we accounted for impacts of hydrodynamic conditions and sediment partitioning on metals distribution within the depth column using joint analyses of Ro and Kx numbers (Figure 9). Due to the high variability of Ro under transient conditions of the channel-delta interface and also unsteady flow during flood propagation, we were not able to obtain significant relationships between hydrodynamic parameters and geochemical coefficients. Nevertheless, Figure 9 clearly demonstrates the increase of the metals in the near-bottom layers at the low values of Ro (dominant suspended mode of transport) and fits the following relationship (Rcor = −0.32, p-value < 0.001, n = 150): The observed pattern can be explained by the absence of coarser fractions in the near-bottom layer due to sand and gravel termination in the upper part of delta [7]. Vice versa, a strong sorting of clay particles under low Ro values conditions became the crucial driver of Kx increase. Hence we can conclude that grain size sorting varied between suspended load and bedload dominated river reaches. In the suspended dominated reach, a clearer sorting of fine particles is observed, which explains the metals' behavior and sorting along the depth profile. In some samples, small portions of sand were also observed in the delta edge, which is mostly explained by its possible delivery from the eroded bank during the flood propagation. Upstream from the delta in the main channel of the Selenga River, the hydrodynamic enrichment of relatively coarse quartz grains near the bottom of the channel can be attributed to observed low values of Kx. Unfortunately, at this point, we did not collect enough samples to prove this statement with a sufficient number of grain size profiles.
The observed patterns of hydrodynamic grain size and metal sorting do not have unique patterns and depend on the hydrodynamic conditions of sediment transport. In similar hydrodynamic conditions with low Rouse number and dominant suspended load, as observed in the case of the Selenga delta, Kx typically has values >1. In the Ganges River, the trace elements in the suspended load As, Ba, Bi, Co, Cr, Cs, Cu, Ga, Rb, V and Zn concentrations were negatively correlated with d84 grain size and were higher in the suspended load compared to bedload samples [6], indicating high numbers of Kx within a vertical profile. Vice versa, in the case of bedload-dominated rivers-the Amazon River and its tributaries [5] and Yangtze River [51]-the Kx value is significantly lower than 1. Here the amount of quartz in the samples is the main factor controlling the concentration of these elements throughout the sampled depth profiles due to hence so-called "quartz dilution effect" [5]. This interrelation between mineralogy and grain size is due to the fact that quartz as a mineral is mostly related to coarser sand sediments (>50 µm or larger), which is proved by empirical evidence worldwide (e.g., [33]). The absence of a significant correlation between Kx and Ro for any of the elements x can be attributed to unsteady flow conditions during flood propagation. Another example was found at the Yellow River flow [51] under conditions of the sand-dominated suspended load. Here due to the possible impact of artificial flood from the reservoir, the extremely high-velocity profiles lead to suspension and diffusion of coarse sand material, which also leads to relatively low values (Kx < 1). The revealed patterns are important for organizing surface water monitoring and assessing the transformation of pollutant fluxes along the river. Existing studies, which rely on single surface water samplings, are associated with large uncertainties. These approaches can be characterized by the incompatibility of datasets obtained for various river reaches. For instance, at the reaches where Kx << 1 or Kx >> 1, surface water sampling will induce high errors in particulate fluxes' estimates. For the Baikal region in particular, the hydroclimatic impact on metals spreading may significantly influence pollutant delivery into the Lake Baikal. Regional climate change has been proven to influence the magnitude and number of extreme flow events [47]. Thus, the impact on sediment storage during these events may be significantly changed compared to previous conditions due to enhanced depositional or erosional processes.
Finally, we conclude that consideration of hydrodynamic conditions is of primary importance for sediment quality monitoring design.

Conclusions
We revealed the sorting impact on erosion products and associated heavy metals by river flow due to interaction between the transient flow and channel bed (near-bottom sediments (bedload) and bottom deposits. The following main conclusions are drawn: 1.
Due to hydrodynamic impacts of varying channel patterns along the river continuum, various patterns of heavy metals are observed along the river course. Depositional channels such as Selenga delta branches act as a better "sorting machine" of heavy metals (30% to 40% vertical distribution is higher in the delta than in the gravel-boulder main channel of the Selenga river).

2.
This increase can be attributed to a higher concentration of the clay-silt fraction near the bed in the deltaic channels due to its heavier mineral structure containing specific microelements. The SSC gradient towards the bottom is further pronounced (Kssc = 1.75 vs. Kssc = 0.99 in the main channel).

3.
This vertical variability is also caused by the pronounced increase of coarse mineral fractions (50-2000 mm) in the near-bottom layer in the depositional delta pattern (K (sand (%)) = 1.75) compared to the unbraided main stem of the Selenga river where sediments are deposited on the bottom, thus having increased content in the near-bottom layers.
Author Contributions: The work was conceived by S.C. The text was mostly written by S.C., V.M. and E.S. V.M. and E.S. were responsible for field data collection and processing. V.E. and L.E. did the hydrochemical data processing. S.C., M.L. and N.K. provided general supervision of the project and editing. All authors have read and agreed to the published version of the manuscript.
Funding: Field studies were supported by the Russian Fund for Basic Research (project17-29-05027) and the RFBR-Russian Geographical Society project (RFBR-17-05-41174-RGS); sampling methodology was developed as a part of RFBR project 18-05-60219. Channel processes analyses were done under the Russian Scientific Foundation project 18-17-00086.