Tracer Study of Mixing and Transport in the Upper Hudson River with Multiple Dams

In October 2001, ～ 0.2?mol of SF6 was injected into the upper Hudson River, a modified natural channel with multiple dams, at Ft. Edward, N.Y. The tracer was monitored for 7 days as it moved ～ 50?km downriver. The longitudinal evolution of the tracer distribution was used to estimate one-dimensional advection (9.0±0.2?km?d?1) and dispersion (17.3±4.0?m2?s?1) along the river axis. Comparison of these results to tracer studies on channels without dams suggests that dams reduce longitudinal dispersion below the value expected in a natural channel with the same discharge. SF6 loss through air–water gas exchange along the river and at two dams (10.7?m combined height) was estimated by observing decay in peak concentration. Losses at dams (approximately 50% per dam) were dominant. The estimated gas exchange at dams was compared to a simple model adapted from those available in literature. Small amounts of tracer were trapped in a canal segment ( ～ 5?km long) that parallels the river, where advection and dispersion were sharply reduced.     

Case Study: Equivalent Widths of the Middle Rio Grande, New Mexico

Successive reaches of the Rio Grande have maintained equivalent channel widths of 50 and 250?m, respectively, over long periods of time. It is hypothesized that alluvial channels adjust bed slope to match the long-term changes in channel width. Analytical relationships show that wider river reaches develop steeper slopes. A modeling approach using daily water and sediment discharges simulates the transient evolution of bed elevation changes. The analytical and numerical models are in very good agreement with the longitudinal profile measurements of the Bosque del Apache reach of the Rio Grande, NM, from 1992 to 1999. The slope of the 50?m wide reach was 50?cm/km and the slope of the 250?m wide reach of the same river increased to 80?cm/km. This unsteady daily transient model compares well with a steady transient solution at a constant discharge close to the mean annual flow. The transient slope adjustments can also be approximated with an exponential model. Accordingly, it takes about 20–25?years for the Rio Grande to achieve about 90% of its slope adjustment.     

3D Unsteady RANS Modeling of Complex Hydraulic Engineering Flows. I: Numerical Model

A general-purpose numerical method is developed for solving the full three-dimensional (3D), incompressible, unsteady Reynolds-averaged Navier-Stokes (URANS) equations in natural river reaches containing complex hydraulic structures at full-scale Reynolds numbers. The method adopts body-fitted, chimera overset grids in conjunction with a grid-embedding strategy to accurately and efficiently discretize arbitrarily complex, multiconnected flow domains. The URANS and turbulence closure equations are discretized using a second-order accurate finite-volume approach. The discrete equations are integrated in time via a dual-time-stepping, artificial compressibility method in conjunction with an efficient coupled, block-implicit, approximate factorization iterative solver. The computer code is parallelized to take full advantage of multiprocessor computer systems so that unsteady solutions on grids with 106 nodes can be obtained within reasonable computational time. The power of the method is demonstrated by applying it to simulate turbulent flow at R ? 107 in a stretch of the Chattahoochee River containing a portion of the actual bridge foundation located near Cornelia, Georgia. It is shown that the method can capture the onset of coherent vortex shedding in the vicinity of the foundation while accounting for the large-scale topographical features of the surrounding river reach.     

Riparian Vegetation Mapping for Hydraulic Roughness Estimation Using Very High Resolution Remote Sensing Data Fusion

For detailed hydraulic modeling, accurate spatial information of riparian vegetation patterns needs to be derived in automatic fashion. We propose a supervised classification for heterogeneous riparian corridors with a low number of spectrally separate classes using data fusion of a Quickbird image and LIDAR data. The approach considers nine land cover classes including three woody riparian species, brush, cultivated areas, grassland, urban infrastructures, bare soil and water. The classical “stacked vector” approach is adopted for data fusion, while the nonparametric weighted feature-extraction method and the pixel-oriented maximum likelihood algorithm are used for feature-reduction and classification purposes, respectively. We test the approach over a 14-km stretch of the Sieve River (Tuscany Region, Italy). A one-dimensional river modeling is applied over the study reach comparing the results of a classification-derived hydraulic roughness map and a traditional ground-based approach. Despite the complex study reach, the classification method produced encouraging accuracies (OKS = 0.77) and represents a useful tool to delineate application domains of flow resistance models suited to different hydrodynamic patterns (e.g., stiff/flexible vegetation). Hydraulic modeling results showed that the remotely derived floodplain roughness parameterization captures the equivalent Manning coefficient over 20 test cross sections with uncertainty distributions described by low mean and standard deviation values.     

Assessing Uncertainty in Mass Balance Calculation of River Nonpoint Source Loads

Several sources of uncertainty are considered in using field data to perform mass balance calculations to estimate nonpoint source (NPS) loads of total dissolved solids (TDS) and selenium (Se) to two reaches along the Lower Arkansas River in Colorado. This approach renders stochastic models of the mass balance equations for each river reach, where the input variables and their associated parameters are treated as random variables described by probability distributions. A data set collected in an intensive field effort conducted over several years is used in developing the models. Monte Carlo simulation solves the stochastic mass balance equations to describe distributions of possible values of the NPS loads. Results indicate that uncertainty in these calculated loads is sizeable. Annual average coefficient of variation (CV) in calculated total NPS TDS loads for sample periods within the two reaches ranges between 0.15 and 2.76, averaging about 1.1. For the Se mass balance along the downstream reach, the average CV in calculated loads is 0.23. The 90% prediction interval width for total NPS TDS load averaged over the 58 sample periods along the upstream reach is 13,691?(kg/day)/km, compared to an overall average mean load of 8,383?(kg/day)/km. It is 8,545?(kg/day)/km averaged over the 61 sample periods downstream, compared to an overall average mean load of 11,183?(kg/day)/km. For the Se load, the overall average 90% prediction interval width is also substantial compared to the overall average mean: 0.028?(kg/day)/km compared to 0.038?(kg/day)/km. Change in stored solute mass within the river over sample periods is found to be a major contributing factor to the calculation of NPS loads. Also, sensitivity analyses are performed that yield information on the relative influence that the degree of uncertainty in each random parameter has on the uncertainty in the calculated solute loads.     

Large-Scale Tracing of Ground Water with Sulfur Hexafluoride

Sulfur hexafluoride (SF6) was injected into a 9 km reach of the Santa Ana River in Orange County, CA, over a period of two weeks. The entire flow of this river, which averaged 2.8 m3 s?1, percolated into the ground in the field area. The tracer was monitored at wells near the river to determine subsurface flow patterns and flow times with an accuracy much greater than could be achieved using numerical simulations of ground-water flow. During the experiment, SF6 effectively tagged 3.7 × 106 m3 of water. The tracer plume was mapped in the subsurface for 18 months and indicates that linear ground-water velocities averaged about 2 km year?1. The tracer reached two wells adjacent to the river (about 200 m away) within three weeks, giving evidence that SF6 was not retarded significantly relative to the ground-water flow. This is in agreement with previous laboratory experiments.     

Soft Bedrock Erosion Modeling with a Two-Dimensional Depth-Averaged Model

Many rivers in Taiwan have steep slopes, are subject to typhoon-induced flood flows, and contain soft bedrock that is exposed at many locations and easily erodible. The occurrence of extensive bedrock erosion has been a major threat to river infrastructure at many locations. Soft bedrock erosion, therefore, is an important process to consider for river projects in Taiwan. In this study, bedrock erosion models are reviewed. A specific model is proposed by combining two existing models incorporating both the hydraulic and abrasive scour mechanisms. The proposed bedrock erosion model is incorporated into a two-dimensional mobile-bed model, and the integrated model is tested by simulating bedrock erosion downstream of the Chi-Chi weir on the Choshui River in Taiwan. A calibration study is performed to determine appropriate values of the model parameters based on two and a half years of measured data. The model is then assessed based on a verification study that compares model predictions of bedrock erosion of the same reach to two additional years of measured data. The bedrock erosion model is found to be suitable for the river reach studied. Further improvement, however, is still necessary, which points to potential future research.     

Effect of Channel Restoration on Flood Wave Attenuation

Stream channel restoration can increase flow storage and energy dissipation of passing flood waves. Elements of restoration design that can enhance attenuation include remeandering, which reduces channel slope and increases channel length relative to the floodplain; restoring channel-floodplain connectivity; and revegetating banks and the floodplain. Reestablishment of floodplain hydraulic function is increasingly a goal of restoration programs, yet the approximate magnitude of possible change to attenuation due to reach-scale restoration remains poorly quantified. We examined the efficacy of channel restoration on flood attenuation using restored reaches and synthetic reaches representing median dimensions of channel restoration projects in North Carolina (e.g., ～ 1?km in length). We applied an industry standard dynamic flood routing model (UNET in HEC-RAS) to route floods in impaired and restored reach models. Floods routed through field-based reach models either exhibited very small increases in attenuation, largely due to assumed increases in floodplain roughness, or a decrease in attenuation. Analysis demonstrated that attenuation of peak discharge is overall most sensitive to channel and valley slope, channel and floodplain roughness, and channel and valley length in decreasing order, but is dependent on flood magnitude. Restoration most impacted floods of intermediate magnitude (between 2- and 50-year return interval), particularly those confined to the channel under the impaired morphology but able to access the floodplain under the restored morphology. Restoration may rehabilitate a channel’s ability to attenuate small to intermediate floods by augmenting flood access to the floodplain, changing channel geometry, and enhancing channel and floodplain roughness over time. However, our study shows that the predominantly small scale of current channel restoration will provide minimally quantifiable enhancement to flood attenuation.     

Case Study: Application of the HEC-6 Model for the Main Stem of the Kankakee River in Illinois

This case study paper presents results on the application of the HEC-6 model to the main stem of the Kankakee River in Illinois, a distance of about 39.3?km. Modeling was performed to develop comprehensive plans for enhancing the aquatic habitats and also to forecast future sedimentation problems if specific management practices are implemented. The paper concentrates on the modeling aspects of this research. The extent of the model was from the Stateline Bridge to Kankakee Dam in Kankakee. The HEC-6 model, originally developed by the Hydrologic Engineering Center (HEC) of the U.S. Army Corps of Engineers, was adapted for this application. The model was run, calibrated, and verified for both the hydraulic and sediment components. The hydraulic component was calibrated through comparison of measured yearly hydrographs with computed values for three gauging stations on the river. The hydrologic component was verified for the same three gauging stations for two yearly hydrographs for 2 additional water years. The sediment component was calibrated with river cross-sectional data collected by the Illinois State Water Survey in 1980 and 1999. The calibrated and verified hydraulic and calibrated sediment components then were used to predict future changes in water surface elevations and thalweg elevations for a 20-year period beginning in 1999, the last date for which river cross-sectional data are available.     

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.     

Robotic Monitoring to Assess Impacts of Zebra Mussels and Assimilative Capacity for a River

Water quality impacts of zebra mussel metabolism over an infested 15?km reach of the Seneca River, N.Y., are documented, based on vertically and temporally detailed robotic monitoring at the reach boundaries during the summer through early fall intervals of 2?years. Substantial reductions over the study reach are documented for dissolved oxygen (DO), pH, fluorometric chlorophyll a, and turbidity, associated with the metabolism of this invader. Violations of New York State water quality standards for DO that would not be resolved by traditional manual monitoring programs were observed. The loss of assimilative capacity caused by the zebra mussel invasion is confounding rehabilitation efforts for a downstream polluted lake that had considered diversion of municipal effluent to the river. The critical role robotic monitoring units would play in an automated control system for an innovative strategy of time-variable river discharge of the effluent is described. Near-real time robotic monitoring provides a more detailed understanding of the impacts of zebra mussels on water quality than traditional less intensive manual measurements.     

Flood Inundation Modeling with an Adaptive Quadtree Grid Shallow Water Equation Solver

Flood risk studies require hydraulic modeling in order to estimate flow depths and other hydraulic variables in the floodplain for a wide range of input conditions. Currently there is a need to improve the computational efficiency of fully two-dimensional numerical models for large-scale flood simulation. This paper describes an adaptive quadtree grid-based shallow water equation solver and demonstrates its capability for flood inundation modeling. Due to the grid dynamically adapting to dominant flow features such as steep water surface gradients and wet-dry fronts, the approach is both efficient and accurate. The quadtree model is applied to a realistic scenario of flood inundation over an urban area of 36?km2, resulting from the flood defenses breaching at Thamesmead on the River Thames, United Kingdom. The results of the simulation are in close agreement with alternative predictions obtained using the commercially available software TUFLOW.     

Dynamic Wave Study of Flow in Tidal Channel System of San Juan River

In this work the complete equations of one-dimensional unsteady flow in open channels in integral form, and compatibility equations at the junctions of a channel network, are solved numerically. Analytical integration in space is used between each pair of consecutive irregular sections of a channel, and the nonprismatic term is expressed in terms of uncoupled functions of the geometry at the sections. The linearized system of equations for each time interval is solved by an elimination method based on a double-sweep algorithm. The model is applied to the estuary of the San Juan River in Venezuela, where oscillating currents by effect of semidiurnal tides take place and the amplitude of the wave at the mouth is amplified toward the inland direction. Alternating drying and filling is simulated by means of slight modifications in the bed geometry of upper river sections. Measured water elevation and flow rates available at two stations are used to calibrate the model, and a very accurate adjustment of the tidal levels observed in the river is obtained.     

Hydraulic Jet Control for River Junction Design of Yuen Long Bypass Floodway, Hong Kong

The Yuen Long Bypass Floodway (YLBF) was designed to collect flows from the Sham Chung River (SCR) and the San Hui Nullah (SHN) and to serve as a diversion channel of the Yuen Long Main Nullah (YLMN). Under a 200-year return period design condition, the floodway was designed (1) to divert a flow of approximately 38?m3/s from the supercritical YLMN flow and (2) to convey a total combined flow of 278?m3/s to downstream within acceptable flood levels. The success of the design depends critically on complicated junction flow interactions that cannot be resolved by 1D unsteady flow models. These features include the supercritical-subcritical flow transition at the San Hui-Floodway (SHN-YLBF) junction and the diversion of part of the supercritical flow from the Main Nullah (YLMN). A laboratory Froude scale physical model was constructed to study water stages and flow characteristics in the floodway and to investigate optimal design arrangements at channel junctions and transitions. This paper summarizes the main features of the unique river junction network, in particular the use of the hydraulic jet principle at the SHN-YLBF junction to lower flood levels. In addition, a numerical flow model is employed to study flow details at the river junctions. The model is based on the general 2D shallow water equations in strong conservation form. The equations are discretized using the total variation diminishing finite-volume method which captures the discontinuity in hydraulic jumps. The numerical model predictions are well supported by the laboratory data, and the theoretical and experimental results offer useful insights for the design of urban flood control schemes under tight space constraints.     

Sensitivity of Field Data Estimates in One-Dimensional Hydraulic Modeling of Channels

A hydraulic simulation study was carried out for Red Hill Creek, an urban stream channel in Ontario, Canada. Over one million simulations were conducted using the HEC-RAS4b model to evaluate the sensitivity of model predictions to field data accuracy, density and estimation techniques and provide guidance toward balancing human resource allocation with model accuracy. Increased cross-sectional discretization and improved estimates of floodplain roughness dominate the accuracy of the results of computed water surface elevations. For a range of field hours available for data collection, a sampling strategy focused on maximizing the number of sparsely detailed cross sections is shown to outperform a sampling strategy using the same number of field hours to sample fewer cross sections with a higher resolution.     

淮河流域水文随机方法的径流图化

Theoretical difficulties for mapping and for estimating river regime characteristics in a large-scale basin remain because of the nature of the variable under study:river flows are related to a specific area, I.e. The drainage basin, and are hierarchically organized in space through the river network with upstream-downstream dependencies. Another limitation is there are not enough gauge stations in developing countries. This presentation aims at de-veloping the hydro-stochastic approach for producing choropleth maps of average annual runoff and computing mean discharge along the main river network for a large-scale basin.The approach applied to mean annual runoff is based on geostatistical interpolation proce-dures coupled with water balance and data uncertainty analyses. It is proved by an applica-tion in the upstream at Bengbu in the Huaihe River Basin, a typical large-scale basin in China.Hydro-stochasitic approach in a first step interpolates to a regular grid net and in a second step the grid values are integrated along rivers. The interpolation scheme includes a con-straint to be able to account for the lateral water balance along the rivers. Grid runoff map with 10 km x 10 km resolution and the discharge map along the river with the 1 km basic length unit are the main results in this study. This kind of statistic approach can be widely used be-cause it avoids the complexity of hydrological models and does not depend on the meteoro-logical data.     

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.     

Dual Discharge Approach to Accessing Assimilative Capacity: Probabilistic Analysis and Management Application

A dual discharge strategy has been proposed for management of the effluent from the Syracuse Metropolitan Treatment Plant (Metro). The approach involves routing the discharge to the Seneca River when assimilative capacity is available there and to Onondaga Lake when it is not. Application of a deterministic modeling approach has demonstrated that the dual discharge strategy is effective in meeting water-quality standards/goals in both the river [dissolved oxygen (DO)] and the lake [total phosphorus (TP)] under summer average conditions of river flow and upstream boundary condition DO. Here, that analysis is extended to include a probabilistic treatment of the impact of natural variability in river flow and DO boundary conditions on the feasibility of this management option. Model simulations, incorporating these key sources of system variability, indicate that the dual discharge strategy will meet the lake management goal for TP ～ 94% of the time, with no attendant violation of river DO standards. Excursions from the lake TP goal, occurring ～ 6% of the time, range from 1–5?μg?L?1, are within the range of uncertainty in indicators applied in identifying trophic status. This novel management option is compared with an in-lake discharge alternative in terms of technical and economic feasibility and public acceptance of resultant water quality. Additional management actions, recommended to accompany implementation of the dual discharge strategy, are discussed.     

Review of Design Guidelines for FRP Confinement of Reinforced Concrete Columns of Noncircular Cross Sections

Current international design guidelines provide predictive design equations for the strengthening of reinforced concrete (RC) columns of both circular and prismatic cross sections by means of fiber-reinforced polymer (FRP) confinement and subjected to pure axial loading. Extensive studies (experimental and analytical) have been conducted on columns with circular cross sections, and limited studies have been conducted on members with noncircular cross sections. In fact, the majority of available research work has been on small-scale, plain concrete specimens. In this review paper, four design guidelines are introduced, and a comparative study is presented. This study is based on the increment of concrete compressive strength and ductility and includes the experimental results from six RC columns of different cross-sectional shapes. The observed outcomes are used to identify and remark upon the limits beyond the ones specifically stated by each of the guides and that reflect the absence of effects not considered in current models. The purpose of this study is to present a constructive critical review of the state-of-the-art design methodologies available for the case of FRP-confined concrete RC columns and to indicate a direction for future developments.