Top elevations of weir structures.
Abstract
As low-head hydraulic structures, instream weirs are built across rivers to control the upstream water surface elevation and the downstream flow conditions. This chapter presents a study of erosion control using instream weirs at downstream of a reservoir; JiJi Weir was built across the longest river in Taiwan, Chuoshui Creek, a mountainous river with steep slopes. Due to the easy-to-be-eroded fine lithology layers of mud, shiver, and sandstones on channel bed, the downstream of JiJi Weir had suffered from severe channel incision and head-cut development problems, which greatly threatens the integrity of the dam. To protect the JiJi Weir and its downstream channel from serious channel erosions, the Water Resources Agency (WRA) of Taiwan proposed erosion control plans that multiple instream weir structures were to be installed along the downstream channel of JiJi Weir. A three-dimensional (3D) numerical model, CCHE3D model with capabilities of simulating bedrock erosions, was used to evaluate those erosion control plans and thus explore for the optimal design.
Keywords
- instream-weir
- bedrock erosion
- channel incision
- erosion control
- three-dimensional model
1. Introduction
JiJi Weir was built across Chuoshui River, the longest river (178.6 km) in Taiwan, originating from the central mountains and flowing into the East China Sea, covering a watershed about 4323 km2. A shallow reservoir was formed due to the construction of the weir. The precipitation and geology dominate the flow and sedimentation processes in the channel and the reservoir. Particularly, during typhoon seasons, floods with high discharges will carry a large amount of sediments flowing through the river channel. Shortly after the weir construction, a pebble-boulder-sized sedimentation layer covering the channel bed was mobilized by the initial clear water releases, and the soft rocks beneath this sedimentation layer were exposed. Within the 6.5 km region downstream of the dam, the channel bed now is featured with erodible lithology layers of mud, shiver, and sandstones, resulted from severe channel incisions and a head-cut development. According to the comparisons of historical longitudinal profiles in Figure 1, the vertical incision of the channel thalweg from 2002 to 2014 is about 10 m (please refer to Figure 2 for the locations of those measured cross sections from CS-106 to CS-116). Looking downstream, Figure 3 shows the incised channel near the JiJi Weir. The white dash line indicates the flat channel bed before the incision. The head-cut was still actively developing and migrating, which has threatened the integrity of the dam (Figure 4).
To prevent this channel from further erosion and protect the JiJi Weir Dam, the Water Resources Agency (WRA) of Taiwan has proposed several erosion control plans [1, 2]. Several weir structures and channel widening were proposed to promote sediment deposition and lower the water surface level for flood protection. This study is to provide alternative plans optimal for erosional control and flood protection by using a 3D numerical model, CCHE3D [3], with capabilities in simulating bedrock erosions with lateral erosions.
The bed rock erosion is closely related to channel lithology, tectonics, and geomorphology properties, contributed significantly by complicated weathering, plucking, abrasion (due to both bed load and suspended load), cavitation, and dissolution processes [4], which is distinguished from river sediment transport and morphologic processes on alluvial beds. Extensive research studies have been carried out for attempting to describe these complex physical mechanisms of bedrock erosions mathematically. For examples, Annandale [5] proposed a conceptual framework to correlate the flow energy to the earth mass erodibility by introducing the erodibility index, which characterizes the capability of earth materials for resisting erosions. Whipple and Tucker [6] considered the stream power to dominate the bedrock erosions. Whipple
From the previous studies, the fluvial impacting factors for the bedrock erosions were identified as follows: bedrock erodibility (strength), stream power, shear stress, sediment supply, and grain size. Accordingly, two types of models [12], namely the hydraulic scour model (stream power-based [5, 13]) and the abrasion scour model [8] have been used in the applications with numerical simulations. Jia
This chapter is to present the development and applications of the CCHE3D bedrock erosion module to study the channel incisions and erosion control problems in the downstream channel of the JiJi Weir. The numerical model was calibrated and validated using field data. Based on the evaluations of the multiple alternatives, an optimal control plan, involving multiple weir constructions with the side channel excavations as well, was proposed to control the head-cut development soft rock and the channel incision.
2. Numerical model
CCHE3D model [3] is a Finite Element Method (FEM) model based on a partially staggered 3D sigma layered mesh system, which consists of multiple layered structured meshes. It was developed for 3D numerical simulation and analysis for free surface turbulent flows in rivers, lakes, reservoirs, and estuaries and associated processes, such as sediment transport [18], morphological changes [17], heat transfer [19], pollutant transport and water quality evaluation [20], etc. Since it can provide more accurate and detailed local flow fields around in-stream structures [18] than 2D depth averaged models, CCHE3D model was selected to evaluate the erosion control plans with multiple weir structures installed in the study channel.
2.1 3D RANS model
The full 3D Reynolds-averaged Navier-Stoker (RANS) equations are solved in CCHE3D model, which are simply listed in the form of index notations for Cartesian coordinates as follows:
where
For surface flows, the free surface kinematic equation is applied:
where
In CCHE3D model, a partially staggered stencil and a velocity correction algorithm are used for solving the momentum and continuity equation. Several turbulence closure schemes including the zero equation models (parabolic, mixing length, and wind-induced) and the
2.2 3D sediment transport model
Sediment transport is one of the most complex and least understood phenomena in nature. In 3D, sediment particles’ movements are highly affected by vertical motion of fluid flows in addition to horizontal movements. This is particularly true in the vicinity of hydraulic structures where the flow impacts on the solid walls of the structures (bridge pier and abutment, for instance). If the structures are very large (dam), the fluid flows in an open channel are forced to change their speed and direction near structures in order to pass through them, the sediment transport capacity due to this impact is adjusted significantly causing localized scouring and deposition over the sediment bed. In CCHE3D model, in addition to general sediment transport capabilities, special sediment transport features, such as local scouring around structures and channel head-cut migration have been developed as well.
In CCHE3D model, the 3D convection-diffusion equation for the suspended sediment is solved as follows:
where
At the free surface, the vertical sediment flux is zero, so the gravity effects
At the bottom, the following condition is applied:
where
Following the non-equilibrium transport approach (
where
2.3 Bedrock erosion model
According to previous studies, the bedrock erodibility (strength), stream power, shear stress, sediment supply, and grain size were identified as important impacting factors on bedrock erosion rate, which lead to two popular bedrock erosion mechanisms, plucking and abrasion, widely used in the numerical models [1, 12, 14, 16]. The first one is corresponding to the so-called stream power-based method that the bedrock erosion rate is considered as a function of the stream power, and the varying shear stress causes the hydraulic scouring on soft bedrock [4, 5, 6, 13], and the other one is the abrasion-based method, which emphasizes on the important role of the sediment supply (both bedload and suspended load) by considering the eroding and shielding effects of sediment on bedrock [7, 8, 11].
Whipple
In the stream power method, the bedrock erosion rate
where
The above stream power method is further improved to take into account lateral erosion by considering the local lateral bed slope. Thus, a factor
where
where
2.4 Coupling of bedrock erosion model and sediment transport model
When the sediment transport is simulated, the boundary condition between the moving sediment particles and the bedrock has to be treated. Sediment particles can deposit over the rock surface and form a deposition layer. A concept of a sediment mixing layer over the softrock surface is adopted. If the thickness of the sediment layer is large, no rock erosion is calculated. If no sediment deposition exists on the bed, the proposed stream power method is used for the rock bed erosion. If the thickness of the sediment layer is within a criterion (mixing layer thickness), the stream-power-induced erosion would be applied at a reduced rate, proportional to the thickness of the mixing layer. The net change rate of the mixing layer is the combined rates of rock erosion and sediment deposition.
There is insufficient knowledge on the interactive and coupling mechanism between the bedrock erosion and the sediment transport. This simple coupling basically considers the shielding effects of the sediments on bedrock, similar to that of the abrasion model, which is concetually reasonable and has been proved and validated in previous studies [7, 8, 11]. However, the actual interactions between the bedrock erosion process and the sediment transport process are much more complicated in nature.
3. Application to downstream of JiJi Weir
In this study, CCHE3D bedrock erosion model was used to simulate the aforementioned channel incision process. The validated model was then used to evaluate erosion control plans. The study domain is a 6.5 km reach from the JiJi Weir to the downstream Minchu Bridge, where there are measured cross sections available: CS-106 to CS-116, and JiJi-30 to JiJi-22 (Figure 2).
3.1 Model setup
Based on the DEM data of 2013 and 2014, several structured meshes covering the study domain from JiJi Weir to Minchu Bridge were generated, which embed the combinations of the planned weir structures and lateral channel excavations.
As shown in Figure 5, zones representing erodible bedrocks, alluvium, structures (non-erodible), have been identified and embedded in simulation domain based on the satellite image of 2014. The erodible bedrock zone resulted from the channel incisions and head-cut development, where both the bedrock erosion and sediment transport simulations are to be conducted. The alluvial zone is covered with pebbles and gravels, and no bedrock erosion simulation will be applied. The non-erodible zone represents the weir structures and other bank/bed protection measures (i.e., concrete blocks) and is assumed non-erodible during the simulations. Therefore, the distribution of non-erodible zones changes with the simulated erosion control plans.
Figure 6 shows the measured bedrock erodibility index distribution in 2014. The whole domain was divided into 11 zones with different erodibility index,
3.2 Model calibration and validation
The site-specific parameters (
In the calibration step, the parameters are adjusted to obtain the best fit of the simulations to the data; while in the validation step, no parameters are changed. Figure 7 shows the hydrographs of these two typhoons, which are considered as the major hydrological events causing significant bed morphological changes in this channel. Due to its high peak discharge (
Through the numerical calibration tests, the parameter set for the stream power method was found to be:
3.3 Erosion control plans
Since 2007, WRA has proposed a few erosion control plans attempting to stop the downstream channel of JiJi Weir from incisions and head-cut development. These erosion control plans included multiple in-stream weir structures, lateral channel excavations, and other engineering measures, such as concrete blocks for bank protections, gabions, etc. With the calibrated and validated parameters, CCHE3D bedrock erosion model was used to evaluate all the erosion control plans with different combinations of control structures [2]. According to the numerical simulations, one plan with three weirs at JiJi-22, 25, and 26 and side excavations from JiJi-27 to CS-112 was confirmed as the most effective in the erosion reduction among all proposed erosion control plans. Based on the numerical simulations of these erosion control plans, optimal design was explored.
In the historical channel profiles along thalweg of the downstream channel from 1998 to 2014 as shown in Figure 1, the bed slope is 0.0412 near the head-cut reach, 0.0057 in the incision reach, and 0.0028 in the transition reach. The head-cut reach and the incision reach are of bare rock channel, actively eroded; the transition reach is sometimes partially covered with sediments, showing alluvial river morphologic and sediment transport features, and thus considered to be more stable. In this study, it is assumed that the channel will be stabilized if the bed slope can be reduced to 0.0028 approximately by installing erosion control weirs in the channel.
Since the 10 m deep head-cut has reached closer to the JiJi Dam, seriously threatening the safety of the dam structure (Figure 3), the protection of the upstream head-cut zone (from JiJi-30 to JiJi-22) is considered as the first priority. The high weir structure at JiJi-25 is capable of significantly reducing the bed rock erosion in the reach from JiJi-26 to JiJi-22. As for the reach from JiJi-30 to JiJi-26, a weir structure is planned at CS-115, and the small reservoir behind this weir structure is designed to slow down the flow and thus reduce the erosion.
For downstream of the head-cut zone, the deep incised main channel (from CS-109 to JiJi-26) is relatively narrow, the water surface elevation in this reach will increase during floods to endanger the embankment of the left bank. Channel lateral excavations are proposed in such a way that the thalweg is kept and the widening lateral excavation (150
To enhance the erosion reduction effects, two additional low-headed weirs are suggested to be installed at CS-111 and CS-113 in this reach to further control the flow and erosion. To prevent the development of a second head-cut between CS-108 and CS-109 (see Figure 1), another low-headed weir structure is proposed as well to be installed at CS-108.5.
The heights of the four weirs at CS-108.5, CS-111, CS-113, and CS-115 are determined in such a way that the small pools formed behind the weirs can approximately protect half of the reach between weirs (slope equal to zero). As illustrated in Figure 10, the bed slope of the reach between Weir-A and B is
Weirs | CS-108.5 | CS-111 | CS-113 | CS-115 | JiJi-26 | JiJi-25 |
---|---|---|---|---|---|---|
Elevation (m) | 155.5 | 163 | 170 | 177 | 181 | 188 |
Figure 11 shows the initial setup for the optimal design. The excavation is denoted in the white polygon area, while six weir strucrures installed at JiJi-25 (CS-116), JiJi-26, CS-115, CS-113, CS-111, and CS-108.5) as non-erodible zones. For the optimal design case simulations, the calibrated parameter sets and the bedrock erodibity index remain the same. The bed property zone is adjusted if the excavation cuts into the alluvium and exposes the soft rock underneath.
In current sediment transport simulations, both the suspended load and the bedload transport are considered with three representative size diameters: 0.0081, 0.03086, and 0.3006
The simulation results of the optimal design were compared with those of the cases without any structures. Figures 12 and 13 show the longitudinal profile of bed changes. In the first 1.2 km reach, where two high weirs at JiJi-25 and JiJi-26 are located, much more depositions were observed in the optimal plan, while the erosions are dominant for the case without structures. Further downstream, in general, the erosions are dominant for both cases, except for small regions upstream of the weir structures, where the water surface is increased, and the velocity is slower . The bed change pattern around all the weir structures is similar: there are depositions at upstream but erosions at downstream, which is expected. Figure 14 shows the bed profile along the thalweg. As can be seen, in the optimal plan, the erosion along the thalweg has been reduced significantly, and the upstream reach with the first two weirs demonstrated obvious deposistion pattern. The two weir structures at JiJi-25 and JiJi-26 do serve the purpose of reducing channel incision and stopping the head-cut development effectively (Figure 15).
Despite the simple coupling method, the coupled bedrock erosion and sediment transport model demonstrated the promoting effects of the proposed optimal design on sediment depositions to protect the channel bed from further eroding.
3.4 Implementation of erosion control plan
The erosion control plan combined with other engineering measures (concrete blocks, gabions, and filling-ups) has been partially or completely implemented. Figure 16 compares the aerial images of JiJi Weir in 2014 and 2021. In this 1.2 km reach of JiJi Weir downstream, a comb weir structure, two drop structures, and some concrete blocks were installed. The head-cut development from 2014 to 2021 seems not very significant, which demonstrated the effectiveness of the erosion control plan.
4. Conclusions
In this study, a simply coupled 3D bedrock erosion and sediment transport model, CCHE3D bedrock model, was developed and applied to the downstream channel of the JiJi Weir, Chushui River, where the serious channel incision and head-cut migration endanger the JiJi Dam. In this simple coupling model, a concept of a sediment mixing layer over the bedrock surface is adopted. The thickness of the sediment layer determines the bed change mechanism as follows: thick sediment layer leads to non-rock erosion calculation, while non-sediment layer indicates the dominant role of the bedrock erosion. Within the mixing layer thickness, the stream-power-induced erosion would be applied at a reduced rate, proportional to the thickness of the mixing layer. The net change rate of the mixing layer is the combined rates of rock erosion and sediment deposition.
The hydrologic data of Typhoon Matmo in 2014 and Typhoon Morakot in 2009 were used to calibrate and validate the site-specific parameters for the stream power method–based bedrock erosion model. This model was also applied to simulate the bed changes with designed channel erosion control structures installed. The objective of the study is to identify an optimal design to alleviate channel incision and stop the head-cut development.
Based on the evaluations of multiple erosion control plans proposed by WRA and the channel evolution analysis, an optimal erosion control plan was identified. Six weir structures at section JiJi-25, JiJi-26, CS-115, CS-113, CS-111, and CS-108.5 and the 150 m-wide lateral excavations from JiJi-27 to CS-108.5 along the right floodplain to widen the channel were proposed. The lateral excavation would keep flows in the thalweg, but divert the water onto the excavated area during floods. According to the longitudinal profiles along the thalweg (Figure 1), the CS-115, CS-113, CS-111, and CS-108.5 were selected as the locations for installing weir structures to control the channel bed slopes. The elevations of those weir structures were determined in such a way that the target weir protects a half of a segment between two neighboring weirs so that the bed slope of the pool between weirs could be reduced by half (Figure 10).
Using the proposed coupled bedrock erosion and sediment transport model, the optimal design and the case without any structures were evaluated and compared. According to the simulated results, the optimal design reduced the channel incision significantly, and the head-cut development was stopped by the two weir structures installed at JiJi-25 and JiJi-26. Sediment depositions were observed not only in the deep channel but also the excavated area. In general, the whole study reach from JiJi Weir to Minchu Bridge still demonstrated erosion pattern except for the first 1.2 km reach at upstream. Without any control structures, the channel would be further deepened, and the development of the head-cut would continue.
Acknowledgments
This work is a result of research sponsored by National Chiao Tung University of Taiwan under Research Agreement No. 300212267A and The University of Mississippi.
It is also a part of research sponsored by the USDA Agriculture Research Service under Specific Research Agreement No. 6060-13000-025-00D (monitored by the USDA-ARS National Sedimentation Laboratory) and The University of Mississippi.
Appendices and nomenclature
C | suspended sediment concentration |
Db | deposition rate (m2s−1) |
Eb | erosion rate (m2s−1) |
Ebb | erosion rate for both bedrock and bank erosion (m2s−1) |
Kh | bedrock erodibility index |
Lb | bedload adaptation length (m) |
qb | bedload transport rate (m2s−1) |
Sb | bed slope factor |
u | x velocity (ms−1) |
v | y velocity (ms−1) |
w | z velocity (ms−1) |
U | depth-averaged velocity of flow (ms−1) |
P | stream power of flow (kW m−2) |
η | water surface (m) |
ν | kinematic viscosity (m2s−1) |
δ | bedload layer thickness (m) |
τ | shear stress (Pa) |
ωs | sediment settling velocity (ms−1) |
σs | Schmidt number |
References
- 1.
Jia Y, Zhang Y. Preliminary Study of Local Scouring Using CCHE3D and Its Application Potential to Softrock Erosion in Choshuixi, Taiwan. National Center for Computational Hydroscience and Engineering. USA: University of Mississippi; 2013. p. 38677 - 2.
Jia Y, Zhang Y. Enhancement and Application of CCHE2D/3D Soft-rock Models to Chuoshui Creek of Taiwan. National Center for Computational Hydroscience and Engineering. USA: University of Mississippi; 2015 - 3.
Jia Y. CCHE3D Technical Manual, National Center for Computational Hydroscience and Engineering. USA: The University of Mississippi; 2013 - 4.
Whipple KX, Hancock GS, Anderson RS. Riverincision into bedrock: Mechanics and relative efficacy of plucking, abrasion, and cavitation. Geological Society of American Bulletin. 2005; 112 (3):490-503 - 5.
Annandale GW. Erodibility. Journal of Hydraulic Research. 1995; 33 (4):471-494 - 6.
Whipple KX, Tucker GE. Dynamics of the stream power river incision model: implications for height limits of mountain ranges, landscape response time scales, and research needs. Journal of Geophysical Research. 2005; 104 (B8):17661-17674 - 7.
Sklar LS, Dietrich WE. Sediment and rock strength controls on river incision into bed rock, Geological Society of America. Geology. 2009; 29 (12):1087-1090 - 8.
Sklar LS, Dietrich WE. A mechanistic model for river incision into bedrock by saltating bed load. Water Resources Research. 2004; W06301 (40):1-21 - 9.
Turowski JM, Hovius N, Meng-Long H, Lague D, Men-Chiang C. Distribution of erosion across bedrock channels. Earth Surface Processes Landforms. 2005; 33 (3):353-363 - 10.
Stock JD, Montgomery DR, Dietrich WE, Sklar L. Field measurements of incision rates following bedrock exposure: Implications for process controls on the long profiles of valleys cut by rivers and debris flows. GSA Bulletin. 2005; 117 (11/12):174-194 - 11.
Lamb MP, Dietrich WE, Sklar LS. A model for fluvial bedrock incision by impacting suspended and bed load sediment. Journal of Geophysics Research. 2008; 113 :F03025 - 12.
Lai Y, Greimann B, Wu K. Soft Bedrock erosion modeling with a two-dimensional depth-averaged model. Journal of Hydraulic Engineering. 2005; 137 (8):804-814 - 13.
Annandale GW. Scour Technology, Mechanics, and Engineering Practice. New York: McGraw Hill; 2006 - 14.
Jia Y, Zhang Y, Wang S-SY. Computational Study of Softrock Erosion in a Mountain River. In: 13th Cross-Straits Symposium of Hydraulic and Water Resources. Taichung, Taiwan; 2009 - 15.
Jia Y, Wang SSY. Numerical model for channel flow and morphological change studies. Journal of Hydraulic Engineering. 1999; 125 (9):924-933 - 16.
Liao C-T, Keh-Chia K, Huang M-W. Development and application of 2-D mobile-bed model with bedrock river evolution mechanism. Journal of Hydro-environment Research. 2013; 8 (3):210-222 - 17.
Jia Y, Zhang Y. Enhancement of CCHE2D/3D Models on Softrock Erosion, Local Mesh Refinement and Their Applications to Taiwan’s Rivers. National Center for Computational Hydroscience and Engineering. USA: University of Mississippi; 2014 - 18.
Jia Y, Scott S, Xu Y, Huang S, Wang SSY. Three-Dimensional numerical simulation and analysis of flows around a submerged weir in a channel bendway. Journal of Hydraulic Engineering. 2005; 131 (8):682-693 - 19.
Chao X, Jia Y, Wang SSY. 3D numerical simulation of turbulent buoyant flow and heat transport in a curved open channel. Journal of Hydraulic Engineering. 2009; 135 (7):554-563 - 20.
Chao X, Jia Y, Shields FD Jr, Wang SSY, Charles M. Three-dimensional numerical modeling of water quality and sediment-associated processes with application to a Mississippi Delta Lake. Journal of Environmental Management. 2009; 91 (7):1456-1466 - 21.
Wu W. CCHE2D Sediment Transport Model, Technical Report No. NCCHE-TR-2001-3, NCCHE, The University of Mississippi