Vertical Distribution of Suspended Sediments above Dense Plants in Water Flow

: Plants in natural water ﬂow can improve water quality by adhering and absorbing the ﬁne suspended sediments. Dense plants usually form an additional permeable bottom boundary for the water ﬂow over it. In the ﬂow layer above dense plants, the ﬂow velocity generally presents a zero-plane-displacement and roughness-height double modiﬁed semi-logarithmic proﬁle. In addition, the second order shear turbulent moment (or the Reynolds stress) are di ﬀ erent from that found in non-vegetated ﬂow. As a result, the turbulent momentum di ﬀ usivity of ﬂow and thus the di ﬀ usivity of sediment will shift, which will cause the vertical proﬁle of suspended sediment and the corresponding Rouse formula deform. A set of physical experiments with three di ﬀ erent diameters of ﬁne suspended sediments was conducted in an indoor water ﬂume. These experiments investigated a new distribution pattern of suspended sediment and the correspondingly deformed Rouse formula in the ﬂow layer over the dense plants. Experimental results showed that above the dense plants, the shear turbulent momentum of ﬂow presented a plant-height modiﬁed negative linear proﬁle, which has been proposed by a previous study, and the vertical distribution of ﬁne suspended sediments presented an equilibrium pattern. Based on the plant-modiﬁed proﬁles of ﬂow velocity and the shear turbulent momentum a new zero-plane and plant-height double modiﬁed Rouse formula were analytically derived. This double-parameter modiﬁed Rouse formula agrees well with the measured proﬁle of suspended sediment concentration experimentally observed in the present study. By adjusting the Prandtl–Schmidt number, i.e., the ratio of sediment di ﬀ usivity to ﬂow di ﬀ usivity, the double-parameter modiﬁed Rouse formula can be applied to submerged dense plant occupied ﬂow.


Introduction
Aquatic plants in natural water flow as are present in a river, lake or reservoir is usually used as a key phytoremediation element in eco-restoration projects because of its many benefits to an aquatic environment [1][2][3]. The beneficial effects of aquatic plants on the deposition, retention and filtration of suspended sediment in the deep part of a flow due to the plants' extensive root system has been widely acknowledged [4][5][6][7]. In the above-plant region, aquatic plants also influence the transport and movement of fine suspended sediment. These fine suspended sediments increase water turbidity, prevent photosynthesis in plants, affect the accessibility of food for larvae, and impact the quality of the aquatic habitat [8,9]. Consequently, understanding of the vertical distribution pattern of suspended sediments in the above-plant region is important for developing our understanding the aquatic eco-system.
In a natural water flow environment, aquatic plants usually grow in dense groups on the flow bed, thereby forming an additional permeable bottom boundary for the flow layer above it. The flow shear, i.e., the vertical gradient of flow velocity (du/dz) presents a remarkable two-layer feature in the within-plant region and the over-plant region [10][11][12][13]. In the within-plant region, the flow shear (du/dz) is approximately zero, as the velocity profile is constant. In the over-plant region, the flow shear (du/dz) no longer follows the traditional logarithmic law of the non-planted bare-bed flow, but is usually modeled by a deformed logarithmic law [14]. For dense submerged plant occupied flow, some studies have acknowledged that the traditional logarithmic law can be modified with the assistance of two parameters: the zero-plane-displacement (d) and the roughness height (z 0 ; Figure 1). Researchers have proposed various expressions for d and z 0 [15][16][17][18]. Stephan and Gutknecht (2002) [14] have given a detailed summary of these expressions. The d and z 0 modified flow shear (du/dz; u is the time-averaged flow velocity in streamwise direction) causes direct changes in the profile of the momentum diffusivity of flow (ε m = −u w /(du/dz) = τ/(ρ·du/dz), u and w are perturbations of flow velocity in streamwise and vertical directions), which will cause changes in the distribution pattern of suspended sediments (e.g., [19]).
Water 2020, 12, x; doi: FOR PEER REVIEW 2 of 13 distribution pattern of suspended sediments in the above-plant region is important for developing our understanding the aquatic eco-system. In a natural water flow environment, aquatic plants usually grow in dense groups on the flow bed, thereby forming an additional permeable bottom boundary for the flow layer above it. The flow shear, i.e., the vertical gradient of flow velocity ( / ) presents a remarkable two-layer feature in the within-plant region and the over-plant region [10][11][12][13]. In the within-plant region, the flow shear ( / ) is approximately zero, as the velocity profile is constant. In the over-plant region, the flow shear ( / ) no longer follows the traditional logarithmic law of the non-planted bare-bed flow, but is usually modeled by a deformed logarithmic law [14]. For dense submerged plant occupied flow, some studies have acknowledged that the traditional logarithmic law can be modified with the assistance of two parameters: the zero-plane-displacement ( ) and the roughness height ( ; Figure. 1). Researchers have proposed various expressions for and [15][16][17][18]. Stephan and Gutknecht (2002) [14] have given a detailed summary of these expressions. The and modified flow shear ( / ; is the time-averaged flow velocity in streamwise direction) causes direct changes in the profile of the momentum diffusivity of flow ( = − /( / ) = /( • / ), ′ and ′ are perturbations of flow velocity in streamwise and vertical directions), which will cause changes in the distribution pattern of suspended sediments (e.g., [19]). represents the time-averaged flow velocity in streamwise direction; represents the vertical coordinate with upward direction positive; and represents the integration constant.
Another important parameter that influences the turbulent diffusivity of flow ( ) and sediment ( ), and thus influences the distribution patter of suspended sediment is the second order shear turbulent momentum (− , which is usually used to represent the Reynolds stress, = − ). It is widely acknowledged that (− ) reaches a maximum value around the top of the plants, and then decreases upward towards the surface of the water [20][21][22]. Based on this feature, Li et al. (2018) [23] assumed that in the flow layer above the submerged plant, the shear stress follows the Vanoni profile ( = ( − ), is the energy slope; is the water depth), and that the flow velocity follows a -double modified logarithmic-law profile. Based on these assumptions, they then derived the parabolic profile of , which resulted in a profile of the -form of the suspended sediment concentration ( ( ) , is the concentration of suspended sediment). However, a further direct formula for suspended sediment concentration ( ) has not been derived due to the complexity of the expression of ( ).
Based on the generally upward-decreasing feature of the profile of , Huai et al. (2019) [24] assumed a zero-value of at the water surface, and then assumed a linear profile of between the water surface and the plant top. They successfully created a random displacement model for the prediction of the profile of the suspended sediment concentration profile ( ). Yang and Choi (2010) Another important parameter that influences the turbulent diffusivity of flow (ε m ) and sediment (ε s ), and thus influences the distribution patter of suspended sediment is the second order shear turbulent momentum (−u w , which is usually used to represent the Reynolds stress, τ = −ρu w ). It is widely acknowledged that (−u w ) reaches a maximum value around the top of the plants, and then decreases upward towards the surface of the water [20][21][22]. Based on this feature, Li et al. (2018) [23] assumed that in the flow layer above the submerged plant, the shear stress follows the Vanoni profile (τ = ρgS(H − z), S is the energy slope; H is the water depth), and that the flow velocity follows a d-z 0 double modified logarithmic-law profile. Based on these assumptions, they then derived the parabolic profile of ε m , which resulted in a profile of the ln-form of the suspended sediment concentration (ln(C s ), C s is the concentration of suspended sediment). However, a further direct formula for suspended sediment concentration (C s ) has not been derived due to the complexity of the expression of ln(C s ).
Based on the generally upward-decreasing feature of the profile of ε m , Huai et al. (2019) [24] assumed a zero-value of ε m at the water surface, and then assumed a linear profile of ε m between the water surface and the plant top. They successfully created a random displacement model for the prediction of the profile of the suspended sediment concentration profile (C s ). Yang

and Choi
Water 2020, 12, 12 3 of 13 (2010) [25] assumed ε m to be parabolic, and derived a parabolic Rouse formula for suspended sediment concentration, which is modified by the plant height.
The previous models and formulae for the prediction of suspended sediment concentration (C s ) are all based on the equilibrium mechanism of the dynamic balance between upward turbulent diffusion and downward gravitational settling, i.e., the convection-diffusion equation. Based on this equilibrium mechanism, Rouse (1937) [26] derived the Rouse formula with the turbulent diffusion being expressed by Fick's first law. Fredsøe and Deigaard (1992) [27] derived the Vanoni formula with the turbulent diffusion being expressed by a mixing length involved term. These two formulae have the same expression. For flow through submergent plants, several different deformed Rouse formulae have been derived under some particular flow, sediment and plant conditions. However, investigation of the profile of suspended sediment in the flow layer above aquatic plants is still limited. Correspondingly, models to predict sediment and plant in this flow layer are also limited.
Fine sediment is commonly observed to be deposited on the foliage of submerged dense plants after a flood or storm. However, the mechanism by which he foliage-deposited sediment diffuses or is transported in the steady high-level flow after a flood has rarely been reported. The objectives of the present study are: (1) to investigate the profile of suspended sediment concentration (C s ) in dense plant occupied flow under a high submergence ratio, with suspended sediment coming from upstream instead of coming upwards from the flow bed; (2) to seek suitable formulae of the profiles of flow shear (du/dz) and the shear turbulent momentum (−u w ); and (3) to derive an analytical formula of the profile of suspended sediment concentration (i.e., a deformed Rouse formula). The results are expected to compensate for the limited data of C s -profile over dense plants, and to provide fundamental guidance for the design of phytoremediation of eco-restoration projects in natural water flow.

Water Flume and Co-Ordinate System
A rectangular water flume was used to conduct a set of indoor physical experiments. The flume was 8 m long, 40 cm wide, and 40 cm high (Figures 2 and 3). To conveniently observe the configuration of plants, the water surface, and the concentration of fine suspended sediments, large pieces of transparent glass were embedded in aluminum frames for two side-walls of the water flume. A 1 m × 1 m × 0.5 m water tank was set up at the upstream end of the flume. The water was supplied from a flat nozzle, and the polymer-sponge mattress was placed on the water surface to minimize the fluctuation of the water surface in the tank ( Figure 3). The water tank was installed on a bracket, so that the bottom of the water tank was 1.0 m above the bottom of the flume ( Figure 3). The water tank was connected to the water flume with an upper-surface-open corridor. The connecting corridor was 38 cm wide, 40 cm high, and had a 35 • tilt with respect to the horizontal. A valve was set at the upstream end of the corridor to adjust flow discharge in the corridor. Flow in the flume was driven by the water head between the water surface in tank and that in the flume. For all the experiments, water depth in the tank was maintained at 0.35 m. The water in the tank was supplied from a pipe. To adjust flow depth (H) and measure flow discharge (Q), a triangular wire was set at the downstream end of the flume ( Figure 3). An outlet pool with an outlet valve was connected to the flume at the downstream edge of the flume (Figures 2 and 3).
valve was set at the upstream end of the corridor to adjust flow discharge in the corridor. Flow in the flume was driven by the water head between the water surface in tank and that in the flume. For all the experiments, water depth in the tank was maintained at 0.35 m. The water in the tank was supplied from a pipe. To adjust flow depth ( ) and measure flow discharge ( ), a triangular wire was set at the downstream end of the flume ( Figure 3). An outlet pool with an outlet valve was connected to the flume at the downstream edge of the flume (Figures 2 and 3). The co-ordinate system is shown in Figures 2 and 3. The origin point in the co-ordinate system was located at the intersection of the upstream edge and the streamwise center line of the flume bottom.

Plant Configurations
The plants were made of plastic. This allowed an imitation of a natural plant, since the material was flexible enough to allow the artificial plants to sway slightly in flow. The plant was designed as triple foliage, with each foliage having a thin plate shape. The plant foliage was 3.2 cm high and 1.2 mm thick ( Figure 4a). In the present study, the time-space double averaged height of swaying plant ℎ = 3.1 cm. Therefore, the submergence ratio of plant /ℎ > 5, which allowed the dense plants to bed area [28] is an important factor used to classify the density of canopy. Nepf (2012) [  The co-ordinate system is shown in Figures 2 and 3. The origin point in the co-ordinate system was located at the intersection of the upstream edge and the streamwise center line of the flume bottom.

Plant Configurations
The plants were made of plastic. This allowed an imitation of a natural plant, since the material was flexible enough to allow the artificial plants to sway slightly in flow. The plant was designed as triple foliage, with each foliage having a thin plate shape. The plant foliage was 3.2 cm high and 1.2 mm thick ( Figure 4a). In the present study, the time-space double averaged height of swaying plant h p = 3.1 cm. Therefore, the submergence ratio of plant H/h p > 5, which allowed the dense plants to form a permeable bottom boundary for flow.
The plants were installed on the flume bed in a staggered distribution pattern (Figure 4b), with the space between two neighbor plants S = 2 cm. The plants occupied the entire flume bed from x = 0 m to the downstream end of flume. The frontal projected area of submerged vegetation per bed area λ f [28] is an important factor used to classify the density of canopy. Nepf (2012) [11] proposed a reference value, λ f = 0.1, to identify the dense canopy needed to generate inflected velocity profile based on experimentally measured data. In the present study λ f = 0.775, which is in the range of dense plants. In the present study, a high submergence ratio, combined with the dense distribution of plants on the bottom of the flow, ensured that the plants formed a permeable boundary. Therefore, flow velocity could be modeled as a modified logarithmic law.
Water 2020, 12, 12 5 of 13 = 0 m to the downstream end of flume. The frontal projected area of submerged vegetation per bed area [28] is an important factor used to classify the density of canopy. Nepf (2012) [11] proposed a reference value, = 0.1, to identify the dense canopy needed to generate inflected velocity profile based on experimentally measured data. In the present study = 0.775, which is in the range of dense plants. In the present study, a high submergence ratio, combined with the dense distribution of plants on the bottom of the flow, ensured that the plants formed a permeable boundary. Therefore, flow velocity could be modeled as a modified logarithmic law.

Bulk Flow Condition and Fine Suspended Sediment Model
Bed slope (S b ) of the flume was adjusted at 0.38% to keep the experimental flow in a quasi-steady condition. During the total four sets of experiments in the present study, the flow discharge (Q) was kept at a constant value of 0.0096 m 3 /s, and the flow depth (H) was 20 cm ( Table 1).
The fine suspended sediment was modeled as spherical balls made of plexiglass. Three different diameters of fine sediments were used: d s = 36, 59, and 80 µm. Three different specific gravities were used for these different diameters of sediment ( Table 2). The concentration profile of the fine suspended sediments was measured at x = 3.5 m (Figure 3). The fine sediment was injected into the experimental flow from the upstream edge of the flume. The center of the injection nozzle was aligned with the center of cross direction of the flume, at a height equal to that of the top of the plants, i.e., at z = 3.2 cm.
Sampling of suspended sediments was taken by a 95 mL sampling bottle over 35 s. Three sampling runs were conducted, with the first run starting from 5 min, the second run from 30 min and the third run from 45 min after the flow surface line became steady. In each sampling run three rigid rods were separately applied in order. Each rod had four sampling bottles attached. The spacing between the bottles was 5 cm for Rod_1 and Rod_3, and 6 cm for Rod_2 ( Figure 5). Thus 12 measurement points along a vertical line were obtained. The volume concentration (g/L) of sediment was obtained by the filtering method. This sampling bottle was calibrated with the standard sampling bottle.

Bulk Flow Condition and Fine Suspended Sediment Model
Bed slope ( ) of the flume was adjusted at 0.38‰ to keep the experimental flow in a quasi-steady condition. During the total four sets of experiments in the present study, the flow discharge ( ) was kept at a constant value of 0.0096 m 3 /s, and the flow depth ( ) was 20 cm ( Table 1).
The fine suspended sediment was modeled as spherical balls made of plexiglass. Three different diameters of fine sediments were used: = 36, 59, and 80 μm. Three different specific gravities were used for these different diameters of sediment ( Table 2). The concentration profile of the fine suspended sediments was measured at = 3.5 m (Figure 3). The fine sediment was injected into the experimental flow from the upstream edge of the flume. The center of the injection nozzle was aligned with the center of cross direction of the flume, at a height equal to that of the top of the plants, i.e., at = 3.2 cm. Sampling of suspended sediments was taken by a 95 mL sampling bottle over 35 s. Three sampling runs were conducted, with the first run starting from 5 min, the second run from 30 min and the third run from 45 min after the flow surface line became steady. In each sampling run three rigid rods were separately applied in order. Each rod had four sampling bottles attached. The spacing between the bottles was 5 cm for Rod_1 and Rod_3, and 6 cm for Rod_2 ( Figure 5). Thus 12 measurement points along a vertical line were obtained. The volume concentration (g/L) of sediment was obtained by the filtering method. This sampling bottle was calibrated with the standard sampling bottle.

Measurement of Flow Velocity and the Observation of Plant Swaying
To obtain time-averaged mean and turbulent flow characteristics of flow, three-dimensional time series of flow velocity was observed by Acoustic Doppler Velocimeter (ADV). Time series of 3-D flow velocities along a vertical profile were measured in the middle section at x = 3.5 m. The time-averaged flow velocities u and w (u = U − u , w = W − w, U and W are the instantaneous flow velocities in streamwise and vertical directions, respectively; u and w are the corresponding velocity perturbations) and the second order turbulent momentum (−u w ) were calculated from the measured time series of velocity. The recording frequency of the ADV was f R = 25 Hz.
The recorded signal was filtered and de-spiked by the acceleration method and the phase-space threshold method [29]. We fed particles into the flow; thus, the spikes in the measured time series were not significant. For all the experimental tests, the filtered and de-spiked valid digital records of each time series were greater than 85%. The correlation coefficient falling in the range of [80, 95], and the ratio of signal to noise (SNR) was in the range of [15.0, 26.3] dB.
A laser emitter was hung over the flume to generate a thin sheet of laser light with a thickness of 1.5 mm in the x − z plane along the longitudinal centerline of the flume (Figure 3). A video camera was fixed outside the glass window of the flume with the lens normal to the center of the laser light sheet (Figure 3). During each testing period, the dynamic movement of the plants was recorded by the camera. Thus, the time-averaged height of the plants was obtained by analyzing the sequential video via the computer.

Diffusivities of Momentum of Flow and Fine Suspended Sediments
The momentum diffusivity (i.e., diffusion coefficient) of flow is defined as ε m ≡ −u w /(du/dz). (1) The diffusivity of fine suspended sediments is defined as where β is the Prandtl-Schmidt number. A β-value of unity is widely used. In applications, the Prandtl-Schmidt number is dependent on flow conditions as well as sediment and plant conditions.

Settling Velocity of Fine Suspended Sediments
The settling velocity of fine sediments (ω s ) is an important parameter for the determination of its distribution pattern. Since the dynamics of fine suspended sediments is influenced by turbulent dissipation, the modified formula of ω s under turbulent flow condition was applied [30]: where ω stokes = (1/18)[(r s − r)/r] gd 2 /ν is the Stokes settling velocity of fine sediment in still water; r is the specific gravity of water, r = 1; v is the kinematic viscosity of water, v = 1.004 × 10 −6 m 2 s −1 at 20 • C; γ is a modification coefficient, γ = 0.034; and k is the Karman constant, the standard value k = 0.4 is used in the present study. The settling velocities of three fine sediments are given in Table 2.
For the three diameters of sediment used in the present study the ratios of u * /ω s are greater than unity ( Table 2), suggesting the fine sediments can be suspended [31,32].
Since the dense plants act as a new permeable boundary for the above flow, the shear velocity (u * ) was calculated from the value of −u w at the top of the plants, which is expressed by: Figure 6 shows that with the dense plants acting as permeable bottom boundary, the over flow velocity follows a semi-logarithmic law modified by d and z 0 :

Vertical Velocity Profile
where k is the von Karman constant (the standard value k = 0.4 was used in the present study); ρ is the density of water; d is the zero-plane displacement; z 0 is the roughness height; and C 0 is a constant. where is the von Karman constant (the standard value = 0.4 was used in the present study); is the density of water; is the zero-plane displacement; is the roughness height; and is a constant. In applications, the values of , , and depend on the specified flow environment and plant configuration. For example, Stephan and Gutknecht (2002) [14] reported that both and take the average value of swaying plant height, that is, /ℎ = 1.0, /ℎ = 1.0, and = 8.5; Nezu and Sanjou (2008) [22] reported that the value of  Figure 7 shows the vertical profile of the second order of turbulent momentum − ′ ′ . Over the canopy, where /ℎ ≥ 1, the value of − ′ ′ reaches a maximum value near the top of the plants, then varies upwards towards the water surface with elevation in a negative linear pattern. In In applications, the values of d, z 0 , and C 0 depend on the specified flow environment and plant configuration. For example, Stephan and Gutknecht (2002) [14] reported that both d and z 0 take the average value of swaying plant height, that is, d/h p = 1.0, z 0 /h p = 1.0, and C 0 = 8.5; Nezu and Sanjou (2008) [22] reported that the value of z 0 /h p falls within the range of [0.15, 0.17] and C 0 = 0 in their experiments.

Second Order Turbulent Momentum
In the present study, d/h p = 1.026, z 0 /h p = 0.15, and C 0 = −0.02. Therefore, the formula for the present study is: Water 2020, 12, 12 8 of 13 Figure 7 shows the vertical profile of the second order of turbulent momentum −u w . Over the canopy, where z/h p ≥ 1, the value of −u w reaches a maximum value near the top of the plants, then varies upwards towards the water surface with elevation z in a negative linear pattern. In general, this pattern agrees with other reports (e.g., [33]). This pattern implies that the dense plants can be considered as rough permeable bottom boundary for the over flow. The experimental data shows a good agreement with a plant-height modified linear formula, which was proposed by Nezu and Sanjou (2008) [22]:

Second Order Turbulent Momentum
It should be noted that near the water surface, the calculated value of −u w from Equation (7) is smaller than the observed data. This is mainly because the fluctuation of the water surface enhances the turbulence, thereby increasing the corresponding turbulent momentum.

Derivation of Turbulent Diffusivity of Flow
Based on Equations (5) and (7), the turbulent momentum diffusion coefficient ( ) can be derived. Taking the derivative of with respect to in Equation (5), yields or: For the convenience of integration, Equation (9) can be re-written as: Re-arranging the right hand of Equation (10) yields: i.e., = * ( ).

Derivation of Turbulent Diffusivity of Flow
Based on Equations (5) and (7), the turbulent momentum diffusion coefficient (ε m ) can be derived. Taking the derivative of u with respect to z in Equation (5), yields or: For the convenience of integration, Equation (9) can be re-written as: Re-arranging the right hand of Equation (10) yields: Water 2020, 12, 12 9 of 13 i.e., Substituting Equations (12) and (7) into ε m ≡ −u w / du dz yields i.e.,

Equilibrium Equation between Turbulent Diffusion and Gravitational Settling of Fine Suspended Sediments
The net flux of suspended sediments per unit area of a plane parallel to flow bed is described by where C s is the volumetric concentration of fine sediments. Substituting of Equation (14) for k m in Equation (15) yields The downward flux of suspended sediments per unit area due to gravity is described by In the dynamic equilibrium state, F eddy + F gravity = 0, that is, Re-arranging Equation (18) in term of dC s /C s yields i.e., Integrating the left-and right-hand sides of Equation (20) yields Since (z − H) < 0, and (z − d) > 0 over the plants, Equation (21) can be re-written as Choosing a reference concentration C a as the integrating constant C as a boundary condition of Equation (22): to make C s = C a at the height level where z = a. Substituting Equation (23) into Equation (22), and then solving Equation (22) yields Equation (24) indicates that the concentration profile of fine suspended sediments is independent of the roughness height z 0 , neither of the constant C in the modified semi-logarithmic law of flow velocity expressed by Equation (5). Equation (24) also indicates that over dense plants, an equilibrium state of fine suspended sediments under the dynamic balance between turbulent diffusion and gravitational settling can be obtained. Comparison between experiments using three diameters of suspended sediments is also shown in Figure 8b. For the smallest diameter (i.e., the least ω s ), although the C a -value was the smallest (Table 2), the C s /C a -value at all higher locations were the largest. For the largest diameter (i.e., the greatest ω s ), the C a was the largest (Table 2), while the C s /C a -value at higher locations was the smallest. For the medium diameter, the concentrations at all height levels lay between the values for the smallest and largest diameters. This phenomenon indicates that given the same flow condition, the effect of turbulent diffusion on the suspension of fine sediments enhanced the decrease of settling velocity of sediment (ω s ), while the gravitational settling effect became weak with the decrease of ω s .  Figure 8b). Comparison between experiments using three diameters of suspended sediments is also shown in Figure 8b. For the smallest diameter (i.e., the least ), although the -value was the smallest (Table 2), the / -value at all higher locations were the largest. For the largest diameter (i.e., the greatest ), the was the largest (Table 2), while the / -value at higher locations was the smallest. For the medium diameter, the concentrations at all height levels lay between the values for the smallest and largest diameters. This phenomenon indicates that given the same flow condition, the effect of turbulent diffusion on the suspension of fine sediments enhanced the decrease of settling velocity of sediment ( ), while the gravitational settling effect became weak with the decrease of .

Discussion
The newly derived Rouse formula, which is modified by the zero-plane displacement (d) and the plant height (h p ) was validated by the observed profile of suspended sediment concentration (C s ) in the present study. This d-h p double modified Rouse formula is particularly suitable for describing dense plants with a high submergence ratio where the d-h p double modified logarithmic law of flow shear (Equation (5)) and the negative linear profile of shear turbulent momentum (Equation (7)) are applicable.
The Prandtl-Schmidt number (β) is influenced by different turbulence structures. In dense plant occupied flow, the characteristics of turbulence, including the turbulent kinetic energy, the skewed velocity and the intermittent features might have influence on the β-value. However, an accurate method for the determination of β-value has not been proposed. Therefore, it needs to be adjusted according to the conditions of the specified application.
In practical applications, the shear flow velocity u * can be estimated by the widely used plant-height modified Vanoni formula u * = gS H − h p [22,23,25,34,35]. For an after-flood flow, the suspended sediments in the above-plant region mainly come from the near-plant-top region, therefore a near-plant-top height level can be chosen for the reference a-value.
In applications, error in the estimation of u * , k, β, ω s , h p , or d will lead to an error in suspension number, i.e., the exponent of the double-parameter modified Rouse formula Equation (24), and thus lead to an error in the estimation of the profile of suspended sediment concentration.
Some previous studies have attempted to make more complete model for the profile of suspended sediments [25]. However, for the dense plant occupied flow, which is usually presented in bio-restoration projects, turbulent diffusion is not the dominant factor for the vertical distribution of suspended sediments due to the direct interference and obstacle effect of plant foliage on sediment. In addition, the movement of sediment in the deep part of plant is limited within the plant foliage or tends to deposited on the flow bed by the retention effect of the plant foliage or roots. Therefore, separate investigation of the over-plant region is required.

Conclusions
A set of indoor experiments was conducted to investigate the vertical distribution of suspended sediments over dense plants in flowing water. The main conclusions are: (1) Over dense plants the velocity profile follows a semi-logarithmic law modified by d and z 0 .
The second order turbulent momentum (−u w ) follows a negative linear function of z and H − h p . Based on these two laws, a new function of turbulent diffusion coefficient (ε m ) was derived.
(2) Based on the dynamic balance between the upward turbulent diffusion and the downward gravitational settling, a d-h p double modified Rouse formula for the vertical distribution profile of fine suspended sediments was derived. The equation agrees well with experimentally observed data in the over-plant region.

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