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

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

Surveyed and modelled one‐year morphodynamics in the braided lower Tana River

Field measurements and morphodynamic simulations were carried out along a 5‐km reach of the sandy, braided, lower Tana River in order to detect temporal and spatial variations in river bed modifications and to determine the relative importance of different magnitude discharges on river bed and braid channel evolution during a time span of one year, i.e. 2008–2009. Fulfilling these aims required testing the morphodynamic model's capability to simulate changes in the braided reach. We performed the simulations using a 2‐D morphodynamic model and different transport equations. The survey showed that more deposition than erosion occurred during 2008–2009. Continuous bed‐load transport and bed elevation changes of ±1 m, and a 70–188‐m downstream migration of the thalweg occurred. Simulation results indicated that, during low water periods, modifications occurred in both the main channel and in other braid channels. Thus, unlike some gravel‐bed rivers, the sandy lower Tana River does not behave like a single‐thread channel at low discharge. However, at higher discharge, i.e. exceeding 497 m³/s, the river channel resembled a single‐thread channel when channel banks confined the flow. Although the spring discharge peaks caused more rapid modifications than slower flows, the cumulative volumetric changes of the low water period were greater. The importance of low water period flows for channel modifications is emphasized. Although the 2‐D model requires further improvements, the results were nevertheless promising for the future use of this approach in braided rivers. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling of meander migration in an incised channel

An updated linear computer model for meandering rivers with incision has been developed. The model simulates the bed topography, flow field, and bank erosion rate in an incised meandering channel. In a scenario where the upstream sediment load decreases （e.g., after dam closure or soil conservation）, alluvial river experiences cross section deepening and slope flattening. The channel migration rate might be affected in two ways： decreased channel slope and steeped bank height. The proposed numerical model combines the traditional one-dimensional （1D） sediment transport model in simulating the channel erosion and the linear model for channel meandering. A non-equilibrium sediment transport model is used to update the channel bed elevation and gradations. A linear meandering model was used to calculate the channel alignment and bank erosion/accretion, which in turn was used by the 1D sediment transport model. In the 1D sediment transport model, the channel bed elevation and gradations are represented in each channel cross section. In the meandering model, the bed elevation and gradations are stored in two dimensional （2D） cells to represent the channel and terrain properties （elevation and gradation）. A new method is proposed to exchange information regarding bed elevations and bed material fractions between 1D river geometry and 2D channel and terrain. The ability of the model is demonstrated using the simulation of the laboratory channel migration of Friedkin in which channel incision occurs at the upstream end.     

Morphodynamic evolution of a lower Mississippi River channel bar after sand mining

In‐channel sand mining by dredge removes large quantities of bed sediment and alters channel morphodynamic processes. While the reach‐scale impacts of dredging are well documented, the effects of the dredged borrow pit on the local flow and sediment transport are poorly understood. These local effects are important because they control the post‐dredge evolution of the borrow pit, setting the pit lifespan and affecting reach‐scale channel morphology. This study documents the observed morphological evolution of a large (1·46 million m³) borrow pit mined on a lateral sandbar in the lower Mississippi River using a time‐series of multibeam bathymetric surveys. During the 2·5 year time‐series, 53% of the initial pit volume infilled with sediment, decreasing pit depth by an average of 0·88 m yr^?1. To explore the controls of the observed infilling, a morphodynamic model (Delft3D) was used to simulate flow and sediment transport within the affected river reach. The model indicated that infilling rates were primarily related to the riverine sediment supply and pit geometry. The pit depth and length influenced the predicted magnitude of the pit bed shear stress relative to its pre‐dredged value, i.e. the bed‐stress reduction ratio (R*), a metric that was correlated with the magnitude and spatial distribution of infilling. A one‐dimensional reduced‐complexity model was derived using predicted sediment supply and R* to simulate patterns of pit infilling. This simplified model of borrow‐pit evolution was able to closely approximate the amount and patterns of sediment deposition during the study period. Additional model experiments indicate that, for a borrow pit of a set volume, creating deep, longitudinally‐shorter borrow pits significantly increased infilling rates relative to elongated pits. Study results provide insight into the resilience of alluvial river channels after a disturbance and the sustainability of sand mining as a sediment source for coastal restoration. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupled hierarchical modeling approach to simulating the geomorphic response of river systems to anthropogenic climate change

Anthropogenic climate change is expected to change the discharge and sediment transport regime of river systems. Because rivers adjust their channels to accommodate their typical inputs of water and sediment, changes in these variables can potentially alter river morphology. In this study, a hierarchical modeling approach was developed and applied to examine potential changes in reach‐averaged bedload transport and spatial patterns of erosion and deposition for three snowmelt‐dominated gravel‐bed rivers in the interior Pacific Northwest. The modeling hierarchy was based on discharge and suspended‐sediment load from a basin‐scale hydrologic model driven by a range of downscaled climate‐change scenarios. In the field, channel morphology and sediment grain‐size data for all three rivers were collected. Changes in reach‐averaged bedload transport were estimated using the Bedload Assessment of Gravel‐bedded Streams (BAGS) software, and the Cellular Automaton Evolutionary Slope and River (CAESAR) model was used to simulate the spatial pattern of erosion and deposition within each reach to infer potential changes in channel geometry and planform. The duration of critical discharge was found to control bedload transport. Changes in channel geometry were simulated for the two higher‐energy river reaches, but no significant morphological changes were found for a lower‐energy reach with steep, cohesive banks. Changes in sediment transport and river morphology resulting from climate change could affect the management of river systems for human and ecological uses. Copyright © 2015 John Wiley & Sons, Ltd.     

Long-term morphological river response to hydrological changes

Hydrologic changes have a great impact on the long-term river morphology. The most common anthropic cause is the construction of dams, which often reduces both the discharge regime and sediment transport, producing a narrowing and degradation of the river bed. In this study we propose a simple, lumped morphodynamic model that describes fluvial cross-section dynamics consequent to changes in discharge and sediment transport induced by external factors. The model provides the temporal dynamics of the river width and bed quote. These dynamics result not to be trivial and can exhibit non monotonic behavior, with aggradations and degradations, and narrowing and widening. The model has been tested on real rivers using data obtained from field studies. The agreement between the outcomes and the data measured in the field works is always satisfactory.     

EFFECTS OF RIVER TRAINING STRUCTURES ON THE SEDIMENT BALANCE OF A SAND - BEARING RIVER

1. INTROOrCTIOXThe upper Tisza in Hungary is a sand bed alluvial river. presenting several river training problems. In order to solve these problems first the laws of free river--bed evolution must be clarified thenthe effects of existing river training s…     

Numerical modeling of bedload and suspended load contributions to morphological evolution of the Seine Estuary (France)

This numerical modeling study (i) assesses the influence of the sediment erosion process on the sediment dynamics and subsequent morphological changes of a mixed-sediment environment, the macrotidal Seine estuary, when non-cohesive particles are dominant within bed mixtures (non-cohesive regime), and (ii) investigates respective contributions of bedload and suspended load in these dynamics. A three dimensional (3D) process-based morphodynamic model was set up and run under realistic forcings (including tide, waves, wind, and river discharge) during a 1-year period. Applying erosion homogeneously to bed sediment in the non-cohesive regime, i.e., average erosion parameters in the erosion law (especially the erodibility parameter, E₀), leads to higher resuspension of fine sediment due to the presence of coarser fractions within mixtures, compared to the case of an independent treatment of erosion for each sediment class. This results in more pronounced horizontal sediment flux (two-fold increase for sand, +30% for mud) and erosion/deposition patterns (up to a two-fold increase in erosion over shoals, generally associated with some coarsening of bed sediment). Compared to observed bathymetric changes, more relevant erosion/deposition patterns are derived from the model when independent resuspension fluxes are considered in the non-cohesive regime. These results suggest that this kind of approach may be more relevant when local grain-size distributions become heterogeneous and multimodal for non-cohesive particles. Bedload transport appears to be a non-dominant but significant contributor to the sediment dynamics of the Seine Estuary mouth. The residual bedload flux represents, on average, between 17 and 38% of the suspended sand flux, its contribution generally increasing when bed sediment becomes coarser (can become dominant at specific locations). The average orientation of residual fluxes and erosion/deposition patterns caused by bedload generally follow those resulting from suspended sediment dynamics. Sediment mass budgets cumulated over the simulated year reveal a relative contribution of bedload to total mass budgets around 25% over large erosion areas of shoals, which can even become higher in sedimentation zones. However, bedload-induced dynamics can locally differ from the dynamics related to suspended load, resulting in specific residual transport, erosion/deposition patterns, and changes in seabed nature.     

2D NUMERICAL SIMULATION OF FLOOD AND FLUVIAL PROCESS IN THE MEANDERING AND ISLAND-BRAIDED MIDDLE YANGTZE RIVER

The characteristics of water flow and sediment transport in a typical meandering and island-braided reach of the middle Yangtze River is investigated using a two-dimensional （2D） mathematical model. The major problems studied in the paper include the carrying capacity for suspended load, the incipient velocity and transport formula of non-uniform sediment, the thickness of the mixed layer on the riverbed, and the partitioning of bed load and suspended load. The model parameters are calibrated using extensive field data. Water surface profiles, distribution of flow velocities, riverbed deformation are verified with site measurements. The model is applied to a meandering and island-braided section of the Wakouzi-Majiazui reach in the middle Yangtze River, which is about 200 km downstream from the Three Gorges Dam, to study the training scheme of the navigation channels. The model predicts the processes of sediment deposition and fiver bed erosion, changes of flow stage and navigation conditions for the first 20 years of impoundment of the Three Gorges Project.     

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…     

Sedimentation of overbank floods in the confined complex channel–floodplain system of the Lower Yellow River,China

The channel boundary conditions along the Lower Yellow River (LYR) have been altered significantly since the 1950s with the continual reinforcement and construction of both main and secondary dykes and river training works. To evaluate how the confined complex channel–floodplain system of the LYR responds to floods, this study presents a detailed investigation of the relationship between the tempo‐spatial distribution of sedimentation/erosion and overbank floods occurred in the LYR. For large overbank floods, we found that when the sediment transport coefficient (ratio of sediment concentration of flow to flow discharge) is less than 0.034, the bankfull channel is subject to significant erosion, whereas the main and secondary floodplains both accumulate sediment. The amount of sediment deposited on the main and secondary floodplains is closely related to the ratio of peak discharge to bankfull discharge, volume of water flowing over the floodplains, and sediment concentration of overbank flow, whereas the degree of erosion in the bankfull channel is related to the amount of sediment deposited on the main and secondary floodplains, water volume, and sediment load in flood season. The significant increase in erosion in the bankfull channel is due to the construction of the main and secondary dykes and river training works, which are largely in a wide and narrow alternated pattern along the LYR such that the water flowing over wider floodplains returns to the channel downstream after it drops sediment. For small overbank floods, the bankfull channel is subject to erosion when the sediment transport coefficient is less than 0.028, whereas the amount of sediment deposited on the secondary floodplain is associated closely with the sediment concentration of flow. Over the entire length of the LYR, the situation of erosion in the bankfull channel and sediment deposition on the main and secondary floodplains occurred mainly in the upper reach of the LYR, in which a channel wandering in planform has been well developed.     

Peak discharge increase in hyperconcentrated floods

The peak of river floods usually decreases in the downstream direction unless it is compensated by freshwater inflow from tributaries. In the Yellow River (China) the opposite is regularly observed, where the peak discharge of river floods increases in the downstream direction (at a rate far exceeding the contribution from tributaries). This flood peak discharge increase is probably related to rapid morphological changes, to a modified bed friction, or to a combination of both. Yet the relative role of these processes is still poorly understood. This paper aims to analyze the relative contribution of bed erosion and friction change to the peak discharge increase, based on available data and a recently developed numerical model. Using this high-resolution, fully coupled morphodynamic model of non-capacity sediment transport, two hyperconcentrated floods characterized by downstream peak discharge increase are numerically reproduced and analyzed in detail. The results reveal that although erosion effects may contribute to the downstream discharge increase (especially in case of extreme erosion), for most cases the increase must be mainly due to a reduction in bed friction during peak discharge conditions. Additionally, based on the concept of channel storage reduction, the effects of decreasing bed friction and (very strong) bed erosion can be integrated in explaining the peak discharge increase.     

Erosion, sediment transportation and accumulation in rivers

The present paper analyses the interrelation between erosion, sediment transportation and accumulation proposed by N. I. Makkaveyev （1908-1983） and its further development in modem studies of river channel processes in Russia. Spatio-temporal linkages between erosion and accumulation are defined considering channel processes at different scales - river longitudinal profile, channel morphological patterns, alluvial bedforms （bars, dunes） and individual sediment particles. Relations between river geomorphic activity, flow transportation capacity and sediment budgets are established （sediment input and output; channel bed erosion and sediment entrainment into flow - termination of sediment transport and its deposition）. Channel planforms, floodplain segments separated by the latter and alluvial channel bedforms are shown to be geomorphic expressions of sediment transport process at different spatial and temporal scales. This paper is dedicated to the 100th anniversary of N. I. Makkaveyev, Professor of the Moscow State University, author of the book ＂River channel and erosion in its basin＂ （1955）. That book is regarded in Russia as the pioneering work which initiated the complex hydrological and geographical studies of channel processes and laid a basis for the theory of unified fluvial erosion-accumulation process.     

Implications of climate change in the twenty‐first century for simulated magnitude and frequency of bed‐material transport in tributaries of the Saint‐Lawrence River

More frequent extreme flood events are likely to occur in many areas in the twenty‐first century due to climate change. The impacts of these changes on sediment transport are examined at the event scale using a 1D morphodynamic model (SEDROUT4‐M) for three tributaries of the Saint‐Lawrence River (Québec, Canada) using daily discharge series generated with a hydrological model (HSAMI) from three global climate models (GCMs). For all tributaries, larger flood events occur in all future scenarios, leading to increases in bed‐material transport rates, number of transport events and number of days in the year where sediment transport occurs. The effective and half‐load discharges increase under all GCM simulations. Differences in flood timing within the tributaries, with a shift of peak annual discharge from the spring towards the winter, compared to the hydrograph of the Saint‐Lawrence River, generate higher sediment transport rates because of increased water surface slope and stream power. Previous research had shown that channel erosion is expected under all GCMs' discharge scenarios. This study shows that, despite lower bed elevations, flood risk is likely to increase as a result of higher flood magnitude, even with falling base level in the Saint‐Lawrence River. Copyright © 2010 John Wiley & Sons, Ltd.     

Sediment mobility and bed armoring in the St Clair River: insights from hydrodynamic modeling

The lake levels in Lake Michigan‐Huron have recently fallen to near historical lows, as has the elevation difference between Lake Michigan‐Huron compared to Lake Erie. This decline in lake levels has the potential to cause detrimental impacts on the lake ecosystems, together with social and economic impacts on communities in the entire Great Lakes region. Results from past work suggest that morphological changes in the St Clair River, which is the only natural outlet for Lake Michigan‐Huron, could be an appreciable factor in the recent trends of lake level decline. A key research question is whether bed erosion within the river has caused an increase in water conveyance, therefore, contributed to the falling lake level. In this paper, a numerical modeling approach with field data is used to investigate the possibility of sediment movement in the St Clair River and assess the likelihood of morphological change under the current flow regime. A two‐dimensional numerical model was used to study flow structure, bed shear stress, and sediment mobility/armoring over a range of flow discharges. Boundary conditions for the numerical model were provided by detailed field measurements that included high‐resolution bathymetry and three‐dimensional flow velocities. The results indicate that, without considering other effects, under the current range of flow conditions, the shear stresses produced by the river flow are too low to transport most of the coarse bed sediment within the reach and are too low to cause substantial bed erosion or bed scour. However, the detailed maps of the bed show mobile bedforms in the upper St Clair River that are indicative of sediment transport. Relatively high shear stresses near a constriction at the upstream end of the river and at channel bends could cause local scour and deposition. Ship‐induced propeller wake erosion also is a likely cause of sediment movement in the entire reach. Other factors that may promote sediment movement, such as ice cover and dredging in the lower river, require further investigation. Copyright © 2012 John Wiley & Sons, Ltd.     

THREE DIMENSIONAL NUMERICAL MODELLING OF FLOW AND SEDIMENT TRANSPORT IN RIVERS

The 3D numerical model, ECOMSED (open source code), was used to simulate flow and sediment transport in rivers. The model has a long history of successful applications to oceanic, coastal and estuarine waters. Improvements in the advection scheme, treatment of river roughness parameterization and shear stress partitioning were necessary to reproduce realistic and comparable results in a river application. To account for the dynamics of the mobile bed boundary, a model for the bed load transport was included in the code. The model reproduced observed secondary currents, bed shear stress distribution and erosion-deposition patterns on a curved channel. The model also successfully predicted the general flow patterns and sediment transport characteristics of a 1-km long reach of the River Klar?lven, located in the north of the county of V?rmland, Sweden.     

Morphological responses in a meandering and island-braided reach of the Middle Yangtze River to the Three Gorges Reservoir impoundment

Lateral migration is an important form of morphological changes on the Middle Yangtze River (MYR), particularly for the lower Jingjiang reach. The Three Gorges Reservoir (TGR) has substantially reduced sediment supply to the downstream river channels since its impoundment in June 2003. The scientific understanding of how decrease of sediment influences the processes of bank erosion and channel adjustments is complex and limited. In the present paper, the morphological responses in a typical meandering and island-braided river segment of the MYR to the filling of the TGR were investigated by a 3-D morphodynamic model. The potential of the 3-D model has been demonstrated by the observed data. The morphological evolutions in the Shishou bend during the first 12 years of the TGR impound-ment were predicted. The effects of the TGR operation on the planform evolutions in the study reach were analyzed based on the simulated results. Sediment load is decreased by 75%due to the early filling of the TGR. The magnitude of bed degradation with less sediment load due to the TGR operation is increased compared with the pre-dam situation. Qualitatively, the overall planform evolution trends in the Shishou bend after the TGR operation are similar to that without the TGR operation. The magnitude of lateral migration has been increased in some part of the channel bend, where the morphological response of the TGR operation exhibits more lateral migration rather than vertical degradation. Scouring at the bank toe enhances bank failure. Decrease of sediment load and weak bank anti-scour ability as well as the significant helical flow can be responsible for intensified bank erosion in the channel bend.     

Estimating bed material fluxes upstream and downstream of a controlled large bifurcation - the Mississippi-Atchafalaya River diversion

Bed material transport at river bifurcations is crucial for channel stability and downstream geomorphic dynamics. However, measurements of bed material transport at bifurcations of large alluvial rivers are difficult to make, and standard estimates based on the assumption of proportional partitioning of flow and bedload transport at bifurcations may be erroneous. In this study, we employed a combined approach based on observed topographic change (erosion/deposition) and bed material transport predicted from a one-dimensional model to investigate bed material fluxes near the engineering-controlled Mississippi-Atchafalaya River diversion, which is of great importance to sediment distribution and delivery to Louisiana's coast. Yang's (1973) sediment transport equation was utilized to estimate daily bed material loads upstream, downstream, and through the diversion during 2004–2013. Bathymetric changes in these channels were assessed with single beam data collected in 2004 and 2013. Results show that over the study period, 24% of the Mississippi River flow was diverted into the Atchafalaya River, while the rest remained in the mainstem Mississippi. Upstream of the diversion, the bed material yield was predicted to be 201 million metric tons (MT), of which approximately 35 MT (i.e., 17%) passed through the bifurcation channel to the Atchafalaya River. The findings from this study reveal that in the mainstem Mississippi, the percentage of bed material diversion (83%) is larger than the percentage of flow diversion (76%); Conversely, the diversion channel receives a disproportionate amount of flow (24%) relative to bed material supply (17%). Consequently, severe bed scouring occurred in the controlled Outflow Channel to the Atchafalaya River, while riverbed aggradation progressed in the mainstem Mississippi downstream of the diversion structures, implying reduced flow capacity and potential risk of a high backwater during megafloods. The study demonstrates that Yang's sediment transport equation provides plausible results of bed material fluxes for a highly complicated large river diversion, and that integration of the sediment transport equation with observed morphological changes in riverbed is a valuable approach to investigate sediment dynamics at controlled river bifurcations.     

Numerical study on the effects of floodplain vegetation on river planform style

The effects of floodplain vegetation on river planform have been investigated for a medium‐sized river using a 2D morphodynamic model with submodels for flow resistance and plant colonization. The flow resistance was divided into a resistance exerted by the soil and a resistance exerted by the plants. In this way it was possible to reproduce both the decrease in bed shear stress, reducing the sediment transport capacity of the flow within the plants, and the increase in hydraulic resistance, reducing the flow velocities. Colonization by plants was obtained by instantaneously assigning vegetation to the areas that became dry at low water stages. This colonization presents a step forward in the modelling of bank accretion. Bank erosion was related to bed degradation at adjacent wet cells. Bank advance and retreat were reproduced as drying and wetting of the computational cells at the channel margins. The model was applied to a hypothetical case with the same characteristics as the Allier River (France). The river was allowed to develop its own geometry starting from a straight, uniform, channel. Different vegetation densities produced different planforms. With bare floodplains, the river always developed a braided planform, even if the discharge was constant and below bankfull. With the highest vegetation density (grass) the flow concentrated in a single channel and formed incipient meanders. Lower vegetation density (pioneer vegetation) led to a transitional planform, with a low degree of braiding and distinguishable incipient meanders. The results comply with flume experiments and field observations reported in the literature.