Modeling Salt Water Intrusion in Tanshui River Estuarine System—Case-Study Contrasting Now and Then

A vertical (laterally integrated) two-dimensional numerical model was applied to study the salt water intrusion in the Tanshui River estuarine system, Taiwan. The river system has experienced dramatic changes in the past half century because of human intervention. The construction of two reservoirs and water diversion in the upper reaches of the river system significantly reduces the freshwater inflow. The land subsidence within the Taipei basin and the enlargement of the river constriction at Kuan-Du have lowered the river bed. Both changes have contributed farther to the intrusion of tidal flow and salt water in the upstream direction. The model was reverified with the earliest available hydrographic data measured in 1977. The overall performance of the model is in reasonable agreement with the field data. The model was then used to investigate the change in salt water intrusion as a result of reservoir construction and bathymetric changes in the river system. The model simulation study reveals that significant salinity increases have resulted from the combined changes. It has been speculated by ecological researchers that the long-term increase in salinity might be the driving force altering the aquatic ecosystem structure in the lower reach of the estuary and the Kuan-Du mangrove swamp, particularly the enlargement of the mangrove area and the disappearance of freshwater marshes. However, concrete proof has not been available since no prototype salinity data were available prior to the reservoir construction. This case study offers the first quantitative estimate of the salinity changes due to human interference in this natural system.     

Compatibility of Reservoir Sediment Flushing and River Protection

In this paper, we propose a system of numerical models for the compatibility assessment of reservoir sediment flushing and protection of downstream river environments. The model system is made up of two simulation models. The first model simulates soil erosion in watershed slopes and sediment transport in the tributary of the reservoir by means of a weighted essentially nonoscillatory (WENO) method, which is conservative and fourth-order accurate in space and time. The second model simulates velocity and suspended solid concentration fields in the reservoirs. This model is based on the three-dimensional (3D) numerical integration of motion and concentration equations, expressed in contravariant form on a generalized boundary-conforming curvilinear coordinate system by using a conservative and higher-order accurate numerical scheme. The proposed system of models is applied to the Pieve di Cadore (Veneto, Italy) reservoir and to its catchment area. By comparing suspended solid concentrations that are discharged through the bottom outlets during flushing operations with suspended solid concentrations in the main river during natural flooding, we perform an assessment of the compatibility between sediment flushing and the protection of the river ecosystem downstream.     

Modeling the Effect of Flow and Sediment Transport on White Sturgeon Spawning Habitat in the Kootenai River, Idaho

Kootenai River white sturgeon spawn in an 18-km reach of the Kootenai River, Id. Since completion of Libby Dam upstream from the spawning reach in 1972, 1974 is the only year with documented significant recruitment of juvenile fish. Where successful in other rivers, white sturgeon spawn over clean coarse material of gravel size or larger. The channel substrate in the current (2008) 18-km spawning reach is composed primarily of sand and some buried gravel; within a few kilometers upstream there is an extended reach of clean gravel, cobble, and bedrock. We used a quasi-three-dimensional flow and sediment-transport model along with the locations of collected sturgeon eggs as a proxy for spawning location from 1994 to 2002 to gain insight into spawning-habitat selection in a reach which is currently unsuitable due to the lack of coarse substrate. Spatial correlations between spawning locations and simulated velocity and depth indicate fish select regions of higher velocity and greater depth within any river cross section to spawn. These regions of high velocity and depth occur in the same locations regardless of the discharge magnitude as modeled over a range of pre- and postdam flow conditions. A flow and sediment-transport simulation shows high discharge, and relatively long-duration flow associated with predam flow events is sufficient to scour the fine sediment overburden, periodically exposing existing lenses of gravel and cobble as lag deposits in the current spawning reach. This is corroborated by video observations of bed surface material following a significant flood event in 2006, which show gravel and cobble present in many locations in the current spawning reach. Thus, both modeling and observations suggest that the relative rarity of extremely high flows in the current regulated flow regime is at least partly responsible for the lack of successful spawning; in the predam flow regime, frequent high flows removed the fine sediment overburden, unveiling coarse material and providing suitable substrate in the current spawning reach.     

Calibration of a Continuous Simulation Fecal Coliform Model Based on Historical Data Analysis

Since 1984, the major water reclamation plants discharging to the Chicago Waterway System (CWS) have not disinfected their effluents. The possible addition of disinfection at these plants is the subject of an ongoing use attainability analysis (UAA). For the UAA, Escherichia coli (E. coli) is used as the indicator of bacterial contamination. However, only a few years of E. coli data are available for the CWS and the treatment plants discharging to the CWS. Thus, it was decided to develop a model based on fecal coliforms for which more data are available and to develop a relation between fecal coliform and E. coli counts for the CWS. A 1:1 relation was found between fecal coliform and E. coli counts in the CWS by Limnotech (2004, written communication) as part of the UAA. In order to evaluate the effects of possible disinfection measures on fecal coliform and related E. coli counts in the CWS, a simple first-order fecal coliform decay model was added to the continuous-simulation flow-water quality model DUFLOW applied to the CWS system. Due to the limited amount (monthly samples) of measured fecal coliform concentration data for the CWS, a reasonable calibration of the model would have been difficult to achieve based on the traditional trial and error method. In this paper, a new concept of model parameter estimation based on historical data analysis and its application to model calibration is presented. The fecal coliform decay rate k was estimated for every reach of the CWS based on analysis of historical data (1990–2003) between each two consecutive sampling locations and the related travel time between these stations. The fecal coliform decay rate then was determined on the basis of many years (14 years, in this case) of monthly fecal coliform samples rather than the few monthly samples taken in a typical calibration period. The results obtained indicate that the calibration process was successful, and a good match between measured and simulated fecal coliform concentrations at almost all locations along the CWS is achieved with one model run for several multiple month periods in 1998, 1999, 2001, and 2002.     

Coupling Bank Stability and Bed Deformation Models to Predict Equilibrium Bed Topography in River Bends

Attempts to include river bank erosion predictors within morphological models are becoming increasingly common, but uncertainty surrounds the procedures used to couple bed deformation and bank erosion submodels in a way that maintains the mass continuity of eroded bank sediment. Herein we present a coupling procedure that comprises two discrete elements. First, immediately following bank failure, the slumped debris comes to rest at the bank toe as a planar surface inclined at an “angle of repose.” Second, we use a fractional transport model to simulate the subsequent erosion and deposition of the failed bank material debris. The method is demonstrated with an example in which an equilibrium bed topography model is combined with a river bank erosion model to predict the morphological response of a river bend to a large flow event.     

Sensitivity Analysis and Predictive Uncertainty Using Inundation Observations for Parameter Estimation in Open-Channel Inverse Problem

The necessity of an improvement of flood risk management has been emphasized by the importance of the damages caused by several disastrous floods in recent years. An obligatory way to better understand and forecast these events is to be able to accurately simulate the water transport in the river. Some parameters embedded in open-channel flow equations are difficult to calibrate in practice and hydraulic models are currently constrained by lack of data. The objective of this study is to test the potential of data typical of remote sensing imagery which may be used together with data assimilation methods to ascertain the value of a set of open-channel model parameters. Indeed, for the modeling of ungauged rivers, it may be necessary to estimate open-channel flow and/or morphological characteristics of the river. The innovation of the proposed methodology is to reconstitute information about the geometry of the river from top sight. A sensitivity analysis using the generalized sensitivity analysis method is carried out, together with an estimation of predictive uncertainty using the generalized likelihood uncertainty estimation, a methodology for calibration and uncertainty estimation of distributed models introduced by Beven and Binley (1992). The study aims at building a framework for parameter optimization in the case of ungauged rivers. The proposed method has been implemented on a reach of the Lèze River, located in the southwest of France.     

Structural Ice Control Alternatives for Middle Mississippi River

The middle Mississippi River, which extends from the mouth of the Missouri River near St. Louis to the confluence with the Ohio River at Cairo, Ill., is a critical navigation route throughout the year. During cold periods, the ice accumulations and ice jams that form on this reach can delay or suspend winter navigation, incurring great costs to industry as well as cities and towns whose economies depend on river commerce. Much of this ice originates in the Missouri River. With the onset of warmer air temperatures, the rapid release of these ice accumulations can result in substantial damage to river structures such as dikes, revetments, and levees. This study analyzed historical data and used numerical hydraulic models to assess the possibility of structural solutions to these ice problems. The study identified reach locations where structural ice control might be possible. A simple computer model then simulated the upstream progression of ice covers on the middle Mississippi to assess the feasibility of various ice control alternatives. An ice retention structure located on the Missouri River near its mouth was found to be the most favorable of the structural ice control options considered.     

Spur Failure in River Engineering

Spurs are river engineering elements used to protect river banks from erosion and to concentrate flow to the river axis. Today, spurs are also employed for promoting environmental conditions along a river bank. These elements are characterized by a large variety of geometrical parameters, of which none is definitely fixed. Based on a preliminary study that identified optimum spur arrangements in a straight river reach the present research project adds to the protection of spurs by riprap. The first two spurs are demonstrated to require a suitable riprap for promoting nearly uniform scour conditions along the entire spur reach. This research investigated the effects of a variety of parameters on spur flow, notably spur length, spur spacing, spur height plus the diameter and the number of riprap rows, along with the main hydraulic and granulometric parameters. Design equations were established based on a large experimental campaign to predict riprap failure in terms of the previous set of variables. The failure modes are described along with a novel failure mechanism. The paper ends with a set of limitations allowing for the application of these results in river engineering.     

Numerical Simulation of Longitudinal and Lateral Channel Deformations in the Braided Reach of the Lower Yellow River

Bank erosion frequently occurs in the Lower Yellow River (LYR), playing an important role in the evolution of this braided river. A two-dimensional (2D) composite model is developed herein that consists of a depth-averaged 2D flow and sediment transport submodel and a bank-erosion submodel. The model incorporates a new technique for updating bank geometry during either degradational or aggradational bed evolution, allowing the two submodels to be closely combined. Using the model, the fluvial processes in the braided reach of the LYR between Huayuankou and Laitongzhai are simulated, and the calculated results generally agree with the field measurements, including the water-surface elevation, variation of water-surface width, and variations of cross-sectional profiles. The calculated average water-surface elevation in the study reach was 0.09?m greater than the observed initial value, and the calculated mean bed elevation for six cross sections was 0.11?m lower than the observed value after 24 days. These errors are attributed to the large variability of flow and sediment transport processes. Sensitivity tests of three groups of parameters are conducted, and these groups of parameters are related to flow and sediment transport, bank erosion, and model application, respectively. Analysis results of parameter sensitivity tests indicate that bank erodibility coefficient and critical shear stress for bank material are sensitive to the simulated bank erosion process. The lateral erosion distance at Huayuankou will increase by 19% as the value of bank erodibility coefficient changes from 0.1 to 0.3, and it will decrease by 57% as the value of critical shear stress for bank increases from 0.6 to 1.2?N/m2. Limited changes of other parameters have relatively small effects on the simulated results for this reach, and the maximum change extent of calculated results is less than 5%. Because the process of sediment transport and bank erosion in the braided reach of the LYR is very complicated, further study is needed to verify the model.     

Optimal Cross-Sectional Spacing in Preissmann Scheme 1D Hydrodynamic Models

Choosing a suitable set of cross sections for the representation of the natural geometry of a river is important for the efficiency of one-dimensional (1D) hydraulic models, but only few guidelines are available for the selection of the most suitable distance between cross sections, depending on the hydraulic problem at hand. This issue is investigated by examining models of a ～ 55?km reach of the River Po, Italy, and a ～ 16?km reach of the River Severn, United Kingdom, for both of which high quality laser scanning altimetry are available. The high-resolution digital terrain models of the two river reaches enabled the construction of a series of hypothetical topographical ground surveys with different spacing between cross sections, which could be used as input to a standard 1D model (UNET). Both historical and synthetic flood events for the two river reaches were simulated, and the results were then analyzed to quantify the accuracy associated with each resolution and to assess how survey resolution impacts the performance of standard 1D models. The study results agree with the available suggestions in the literature and provide useful guidelines for 1D hydrodynamic modeling.     

Multicrack Detection on Semirigidly Connected Beams Utilizing Dynamic Data

The problem of crack detection has been studied by many researchers, and many methods of approaching the problem have been developed. To quantify the crack extent, most methods follow the model updating approach. This approach treats the crack location and extent as model parameters, which are then identified by minimizing the discrepancy between the modeled and the measured dynamic responses. Most methods following this approach focus on the detection of a single crack or multicracks in situations in which the number of cracks is known. The main objective of this paper is to address the crack detection problem in a general situation in which the number of cracks is not known in advance. The crack detection methodology proposed in this paper consists of two phases. In the first phase, different classes of models are employed to model the beam with different numbers of cracks, and the Bayesian model class selection method is then employed to identify the most plausible class of models based on the set of measured dynamic data in order to identify the number of cracks on the beam. In the second phase, the posterior (updated) probability density function of the crack locations and the corresponding extents is calculated using the Bayesian statistical framework. As a result, the uncertainties that may have been introduced by measurement noise and modeling error can be explicitly dealt with. The methodology proposed herein has been verified by and demonstrated through a comprehensive series of numerical case studies, in which noisy data were generated by a Bernoulli–Euler beam with semirigid connections. The results of these studies show that the proposed methodology can correctly identify the number of cracks even when the crack extent is small. The effects of measurement noise, modeling error, and the complexity of the class of identification model on the crack detection results have also been studied and are discussed in this paper.     

Case Study: Two-Dimensional Model Simulation of Channel Migration Processes in West Jordan River, Utah

The paper presents the application of a two-dimensional depth-averaged numerical model to simulate the lateral migration processes of a meandering reach in the West Jordan River in the state of Utah. A new bank erosion model was developed and then integrated with a depth-averaged two-dimensional hydrodynamic model. The rate of bank erosion is determined by bed degradation, lateral erosion, and bank failure. Because bank material in the West Jordan River is stratified with layers of cohesive and noncohesive materials, a specific bank erosion model was developed to consider stratified layers in the bank surface. This bank erosion model distinguishes itself from other models by relating bank erosion rate with not only flow but also sediment transport near the bank. Additionally, bank height, slope, and thickness of two layers in the bank surface were considered when calculating the rate of bank erosion. The developed model was then applied to simulate the processes of meandering migration in the study reach from 1981 to 1992. Historical real-time hydrographic data, as well as field survey data of channel geometry and bed and bank materials, were used as the input data. Simulated cross-sectional geometries after this 12-year period agreed with field measurements, and the R2 value for predicting thalweg elevation and bank shift are 0.881 and 0.706, respectively.     

Simplified Model for Assessing Meander Bend Migration Rates

This paper describes a simple analytical model for assessing meander bend migration rates. The new model is innovative in that it can be applied directly using output data from one-dimensional hydraulic models. The method is designed as a first approximation that can be applied rapidly to large reaches of river systems, and should not be considered as a replacement for more detailed bank stability assessments. The model is based on a conversion factor for one-dimensional shear stress. Comparison of the conversion factor with existing data sets is presented. Converted shear stresses are translated to bank erosion volumes, which in turn are related to meander bend migration rates, using bank geometry and an appropriate sediment transport law. Results are presented for an evaluation of meander bend migration rates for a restoration feasibility study for a project on the Williamson River, Oregon.     

Influence of Stratification and Shoreline Erosion on Reservoir Sedimentation Patterns

Sedimentation in the main pool of a deep (maximum depth: 50?m), 227?km2 hydropower reservoir was modeled using a three-dimensional numerical model of hydrodynamics and sedimentation for different wind, inflow, and outflow conditions. Short-term velocity measurements made in the reservoir were used to validate some aspects of the hydrodynamic model. The effects of thermal stratification on sedimentation patterns were investigated, since the reservoir is periodically strongly stratified. Stratification alters velocity profiles and thus affects sedimentation in the reservoir. Sedimentation of reservoirs is often modeled considering only the deposition of sediments delivered by tributaries. However, the sediments eroding from the shorelines can contribute significantly to sedimentation if the shorelines of the reservoir erode at sufficiently high rates or if sediment delivery via tributary inflow is small. Thus, shoreline erosion rates for a reservoir were quantified based on measured fetch, parameterized beach profile shape, and measured wind vectors, and the eroded sediments treated as a source within the sedimentation modeling scheme. The methodology for the prediction of shoreline erosion was calibrated and validated using digital aerial photos of the reservoir taken in different years and indicated approximately 1?m/year of shoreline retreat for several locations. This study revealed likely zones of sediment deposition in a thermally stratified reservoir and presented a methodology for integration of shoreline erosion into sedimentation studies that can be used in any reservoir.     

Case Study of Fluvial Modeling of River Responses to Dam Removal

A plan was made to remove Matilija Dam on the Ventura River. With dam removal, the delta in the reservoir and the downstream channel were expected to undergo major changes in morphology. The FLUVIAL-12 model was employed to simulate reservoir and river channel responses after dam removal. As a first step, the model was calibrated using the Ventura River data to establish its validity. In calibration, the model was used to simulate the fluvial processes starting from the time of dam completion. The simulated sediment deposition above the Matilija Dam matches closely with the deposits measured by the U.S. Bureau of Reclamation. A large amount of sediment was stored in the reservoir; some of the stored sediment would be transported downstream after dam removal. An important consequence of dam removal is the major increase of sediment release to the river channel downstream. The sediment supply to the downstream reach is not only from the stored sediment in the reservoir but also from natural sediment inflow from the upstream watershed. Therefore, sediment supply to the downstream reach will exceed the natural sediment flow before the dam presence. This situation tends to overload the downstream reach with sediment, resulting in excessive deposition. The amount of sediment release from the area above the removed dam is closely related to the changes in reservoir morphology. It is necessary to model changes in the channel bed profile and channel width during erosion in order to determine the amount of sediment removal. The amount of sediment release may not be simulated using an erodible-bed model but it may be determined using an erodible boundary model.     

Resolving Disputes over Reservoir-River Operation

The main objective of this work is to allocate water with acceptable quality from a reservoir-river system to different water users and to determine the maximum acceptable pollution load that each user discharges to the river. To achieve this, a bargaining model is developed based on a system dynamics approach, and it is compared with an alternative model based on the Nash bargaining theory. The quality of water released from the reservoir is assessed considering selective withdrawal schemes, obtained from schemes based on operating rules derived from an artificial neural network (ANN) and a K-nearest neighborhood (NN) model. A river water-quality simulation model is also linked with the two bargaining models. These models are then applied to the Karkheh reservoir-river system and using a 50-year monthly time series of qualitative and quantitative data of this system. Finally, the performance criteria for the results of the models are determined and compared. This study shows the significant value of the system dynamics approach in resolving disputes over water allocation among stakeholders/water users of reservoir-river systems considering the water-quality aspects.     

Field Measurements and Simulation of Bridge Scour Depth Variations during Floods

An understanding of bridge scour mechanisms during floods in a fluvial river is very important for cost-effective bridge foundation design. Reliable bridge scour data for flood events are limited. In this study, field experiments were performed at the Si-Lo Bridge in the lower Cho-Shui River, the longest river in Taiwan, to collect scour-depth data using a sliding magnetic collar, a steel rod, and a numbered-brick column. By separating each scour component, a methodology for simulating the temporal variations of the total scour depth under unsteady flow conditions is proposed. The proposed total-scour model integrates three scour components, namely general scour, contraction scour, and local scour. The collected field data, comprising both general scour and total scour depths, are used to validate the applicability of the proposed model. Based on the peak flow discharges during floods, a comparison of the local scour depths calculated using several commonly used equilibrium local scour formulas indicates that most equations may overestimate the local scour depth.     

Sediment Exchange between a River and Its Groyne Fields: Mobile-Bed Experiment

Experiments have been carried out in a mobile-bed laboratory flume in order to study the sediment exchange process between the main channel and the groyne fields. The flume represented half the width of a schematized river reach with a series of groynes. The experiment was designed to represent typical dimensions of the Dutch River Waal at a geometrical scale of 1:100. The conditions were set to guarantee bed load as well as suspended load sediment transport. Conditions with submerged and emerged groynes were investigated. In addition to traditional measurements, viz., bed-level changes, suspended sediment concentrations, and flow velocities, bed-form propagation was measured in two dimensions using a the particle image velocimetry technique. The results were analyzed with focus on sediment exchange mechanisms and sediment transport patterns. The results demonstrate that under all flow conditions there is a net import of sediment into the groyne fields. The prevailing transport mechanisms vary with the flow stage: if the groynes are emerged it is mainly advection by the primary circulation cell, whereas if the groynes are submerged it is rather residual advection by large-scale coherent flow structures (in a straight reach). Additional entrainment of sediment by enhanced turbulence complicates the erosion/deposition patterns.     

Development of Refined BOD and DO Models for Highly Polluted Kali River in India

Most commonly used river water quality models for biochemical oxygen demand (BOD) and dissolved oxygen (DO) simulations are mainly based on advection, decay, settling, and loading functions. Using these concepts, refined river water quality models for BOD and DO simulations are developed in the present work considering a large number of physically based parameters and input variables. The refined models developed can be transformed to some of the commonly used river water quality models, if physically based parameters and input variables are omitted or removed. To test the applicability of the refined models developed and commonly used models, a total of 732 water quality and flow data sets are collected during March 1999–February 2000 from 22 sampling stations of the River Kali in India. River Kali is a highly polluted river in India and receives continuous inflow of untreated point source pollution from municipal and industrial wastes and nonpoint source pollution from agricultural areas. Newton–Raphson technique is used to optimize the model parameters during calibration and the performance of different models are evaluated using error estimation, viz. standard error and mean multiplicative error, and correlation statistics (r2). The results indicate that the BOD–DO models proposed by Camp in 1963 provide better results in comparison to other commonly used models. Moreover, the refined models developed for BOD and DO simulations minimize error estimates and improve correlation between observed and computed BOD and DO values of River Kali.