Modeling the impact of dam removal on channel evolution and sediment delivery in a multiple dam setting

Dam removal can generate geomorphic disturbances, including channel bed and bank erosion and associated abrupt/pulsed release and downstream transfer of reservoir sediment, but the type and rate of geomorphic response often are hard to predict. The situation gets even more complex in systems which have been impacted by multiple dams and a long and complex engineering history. In previous studies one-dimensional (1-D) models were used to predict aspects of post-removal channel change. However, these models do not consider two-dimensional (2-D) effects of dam removal such as bank erosion processes and lateral migration. In the current study the impacts of multiple dams and their removal on channel evolution and sediment delivery were modeled by using a 2-D landscape evolution model (CAESAR-Lisflood) focusing on the following aspects: patterns, rates, and processes of geomorphic change and associated sediment delivery on annual to decadal timescales. The current modeling study revealed that geomorphic response to dam removal (i.e., channel evolution and associated rates of sediment delivery) in multiple dam settings is variable and complex in space and time. Complexity in geomorphic system response is related to differences in dam size, the proximity of upstream dams, related buffering effects and associated rates of upstream sediment supply, and emerging feedback processes as well as to the presence of channel stabilization measures. Modeled types and rates of geomorphic adjustment, using the 2-D landscape evolution model CAESAR-Lisflood, are similar to those reported in previous studies. Moreover, the use of a 2-D method showed some advantages compared to 1-D models, generating spatially varying patterns of erosion and deposition before and after dam removal that provide morphologies that are more readily comparable to field data as well as features like the lateral re-working of past reservoir deposits which further enables the maintenance of sediment delivery downstream.     

How sediment composition and flow rate influence downstream channel morphodynamics upon dam removal

The process of dam removal establishes the channel morphology that is later adjusted by high-flow events. Generalities about process responses have been hypothesized, but broad applicability and details remain a research need. We completed laboratory experiments focused on understanding how processes occurring immediately after a sediment release upon dam removal or failure affect the downstream channel bed. Flume experiments tested three sediment mixtures at high and low flow rates. We measured changes in impounded sediment volume, downstream bed surface, and rates of deposition and erosion as the downstream bed adjusted. Results quantified the process responses and connected changes in downstream channel morphology to sediment composition, temporal variability in impounded sediment erosion, and spatial and temporal rates of bedload transport. Within gravel and sand sediments, the process response depended on sediment mobility. Dam removals at low flows created partial mobility with sands transporting as ripples over the gravel bed. In total, 37% of the reservoir eroded, and half the eroded sediment remained in the downstream reach. High flows generated full bed mobility, eroding sands and gravels into and through the downstream reach as 38% of the reservoir eroded. Although some sediment deposited, there was net erosion from the reach as a new, narrower channel eroded through the deposit. When silt was part of the sediment, the process response depended on how the flow rate influenced reservoir erosion rates. At low flows, reservoir erosion rates were initially low and the sediment partially exposed. The reduced sediment supply led to downstream bed erosion. Once reservoir erosion rates increased, sediment deposited downstream and a new channel eroded into the deposits. At high flows, eroded sediment temporarily deposited evenly over the downstream channel before eroding both the deposits and channel bed. At low flows, reservoir erosion was 17–18%, while at the high flow it was 31–41%.     

Simulating bed evolution following the Barlin Dam (Taiwan,China) failure with implications for sediment dynamics modeling of dam removal

The failure of the Barlin Dam in Taiwan, China offers an important case study for evaluating concepts in modeling the rapid erosion and channel recovery following intentional and unplanned dam removals. We present a modeling effort that applied a 1D and quasi-2D uncoupled hydraulics and sediment model (NETSTARS) to evaluate how discretization and parameterization influence thefitofbed elevationpredic-tions to observations following dam failure. Our analysis evaluated the model sensitivity to sediment transport function, active layer thickness, and number of stream tubes used to define the cross-section. Results indicate that a) the model is more sensitive to active layer thickness and sediment transport function than to the number of stream tubes, b) development of dam removal models are likely to benefit from varying the active layer thickness in time, and c) increased lateral discretization does not appear to improve model fit in the steep and rapidly changing river environment at our site. We conclude with discussion on differences between, identifying the need for, and general use of 1D, quasi-2D, and fully 2D models in dam removal and failure analysis.     

DAM IMPACTS ON AND RESTORATION OF AN ALLUVIAL RIVER-RIO GRANDE, NEW MEXICO

1 INTRODUCTION The construction of more than 75,000 dams and reservoirs on rivers in the United States (Graf, 1999) has resulted in alteration of the hydrology, geometry, and sediment flow in many of the river channels downstream of dams. Additionally, hydrologic and geomorphic impacts lead to changes in the physical habitat affecting both the flora and fauna of the riparian and aquatic environments. Legislation for protection of endangered species as well as heightened interest in ma…     

MITIGATION OF SEDIMENTATION HAZARDS DOWNSTREAM FROM RESERVOIRS

1 INTRODUCTION Increasing attention is being given to sedimentation hazards downstream from reservoirs as dams built during the past century accumulate progressively greater volumes of sediment. The sediment storage both decreases reservoir capacity and operating efficiency of the dam, and creates a 搒ediment-shadow?downstream where sediment-starved flows commonly erode channel boundaries and create long-term channel instabilities. Numerous studies have documented downstream channel change…     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.     

A decadal-scale numerical model for wandering,cobble-bedded rivers subject to disturbance

Alluvial rivers are composed of self-formed channels which are sensitive to disturbances in their flow and sediment-supply regimes. Regime changes commonly occur over decadal and longer timescales and can be caused by anthropogenic alterations such as dam construction and removal. Advances in numerical modeling have increased our ability to explore geomorphic adjustments over long timescales; however, many models designed to be run for decades or longer assume that banks are immovable or that channel width is constant. Since river channels often respond to disturbance by adjusting their geometry, this is a significant shortcoming. To investigate the impact of long-term sediment supply alterations on channel geometry and stability, we have adapted MAST-1D, a reach-scale bed evolution model, to incorporate functions for bank erosion, vegetation encroachment, and local avulsions. The model is designed for medium-large, coarse multithreaded rivers and can be run over long (decades–centuries) timescales. Bank erosion is a function of the mobility and transport capacity for structurally-important grains which protect the bank toe. Vegetation growth is proportional to point bar width and occurs during conditions of low shear stress. Local avulsions occur when aggradation causes channel depth to drop below a threshold. We apply the model to the Elwha River in Washington, USA with the goal of investigating if and when the river recovers from dam emplacement and removal. The Elwha was dammed for nearly 100 years, and then two dams were removed, releasing a large pulse of sediment. We have modeled the set of reaches between the two dams. Our simulations suggest that channel response to dam emplacement occurs gradually over several decades but that the channel recovers to near pre-dam conditions within about a decade following the removal. The dams leave a lasting legacy on the floodplain, which does not completely recover, even after two centuries. © 2019 John Wiley & Sons, Ltd.     

Variability of the useful life of reservoirs in tropical locations: A case study from the Burdekin Falls Dam,Australia

Large dam construction for irrigation, hydropower, water supply, and flood control in tropical to sub-tropical areas increased markedly after 1950. Many of the dams built during this period have filled with sediment and no longer perform to their original design specifications; in many cases, forecast dam life was greatly overestimated. This study investigates the useful life of Burdekin Falls Dam (BFD), northeastern Australia and compares the findings with other tropical reservoirs. Using two independent methods it is estimated that between 61 and 65 million m³ of sediment has been deposited in the reservoir over the 24 years of operation through 2011. This sediment volume equates to an average of 0.15% of capacity lost per year since construction was completed. If current sediment loads/climate regimes persist, reservoir capacity will be reduced by 50% after 345 years. However, the useful life of the BFD reduces to just 276 years when drawdown data are considered; these data show reservoir use would be affected once 40% of storage was filled with sediment, with a 60% drawdown return period of 1 in 15 years. When compared to similar large tropical to sub-tropical reservoirs, the BFD has a slightly longer reservoir useful life than dams in India and a much longer half-life than for both similar-sized and larger dams in China, Brazil, and Iran. Properties of the BFD that promote a longer useful life include a lower trap efficiency, relatively low annual sediment load delivered to the reservoir, limited sediment deposition behind the dam wall (and uniform distribution of deposited sediment), and the export of highly turbid annual floodwaters before settling and deposition of any remaining sediment within the reservoir.     

Morphological responses and sediment processes following a typhoon‐induced dam failure,Dahan River,Taiwan

The rates and styles of channel adjustments following an abrupt and voluminous sediment pulse are investigated in the context of site and valley characteristics and time‐varying sediment transport regimes. Approximately 10.5 x 10⁶ m³ of stored gravel and sand was exposed when Barlin Dam failed during Typhoon WeiPa in 2007. The dam was located on the Dahan River, Taiwan, a system characterized by steep river gradients, typhoon‐ and monsoon‐driven hydrology, high, episodic sediment supply, and highly variable hydraulic conditions. Topography, bulk sediment samples, aerial photos, and simulated hydraulic conditions are analyzed to investigate temporal and spatial patterns in morphology and likely sediment transport regimes. Results document the rapid response of the reservoir and downstream channel, which occurred primarily through incision and adjustment of channel gradient. Hydraulic simulations illustrate how the dominant sediment transport regime likely varies between study periods with sediment yield and caliber and with the frequency and duration of high flows. Collectively, results indicate that information on variability in sediment transport regime, valley configuration, and distance from the dam is needed to explain the rate and pattern of morphological changes across study periods. Copyright © 2013 John Wiley & Sons, Ltd.     

Geomorphic monitoring and response to two dam removals: rivers Urumea and Leitzaran (Basque Country,Spain)

Dam removal has been demonstrated to be one of the most frequent and effective fluvial restoration actions but at most dam removals, especially of small dams, there has been little geomorphological monitoring. The results of the geomorphological monitoring implemented in two dams in the rivers Urumea and Leitzaran (northern Spain) are presented. The one from the River Urumea, originally 3.5 m high and impounding 500 m of river course, was removed instantaneously whereas that in the River Leitzaran, 12.5 m high, and impounding 1500 m of river course, is in its second phase of a four‐stage removal process. Changes in channel morphology, sediment size and mobility and river bed morphologies were assessed. The monitoring included several different techniques: topographical measurements of the channel, terrestrial laser scanner measurements of river bed and bars, sediment grain size and transport; all of them repeated in four (May, August, November 2011 and May 2012) and five (July and September 2013, April and August 2014 and June 2015) fieldwork campaigns in the River Urumea and River Leitzaran, respectively. Geomorphic responses of both dam removals are presented, and compared. Morphological channel adjustments occurred mainly shortly after dam removals, but with differences among the one removed instantaneously, that was immediate, whereas that conducted by stages took longer. Degradational processes were observed upstream of both dams (up to 1.2 m and 4 m in the River Urumea and River Leitzaran, respectively), but also aggradational processes (pool filling), upstream of Inturia Dam (2.85 m at least). Less evident aggradational processes were observed downstream of the dams (up to 0.37 m and 0.50 m in the River Urumea and River Leitzaran, respectively). Flood events, especially a 100 year flood registered during the monitoring period of Mendaraz Dam removal, reactivated geomorphological processes as incision and bank erosion, whereas longitudinal profile recovery, grain‐size sorting and upstream erosion took longer. Copyright © 2016 John Wiley & Sons, Ltd.     

Evolution of river planforms downstream of dams: Effect of dam construction or earlier human‐induced changes?

When studying the evolution of landscape, it is difficult to discriminate the influence of anthropogenic from natural causes, or recognise changes caused by different sources of human action. This is especially challenging when the influence of certain sources is overprinted. For instance, although dam closure is the most common method of altering river courses, dam construction is often preceded by hydro‐technical works such as channel straightening, embankment construction or sediment mining. Both dam construction and the hydro‐technical works that precede dam closure can result in changes in the balance between sediment supply and transport capacity, and often, changes in river planform. The main objective of this study was to verify whether the works preceding dam closure are an important driver of river planform changes on the lower Drava River (Hungary). The case study is based on geological and geophysical surveys, as well as the analysis of historical maps covering an anabranching, 23 km long valley section. We show that channel straightening conducted prior to dam closure resulted in a transition from a meandering to sinuous planform with channel bars. Dam construction itself then caused enhanced incision, exposure of bar surfaces, vegetation encroachment and the formation of an anabranching planform. Based on this study, we developed models of alluvial island and channel planform evolution downstream of dams. Dam construction enhances channel incision, narrowing, and the reduction of flow caused by earlier hydro‐technical works. Many rivers downstream of dams experience episodes of anabranching or wandering, with a multi‐thread pattern replacing sinuous, braided and meandering courses. When incision continues, river patterns evolve from anabranching to sinuous via the attachment of alluvial islands to floodplains. However, the timing and sequence of these changes depend on hydrological and sediment supply regimes, geomorphic settings and anthropogenic actions accompanying dam construction. Copyright © 2018 John Wiley & Sons, Ltd.     

A semi-analytical method for predicting the outflow hydrograph due to dam-break in natural valleys

The paper presents a semi-analytical method for predicting the flow rate hydrograph due to a hypothetical sudden and total dam failure in a natural valley. The method generalizes the approach proposed by Hunt for the dam-break problem in a rectangular frictionless sloping channel to a valley with a cross-section area expressed by a power-law function of water depth, in order to take into account the most common shapes of natural valleys. The parameters of the deriving model can be set by exploiting data usually available concerning the dam section geometry and the reservoir storage-depth curve. The application of the technique to three different reservoirs is discussed. The results show that the flow rate hydrographs obtained at the dam site agree with the ones calculated by means of a finite volume numerical code based on two-dimensional shallow water equations. The method requires moderate computational and data collecting effort, so it can be regarded as a useful alternative to other procedures commonly adopted in the practice.     

The role of landscape setting in minimizing hydrogeomorphic impacts of flow regulation

The Tongariro Power Development Scheme (TPDS) is used to regulate flow in the headwaters of the largest catchment on the North Island of New Zealand (the Waikato). Two small dams, the Rangipo Dam and the Poutu Intake Dam, were constructed in 1973 and 1983. The flow regime of the river is managed to divert freshes into the power scheme, but allows flows larger than 100 m³ s^?1 to be released, to rework and transport sediment through the catchment. Analysis of aerial photos and maps spanning 1928 to 2007, alongside field measurements, show that there have been few hydrogeomorphic adjustments since dam construction. This includes limited changes to channel geometry, channel planform and bed material organization immediately downstream of the dams. In addition, offsite effects are minimal, both 500 m downstream of each dam, and in the more sensitive, less confined reaches in the lower catchment (11 km downstream of the Poutu Intake dam). The limited changes can be attributed to the locations of the dams within reaches characterised by bedrock gorges and confined within terraces. These locations act to flush sediments and impose margins that allow minimal adjustment of the channel. Bed material within this reach is characterised by the presence of a boulder lag. This is sourced from long-term incision into lahar deposits, and acts to limit the rate of incision, creating a steep and stable base upon which active fractions are transported. Just as importantly, significant storage in the low-relief volcanic plateau located in the upper catchment acts to disconnect and store the high sediment yields generated by active volcanic cones in the western sub-catchment upstream of the dams. This limits the rate of sediment supply to regulated reaches. Findings from this study show that analysis of reach-scale controls is essential in framing dam site locations in relation to the distribution of reaches and landscape units across the catchment. In this instance, tributary inputs downstream of the dams do not replenish the sediment and flow removed at the dam locations, as has been observed in other regulated systems. Rather, the river itself is resilient to change and flow variability is well managed allowing geomorphically effective floods to occur. Landscape setting is a key consideration in determining the hydrogeomorphic impact of flow regulation.     

Optimization of the Sefid-Roud Dam desiltation process using a sophisticated one-dimensional numerical model

Although water and soil conservation activities reduce reservoir sedimentation, it is inevitable that reservoirs fed by rivers transporting high amounts of sediment will experience sedimentation. The Ghezel-Ozan and Shah-Roud rivers, which flow to the Sefld-Roud reservoir dam, are both highly sediment-laden and transport significant amounts of sediment in both bed load and suspended load forms to the reservoir. Hence, it seems that the only practical way to remove the sediment from the reservoir is to flush it out using the Chasse method. In the present paper, field measurements of Chasse operation characteristics taken in previous years are presented, and a numerical model that simulates this process is introduced. After calibrating the model using field measured data, the calculated results （for reservoir pressure flushing and released sediment volume） of the numerical model were compared with other measured data for the same Chasse operation and the results agree well. Finally, using the numerical simulation results, the best approaches to ensure highly effective flushing while conserving reservoir water are presented （at least for the Sefid-Roud dam）. The operation of the bottom outlet gates, the shape of the output hydrograph, and the reservoir water level variation during flushing were optimized. In addition, the numerical model and related parameters, which need to be calibrated, are discussed.     

A non-reflecting boundary condition for the finite element modeling of infinite reservoir with layered sediment

The design of seismic resistant concrete gravity dam necessitates accurate determination of hydrodynamic pressure developed in the adjacent reservoir. The hydrodynamic pressure developed on structure is dependent on the physical characteristics of the boundaries surrounding the reservoir including reservoir bottom. The sedimentary material in the reservoir bottom absorbs energy at the bottom, which will affect the hydrodynamic pressure at the upstream face of the dam. The fundamental parameter characterizing the effect of absorption of hydrodynamic pressure waves at the reservoir bottom due to sediment is the reflection coefficient. The wave reflection coefficient is determined from parameters based on sediment layer thickness, its material properties and excitation frequencies. An analytical or a closed-form solution cannot account for the arbitrary geometry of the dam or reservoir bed profile. This problem can be efficiently tackled with finite element technique. The need for an accurate truncation boundary is felt to reduce the computational domain of the unbounded reservoir system. An efficient truncation boundary condition (TBC) which accounts for the reservoir bottom effect is proposed for the finite element analysis of infinite reservoir. The results show the efficiency of the proposed truncation boundary condition.     

NUMERICAL SIMULATION OF SEDIMENT TRANSPORT IN ALLUVIAL RIVER WITH FLOODPLAINS

1 INTRODUCTION The transport of sediment in rivers with active floodplains is a two-dimensional process because the main channel and the floodplain can have very different transport capacities. Therefore, two-dimensional (2D) models are often used to simulate the streamwise and transverse variations of sediment erosion and deposition. Many 2D numerical models have been presented to simulate sediment transport in floodplains (James, 1985; Pizzuto, 1987; Howard, 1992; Nicholas and Walli…     

Upstream morphological effects of a sequential check dam adjustment process

Sediment flushing and the morphological responses to the procedure of check dam removal are still unclear. Following laboratory experiments that revealed three stages (deepening, widening, and volume release) of check dam adjustment, a check dam built in 2007 at Landao Creek in central Taiwan was adjusted in 2015 by removing central bars and cutting 2.5 m from the middle two piers (stage 1 + 2), with the purpose of regulating sediment transfer and keeping the thalweg at the center of the channel, while also preventing hill slope toe erosion. In 2019, four central piers were removed (stage 3) to increase the volume of sediment released. Annual surveys were conducted after the initial adjustment in 2015 through to 2020 using unmanned aerial vehicles (UAVs). The check dam adjustments revealed that the channel had narrowed and stabilized as indicated by regenerating riparian vegetation. Additionally, distinct terraces had formed on the hill slope toes of the creek channel in proximity to the check dam. The meander upstream weakened following the dam adjustments. This study combining laboratory experiments with actual field observation contributed immensely to check dam decommissioning. Additionally, this study illustrated how an adjustable check dam may aid regulation of sediment transport and thereby sediment balance. It can be adjusted accordingly based on the prevailing channel condition.     

考虑河道整治工程影响的一维水沙数学模型及其应用

在冲积河流上,各类河道整治工程的实施一定程度上会影响河床的演变过程.本文通过改进现有的一维水沙数学模型,重点研究护岸及护滩(底)这类限制河床进一步冲刷的整治工程对水沙输移及河床冲淤变形的影响.首先对固定断面的各节点采用特定的代码进行标记,以此区分河漫滩、有或无整治工程的主槽区域.然后对悬沙输移及河床冲淤变形模块进行改进:当断面发生淤积时,其形态调整不受整治工程的影响,淤积量将在整个断面上进行分配;当发生冲刷时,断面形态调整仅发生在未实施工程的位置或者受工程限制但形成了一定厚度淤积层的区域;在实施了整治工程且无法提供沙源的区域,河床冲刷则不会发生.最后长江中游荆江段2016年的模拟结果表明:考虑整治工程情况时,改进模型计算的河道冲刷量偏小且与实测值更为接近;河槽断面形态与实测结果也更为吻合.     

INFLUENCE OF THE UPSTREAM RESERVOIRS ON THE ADJUSTMENT OF DOWNSTREAM ALLUVIAL CHANNEL

I. INTRODUCTIONThe Yellow River is a heavily sediment--laden river. The sediment load of the Yellow River ranks the first in the world while its annual runoff is only of medium size. Toharness the river, it is necessary to build reservoirs for regulating runoff to meet the demands of economic development. Since the founding of PRC in 1949, I S4 large and medium--sized reservoirs have been constructed on the main stem and the tributaries with atotal storage capacity of 84.5 billion m3.…     

Real-time probabilistic sediment concentration forecasting using integrated dynamic network and error distribution heterogeneity

Sediment forecasting at a dam site is important for the operation and management of water and sediment in a reservoir. However, the forecast results generally have some uncertainties, which may hinder the operation of the dam. In this study, a real-time sediment concentration probabilistic forecasting model is proposed based on a dynamic network model. Under this framework, the Elman neural network (ENN) and nonlinear auto-regressive with exogenous inputs (NARX) neural network models were established for sediment concentration forecasting with different lead times. A hybrid algorithm, which combined the Levenberg–Marquardt algorithm and real-time recurrent learning, was used to train the model. Using the aforementioned method, the sediment concentration was forecast for at the Sanmenxia Dam, China, and, subsequently, the forecast results were evaluated. Among the selected lead time, the results at 5 h exhibited the highest accuracy and practical significance. Compared with the ENN model, the sediment concentration peak error using the NARX neural network was reduced by 4.5%, and the sediment yield error was reduced by 0.043%. Therefore, the NARX neural network was selected as the deterministic sediment forecasting model. Additionally, the probability density function of the sediment concentration was derived based on the heterogeneity of the error distribution, and the sediment concentration interval, with different confidence levels, expected values, and median values, was forecast. The Nash–Sutcliffe coefficient of efficiency for the sediment concentration, as forecasted based on the median value, was the highest (0.04 higher than that using a deterministic model), whereas the error of the sediment concentration peak and sediment yield remained unaltered. These results indicated the accuracy and superiority of the proposed real-time sediment probabilistic forecasting hybrid model.