Engineered rivers are often prone to channel bed incision. This decreases the channel-floodplain connection, hampers navigation where nonerodible reaches increasingly protrude from the bed, and can destabilize structures. Here we inventorize causes and characteristics of channel incision measures. We elaborate on how channel bed incision is a transient channel response toward a new equilibrium channel state. Causes of incision comprise base level fall, channel narrowing (e.g., due to river training), channel shortening (bend cut-offs), an increased channel-forming discharge (e.g. due to climate change), and a decrease (or fining or coarsening) of the sediment flux from the upstream part of the basin. Finally, we discuss two measures that may mitigate channel bed incision: sediment nourishments and longitudinal training walls.

@en

Channel adjustment in engineered rivers is often associated with channel bed incision (e.g., Chowdhury et al., 2023, Czapiga et al., 2022a, 2022b, Ylla Arbós et al., 2021). Channel bed incision reduces the stability of in-river structures, exposes river-crossing cables and pipelines, and the spatial variability of channel bed incision due to less erodible reaches creates shipping bottlenecks. Various measures have been implemented to cope with these issues. They range from sediment nourishments to erosion control structures (e.g., Habersack and Piégay, 2007). Our objective is to assess the potential of floodplain lowering and sediment nourishments in mitigating large-scale channel bed incision in engineered rivers affected by climate change, considering a spatial scale of hundreds of kilometres. Our domain of interest is the Rhine River between Bonn, Germany, to Gorinchem, Netherlands. This reach has been extensively channelized during the 18th-20th centuries for improved navigation and flood protection (e.g., Ylla Arbós et al., 2021).@en

Climate change is responsible for global shifts in precipitation patterns and an overall in-crease in global temperatures. The transi-tions are anticipated to modify the river hydro-graph and sea level. The changes to the hy-drograph are also likely to influence sediment flux. These alterations imply shifts in both up-stream and downstream boundaries for river bifurcations. However, the resulting bifurca-tion response remains uncertain and warrants further investigation. Our objective is to un-derstand the extent of large-scale and long-term response of river bifurcations to climate change. We take the Upper Dutch Rhine bifur-cation region as our case study and develop a 1D hydro-morphodynamic model representing the system to achieve this goal.@en

River bifurcations divide the water and sediment over two downstream branches or bifurcates. As the changing climate adjusts the boundary conditions (i.e., base level, hydrograph, and sediment flux) for bifurcations, it will affect their flow and sediment partitioning over the bifurcates. Our objective is to provide insight into the response of a bifurcation to sea level rise (SLR). To this end, we compare the response of an idealized bifurcation in an engineered river (i.e., with a fixed planform and width) to SLR to the one of a single channel, using a one-dimensional numerical model system.@en

Climate change puts pressure on river systems, as it increasingly alters the river controls. Engineered rivers with a fixed planform respond to climate change and human intervention by adjusting the channel slope and bed surface grain size distribution. This response often consists of channel bed incision, over hundreds of kilometres, and during decades to centuries, resulting in serious disruptions of inland navigation, increased flood risk, and ecological degradation. Here we investigate how the lower Rhine River (Bonn, Germany – Vuren, the Netherlands, including the Pannerden bifurcation) continues to adjust to channelization measures of the 19th century (Ylla Arbós et al., 2021), and responds to different climate scenarios of control change over the 21st century, using a schematized one-dimensional numerical model for mixed size sediment.@en

Climate change is causing the global sea level to rise. Research and discussion of the effects of sea level rise are often focused along coastlines. However, the effects of higher water level and changing morphodynamics can reach far inland via rivers. This study uses a one-dimensional numerical model to analyze bed level response to sea level rise. The model simulates 100 years of steady sea level rise on a 1000 km long, fixed width channel. Sea level rise creates a backwater curve which grows in the upstream and vertical directions. The transient response of the channel bed is an aggradation wave the grows in the upstream and vertical directions. We studied different cases which vary in sediment flux, flow discharge, grain size, and rate of sea level rise and found changes in the rate of growth of the aggradation in both directions. All runs start in an equilibrium state and run for 100 years.

This study finds a close relationship between the equilibrium slope and depth of the channel and the shape of the backwater curve. This relationship drives the aggradation patterns. For example, for the same amount of sea level rise, a flat, deep channel has a longer backwater curve with a smaller relative increase in depth than that of a steeper, shallower channel. The backwater curve drives the aggradation patters, such that the flat, deep channel then has aggradation over a longer reach, and a smaller increase in bed level than the steeper shallower channel.

With the cases modeled in this study, three general trends in aggradation rates emerge: (1) an aggradation mound that grows quickly upstream, with slower increase in bed level as found in flatter, deeper channels; (2) a faster increase in bed level with slower upstream growth, found in cases with steeper slope and shallower depths; and (3) faster growth in both bed level and upstream direction caused by an increased rate of sea level rise.

Since natural channels are often complex with multiple sources of discharge inputs, a tributary case is also
included. Starting from equilibrium state and applying sea level rise to the downstream boundary, the model
shows aggradation waves in the regions downstream and upstream of the confluence, starting from the downstream boundary and the confluence. There is also degradation just downstream of the confluence. The scour hole grows at first, in depth and the downstream direction then reduces. The aggradation moving upstream intersects with the degradation moving downstream and fills in the scour hole. In some cases, we see the scour fills in within the 100 year time frame, resulting in a net increase in bed level. In a tributary system, the risk of scour is greatest for tributaries with high flow discharge or low sediment flux.

The transient response of the bed level to sea level rise is aggradation. As the water level continues to rise, the bed level is expected to do the same, but at a lower rate. In this study, the nominal rate of sea level rise is 10 mm/yr. The fastest rate of bed level rise in the model results is less than 6mm/yr, creating an ever increasing water depth. This is beneficial for shipping and navigation in the channel, which would not require dredging of the aggradated material. However, the reduction in bankfull volume with rising water levels is dangerous for flood control.

In the last few years, the climate crisis has been accelerating at a dizzying pace and poses an emergency threat to our planet. Rainfalls have transformed into intense downpours, and flash flooding, combined with the sea level rise, leads to a higher risk of inundation of densely populated coastal cities. Moreover, the melting of the permafrost can lead to significant landslides, triggering the generation of mega-tsunami waves. The latter can have a catastrophic impact, not only on the infrastructure but also on human life. Man-induced climate change is responsible for the increase in frequency and intensity of extreme natural events, such as tsunamis, floods, and storm surges. Recent research indicates that almost one-fourth of the world population lives at high-risk locations to at least 0.15m of inundation depths with a return period of 1 in 100 years. Therefore,
the urge of implementing protection measures against unsteady flows is imperative.

Undoubtedly, the involvement of engineers can play a pivotal role in order to analyze these flows in the built environment and provide sufficient coastal and building plans to ensure safety and reduce reconstruction
costs. The behavior of the unsteady flow around a structure is not a well-understood topic and results in a lack of accuracy and reliability. More insights are required into the fluid-structure interactions to come up
with a safe building design.

In the present study, unsteady flows are generated using the dam-break technique in line with previous research. The Thesis aims to model, validate, and implement a simplified approach to analyze the complexity of the hydrodynamic behavior of unsteady flows around impervious buildings, with different orientations
and blockage ratios. To do so, the research introduces a numerical simulation method of a dam-break wave, using the two-dimensional, and non-rotational shallow water equations. The Galerkin finite-element model is applied for the discretization of the solution on a limited domain.

Initially, the flow of the dam-break wave was validated in the absence of the structure, using a dam-break experimental work for comparison. Then, in order to insert a structure in the domain, a second experiment with a structure is used, which generates tsunami-like waves using the vertical release technique. The gate, which represents the dam, is located at x=0 and opens instantaneously at t=0s. Behind the gate, the reservoir maintains an amount of water that flows, after the opening of the gate, into the channel generating the shock
wave. At the channel downstream of the gate, initial water levels are considered and the behavior of a bore propagation is simulated. The building is located on the downstream side and different impervious building configurations were studied in terms of orientation and shape to investigate the impact of the unsteady flow.
To analyze the complex hydrodynamic processes of flooding, four fixed points around the structure are set to measure the action of the bore on different impervious building configurations. The water elevations and the averaged velocity profiles in time were derived at each point. Moreover, the horizontal forces in the x and y directions are calculated by the model, integrating the stresses over the wet surface of the buildings.

The general behavior of the fluid-structure interactions is captured well, especially upstream of the building, and insights are gained regarding the behavior of the fluid around the structure and the parameters of
influence. Results showed that orientation changes completely the impact on the building configurations. The separation of the flow and the blockage ratio are the main parameters that are influenced by the angle of rotation and change the behavior of the loading process at the initial impulsive phase, when the wave arrives at the structure, and at the hydrodynamic phase, where the flow has a quasi-steady behavior. Overall, a good agreement is achieved with the experimental data, although the numerical model overestimates the loads
acting on the different building configurations. Results proved that the orientation of the building with respect
to the flow facilitates the flow around the structure and contributes to lower water levels, and to a better distribution of the horizontal loads on the surface. The best results were achieved for an angle of rotation of 45°, where symmetrical separation of the flow is also playing an important role.

Sediment transport capacity and supply of sediment to a river channel increase significantly during peak flow events. Here we study how a river bifurcation system (partitioning water and sediment over its downstream branches) responds to peak flow events. We focus on the Pannerdense Kop bifurcation in the Dutch Rhine River, an engineered system where planform and channel width are fixed. We analyze water discharge and bed level data measured over the last century. We observe rapid aggradation in one of the branches (Pannerden Channel) following the peak flow events of 1993 and 1995, and little to no bed level change in the other branch (Waal). Prior to the event, both branches eroded, and the upstream part of the Pannerden Channel had a greater erosion rate than the Waal. After the 1993 and 1995 peak flow events, the erosion in the upstream part of the Pannerden Channel slowed significantly, whereas the upstream part of the Waal branch continued to erode (though at a smaller pace than before the peak flow events). This differential erosion has resulted in a gradual increase of water discharge toward the Waal branch. Interestingly, the bifurcation system does not appear to respond equally to all peak flow events. We hypothesize that the bifurcation response to the 1993 and 1995 peak flows differs from previous peak flows because of the sequence of the two events. Between the 1993 and 1995 events, the system may not have had sufficient time to disperse the sediment deposited at the upstream end of the Pannerden channel. Another reason for the response to the 1993 and 1995 peak flows to differ from previous events may be that the channel bed surface within the region of interest has coarsened significantly. This study illustrates the importance of peak flows regarding bifurcation dynamics, and further research is focused on the interaction between bifurcation dynamics and the dynamics of the larger-scale system.@en

While most of the world's large rivers are heavily engineered, channel response to engineering measures on decadal to century and several 100 km scales is scarcely documented. We investigate the response of the Lower Rhine River (Germany-Netherlands) to engineering measures, in terms of channel slope and bed surface grain size. Field data show domain-wide incision, primarily associated with extensive channel narrowing. Remarkably, the channel slope has increased in the upstream end, which is uncommon under degradational conditions. We attribute the observed response to two competing mechanisms: bedrock at the upstream boundary increases the channel slope over the upstream part of the alluvial reach to compensate for the reduction of net annual sediment mobility, and extensive channel narrowing reduces the equilibrium slope. Another striking feature is the advance and flattening of the gravel-sand transition, suggesting its gradual fading due to an increasingly reduced slope difference between the gravel and sand reaches.

@en

The objective of this thesis is to assess the effects of climate change on the initial, transient and equilibrium response of mixed-sediment river. Climate change will cause changes to hydrographs, which in turn affects sediment transport capacity and thereby the river bed profile. In order to create a better understanding of the processes in play we analyse a theoretical situation using the Lower Rhine as reference. The river Rhine is a heavily engineered river. The exact effects of climate change are therefore difficult to predict, as human intervention is significant. A simplified version of the Rhine can however give an insight into river response to climate change.

First an analyses of historical discharge and future discharge of the Lower Rhine was made, which gave an insight in the changes to the hydrograph that can be expected in the future. Next a model was made representing a theoretical river reach. For this river reach the upstream hydrodynamic boundary is varied and the results of these variations are analysed. Using the results of this model and the hydrograph expectations it was discussed what can be expected to happen to a river bed under these circumstances.

Sediment nourishments have become an increasingly attractive alternative to deal with continuous bed degradation problems in the river system. When supplied to the river, sediment nourishments induce a sediment wave that propagates through the system causing changes in bed level and bed surface texture. In the present study we use a one-dimensional numerical model to analyse how the propagation of these sediment waves varies under different conditions. We find that the propagation of the sediment wave tends to be predominantly dispersive for large spatial variations of the flow, which is the dominant case when the height of the wave is large in comparison to the flow depth, or when the sediment is dispersed over a large enough reach to lead to a backwater effect. In artificial sediment nourishments the height of the sediment wave is usually limited to prevent creating an additional obstacle to navigation, and it is then that the temporal variations in the flow, due to high discharge events for example, determine the propagation of the wave. We study the morphodynamic response in a mixed-size sediment river system to a sediment nourishment, and identify the physical mechanisms that impact the propagation of mixed-sized sediment waves. In order to prevent it from being easily flushed downstream, the grain size distribution of the nourished sediment is typically selected to be coarser than the river bed. This reduces sediment mobility and enhances sediment deposition over the nourished reach, which in turn causes a deficit in sediment supply to the downstream reach and leads to the formation of a degradational wave. The propagation celerity of the degradational wave is found to be significantly faster than that of the sediment wave, which means that the effects of the nourishments in a river system can be observed significantly further downstream than from where the front of the sediment wave is found. We find that the celerities of both the sediment wave and the degradational wave are predominantly affected by the adjustments in surface grain size in the system. These adjustments migrate in the river system as a fining wave, which both accelerates the sediment wave, and decelerates the degradational wave. We learn that higher content of the finer fractions in the nourished sediment accelerate the evolution of the wave and lead to a predominantly translational behaviour. It also affects the magnitude of the incision depth, the closer the composition of the nourished sediment is to that of the bed surface, the less severe the scour. And even though it does not affect the propagation celerity of the degradational wave, it does impact how far downstream the wave travels, affecting a larger river reach the coarser the grain size composition of the sediment wave is. Finally, we analyse the results from a nourishment pilot project carried out in the Bovenrijn from 2016-2019. We find that the propagation of the sediment wave associated with the nourishment is only observed during high discharge events. The wave is predominantly dispersive, and only during a prolonged extreme discharge event did the wave show slight translational behaviour. We observe that lateral sorting mechanisms caused the tracer sediment to have a different trajectory from the sediment wave, which highlights the limitations of using a one-dimensional model. We also observe that the changes in bed level are related to the changes in surface grain size, with sediment deposition over the nourished reach and bed degradation just downstream.

The Rhine reach comprising the Niederrhein and the Upper Rhine Delta is characterized by a long history of engineering interventions. The area is densely populated and the Rhine, the main inland waterway in Europe, justifying the need for such measures. Interventions have had a strong impact on morphodynamic channel characteristics. We use bed level data collected since 1926 to assess large scale bed level change in the Upper Rhine Delta and Niederrhein. The field data show various trends both at the large scale and locally. A reduction of the main channel slope is observed in the Waal and the IJssel branches of the Rhine Delta, achieved through bed degradation. In the Pannerdensch Kanaal-Nederrijn-Lek branch, the morphoynamics are controlled by three weirs, downstream of which degradation has developed. The downstream half of the Niederrhein has degraded, resulting in a slope increase. Most of the domain has remained relatively stable for at least 20 years. The observed trends provide insight on the river response to interventions, which helps to better understand degradation processes, in model calibration, and in anticipating future change of channel geometry.@en