Evaluating Regime Change of Sediment Transport in the Jingjiang River Reach , Yangtze River , China

The sediment regime in the Jingjiang river reach of the middle Yangtze River has been significantly changed from quasi-equilibrium to unsaturated since the impoundment of the Three Gorges Dam (TGD). Vertical profiles of suspended sediment concentration (SSC) and sediment flux can be adopted to evaluate the sediment regime at the local and reach scale, respectively. However, the connection between the vertical concentration profiles and the hydrologic conditions of the sub-saturated channel has rarely been examined based on field data. Thus, vertical concentration data at three hydrological stations in the reach (Zhicheng, Shashi, and Jianli) are collected. Analyses show that the near-bed concentration (within 10% of water depth from the riverbed) may reach up to 15 times that of the vertical average concentration. By comparing the fractions of the suspended sediment and bed material before and after TGD operation, the geomorphic condition under which the distinct large near-bed concentrations occur have been examined. Based on daily discharge-sediment hydrographs, the reach scale sediment regime and availability of sediment sources are analyzed. In total, remarkable large near-bed concentrations may respond to the combination of wide grading suspended particles and bed material. Finally, several future challenges caused by the anomalous vertical concentration profiles in the unsaturated reach are discussed. This indicates that more detailed measurements or new measuring technologies may help us to provide accurate measurements, while a fractional dispersion equation may help us in describing. The present study aims to gain new insights into regime change of sediment suspension in the river reaches downstream of a very large reservoir.


Introduction
Characteristics of the sediment regime are one of the most important factors that control the processes in the transportation and sedimentation zones of a fluvial system [1].Changes in the sediment regime cause not only varied fluvial evolution processes but also the inner physical processes of sediment suspension.Using a reach scale, the sediment regime of a river reach can be expressed as the suspended sediment load (sediment flux).Regarding the physical processes of suspension on a local scale, they may be related to the vertical profiles of suspended sediment concentration (SSC), which are essential for estimating sediment flux, and also numerical simulation [2].Thus, a large number of researchers have paid, and are paying, attention to the vertical concentration profiles.During the early stage of physical or analytical analyses, many studies have investigated steady and equilibrium sediment transport [3].For equilibrium sediment transport, traditional advection-diffusion equation (ADE) theory and the improvements (including sediment turbulent diffusion coefficient, vertical flow velocity distribution, and suspension index) are widely used [4,5].For calculating profiles of SSC, numerical models produce good results for uniform sediments in equilibrium situations [6,7].However, this is not the case for non-uniform sediments in equilibrium situations, not to mention non-uniform sediments in non-equilibrium situations [6,7].For instance, Nicholas et al. [8] compared observed flow widths with results from a theoretical model developed for non-equilibrium (aggradation) conditions.Analysis results illustrated that models derived for equilibrium conditions may have limited utility in non-equilibrium situations, despite their widespread use to date [8].
In the non-equilibrium condition, the deviations of the vertical profiles of SSC from the local equilibrium state are different from those of the equilibrium state [9,10].For hydraulic engineers, the vertical concentration profiles are essential for estimating erosion and deposition using the sediment balanced method, and also for estimating sediment-associated contaminants [11,12].For geomorphologies, non-equilibrium suspended sediment transport (SST) is one of the key processes driving the channel evolution by erosion or deposition [13].Thus, various approaches have been employed in non-uniform sediments in non-equilibrium sediment transport.However, most of these analyses are based on flume data or theoretical analysis [14,15].The connection between the hydrologic condition of the non-equilibrium channel and the vertical concentration profiles has rarely been examined based on field data [16].
The regime of sediment transport in the Jingjiang river reach (JJRR, about 102 km downstream of the Three Gorges Dam (TGD)) has been significantly changed from quasi-equilibrium to sub-saturation, as the impoundment of the TGD.More and more researchers are focusing on the changed sediment regime of the sub-saturated channel in the reach scale.However, field observation has revealed frequent occurrences of remarkable large concentrations in the near-bed zone by extending the vertical measuring to the near-bed zone (i.e., with a distance of 0.1 m from the riverbed).Thus, it is also necessary to analyze the changed sediment regime in the local scale.The contribution of the remarkable, large, near-bed concentration on the sediment flux in the unsaturated channel has been analyzed [17].
Here, this paper aims to evaluate the changed sediment regime of the JJRR from both a local-scale view (vertical concentration profiles) and a reach-scale view (sediment flux).It also helps one to understand the connection between the concentration profiles (local-scale) and the hydrologic conditions of the sub-saturated channel.
Firstly, the vertical profiles of SSC are analyzed by their characteristics, hydrodynamic conditions, and geomorphic conditions to exhibit the changed sediment regime at the local scale.Secondly, the temporal variation of sediment flux is analyzed with the cumulative anomaly method to exhibit the changed sediment regime at the reach scale.Thirdly, the relationship between SSC and discharge is analyzed with the method of the sediment rating curve (SRC), which is used to investigate the conditions for the occurrence of the remarkable large near-bed concentration.Finally, the challenges that may be caused by the remarkable, large, near-bed concentration are discussed.

The Study Reach
The length of the Yangtze River (YR) is 6.3 × 10 3 km, and the drainage area is 1.8 × 10 6 km 2 (Figure 1).The TGD, located at the end of the upper reach of the YR, controls a drainage area of 1.0 × 10 6 km 2 .The dam is 185 m high, and storage capacity of the reservoir is 3.93 × 10 4 Hm 3 .The main purposes of the project are flood control, power generation, and navigation.
The JJRR (between Zhicheng and Chenglingji stations) is about 60 km downstream from the Yichang hydrological station.There are five hydrological stations in the JJRR: Zhicheng, Shashi (Jing 45), Xinchang (Jing 84), Jianli, and Chenglingji (Figure 1b).Moreover, there are three distributary channels on the reach through which YR delivers water and sediment into Dongting Lake, Songzikou, Taipingkou, and Ouchikou.The river reach can be divided by Ouchikou into the upper and lower sub-reaches (Figure 1b).The lengths of these two sub-reaches are 172 km and 176 km, respectively.
Water 2018, 10, x FOR PEER REVIEW 3 of 21 Taipingkou, and Ouchikou.The river reach can be divided by Ouchikou into the upper and lower sub-reaches (Figure 1b).The lengths of these two sub-reaches are 172 km and 176 km, respectively.

Data
Both daily data and field surveyed data were used in the present study.Daily discharge and suspended sediment load (SSL) are measured and published by the Yangtze River Water Resources Commission (YRWRC).Data at three stations (Zhicheng, Shashi, and Jianli) in 1987-2014 were obtained from yearbooks.Distances of these three stations downstream from the Yichang station are 60 km, 157 km, and 300 km, respectively.
Field surveyed data were provided by the Jingjiang Hydrology and Water Resources Surveying Bureau (JHWRSB), YRWRC.Data include vertical concentration, sediment gradation, and corresponding velocities.Table 1 manifests the years of measurement and the numbers of vertical profiles.Profiles measured before and after 2010 are different.Before 2010, there are five measuring points in each vertical profile.The normalized depths are approximately 1.0H, 0.8H, 0.4H, 0.2H, and 0.1H (H is the water depth of the vertical line, in meters), respectively.After 2010, there are seven measuring points in each vertical profile.The normalized depths of the upper five points are approximately 1.0H, 0.8H, 0.4H, 0.2H, and 0.1H, respectively.The other two points were measured in the near-bed region with constant distances from the riverbed of 0.5 m and 0.1 m.

Data
Both daily data and field surveyed data were used in the present study.Daily discharge and suspended sediment load (SSL) are measured and published by the Yangtze River Water Resources Commission (YRWRC).Data at three stations (Zhicheng, Shashi, and Jianli) in 1987-2014 were obtained from yearbooks.Distances of these three stations downstream from the Yichang station are 60 km, 157 km, and 300 km, respectively.
Field surveyed data were provided by the Jingjiang Hydrology and Water Resources Surveying Bureau (JHWRSB), YRWRC.Data include vertical concentration, sediment gradation, and corresponding velocities.Table 1 manifests the years of measurement and the numbers of vertical profiles.Profiles measured before and after 2010 are different.Before 2010, there are five measuring points in each vertical profile.The normalized depths are approximately 1.0H, 0.8H, 0.4H, 0.2H, and 0.1H (H is the water depth of the vertical line, in meters), respectively.After 2010, there are seven measuring points in each vertical profile.The normalized depths of the upper five points are approximately 1.0H, 0.8H, 0.4H, 0.2H, and 0.1H, respectively.The other two points were measured in the near-bed region with constant distances from the riverbed of 0.5 m and 0.1 m.For profiles measured at Shashi station, measured flow depths range from 2.09 m to 19.7 m, and the normalized depths of these two near-bed points are approximately (0.02-0.2)H.For profiles measured at Jianli station, measured flow depths range from 2.13 m to 21.2 m, and the normalized depths of these two near-bed points are approximately (0.0047-0.047)H.1986, 1998, 2002, 2010, 2011, 2012, 2013 180 Equipment adopted in the field survey includes global navigation satellite system (GNSS, antenna, and receiver), total station, vessel-mounted acoustic Doppler current profilers (ADCP), digital level, GPS, and laser particle size analyzer (LPSA) [17].All of this equipment has been strictly examined by professional surveyors.According to the technical manual by JHWRSB, several technical standards and criteria are illustrated, as follows: (1) Water depth is measured with ADCP, and verified with fish lead.
(2) ADCP is applied to measure the flow twice.The deviation between each measured discharge and the average discharge should be less than ±5%.Otherwise, the data should be re-measured.(3) Vertical lines for measuring the velocity are positioned with a real-time kinematic (RTK) GNSS.(4) Suspended sediments are sampled with samplers in the field, and its grain-size distributions are analyzed by the sieving method.All the sampling procedures should be finished within dozens of minutes, and sieving analyses should be finished following the manual operation.For all the five points measured before 2010 and the upper five points measured after 2010, the sampler shown in Figure 2a is adopted.For the two near-bed points measured after 2010, bottom-touched automatic-closing samplers shown in Figure 2b,c are adopted.With this kind of sampler (Figure 2b,c), the distances of the two near-bed points from the riverbed can be ensured.However, the potential disturbance on the riverbed when settling the sampler may still be questionable.Therefore, except a slowed down settling velocity of the sampler, the near-bed concentrations may be double-checked after sampling and grain-size analyses.The assumption is that the measured concentration is reliable only if the vertical concentration profiles of particles finer than 0.062 mm have no abrupt changing point in the near-bed region [17].Thus, bed materials can be distinguished from suspended sediments, which may be caused by disturbance during the sampling processes.

Methodology
Methods of cumulative anomaly and SRCs are applied to investigate the changed sediment regime of un-saturated channel in reach scales.ADEs and the sediment-balance method are adopted to evaluate new challenges for the non-saturated channel.

Cumulative Anomaly
The cumulative anomaly has been widely applied to detect the trends and stages of time series [18].For time series, the cumulative anomaly (X t ) for data point x i can be expressed as: in which x i is the value of the time series; x m is the mean value of the series x i ; and n is the length of the time series.The mass curve of cumulative anomaly values describes the change process of a particular parameter by comparing it with the mean value of that parameter over the whole study period [18].For instance, if the observed data of a given year is greater than the overall mean, the anomaly values will be positive and the mass curve will rise.

SRCs
The SRC is defined as the statistical relationship between SSC (kg•m −3 ) and discharge.SRCs are normally used to describe the flow-sediment relationships in river systems for various purposes [19].The method is often applied to estimate the SSC or SSL in ungauged regions [20].
The power function of SRC is usually expressed as one of the following two formats: , and e and f are the sediment rating coefficient and exponent, respectively.Linear regression, i.e., Equation (3), is often used to derive the values of the rating coefficient e and exponent f from empirical data (e.g., [21]).The e-parameter contains information on the conversion of Q into SSC, namely, the erosion severity index [22].A high value of the e-parameter indicates easily-erodible materials and high loads of transported materials [20].The exponent f -coefficient corresponds to the erosive power of the river and the influence on the sediment supplement.High values of the f -coefficient indicate a slight increase in discharge would significantly enhance the erosive power of the river [20].Moreover, the f -coefficient is also related to climate, channel morphology, the grain-size distribution of sediments, and erodibility within the river basin [22].

ADEs
Three equations are adopted to evaluate the challenge in simulating the vertical concentration profiles in the unsaturated channel: the Rouse Equation, the Han Equation, and the fADE.These three equations may represent typical concentration profiles in equilibrium transport, non-equilibrium transport, and non-Fickion suspension in the non-equilibrium condition, respectively.
By assuming Fick's first law for sediment diffusion in turbulence, the diffusion theory of sediment suspension can be expressed with traditional ADE model.The equation for equilibrium sediment transport can be expressed as: in which S is the sediment concentration (W•L −3 ), y is the vertical coordinate (L), ω is the sediment settling velocity (L•T −1 ), and ε sy is the sediment turbulent diffusion coefficient in the y direction (L•T −1 ).By replacing ε sy with the fluid eddy viscosity and assuming ε m = κu * (1 − y/H)y based on the Karman-Prandtl logarithmic velocity profile, Rouse et al. [4] obtained the analytical vertical concentration profile: (5) in which S a is a reference concentration at a given height above the river bed (W•L −3 ), α is the reference height (L), H is the flow depth (L), κ is von Karman's constant, and u * is the shear velocity (L•T −1 ).The frictional flow velocity can be calculated as u * = τ b /ρ, in which τ b is the near-bed shear stress and ρ is the water density.The Rouse equation, Equation (5), has been widely used for decades.However, its drawbacks are obvious: the concentration is calculated as zero at the water surface and infinity at the river bed.These drawbacks are caused by implicit assumptions, equilibrium sediment transport, and Fick's first law.Extensive efforts have been put forward to improve the equation [10,[23][24][25].
The assumption of equilibrium sediment transport limits its application in the field in which the sediment regime is far from equilibrium [26].Among others, Han et al. [9] have modified Equation (4) to meet the non-equilibrium condition: in which q s is the net flux due to the imbalance between downward sediment settling and upward turbulent dispersion.Han et al. [9] obtained the solution of Equation (6): in which c is non-equilibrium coefficient defined as a function of the degree of saturation c = According to Fick's First Law, particles may make "local" jumps, as most particle jumps induced by turbulence are constrained to a small distance (∆y) in a given ∆t according to the Central Limit Theorem (CLT), and ∆y is characterized by the length of Representative Elementary Volume (REV) [27,28].However, modern observations of coherent structures in turbulence have proven that high-speed currents may intermittently sweep the bed, carry much sediment, and eject directly into the upper part of the water column [27,29].The non-Fickion suspension indicates that particles may make many "nonlocal" jumps in burst-like suspension events [27].The non-Fickion suspension in non-equilibrium sediment suspension can be described by a fractional advection-dispersion equation (fADE) [27]: in which α is the order of the fractional derivative (0 < α ≤ 1); ε sy (L α •T −1 ) is the depth-averaged diffusivity expressed as follows: We assume: in which c 0 is called the non-equilibrium coefficient, and c 0 is a function of (S/S * ): in which S * is the sediment transport capacity of flow; if S = S * , then c 0 = q s = 0, which means the flow reaches equilibrium; if S < S * , then c 0 > 0 and q s < 0, which means the flow is unsaturated with sediment; and if S > S * , then c 0 < 0 and q s > 0, which means flow is supper-saturated with sediment.
The analytical solution of Equation ( 8) is obtained: in which E a () is called Mittag-Leffler function (y ≤ a ≤ H), which can be calculated by series expansion or the MATLAB open-source code.

Contribution of Near-Bed Concentrations on SST Rate
The SST rate can be estimated with measured vertical concentration profiles, and the deposition/ erosion amount of the channel can be estimated by the sediment balance method.Changed vertical concentration profiles in the unsaturated channel lead to new challenges in field surveying and estimating the SST rate.The field survey should be conducted according to its new characteristic (as extended measuring of near-bed zone by JHWRSB).The effect of the changed vertical concentration profiles on the SST rate can be estimated by evaluating the contribution of two near-bed concentrations on the SST rate: in which Q s (7) and Q s(5) are cross-sectional SST rates by seven points and, artificially, by five points, respectively.The cross-sectional sediment flux can be estimated with all measured vertical profiles.The delta-shaped area between the left bank and the first vertical line from the left bank is also considered [17].The delta-shaped area between the right bank and the last vertical line from the left bank is not considered due to the difficulty of identifying the right bank [17].

Change of the Sediment Regime
Vertical concentration profiles and SST rates are analyzed to evaluate the sediment regime of local and reach scales, respectively.Then, sediment rating curves are analyzed to investigate the availability of sediment for suspension from the riverbed.

Vertical Profiles of SSC
3.1.1.Remarkable Large Concentrations in the Near-Bed Zone e depth-averaged concentration S avg (kg•m −3 ) is calculated as S avg = 1 H S i dy, in which S i is the measured concentration at each point in a vertical line (kg•m −3 ) and y is the vertical axis of the Cartesian coordinate system.The coordinate origin is located on the riverbed.To present the deviation of the near-bed concentration from the average, the relative concentration is defined as S a /S avg , in which S a is the measured concentration at the reference height.The reference height is the minimum distance of measuring points from the river bed.The reference height a equals 0.4-0.5 m for profiles measured before 2010, while the value is 0.1 m for profiles after 2010.
Figure 3 shows the typical vertical profiles measured at Shashi station before and after the impoundment of TGD.The relative concentrations (S a /S avg ) of these four vertical profiles are 2.47, 4.1, 2.10, and 1.52, respectively.Thus, remarkable large concentrations can be observed in the near-bed zone.The remarkable large near-bed concentration is defined by comparison to the vertical average concentration.This can be observed in vertical concentration profiles with large or small measured near-bed concentrations (3.01 kg•m −3 for Figure 3a and 0.041 kg•m −3 for Figure 3b).
The statistic characteristics of hundreds of vertical profiles are summarized in Table 2. Before TGD operation, the maximum relative concentrations (S a /S avg ) are 2.76 (Zhicheng), 2.30 (Shashi), and 3.62 (Jianli), respectively.After TGD operation, the maximum relative concentrations (S a /S avg ) are 15.05 (Shashi), and 14.45 (Jianli), respectively.The average relative concentrations before TGD are 1.35 (Shashi) and 1.57 (Jianli), respectively.These values jump to 3.22 (Shashi) and 4.89 (Jianli) after TGD, respectively.Thus, both the average and maximum values of relative concentrations (S a /S avg ) increase significantly after the impoundment of the TGD.
Relatively high concentrations near the bed region have also been pronounced in the estuary of YR, Qiantangjiang River, and the Taizhou sea area [30][31][32][33].The near-bed concentration of the S-type vertical distribution is also relatively large [34].However, the relative concentrations are relatively small (shown in Table 2).The maximum relative concentration (S a /S avg ) is approximately 4.3, and the minimum relative height of measured points from river bed is approximately 0.02.In total, the tailing phenomena in the unsaturated channel are much more apparent.
Water 2018, 10, x FOR PEER REVIEW 8 of 21

Change of the Sediment Regime
Vertical concentration profiles and SST rates are analyzed to evaluate the sediment regime of local and reach scales, respectively.Then, sediment rating curves are analyzed to investigate the availability of sediment for suspension from the riverbed.

Remarkable Large Concentrations in the Near-Bed Zone
The depth-averaged concentration Savg (kg•m −3 ) is calculated as = , in which Si is the measured concentration at each point in a vertical line (kg•m −3 ) and y is the vertical axis of the Cartesian coordinate system.The coordinate origin is located on the riverbed.To present the deviation of the near-bed concentration from the average, the relative concentration is defined as Sa/Savg, in which Sa is the measured concentration at the reference height.The reference height is the minimum distance of measuring points from the river bed.The reference height a equals 0.4-0.5 m for profiles measured before 2010, while the value is 0.1 m for profiles after 2010. Figure 3 shows the typical vertical profiles measured at Shashi station before and after the impoundment of TGD.The relative concentrations (Sa/Savg) of these four vertical profiles are 2.47, 4.1, 2.10, and 1.52, respectively.Thus, remarkable large concentrations can be observed in the near-bed zone.The remarkable large near-bed concentration is defined by comparison to the vertical average concentration.This can be observed in vertical concentration profiles with large or small measured near-bed concentrations (3.01 kg•m −3 for Figure 3a and 0.041 kg•m −3 for Figure 3b).
Relatively high concentrations near the bed region have also been pronounced in the estuary of YR, Qiantangjiang River, and the Taizhou sea area [30][31][32][33].The near-bed concentration of the S-type vertical distribution is also relatively large [34].However, the relative concentrations are relatively small (shown in Table 2).The maximum relative concentration (Sa/Savg) is approximately 4.3, and the minimum relative height of measured points from river bed is approximately 0.02.In total, the tailing phenomena in the unsaturated channel are much more apparent.The coarsening of suspended sediment can be revealed by mean diameters of vertical profiles (shown in Table 2).For vertical lines measured at Jianli station, the average mean diameters after TGD range from 0.0087 mm to about 0.22 mm, while the values before TGD may range from 0.0037 mm to 0.07 mm.For vertical lines measured at Shashi station, the average mean diameters after TGD change from 0.0058-0.027mm to 0.0098-0.224mm.Thus, the coarsening process of suspended sediment after TGD operation is apparent, especially at Shashi station.
Both of the near-bed SSC and sediment flux are connected with hydrodynamic conditions [35].Thus, hydrodynamic parameters (depth, velocity, and vertical average concentration) are also summarized in Table 2.The values of ω/u * vary dramatically (shown in Table 2).Remarkable large values of ω/u * indicate that relatively coarser particles may be suspended with a small bed shear velocity.Tang et al. [32] pointed out that the near-bed concentration may be more apparent when (ω/κu * ) > 5. Therefore, larger values of log(ω/u * ) in non-equilibrium channels contribute to remarkable large concentrations in the near-bed zone.
Normally, concentration at the water surface should be smaller than that at the riverbed.This means that the concentration gradient (estimated with d(Si − Si−1)/d(yi − yi−1)) should be negative.However, the occurrence of vertical concentration profiles with positive gradients increases (Table 3).For data measured before TGD, about 20% (Zhicheng, Shashi, and Jianli stations) vertical concentration profiles are characterized with positive gradients.For data measured after TGD, about 28-47% (Shashi and Jianli) of vertical concentration profiles are characterized with positive gradients (Table 3).The differential equation for vertical SSC distribution of non-uniform particles indicates that the concentration of finer particles group near the surface may be greater than that near the bed with wide grading non-uniformed suspended sediment [36].Moreover, concentration gradients may be positive when the diameter of coarser particles is about eight times that of finer particles [36].This indicates wide grading sediment particles in the non-equilibrium channel contribute to the remarkable large near-bed concentration.The coarsening of suspended sediment can be revealed by mean diameters of vertical profiles (shown in Table 2).For vertical lines measured at Jianli station, the average mean diameters after TGD range from 0.0087 mm to about 0.22 mm, while the values before TGD may range from 0.0037 mm to 0.07 mm.For vertical lines measured at Shashi station, the average mean diameters after TGD change from 0.0058-0.027mm to 0.0098-0.224mm.Thus, the coarsening process of suspended sediment after TGD operation is apparent, especially at Shashi station.
Both of the near-bed SSC and sediment flux are connected with hydrodynamic conditions [35].Thus, hydrodynamic parameters (depth, velocity, and vertical average concentration) are also summarized in Table 2.The values of ω/u * vary dramatically (shown in Table 2).Remarkable large values of ω/u * indicate that relatively coarser particles may be suspended with a small bed shear velocity.Tang et al. [32] pointed out that the near-bed concentration may be more apparent when (ω/κu * ) > 5. Therefore, larger values of log(ω/u * ) in non-equilibrium channels contribute to remarkable large concentrations in the near-bed zone.
Normally, concentration at the water surface should be smaller than that at the riverbed.This means that the concentration gradient (estimated with d(S i − S i−1 )/d(y i − y i−1 )) should be negative.However, the occurrence of vertical concentration profiles with positive gradients increases (Table 3).For data measured before TGD, about 20% (Zhicheng, Shashi, and Jianli stations) vertical concentration profiles are characterized with positive gradients.For data measured after TGD, about 28-47% (Shashi and Jianli) of vertical concentration profiles are characterized with positive gradients (Table 3).The differential equation for vertical SSC distribution of non-uniform particles indicates that the concentration of finer particles group near the surface may be greater than that near the bed with wide grading non-uniformed suspended sediment [36].Moreover, concentration gradients may be positive when the diameter of coarser particles is about eight times that of finer particles [36].This indicates wide grading sediment particles in the non-equilibrium channel contribute to the remarkable large near-bed concentration.Relatively large concentrations in the near-bed zone in estuary areas can be observed during tide periods and non-flood seasons [30].Thus, the hydrologic condition (suspended sediment and bed materials) under which the remarkable large concentrations in the near-bed zone may occur is also analyzed.
With the operation of TGD, both the suspended and bed materials become much coarser, except for the dramatic decrease of the annual suspended sediment load (SSL) [37].The temporal variation of suspended sediment size in the downstream channel may exhibit some complex processes [38].This means that the saturation recovery process of fine sediments that responds to changing hydraulic conditions is different from that of the coarser one.Additionally, under non-uniform and non-equilibrium conditions, suspended sediment dynamics require a relatively large distance or time to approach equilibrium [13].Thus, the variation trend of non-uniform sediment varies in different river reaches.
The mean diameter of SSC measured at Yichang station decreases from 0.009 mm (before impoundment) to 0.005 mm (2003-2008) [37].Size distributions of suspended sediment and bed materials in the JJRR are shown in Figures 4 and 5.For Zhicheng station, the proportion of grouped particles with diameters of 0.008-0.062mm increases, while the proportions of finer (<0.008 mm) and coarser (0.062-0.5 mm) groups decreases.For Shashi station, the non-uniform sediment becomes a kind of wide grading (as shown in Figure 4b), as the proportion of grouped particles with diameters finer than 0.008 mm decreases and the proportion of grouped particles with diameters coarser than 0.062 mm increases.The mean diameter measured at Jianli station decreases from 0.009 mm (before impoundment) to 0.045 mm (2003)(2004)(2005)(2006)(2007)(2008) [37].This indicates that the coarser particles (>0.062 mm) around Shashi station have recovered by bed erosion.When flow transports the reach around Jianli station, the change of coarser particles (>0.062 mm) is much more apparent.Thus, wide grading suspended sediment at Shashi and Jianli stations may contribute to the remarkable large concentration in the near-bed zone.Relatively large concentrations in the near-bed zone in estuary areas can be observed during tide periods and non-flood seasons [30].Thus, the hydrologic condition (suspended sediment and bed materials) under which the remarkable large concentrations in the near-bed zone may occur is also analyzed.
With the operation of TGD, both the suspended and bed materials become much coarser, except for the dramatic decrease of the annual suspended sediment load (SSL) [37].The temporal variation of suspended sediment size in the downstream channel may exhibit some complex processes [38].This means that the saturation recovery process of fine sediments that responds to changing hydraulic conditions is different from that of the coarser one.Additionally, under non-uniform and non-equilibrium conditions, suspended sediment dynamics require a relatively large distance or time to approach equilibrium [13].Thus, the variation trend of non-uniform sediment varies in different river reaches.
The mean diameter of SSC measured at Yichang station decreases from 0.009 mm (before impoundment) to 0.005 mm (2003-2008) [37].Size distributions of suspended sediment and bed materials in the JJRR are shown in Figures 4 and 5.For Zhicheng station, the proportion of grouped particles with diameters of 0.008-0.062mm increases, while the proportions of finer (<0.008 mm) and coarser (0.062-0.5 mm) groups decreases.For Shashi station, the non-uniform sediment becomes a kind of wide grading (as shown in Figure 4b), as the proportion of grouped particles with diameters finer than 0.008 mm decreases and the proportion of grouped particles with diameters coarser than 0.062 mm increases.The mean diameter measured at Jianli station decreases from 0.009 mm (before impoundment) to 0.045 mm (2003-2008) [37].This indicates that the coarser particles (>0.062 mm) around Shashi station have recovered by bed erosion.When flow transports the reach around Jianli station, the change of coarser particles (>0.062 mm) is much more apparent.Thus, wide grading suspended sediment at Shashi and Jianli stations may contribute to the remarkable large concentration in the near-bed zone.The variations of volume fractions of bed materials indicate a remarkable coarsening process in the JJRR, especially the river reach of Zhicheng-Shashi (Figure 5).Zhang [30] pointed out that concentrations in the near-bed zone are more apparent with larger values of relative particle size (Di/Dm, particle-size over averaged particle size) of non-uniform sediment.Thus, the lifting and suspension of wide grading bed materials may also contribute to the remarkable large concentration in the near-bed region of the non-equilibrium channel.The variations of volume fractions of bed materials indicate a remarkable coarsening process in the JJRR, especially the river reach of Zhicheng-Shashi (Figure 5).Zhang [30] pointed out that concentrations in the near-bed zone are more apparent with larger values of relative particle size (D i /D m , particle-size over averaged particle size) of non-uniform sediment.Thus, the lifting and suspension of wide grading bed materials may also contribute to the remarkable large concentration in the near-bed region of the non-equilibrium channel.The variations of volume fractions of bed materials indicate a remarkable coarsening process in the JJRR, especially the river reach of Zhicheng-Shashi (Figure 5).Zhang [30] pointed out that concentrations in the near-bed zone are more apparent with larger values of relative particle size (Di/Dm, particle-size over averaged particle size) of non-uniform sediment.Thus, the lifting and suspension of wide grading bed materials may also contribute to the remarkable large concentration in the near-bed region of the non-equilibrium channel.

Temporal Change of SSL
At the reach scale, the changed sediment regime can be described with temporal and spatial variations of SSL.The temporal and spatial variations of water runoff have also been analyzed to make a comparison.The temporal variation of water runoff and SSL show a decreasing trend, especially the SSL (Figure 6).
Water 2018, 10, x FOR PEER REVIEW 13 of 21

Temporal Change of SSL
At the reach scale, the changed sediment regime can be described with temporal and spatial variations of SSL.The temporal and spatial variations of water runoff have also been analyzed to make a comparison.The temporal variation of water runoff and SSL show a decreasing trend, especially the SSL (Figure 6).

Temporal Change of SSL
At the reach scale, the changed sediment regime can be described with temporal and spatial variations of SSL.The temporal and spatial variations of water runoff have also been analyzed to make a comparison.The temporal variation of water runoff and SSL show a decreasing trend, especially the SSL (Figure 6).

Availability of Sediment Sources
The availability of weathered sediments in the unsaturated channel can be revealed by rating parameters.Thus, the SSC-Q rating curves with data measured before and after TGD operation at the JJRR is adopted to investigate the hydrological condition.
According to data measured in 2013, values of the f-parameter are 2.27 (Shashi) and 1.67 (Jianli), respectively, and values of the log(e)-parameter are −7.70 (Shashi) and −10.70 (Jianli), respectively (Figure 8).The long-term average values by Yang et al. [39] are illustrated as data before dam operation (Figure 8).The long-term average values of the log(e)-parameter and the f-parameter for Xinchang station are −1.92 and 1.17, respectively, and the values for Jianli station are −0.97 and 0.96, respectively [39].Historical data adopted in the analysis of Yang et al. [39] were measured at four stations: Yichang , Xinchang (1956Xinchang ( -1984)), Jianli (1954Jianli ( -1986)), and Chenglingji .According to Figures 6 and 7, parameters by Yang et al. [39] may represent the hydrological conditions during the period before the changing point.Thus, with the operation of TGD, the values of the f-parameter increase while the values of the log(e)-parameter decrease.Instantaneous SST rates are not only a function of the transport capacity of a river, but also of sediment availability [22].For sediment rating parameters estimated through regression analysis, the e-coefficient is an erosion severity index in the river channel and is associated with the availability of weathered sediments in the basin area [22].A high e-value indicates that this area is characterized with easily erodible materials and high loads of materials transported by runoffs [20].The exponent f-coefficient corresponds to the erosive power of the river and the influence on the sediment supply from the entire basin surface.High values indicate a considerable increase in erosive power and sediment-carrying capacity with an increase in river discharge [20].Thus, a slight increase in discharge enhances the erosive power of the river significantly.For the JJRR, the increased value of the f-parameter indicates a considerable increase of erosive power and sediment-carrying capacity  Thus, after 2005, both the runoff and SSL are less than the long-term average value.

Availability of Sediment Sources
The availability of weathered sediments in the unsaturated channel can be revealed by rating parameters.Thus, the SSC-Q rating curves with data measured before and after TGD operation at the JJRR is adopted to investigate the hydrological condition.
According to data measured in 2013, values of the f -parameter are 2.27 (Shashi) and 1.67 (Jianli), respectively, and values of the log(e)-parameter are −7.70 (Shashi) and −10.70 (Jianli), respectively (Figure 8).The long-term average values by Yang et al. [39] are illustrated as data before dam operation (Figure 8).The long-term average values of the log(e)-parameter and the f -parameter for Xinchang station are −1.92 and 1.17, respectively, and the values for Jianli station are −0.97 and 0.96, respectively [39].Historical data adopted in the analysis of Yang et al. [39] were measured at four stations: Yichang , Xinchang (1956Xinchang ( -1984)), Jianli (1954Jianli ( -1986)), and Chenglingji .According to Figures 6 and 7, parameters by Yang et al. [39] may represent the hydrological conditions during the period before the changing point.Thus, with the operation of TGD, the values of the f -parameter increase while the values of the log(e)-parameter decrease.

Availability of Sediment Sources
The availability of weathered sediments in the unsaturated channel can be revealed by rating parameters.Thus, the SSC-Q rating curves with data measured before and after TGD operation at the JJRR is adopted to investigate the hydrological condition.
According to data measured in 2013, values of the f-parameter are 2.27 (Shashi) and 1.67 (Jianli), respectively, and values of the log(e)-parameter are −7.70 (Shashi) and −10.70 (Jianli), respectively (Figure 8).The long-term average values by Yang et al. [39] are illustrated as data before dam operation (Figure 8).The long-term average values of the log(e)-parameter and the f-parameter for Xinchang station are −1.92 and 1.17, respectively, and the values for Jianli station are −0.97 and 0.96, respectively [39].Historical data adopted in the analysis of Yang et al. [39] were measured at four stations: Yichang , Xinchang (1956Xinchang ( -1984)), Jianli (1954Jianli ( -1986)), and Chenglingji .According to Figures 6 and 7, parameters by Yang et al. [39] may represent the hydrological conditions during the period before the changing point.Thus, with the operation of TGD, the values of the f-parameter increase while the values of the log(e)-parameter decrease.Instantaneous SST rates are not only a function of the transport capacity of a river, but also of sediment availability [22].For sediment rating parameters estimated through regression analysis, the e-coefficient is an erosion severity index in the river channel and is associated with the availability of weathered sediments in the basin area [22].A high e-value indicates that this area is characterized with easily erodible materials and high loads of materials transported by runoffs [20].The exponent f-coefficient corresponds to the erosive power of the river and the influence on the sediment supply from the entire basin surface.High values indicate a considerable increase in erosive power and sediment-carrying capacity with an increase in river discharge [20].Thus, a slight increase in discharge enhances the erosive power of the river significantly.For the JJRR, the increased value of the f-parameter indicates a considerable increase of erosive power and sediment-carrying capacity  Instantaneous SST rates are not only a function of the transport capacity of a river, but also of sediment availability [22].For sediment rating parameters estimated through regression analysis, the e-coefficient is an erosion severity index in the river channel and is associated with the availability of weathered sediments in the basin area [22].A high e-value indicates that this area is characterized with easily erodible materials and high loads of materials transported by runoffs [20].The exponent f -coefficient corresponds to the erosive power of the river and the influence on the sediment supply from the entire basin surface.High values indicate a considerable increase in erosive power and sediment-carrying capacity with an increase in river discharge [20].Thus, a slight increase in discharge enhances the erosive power of the river significantly.For the JJRR, the increased value of the f -parameter indicates a considerable increase of erosive power and sediment-carrying capacity with an increase in river discharge.A decreased value of the log(e)-parameter indicates a decreased availability of suspended sediments in the river channel.
The saturation recovery process in the JJRR can be described with Figures 4 and 5, while the changed interactions between suspended sediment and bed material can be explained by the changed values of rating curves' coefficients (Figure 8).
For river reaches around Zhicheng station, coarsening of bed material is obvious.As to its limited distance from the TGD dam (60 km), the saturation recovery process of suspended sediment is limited.
For river reaches around Shashi station, coarsening of bed material is also obvious.Finer-grouped suspended sediments (particles with diameter < 0.008 mm) are transported from the upper river reach instead of local bed materials.Coarser suspended particles (with diameter of 0.062-0.05mm) have recovered, which are eroded from local bed materials as an increased value of the sediment rating exponent (f -parameter).
For river reaches around Jianli station, the coarsening of the bed material is limited.The coarsening of suspended sediment is obvious, especially after flooding (Figure 9).Coarser suspended particles (with diameter of 0.062-0.05mm) may come from the erosion of the local river bed and transportation from upper reaches by an increased value of the f -parameter.
The variability in the relationship between the sediment concentration and water discharge (namely hysteretic patterns) has also been used to explain the variability of sediment sources from one flood to another [40].Whitaker et al. [41] compared the estimated sediment yield for an un-sampled flood event, a single-event surrogate SSC-Q rating curve, or a long-term SSC-Q rating curve.Analysis revealed that the long-term SSC-Q rating curve estimates approximately 37% higher than the single-event SSC-Q estimates, which indicates a high degree of uncertainty.Analysis shows a high variability in the rating curve between similar floods [41].This is the reason for erosion during the flood season.
Water 2018, 10, x FOR PEER REVIEW 15 of 21 with an increase in river discharge.A decreased value of the log(e)-parameter indicates a decreased availability of suspended sediments in the river channel.
The saturation recovery process in the JJRR can be described with Figures 4 and 5, while the changed interactions between suspended sediment and bed material can be explained by the changed values of rating curves' coefficients (Figure 8).
For river reaches around Zhicheng station, coarsening of bed material is obvious.As to its limited distance from the TGD dam (60 km), the saturation recovery process of suspended sediment is limited.
For river reaches around Shashi station, coarsening of bed material is also obvious.Finergrouped suspended sediments (particles with diameter < 0.008 mm) are transported from the upper river reach instead of local bed materials.Coarser suspended particles (with diameter of 0.062-0.05mm) have recovered, which are eroded from local bed materials as an increased value of the sediment rating exponent (f-parameter).
For river reaches around Jianli station, the coarsening of the bed material is limited.The coarsening of suspended sediment is obvious, especially after flooding (Figure 9).Coarser suspended particles (with diameter of 0.062-0.05mm) may come from the erosion of the local river bed and transportation from upper reaches by an increased value of the f-parameter.
The variability in the relationship between the sediment concentration and water discharge (namely hysteretic patterns) has also been used to explain the variability of sediment sources from one flood to another [40].Whitaker et al. [41] compared the estimated sediment yield for an unsampled flood event, a single-event surrogate SSC-Q rating curve, or a long-term SSC-Q rating curve.Analysis revealed that the long-term SSC-Q rating curve estimates approximately 37% higher than the single-event SSC-Q estimates, which indicates a high degree of uncertainty.Analysis shows a high variability in the rating curve between similar floods [41].This is the reason for erosion during the flood season.

Accurate Measuring in Field Surveys
Accurate measurement of vertical concentration profiles is important for investigating the vertical distribution of SSC and the accurate estimation of sediment transport.For the field survey in the JJRR, the two near-bed points are measured with constant distances of 0.5 m and 0.1 m from the bed, respectively.However, the small distances may be doubtful due to difficulties in distinguishing the bed material and suspended sediment, especially when the bed is disturbed during settling sampling.Equipment and criteria adopted in the field survey have been described in detail in the Methodology section, and their availability and methods of improvement with regard to accuracy in the future are discussed.
The availabilities of ADCP and GNSS have been certified by field surveys conducted by other researchers.The error of hydrologic measurement with ADCP is ±1%, and the random uncertainty is

Accurate Measuring in Field Surveys
Accurate measurement of vertical concentration profiles is important for investigating the vertical distribution of SSC and the accurate estimation of sediment transport.For the field survey in the JJRR, the two near-bed points are measured with constant distances of 0.5 m and 0.1 m from the bed, respectively.However, the small distances may be doubtful due to difficulties in distinguishing the bed material and suspended sediment, especially when the bed is disturbed during settling sampling.
Equipment and criteria adopted in the field survey have been described in detail in the Methodology section, and their availability and methods of improvement with regard to accuracy in the future are discussed.
The availabilities of ADCP and GNSS have been certified by field surveys conducted by other researchers.The error of hydrologic measurement with ADCP is ±1%, and the random uncertainty is 2-5% for measuring cross-sections in the middle and lower YR channel [42].The integrated system (GPS and ADCP) delivers positioning with <15 cm accuracy vertically and <30 cm horizontally [42].For GNSS, the differences between coordinates of the ground control points estimated using RTK-GNSS, and those computed from the classical topographic measurements during processing of the photogrammetric blocks have been analyzed by Pagliari et al. [43].It was pointed out that an RTK-GNSS survey may be sufficient to reach the requested tolerance [43].GNSS receivers have been adopted in field surveys in recent years, as their main advantages are of being faster, cheaper, and easier to use during the surveying phase than classical topographic survey [43,44].For sediment grain-size analysis, the results by LPSA and sieving methods are comparable when particles are spherical [45].Sieving methods for coarser particles are more acceptable than for finer ones [45].Thus, the sieving method adopted by JHWRSB for fine suspended particles needs further verification.
In total, various efforts have being applied to measure the SSC profiles in field surveying, including methods and instrumentation [46][47][48].Some new field instrumentation, including magnetic tracers, drone-based or plane-based topography surveys, and video analysis by spectral cameras can also be used to measure sediment transport and morphology evolution in future field surveys [49].Except for these new instruments, already existing instruments can also be used with more details.For instance, ADCPs are increasingly used to measure the three-dimensional velocity distribution and sediment transport by measuring the Doppler shift in the backscattered signals from an array of acoustic beams [49,50].According to the backscatter intensity of ADCP beams, SSC can also be estimated with inversion [51].

Describing Vertical Concentration Profiles
The vertical distribution of SSC needs to be described accurately for an accurate simulation of channel migration, i.e., a one-dimensional simulation for the transport of sediment mixture in non-equilibrium conditions [52].Thus, finding a proper way to describe the vertical concentration profiles of the unsaturated channel is one of the key problems for hydraulic researchers [53].
Measured sediment profiles are compared with calculations by three equations, i.e., the Rouse Equation (Equation ( 5)), the Han Equation (Equation ( 7)), and the new fADE (Equation (12)).Figure 10 indicates the Rouse Equation is not applicable in the non-equilibrium reach.Compared to the Rouse Equation, the Han Equation for non-equilibrium, sediment-laden flow performs better in the upper part of a vertical but underestimates the concentration near the river bed, i.e., within 10% of the water depth.The fADE is able to characterize the dynamics of the non-equilibrium sediment suspension in "starving" river reaches.
No doubt, Equation (7) for non-equilibrium conditions can offer a better description of sediment profiles in river reaches with unsaturated or oversaturated flows.The inaccurate prediction is partly due to our poor understanding of the mechanism of sediment suspension, i.e., the spatially-stochastic transport behavior in sediment dispersion.For the changed sediment regime of the JJRR, more efforts are required to improve the description of vertical profiles of SSC, especially the physical processes.Additionally, "nonlocal" jumps in burst-like suspension events indicate that bursting eddies may sweep particles from the water bottom and directly eject them to the upper part of water bodies, which may lead to ecological-related problems.The ecology consequences include water quality of the riverine ecosystems, as nutrients and pollutants are often closely connected with sediments [49].The three outlets (Songzikou, Taipingkou, and Ouchikou) have a special setting of the surface water system between YR and Dongting Lake, and the three diversion channels could interact directly with the lake.Thus, it may also have a profound impact on the biological connection between the main stream and the downstream lakes of Poyang and Dongting.
equilibrium conditions [52].Thus, finding a proper way to describe the vertical concentration profiles of the unsaturated channel is one of the key problems for hydraulic researchers [53].
Measured sediment profiles are compared with calculations by three equations, i.e., the Rouse Equation (Equation ( 5)), the Han Equation (Equation ( 7)), and the new fADE (Equation ( 12)). Figure 10 indicates the Rouse Equation is not applicable in the non-equilibrium reach.Compared to the Rouse Equation, the Han Equation for non-equilibrium, sediment-laden flow performs better in the upper part of a vertical but underestimates the concentration near the river bed, i.e., within 10% of the water depth.The fADE is able to characterize the dynamics of the non-equilibrium sediment suspension in "starving" river reaches.No doubt, Equation (7) for non-equilibrium conditions can offer a better description of sediment profiles in river reaches with unsaturated or oversaturated flows.The inaccurate prediction is partly due to our poor understanding of the mechanism of sediment suspension, i.e., the spatially-stochastic transport behavior in sediment dispersion.For the changed sediment regime of the JJRR, more efforts are required to improve the description of vertical profiles of SSC, especially the physical processes.Additionally, "nonlocal" jumps in burst-like suspension events indicate that bursting eddies may sweep particles from the water bottom and directly eject them to the upper part of water bodies, which may lead to ecological-related problems.The ecology consequences include water quality of the riverine ecosystems, as nutrients and pollutants are often closely connected with sediments [49].The three outlets (Songzikou, Taipingkou, and Ouchikou) have a special setting of the surface water system between YR and Dongting Lake, and the three diversion channels could interact directly with the lake.Thus, it may also have a profound impact on the biological connection between the main stream and the downstream lakes of Poyang and Dongting.

Estimating SST Rate
Measuring and describing the vertical SSC profiles accurately is also imperative to the quantification of the SST rate [2].If near-bed concentrations are not measured or simulated correctly, the estimation of the SST rate may also be influenced.
Figure 11 shows the estimated contributions by two near-bed points.Both of these two extremely large ratios (larger than 60%) occurred at Jianli station in 1996.For these two sections, the estimated sectional average velocities are different from measured data with a large deviation.This indicates that the assumption for Equation ( 13) is not adoptable in these two cross-sections.Thus, these two sections are excluded from discussion. Figure 11 shows that, if the near-bed concentrations are not included, the sediment flux may be underestimated by approximately 23.5% (Zhicheng), 9.35% (Shashi), and 18.68% (Jianli), respectively.Thus, for the unsaturated JJRR, sediment concentration within 10% of the water depth from the riverbed cannot be ignored.Otherwise, it will lead to underestimated results when calculating the erosion/deposition of the JJRR by the sediment-

Estimating SST Rate
Measuring and describing the vertical SSC profiles accurately is also imperative to the quantification of the SST rate [2].If near-bed concentrations are not measured or simulated correctly, the estimation of the SST rate may also be influenced.
Figure 11 shows the estimated contributions by two near-bed points.Both of these two extremely large ratios (larger than 60%) occurred at Jianli station in 1996.For these two sections, the estimated sectional average velocities are different from measured data with a large deviation.This indicates that the assumption for Equation ( 13) is not adoptable in these two cross-sections.Thus, these two sections are excluded from discussion. Figure 11 shows that, if the near-bed concentrations are not included, the sediment flux may be underestimated by approximately 23.5% (Zhicheng), 9.35% (Shashi), and 18.68% (Jianli), respectively.Thus, for the unsaturated JJRR, sediment concentration within 10% of the water depth from the riverbed cannot be ignored.Otherwise, it will lead to underestimated results when calculating the erosion/deposition of the JJRR by the sediment-transport balance method and volume method.Additionally, an underestimated SST rate at Zhicheng station may also lead to overestimated sedimentation in the TGR, which is important for the water storage capacity of the reservoir.The reservoir's functions of water supply, energy production, navigation, flood control, and increased maintenance costs are correlated with the reservoir's water storage capacity [49].Thus, underestimated sedimentation at Zhicheng station may also affect the reservoir's waterway systems, hydraulic schemes, intakes and outlets, and flow regulation [49].Fortunately, the influence of near-bed concentration on estimating sediment flux can also be verified by comparing the estimated deposition with the sediment-balance method and the geometric method.
Based on the hydrology data and geometry data measured in 2002-2008, contributing factors that may cause deviation during estimating deposition were distinguished [54,55].Yuan et al. [53] pronounced that the contribution ratio by near-bed concretion varies in different hydrology years and river reaches (Table 4).For a certain hydrological station, the contribution ratios may also vary with the inflow condition.In total, they partially verified the reliability of measured vertical concentration profiles by comparing contribution ratios.
When suspension index A = ω/(βku * ) = 2.0, suspended sediment may concentrate in areas with a distance of 0.2 times of the water depth from the riverbed, in which β is a coefficient for non-uniform sediment [30].The near-bed concentration may affect the estimated sediment load only if A > 0.15 [30].If A = 2.0, the ratio of near-bed transportation over the total transportation is approximately 30% [30].The values of ω/u * with data measured after TGD in Table 2 are much larger than the criteria by Zhang [30], and a large contribution of the near-bed transportation over the total transportation in Figure 11 can be certified.The ratio of the near-bed transportation over the total transportation is approximately 25% when Di/Dm > 3, in which Di/Dm is the particle size the over average size [30].The wide grading distribution of SSC after TGD (as shown in Figures 4 and 5) also indicates a large contribution to the total transportation.Thus, the potential doubt that refers to difficulties in distinguishing bed material and suspended sediment can partially be certified.Moreover, the impact of a remarkable large near-bed concentration must be given more attention in the future, especially at river reaches around Shashi station.In total, fluvial and dynamic characteristics of distinct, large concentrations in near-bed regions may affect sediment flux, erosion/deposition estimation, sediment rating curve, and the channel geometry followed.This brings river engineers a lot of new problems, i.e., degradation of navigation, and dramatic bank erosion [56].Additionally, the estimation of erosion in the JJRR has a profound impact on the hydrologic connection between the middle YR and Dongting Lake [57].Fortunately, the influence of near-bed concentration on estimating sediment flux can also be verified by comparing the estimated deposition with the sediment-balance method and the geometric method.
Based on the hydrology data and geometry data measured in 2002-2008, contributing factors that may cause deviation during estimating deposition were distinguished [54,55].Yuan et al. [53] pronounced that the contribution ratio by near-bed concretion varies in different hydrology years and river reaches (Table 4).For a certain hydrological station, the contribution ratios may also vary with the inflow condition.In total, they partially verified the reliability of measured vertical concentration profiles by comparing contribution ratios.
When suspension index A = ω/(βku * ) = 2.0, suspended sediment may concentrate in areas with a distance of 0.2 times of the water depth from the riverbed, in which β is a coefficient for non-uniform sediment [30].The near-bed concentration may affect the estimated sediment load only if A > 0.15 [30].If A = 2.0, the ratio of near-bed transportation over the total transportation is approximately 30% [30].The values of ω/u * with data measured after TGD in Table 2 are much larger than the criteria by Zhang [30], and a large contribution of the near-bed transportation over the total transportation in Figure 11 can be certified.The ratio of the near-bed transportation over the total transportation is approximately 25% when D i /D m > 3, in which D i /D m is the particle size the over average size [30].The wide grading distribution of SSC after TGD (as shown in Figures 4 and 5) also indicates a large contribution to the total transportation.Thus, the potential doubt that refers to difficulties in distinguishing bed material and suspended sediment can partially be certified.Moreover, the impact of a remarkable large near-bed concentration must be given more attention in the future, especially at river reaches around Shashi station.In total, fluvial and dynamic characteristics of distinct, large concentrations in near-bed regions may affect sediment flux, erosion/deposition estimation, sediment rating curve, and the channel geometry followed.This brings river engineers a lot of new problems, i.e., degradation of navigation, and dramatic bank erosion [56].Additionally, the estimation of erosion in the JJRR has a profound impact on the hydrologic connection between the middle YR and Dongting Lake [57].

Conclusions
Sediment regime of the JJRR has changed dramatically with the operation of TGD.Based on daily and vertical concentration data at the three hydrological stations (Zhicheng, Shashi, and Jianli) in the JJRR, the sediment regime and its potential challenges are analyzed in this study.The conclusions can be summarized as follows: (1) In the unsaturated JJRR, measurements have revealed anomalous vertical profiles of SSC, as the near-bed concentrations normalized with the vertical average concentration are dramatically larger than that of the pre-equilibrium channel.The near-bed concentration (within 10% of the water depth from the river bed) may reach up to 15 times that of the vertical average concentration in the non-equilibrium channel.(2) In the unsaturated JJRR, the combination of wide grading suspended sediment and coarsened bed materials in non-equilibrium channel contribute to a remarkably large concentration in the near-bed zone.For the river reach around Shashi station, the remarkably large near-bed concentration is more apparent.(3) More detailed measurements or new measuring technologies may assist in the accurate measurement of vertical concentration profiles.A fractional dispersion equation may help to provide accurate descriptions.The outcomes can provide useful information for modeling the morphologic change of the JJRR.

Figure 3 .
Figure 3.Typical measured vertical profiles of SSC and vertical averaged concentration in the JJRR, measured at Shashi hydrological station before (c,d) and after 2010 (a,b).

Figure 3 .
Figure 3.Typical measured vertical profiles of SSC and vertical averaged concentration in the JJRR, measured at Shashi hydrological station before (c,d) and after 2010 (a,b).

Figure 6 .
Figure 6.Temporal variation of water runoff and SSL in the JJRR.

Figure 7 Figure 6 .
Figure 7 shows the mass curves of cumulative anomaly of annual runoff and SSL at the three stations (1987-2014).It shows that the water runoff and SSL show an increasing trend before 2000-2001 and then a decreasing trend after 2000-2001.Three Gorges Reservoir (TGR) was launched to construct in 1994, and TGD was completed and firstly impounded water in June 2003.The TGR started to operate regularly in June, 2006 and became fully operation in 2008.The other changing points of water runoff are 1993, 2005, and 2010, which are in accordance with the operation of TDR.The changing point of 1997 may be caused by the great 1998 flood, which was the highest in recorded history.Thus, after 2005, both the runoff and SSL are less than the long-term average value.

Figure 7
Figure 7 shows the mass curves of cumulative anomaly of annual runoff and SSL at the three stations (1987-2014).It shows that the water runoff and SSL show an increasing trend before 2000-2001 and then a decreasing trend after 2000-2001.Three Gorges Reservoir (TGR) was launched to construct in 1994, and TGD was completed and firstly impounded water in June 2003.The TGR started to operate regularly in June, 2006 and became fully operation in 2008.The other changing points of water runoff are 1993, 2005, and 2010, which are in accordance with the operation of TDR.The changing point of 1997 may be caused by the great 1998 flood, which was the highest in recorded history.

Figure 6 .
Figure 6.Temporal variation of water runoff and SSL in the JJRR.

Figure 7 Figure 7 .
Figure 7 shows the mass curves of cumulative anomaly of annual runoff and SSL at the three stations (1987-2014).It shows that the water runoff and SSL show an increasing trend before 2000-2001 and then a decreasing trend after 2000-2001.Three Gorges Reservoir (TGR) was launched to construct in 1994, and TGD was completed and firstly impounded water in June 2003.The TGR started to operate regularly in June, 2006 and became fully operation in 2008.The other changing points of water runoff are 1993, 2005, and 2010, which are in accordance with the operation of TDR.The changing point of 1997 may be caused by the great 1998 flood, which was the highest in recorded history.Thus, after 2005, both the runoff and SSL are less than the long-term average value.

Figure 8 .
Figure 8. Coefficients and exponents of sediment rating curves (a) for log(e) and (b) for f; solid lines are data drawn from Yang et al. [39].

Figure 8 .
Figure 8. Coefficients and exponents of sediment rating curves (a) for log(e) and (b) for f; solid lines are data drawn from Yang et al. [39].

Figure 8 .
Figure 8. Coefficients and exponents of sediment rating curves (a) for log(e) and (b) for f ; solid lines are data drawn from Yang et al. [39].

Figure 9 .
Figure 9. Cumulative volume fractions of suspended sediment at 2016: (a) measured in July 2016 and (b) data measured in September 2016.

Figure 9 .
Figure 9. Cumulative volume fractions of suspended sediment at 2016: (a) measured in July 2016 and (b) data measured in September 2016.

Figure 10 .
Figure 10.Comparison between measured and calculated vertical profiles of suspended sediment concentration.Data in (a) were measured on 10 August 2011 at Shashi station, data in (b) were measured on 26 September 2013 at Shashi station, and data in (c) were measured on 27 November 2017 at Jianli station.

Figure 10 .
Figure 10.Comparison between measured and calculated vertical profiles of suspended sediment concentration.Data in (a) were measured on 10 August 2011 at Shashi station, data in (b) were measured on 26 September 2013 at Shashi station, and data in (c) were measured on 27 November 2017 at Jianli station.

Figure 11 .
Figure 11.Contribution of the near-bed load on sediment flux.

Figure 11 .
Figure 11.Contribution of the near-bed load on sediment flux.

Table 2 .
Statistical characteristics of vertical profiles.

Table 3 .
Characteristics of vertical profiles with positive and negative gradations at Shashi and Jianli hydrological stations.

Table 3 .
Characteristics of vertical profiles with positive and negative gradations at Shashi and Jianli hydrological stations.

Table 4 .
Contribution by various factors on sediment flux 1 .

Table 4 .
Contribution by various factors on sediment flux 1 .Data drawn from Yuan et al. [53] and Duan et al. [54]; 2 total sediment flux is estimated with the geometry method.