Dual Behavior of Suspended Sediment Concentration in a Regulated River ()
1. Introduction
Better understanding on the transport of suspended sediment from rivers to the ocean has remained a challenge for an increasing number of reasons. In Japan, the overall erosion rate of the coastline was 7.2 × 105 m2/year before 1980, but sharply increased to 1.6 × 106 m2/year after 1980 due largely to the reduction of suspended sediment supply. In the United States, Environmental Protection Agency (EPA) identifies suspended sediment as the single most widespread pollutant affecting the beneficial uses of the Nation’s rivers and streams [1]. High concentration of suspended sediment may clog fish gills, either killing them or reducing their growth rate. When the water slows down, the suspended sediment settles out and drops to the bottom, a process called siltation. This causes the water to be clear, but as the sediment settles it may change the bottom smothering bottom-dwelling organisms, covering breeding areas, and smothering eggs. Therefore, sound management of suspended sediment transport is a multi-purpose undertaking.
The relationship between river discharge and suspended sediment concentration is generally expressed by a rating curve, which is described by a power function (Walling [2], Sickle and Beschta [326], Steegen et al. [4], Morehead et al. [5]).
(1)
in which SSC is the suspended sediment concentration (g/m3), Q is the flow discharge (m3/s) and a, b are empirically determined regression coefficient. However, various studies (Klein [6], Seeger et al. [4], Williams [7]) showed that the relations between SSC and Q are highly variable either within or between events. Resulting SSC-Q hysteresis were widely examined at the event scale to interpret geomorphic processes and outline the spatial distribution of sediment sources (Crawford [8], Lenzi and Marchi [9], Jansson [10], Krueger [11], Duvert [12]). Based on these studies, the SSC–Q relationships can be mainly classified into three classes. First class: peaks of SSC and discharge arrive simultaneously; Second class: the SSC peak arrives earlier than the discharge peak; Third class: the SSC peak arrives later than the discharge peak. The mechanism behind the second class may be the remobilization of in-channel sediment deposits or adjacent extra-channel sources where sediment is transported a short distance by runoff (Williams [7], Jansson [10]). The third class may be related the arrival of remote sediment from external channel sources during long duration rainfall events or correspond to the absence of in-channel sediment delivery (Brasington and Richards [13], Lenzi and Marchi [9], Orwin [14]). However, little study has been done on the suspended sediment hysteresis in relation to sediment carrying capacity and channel morphology.
In this study, the patterns of SSC-Q relation during three floods of different magnitude in a regulated gravel river were studied. It is aimed to better understand the hysteresis in relation to flow and channel cross-sectional characteristics and from the perspective of channel sediment carrying capacity.
2. Materials
2.1. Study Site
The Tenryu River (Figure 1) originates in the central Honshu Mountains. From its source at Lake Suwa in Okaya, Nagano Prefecture, it extends 213 km, grazing Aichi Prefecture en route to its mouth at the boundary between the cities of Hamamatsu and Iwata in Shizuoka Prefecture. It drains 5050 km² into the Enshunada Coast. In length, it is the 9th longest river in Japan. There are five dams along the mainstream; Yasuoka Dam, Hiraoka Dam, Sakuma Dam, Akiha Dam and Funagira Dam. The Lower Tenryu River is the part from the Kashima site to the river mouth, approximately 25 km long. Due to dam constructions in the upstream of the river, particularly after the completion of the Sakuma Dam in 1956, the
sediment supply from the upstream was very much reduced, resulting in significant morphological change in the downstream reaches and beach erosion in the Enshunada Coast as well (Nagashima [15], Liu [16]). By analyzing the riverbed aggradations data in the reservoir formed by the Sakuma Dam, it was found that the total amount of the sediment trapped in the reservoir during the period from 1956 to 2000 was about 1.2 × 108 m3. Thus, the annual sediment trapping rate by the Sakuma Dam was estimated to be about 2.6 × 106 m3/yr, which virtually prevented all bed material sediment loads from transporting to downstream reaches. As a result, the active remaining source of beach-forming sediment supply to the coast after the dam construction is the Keta River, which is a tributary of the Tenryu River and located downstream of all sediment-trapping dams but upstream of the Kashima site.
2.2. Data
Two large floods occurred in the river in 2006 and 2007, respectively. Hamamatsu River Office conducted fieldobservation of SSC and flow discharge during the two floods. The depth-averaged SSC data, supplied to authors from the river administration office, were used in the analysis of this study. We also collected from the river administration office channel cross-sectional data and continuous water level records for the calculation of flow velocity and depth. In 2008, we installed a back scattering turbidity meter (Compact-HTW) 1m above the river bed at the Kakedsuka gauging station (3.3 km from the river mouth) for continuous monitoring. During the summer of 2008, there was a small flood and the corresponding turbidity data were recorded and converted to SSC according to the calibrated relationship between turbidity and SSC compiled through pre-field lab test.
2.3. Methods
First of all, flood frequency analysis was conducted for the Lower Tenryu River to determine the return periods of the three floods in consideration. The Log-Pearson Type III distribution was used in the analysis.
To explain the difference among the three cases in terms of carrying capacity, we adopted the energy balance theories of Velikanov [17], Bagnold’s [18] and Zhang ([19,20]), which assumed that (1) the amount of energy supplied by the fluid equals frictional energy losses and that required to keep sediment in suspension; (2) the energy supplied by sediment equals frictional energy losses by the sediment.
By dimensional analysis, the difference in energy dissipation rate between clear water and sediment-laden flow could be expressed as
(2)
where A = channel cross-sectional area; g, gs = specific weights of water and sediment; w = settling velocity; Cv = suspended sediment concentration by volume, k1 and a = coefficients.
By definition,
(3)
where J and Js are energy slopes of clear water and sediment-laden flow, respectively. They can be expressed as following
(4)
(5)
where l and ls are friction coefficient without and with suspended sediment, respectively. Mathematical manipulation of Equations (2-5) led to
(6)
Since the difference between l and ls is due to the presence of suspended sediment, (l - ls) in Equation (6) could be assumed to be a function of sediment concentration as
(7)
Substitution of Equation (7) into Equation (6) and mathematical manipulation led to
(8)
where C* can be considered as the saturation concentration of suspended sediment (kg/m3) under the velocity U; R is the hydraulic radius; K and m are empirical parameters. Following the work of Zhang ([19,20]), which was based on flume experiments and field observation, the values of K and m were assigned to be 0.07 and 1.4 in this study. Then, the transport capacity of suspended load was estimated in the present study as
(9)
To make use of this formula, the settling velocity must be determined. In the autumn of 2006, we surveyed surficial sediment size distributions in a vegetated zone on the floodplain 14 km from the river mouth. The type of vegetation was hard tree with an average height of 5 m. After the flood of 2007, we survey surficial sediment size distributions on a vegetated sandbar 5 km from the river mouth. Figure 2 shows the measured sediment size distributions on the vegetated flood plain and on the vegetated sandbar as well. Because trees trap suspended sediment, the surface deposits on the flood plain and sand bar reflect what was in suspension during the last
Figure 2. Size distribution of surficial deposits.
flood. It can be seen from Figure 2 that surface deposits were mainly in the range of 0.1 mm – 1 mm in vegetated zones. Therefore, the mean particle size of suspended sediment (Dm) was taken to be 0.3 mm. The settling velocity was then computed according to Dietrich [21].
(10)
To characterize the three floods in terms of their carrying capacities, the saturation concentrations at the peak discharges of the three floods were computed by Equation (9).
3. Results and Discussions
Figures 3-5 show flow discharges and SSC for the three target flood events. The peak flows were approximately 8000, 6000 and 1000 m3/s, respectively, while the maximum concentrations of suspended sediment were approximately 1400, 1200 and 110 g/m3, respectively. The hydrograph of the 2007 flood was relatively symmetrical inserting suggesting no hysteresis between discharge and water level within this flood event.
A four parameter Pearson Type-Ill distribution given below has been used as the form of inflow hydrograph for testing flood routing methodologies by various investigators (Ponce and Theurer, [22]; Garbrecht and Brunner [23]).
(11)
where Ib is the initial steady flow, Ip is the peak flow, t is the time to peak, and c is the skewness factor. The symmetry exhibited in the 2007 flood can be generated by
Figure 3. Flow discharge and suspended sediment concentration at Kashima station on July 15, 2007.
Figure 4. Flow discharge and suspended sediment concentration at Kashima station on July 19, 2006.
Figure 5. Flow discharge and suspended sediment concentration at Kakedsuka station on July 4, 2008.
Equation (11) if the value of skewness factor is given less than 1.15.
According to the flood frequency curve (Figure 6), the flood in 2006 was a 10-year flood and the 2007 flood was a 30-year flood. Therefore, they may be ranked as large and very large, respectively. For the 30-year flood, the variation of SSC followed the hydrograph well having its peak slightly ahead of the peak flow. For the 10-year and small floods, however, the peaks of SSC lagged behind the peaks of discharge. By using surveyed cross-sectional data, the maximum cross-sectional average velocities at the Kashima station for the 2006 and
Figure 6. Flood frequency in the Lower Tenryu River.
2007 floods were calculated to be 4.0 m/s and 4.9 m/s, respectively. The maximum velocity at the Kakedsuka station during the small flood of 2008 was 0.8 m/s. Thus, the saturation levels of SSC at monitoring sites for the three floods were estimated by Equation (9) to be approximately 2700 g/m3, 1300 g/m3 and 100 g/m3, respectively. It is clear that the observed maximum suspended concentrations of the 10-year flood in 2006 and the small flood in 2008 were around their thresholds. However, during the 30-year flood, SSC was just 55% of the carrying capacity. Following the concept of Julien [24], the sediment transport process can be described as being transport-limited mode in the 2006 and 2008 floods, but a supply-limited mode in the 2007 flood.
A characteristic of the Kashima site is that there is a point bar on the left side of the channel (Figure 7) so that it resembled a two-stage channel on the left side. Figure 8 shows the water level variations for the 2006 and 2007 floods, respectively. Compared to the 2006 flood, the 2007 flood had steeper rising limb. Because of the channel cross-sectional configuration and the characteristics of water level variation, the 2006 flood inundated the point bar shallowly. In a two-stage channel, horizontal vortices around the interface between the main channel and the flood plain, lateral vortices above the flood plain and longitudinal vortices in the main channel may be generated as illustrated in Figure 9 (Ishigaki and Imanoto [25], Knight and Shiono [26]). The experiment by Ikeda et al. [27] showed that the large-scale vortices appeared when the water depth on the floodplain is relatively shallow. Such large-scale vortices can cause re-suspension of sediment deposited on the flood plain prior to a flood (Ikeda et al. [28]). In the 2006 flood, because the water level stopped rising further after shallowly inundated the point bar, coherent vortices must be present and dominant. Figure 10 shows the relationship between measured velocity and water level during this flood. As it is revealed, the maximum cross-sectional average velocity appeared as the water just overflowed the point bar, resulting in a time lag between maximum velocity and highest water level. Since there was no de
Figure 7. Point bar at the Kashima site.
Figure 8. Water level variations in the 2006 and 2007 flood events.
Figure 9 Flow structure in a two-stage channel (modified from Knight and Shiono)
layed arrival of sediment from tributary in the 2006 flood, the hysteresis can be considered as a mechanism by which the peak flow velocity entrained in-channel sediment deposits and the enhanced cross-sectional mixing due to channel configuration sustained the suspended sediment concentration during the falling limb of the flood.
During the 2007 flood, the water level rose quickly through the level of point bar and continued to rise to a depth of 2.8 m over the point bar. Therefore, the time lag between maximum velocity and highest water level was negligible and the lateral momentum exchange was insignificant in the 2007 flood. In other words, the flood in 2007 behaved as like a flood in a single section channel.
According to Equation (9), at the time when the water level reaches its highest, the flow would be able to carry more sediment if the maximum velocity coincided. Therefore, the time lag between maximum velocity and highest water level reduces the flow capacity to transport suspended sediment.
According to the Kleitz-Seddon formula of flood propagation, flood-wave celerity is greater than flow velocity. Then, one may conclude that the peak of suspended sediment concentration should arrive at a gauging station after the flood-wave peak has passed because suspended sediment is transported at stream flow velocity. The small flood in 2008 was just the case. However, the 2006 flood was a different case. The reason is that the quasi-steady assumption used in deriving the Kleitz-Seddon formula did not hold in large floods. When the effect of water surface variation is considered, differentiating the Manning formula of flow velocity with time led to
(12)
where I is the water surface slope (–dH/dt). At the time of highest water level, dR/dt = 0; –dH/dt < 0 at upstream and –dH/dt > 0 at downstream, consequently,
(13)
This implies that the peak of flow velocity is ahead of highest water level when the variation of water surface slope cannot be neglected.
By comparing flood discharge data with the morphological change of Enshuda Coast, it was found that the river discharge below 6.0 × 103 m3/s supplied little sediment to the coast (Sato [29]). This confirms the finding that a 10-year flood in the Lower Tenryu River is not sufficiently effective in generating sediment flux from the river to the coast.
Since very large floods did not occur frequently in reality, the strategy to increase sediment delivery would have to be placed on enhancing the carrying capacity of median floods through channel manipulation.
4. Conclusions
In the Lower Tenryu River, suspended sediment was transported in different modes according to the magni-
Figure 10. Time lag between maximum velocity and highest water level.
tudes of flood. In a 30-year flood, the suspended sediment transport was supply-limited. However, during a 10-year flood and a small frequent flood, the mechanism was found to be transport-limited. In the supply-limited case, the variation of SSC followed the hydrograph well. In the transport-limited case, however, there was a time lag between the peak discharge and maximum concentration. It was further clarified that the cause of hysteresis in the 10-year flood was related to a point bar, while the time lag in the small flood can be attributed to the fact that flood celerity is greater than flow velocity. Besides, it was revealed that the symmetry of the 30-year flood hydrograph could be reproduced by a four parameter Pearson Type-Ill distribution with a value of less than 1.15 for skewness.
The findings of the present study indicate that the occurrence of hysteresis in the relationship between sediment concentration and flow discharge is not just a matter of sediment source but related to flood magnitude and affected by channel morphology. It also serves as a general warning that sediment-discharge rating curves must be used with caution; a relation between suspended sediment concentration and flow discharge obtained from one storm event may not be translatable to events of different magnitudes. Since most of the annual sediment transport usually takes place during a few events (Kronvang et al. [30]), thus it is very important to develop better understanding on the short-term dynamics of suspended sediment and refine the rating curve approach in order to have an accurate estimate of annual suspended sediment loading and discuss management strategy. In the case of the Lower Tenryu River, the suspended sediment hysteresis in a 10-year flood implies that a new strategy must be sought to increase sediment delivery for the purpose of beach protection along the Enshunada Coast.
5. Acknowledgements
This work is part of the Enshunada project supported by Japan Science and Technology Agency. The authors are grateful to Hamamatsu River and Road Administration Office for the supply of valuable field measurement data.