Stage–discharge estimation in straight compound channels using isovel contours

For various hydrological applications such as flood control projects, a knowledge of stage–discharge relationship is of particular interest to river engineers. Stage–discharge curves in compound channels cannot be easily predicted in comparison with single channels due to their 3D characteristics of flow. In this paper, the concept of cross‐sectional isovel contours is used for estimation of stage–discharge curves in compound channels. The multivariate Newton's method is applied to the difference between the observed and estimated data to optimize the exponent values of the governing parameters. The accuracy of the proposed model is tested successfully against available experimental results, which are taken from the Flood Channel Facility (FCF) laboratory. Then the results are compared with the Single and Divided Channel Methods (SCM and DCM, respectively), the Weighted Divided Channel Method (WDCM), the exchange discharge method (EDM), and the Coherence Method (COHM). The average values of mean absolute percentage error (MAPE) and normalized root mean square error (NRMSE) in discharge estimation based on each referenced section at any level for 6 sections of the experimental cases are within 3.1% and 0.023, respectively. The biggest advantage of the proposed method is its inherent simplicity, which does not need any calibration.     

黄河径流量的历史演变规律及成因

基于黄河上、中和下游的径流及气候资料,对径流的年代际变化规律及与气候变化的关系进行了分析.结果表明：黄河流域的径流均存在显著的年代际变化趋势,径流的显著特征是从20世纪80年代开始的减少趋势,但并未达到历史的最低,径流减少的趋势在下游比上游更显著,而这种变化趋势与流域的气候变化趋势基本一致,说明在年代际尺度上,径流的变化主要受气候的控制;在不同季节,这种关系有明显差异,如在冬季两者的变化趋势有较大差异.分析还发现,近年来流域地表的干化是流域径流减少的原因,气温的升高更加剧了流域地表干化.     

Response of salinity distribution around the Yellow River mouth to abrupt changes in river discharge

To investigate how salinity changes with abrupt increases and decreases in river discharge, three surveys were conducted along six sections around the Yellow River mouth before, during and after a water regulation event during which the river discharge was increased from ∼200 to >3000 m³ s⁻¹ for the first 3 days, was maintained at >3000 m³ s⁻¹ for the next 9 days and was decreased to <1000 m³ s⁻¹ for the final 4 days. The mean salinity in the Yellow River estuary area during the event varied ∼1.21, which is much larger than its seasonal variation (∼0.50) and interannual variation (∼0.05). Before the event, a small plume was observed near the river mouth. During the event, the plume extended over 24 km offshore in the surface layer in the direction of river water outflow. After the event, the plume diminished in size but remained larger than before the event. The downstream propagation of the plume (as in a Kelvin wave sense) was apparent in the bottom layer during the second survey and in both the surface and bottom layers during the third survey. The plume sizes predicted by the formulas from theoretical studies are larger than those we observed, indicating that factors neglected by theoretical studies such as the temporal variation in river discharge and vertical mixing in the sea could be very important for plume evolution. In addition to the horizontal variation of the plume, we also observed the penetration of freshwater from the surface layer into the bottom layer. A comparison of two vertical processes, wind mixing and tidal mixing, suggests that the impact of wind mixing may be comparable with that of tidal mixing in the area close to the river mouth and may be dominant over offshore areas. The change in Kelvin number indicates an alteration of plume dynamics due to the abrupt change in river discharge during the water regulation event.     

A revised procedure for discharge measurement by means of the salt dilution method

Discharge measurement by means of injection of a NaCl-solution and integration of the electrical conductivity as a function of time is a traditional and well-documented method for use in turbulent streams. The ease of the chemical measurements permits results to be calculated in the field. This work has shown that dry fine-grained salt may substitute the salt solution. Correction factors for water temperature and background conductivity are derived from experiments. A simple procedure for calculation of the discharge is described.     

Estimation of debris flow discharge coefficient considering sediment concentration

A sabo dam has a purpose to block the path of debris flow. However, when overflow occurs, a sabo dam works as a weir, a vertical obstruction, where the fluid must flow over. Many empirical formulas and discharge coefficients for weirs relating flow depth to discharge have been proposed to calculate overflow discharges. However, only a few studies about overflow discharge coefficients are available in the case of debris flow. In this paper, experiments and numerical simulations were done to estimate debris flow discharge coefficients by considering the sediment concentration. In the numerical simulation, a complete overflow equation and a free overfall equation were implemented to calculate debris overflow discharges at a sabo dam. To determine the discharge coefficients for each equation, single factor regression analysis was used. Laboratory experiments were done to calibrate and to compare with the simulation. Study results showed that the discharge coefficients increase as the sediment concentration increases. This finding suggests debris flow discharge coefficients are derived to calculate the debris overflow discharges at a sabo dam.     

A fast method of flood discharge estimation

Discharge, especially during flood periods, is among the most important information necessary for flood control, water resources planning and management. Owing to the high flood velocities, flood discharge usually cannot be measured efficiently by conventional methods, which explains why records of flood discharge are scarce or do not exist for the watersheds in Taiwan. A fast method of flood discharge estimation is presented. The greatest advantage of the proposed method is its application to estimate flood discharge that cannot be measured by conventional methods. It has as its basis the regularity of open‐channel flows, i.e. that nature maintains a constant ratio of mean to maximum velocities at a given channel section by adjusting the velocity distribution and the channel geometry. The maximum velocity at a given section can be determined easily over a single vertical profile, which tends to remain invariant with time and discharge, and can be converted to the mean velocity of the entire cross‐section by multying by the constant ratio. Therefore the mean velocity is a common multiple of maximum velocity and the mean/maximum velocity ratio. The channel cross‐sectional area can be determined from the gauge height, the water depth at the y‐axis or the product of the channel width multiplied by the water depth at the y‐axis. Then the most commonly used method, i.e. the velocity–area method, which determines discharge as the product of the cross‐sectional area multiplied by mean velocity, is applied to estimate the flood discharge. Only a few velocity measurements on the y‐axis are necessary to estimate flood discharge. Moreover the location of the y‐axis will not vary with time and water stage. Once the relationship of mean and maximum velocities is established, the flood estimation can be determined efficiently. This method avoids exposure to hazardous environments and sharply reduces the measurement time and cost. The method can be applied in both high and low flows in rivers. Available laboratory flume and stream‐flow data are used to illustrate accuracy and reliability, and results show that this method can quickly and accurately estimate flood discharges. Copyright © 2004 John Wiley & Sons, Ltd.     

Countermeasure of river bend scour using a combination of submerged vanes and riprap

A series of laboratory flume experiments were done in a large-scale 180° bend with non-cohesive sediment to find optimal or effective protection works at a bend. Detailed study of the scour and flow field dynamics with and without protection works was done. Spatially dense, high frequency velocity data were collected and analyzed to describe the pattern and magnitude of three-dimensional(3 D)velocity throughout the bend. Characterizing the role of flow field dynamics on the pattern of deposition and erosion through experimental measurements provided valuable data about how such flow features contribute to scour and about the performance of the protection works. From the experimental results, it is revealed that for a perennial river it is not possible to protect from scour either with riprap or with submerged vanes alone. Protection from scour at a bend can be achieved with proper combination of these two works. First, submerged vanes can protect the toe, and, second, riprap can protect the upper part of the slope if it is not damaged through toe erosion. The experiments convincingly demonstrate the efficiency of this bank protection technique.     

Increasing discharge from the Mackenzie River system to the Arctic Ocean

The Mackenzie River, Canada's longest and largest river system, provides the greatest Western Hemisphere discharge to the Arctic Ocean. Recent reports of declining flows have prompted concern because (1) this influences Arctic Ocean salinity, stratification and polar ice; (2) a major tributary, the Peace River, has large hydroelectric projects, and further dams are proposed; and (3) the system includes the extensive and biodiverse Peace–Athabasca, Slave and Mackenzie deltas. To assess hydrological trends over the past century that could reflect climate change, we analysed historic patterns of river discharges. We expanded the data series by infilling for short gaps, calculating annual discharges from early summer‐only records (typical r² > 0.9), coordinating data from sequential hydrometric gauges (requiring r² > 0.8) and advancing the data to 2013. For trend detection, Pearson correlation provided similar outcomes to non‐parametric Kendall's τ and Spearman's ρ tests. There was no overall pattern for annual flows of the most southerly Athabasca River (1913–2013), while the adjacent, regulated Peace River displayed increasing flows (1916–2013, p < 0.05). These rivers combine to form the Slave River, which did not display an overall trend (1917–2013). The more northerly, free‐flowing Liard River is the largest tributary and displayed increasing annual flows (1944–2013, p < 0.01, ~3.5% per decade) because of increasing winter, spring, and summer flows, and annual maximum and minimum flows also increased. Following from the tributary contributions, the Mackenzie River flows gradually increased (Fort Simpson 1939–2013, p < 0.05, ~1.5% per decade), but the interannual patterns for the Liard and other rivers were correlated with the Pacific Decadal Oscillation, complicating the pattern. This conclusion of increasing river flows to the Arctic Ocean contrasts with some prior reports, based on shorter time series. The observed flow increase is consistent with increasing discharges of the large Eurasian Arctic drainages, suggesting a common northern response to climate change. Analyses of historic trends are strengthened with lengthening records, and with the Pacific Decadal Oscillation influence, we recommend century‐long records for northern rivers. Copyright © 2016 John Wiley & Sons, Ltd.     

Mesoscalic estimation of nitrogen discharge via drainage systems

A complex approach has been developed for estimating mesoscalic nitrogen discharges via drainage systems using spatial information about land use, drainage areas, nitrogen balances and soil and site characteristics. Determining the total drainage area involves certain difficulties for larger areas, as on the one hand, the available databases are incomplete, and on the other hand the localisation and digitalisation of large subsurface drainage areas is a very time-consuming process. Knowledge of the history and causes of drainage systems in landscapes is required. To solve this problem a method has been developed to calculate the drainage areas for large catchments. In order to obtain a complete data set of subsurface drainage areas, representative areas were selected to enable the proportion of subsurface drainage area to be determined for various soil and site characteristics. These proportions were extrapolated to the entire area and the approach tested in the Mulde River Catchment Area in Germany.The rate of drained arable land is about 25.2% of the total area, which can be broken down into grassland (19.0%) and arable land (27.4%). The results differ for sandy soils with up to 8% drained areas and 57.8% for stagnant soils. This shows that the proportion of drained land is highly dependent on the nature of the soil in the catchment area, which has profound implications for approaches to nitrogen modelling.Average nitrogen discharge for the whole catchment area via drainage water was 33 kg ha⁻¹ yr⁻¹ in the 1980s and 10 kg ha⁻¹ yr⁻¹ in the 1990s. The nitrogen discharge varies from one soil type to another: in regions with sandy substrate (11,900 ha) discharge was 34 kg ha⁻¹ yr⁻¹ in the 1980s (14 kg ha⁻¹ yr⁻¹ in the 1990s), while in areas with loess lessivé soils (89,200 ha) it was about 26 kg ha⁻¹ yr⁻¹ in the 1980s (9 kg ha⁻¹ yr⁻¹ in the 1990s). The reduction can be explained by the complete change in farming strategy since the demise of the former German Democratic Republic (GDR).The approach shown is well suited to future model approaches on a regional scale. By creating and integrating new data sets derived from modern GIS operations the approach reduces the uncertainty of water and nitrogen modelling. This gives us a better understanding of nitrogen discharges into surface and groundwater and temporal discharge dynamics. The discharge data are highly valuable to predict environmental protection measurements for streams, lakes, coastal waters and groundwaters.     

河北省黄壁庄水库坝基渗漏单井示踪试验研究

以河北省黄壁庄水库副坝坝基的渗漏塌坑为例，利用“九五”国家重点科技成果推广项目“智能化地下水动态参数测量仪”在天然流场下，测量出水库自然渗漏量地下水的渗透流速流向、垂直向流速流向，每米含水层的渗透系数等工程必须的水文地质参数。确定性地诊断和检测出黄壁水库坝基渗漏隐患出现的具体位置和高程，从而为大坝的渗漏处理和安全处理提供详细的现场测量数据。     

Efficient methods of discharge measurements in rivers and streams based on the probability concept

Regularities exist in fluid flows and can be represented by a set of constants. These constants are functions of the parameter of a probability distribution that exhibits resilience and stability under various flow conditions. Together, these regularities form a network and interact with each other, such that if one is known then the others can be determined from it. The regularities and their network explain the various fluid‐flow phenomena and can be used in analysis of rivers and streams. For example, they can be used as the basis to develop simple and efficient methods of discharge measurements as presented herein, which only require velocity sampling at a single point on a water surface or a few points on a single vertical. Because of their simplicity and the short time requirement, these methods can be easily automated for collecting discharge data in unsteady, high flows that are badly needed for real‐time flow forecasting and design of flood control structures, and for advancing the fundamental, scientific knowledge in hydrology. Copyright © 2005 John Wiley & Sons, Ltd.     

Synoptic measurements of episodic offshore flow events in the central mid-Atlantic Bight

Sub-inertial across-shelf flows over the central mid-Atlantic Bight were examined using long-range CODAR SeaSonde type HF-radar surface velocities for the period August 2002 through February 2004. Across-shelf flow in this region was found to be episodic, with offshore flow occurring six times more often than onshore flow. Several recurring spatial patterns in offshore flow are identified, and two of these patterns emerge as statistically significant using empirical orthogonal functions. The most common type of offshore flow was a shelf-wide flow with events occurring throughout the year. The most energetic and well-defined shelf-wide flow events occurred during October through April, when the water column was less stratified. Other offshore flow patterns were more localized, with a flow at the bend or ‘point’ in the New Jersey coast being most common and most energetic. The offshore flow events identified here are capable of driving significant offshore transport, with more than one-third of offshore flow episodes having potential transport distances of greater than 14 km. In addition, clusters of multiple offshore flow events are shown to be capable of transporting near-surface material across much of the continental shelf.     

特征差分法在水位推流上的应用

本文阐述了在符合明渠非恒定流一维定床条件的水文测站,应用变化的特征差分格式,将两断面实测水位过程线推算为流量过程线的原理、方法和应用实例,对理论推算的流量过程和实测流量过程进行了比较分析。从泗洪站和运河站的成果看出,推流和实测过程线基本吻合,精度较高。因此可以明显看出,应用本法推流,不但安全可靠,具有节约人力、物力等经济效益,还可提高测验精度、防止高洪期流量的错测、漏测,保证洪水资料的连续和完整。在湖泊水量平衡计算中,同样可以借助两组水尺获得的水位资料进行推流,无需设置专用流量站。     

Indirect estimation of ungauged peak discharges in a bedrock channel with reference to design discharge selection

Geomorphological evidence and recent trash lines were used as stage indicators in a step-backwater computer model of high discharges through an ungauged bedrock channel. The simulation indicated that the peak discharge from the 26.7 m² catchment was close to 150m³s^?1 during the passage of Hurricane Charlie in August 1986. This estimate can be compared with an estimate of 130–160 m³s^?1 obtained using the Flood Studies Report (FSR) unit hydrograph methodology. Other palaeostage marks indicate that higher stages have occurred at an earlier time associated with a discharge of 200 m³s^?1. However, consideration of both the geometry of a plunge pool and transport criteria for bedrock blocks in the channel indicates that floods since 1986 have not exceeded 150 m³s^?1. Given that the estimated probable maximum flood (PMF) calculated from revised FSR procedure is at least 240 m³s^?1, it is concluded that compelling evidence for floods equal to the PMF is lacking. Taking into consideration the uncertainty of the discharge estimation, the 1986 flood computed using field evidence has a minimum return period of 100 years using the FSR procedure. This may be compared with a return period for the same event in the neighbouring gauged River Greta of > 100 years and a rainfall return period of 190 years. In as much as discharges of similar order to FSR estimates are indicated, it is concluded (a) that regional geomorphological evidence and flood simulation within ungauged catchments may be useful as a verification for hydrological estimates of recent widespread flood magnitude and (b) that palaeohydraulic computation can be useful in determining the magnitude of the local maximum [historic] flood when determining design discharges for hydraulic structures within specific catchments.     

一维河网不稳定流模型在黄河下游枯季水量调度中的应用

随着黄河下游工农业生产的不断发展,人民群众生活质量不断提高,黄河水资源供需矛盾越来越突出.确保黄河不断流必然面临着越来越大的压力,必须提高水量调度的快速反映能力.本文讨论了一维河网不稳定流数值计算模型在黄河下游枯水演进中的数学基础、 实现方法和在此基础上建立的枯季水量调度软件的应用.从对八十年代以来的枯水年份的模拟情况来看,流量拟合结果均能达到乙级以上,这对于水资源短缺的黄河下游的引黄调度具有重要的参考价值.     

Analytic stage–discharge formulae for flow in straight trapezoidal open channels

Analytic stage–discharge formulae are derived for flow in straight trapezoidal channels, based on the 2D analytic velocity distribution in open channels given by Shiono and Knight [Shiono K, Knight DW, Turbulent open-channel flows with variable depth across the channel. J Fluid Mech 1991;222:617–46]. A simple hand-calculation method is provided. Legendre incomplete elliptic integrals of the first and second kinds and a binomial series expansion are used in the derivation of these analytic formulae, together with physically based hydraulic parameters, such as local friction factor (f), dimensionless eddy viscosity (λ) and secondary flow (Γ). The stage–discharge results obtained from the formulae are shown to be in good agreement with experimental data, as are the corresponding analytic velocity and boundary shear stress distributions. The influences of f, λ and Γ on the stage–discharge relationship are also discussed.     

The impact of single natural vegetation elements on flow characteristics

In this experimental study, measurements were conducted to explore the impacts of different forms of individual natural vegetative elements within the flow domain on velocity and turbulence characteristics. All the experiments were performed in a flume measuring 26 m in length, 0·98 m in width and 0·85 m in depth, and real tree saplings were utilized to represent the vegetation element. In order to analyse this commonly observed nature phenomenon in floodplains, trees with wide trunks were classified into three groups on the basis of their volume versus height relation. Throughout the velocity measurements three acoustic Doppler velocimeters were employed. Time‐averaged velocity, streamwise and vertical turbulence intensities and turbulent kinetic energy parameters were examined. Additionally, a formulation that gives the velocity profile at a certain distance downstream of vegetation was introduced and the validity of the proposed formulation was checked with experimental data. It is seen that despite their porous structures, the presence of vegetation considerably disturbs the flow field and dissipates a remarkable amount of energy by turbulence. Copyright © 2008 John Wiley & Sons, Ltd.     

Time-mean structure of secondary flows in open channel with longitudinal bedforms

Secondary motions are commonly present in open channel flows. This study aims to investigate, both experimentally and analytically, time-mean characteristics of cellular secondary flows generated by longitudinal bedforms. Experiments were conducted in a tilting, rectangular flume with six different longitudinal bedforms, including alternate bed strips with different roughness heights and bed ridges of wavy and rectangular shapes. Various flows were sampled using a two-dimensional laser Doppler anemometer (LDA) and a one-dimensional ultrasonic Doppler velocimeter (UDV). Experimental results demonstrate secondary flows appearing basically in cellular fashion over the modeled longitudinal bedforms. It is also shown that the cellular structures can be described analytically with kinematic considerations. The discrepancies between theoretical and measurement results are discussed. An empirical relationship between maximum vertical velocity and bed configuration is finally proposed based on the experimental data.     

Suspended load in flows on erodible bed

Steady state suspended-load of sediment transported in flow over erodible beds usually is treated by the advection-diffusion approach, though in recent years, it is being treated as a two-phase flow phenomenon incorporating kinetics of sediment particles. Among the advection-diffusion approaches, Rouse＇s equation is the well-known, although a number of researchers in later periods have attempted to improve it by modifying the mixing length concept taking into account other aspects. In this paper, the advection-diffusion approach and associated logarithmic law of flow velocity are revisited. It is concluded from the logarithmic law that the Reynolds shear stress is a linear function of height above the bed, which reduces to bed shear stress in the case of a long horizontal channel. As a consequence, it is shown that the volumetric concentration of sediment is best approximated by the sum of two power laws of height above the bed. An equation is derived for the suspended-load transport rate in terms of elementary functions.     

Impact of wastewater discharge on the channel morphology of ephemeral streams

The impact of wastewater flow on the channel bed morphology was evaluated in four ephemeral streams in Israel and the Palestinian Territories: Nahal Og, Nahal Kidron, Nahal Qeult and Nahal Hebron. Channel changes before, during and after the halting of wastewater flow were monitored. The wastewater flow causes a shift from a dry ephemeral channel with intermittent floods to a continuous flow pattern similar to that of humid areas. Within a few months, nutrient‐rich wastewater flow leads to rapid development of vegetation along channel and bars. The colonization of part of the active channel by vegetation increases flow resistance as well as bank and bed stability, and limits sediment availability from bars and other sediment stores along the channels. In some cases the established vegetation covers the entire channel width and halts the transport of bed material along the channel. During low and medium size flood events, bars remain stable and the vegetation intact. Extreme events destroy the vegetation and activate the bars. The wastewater flow results in the development of new small bars, which are usually destroyed by flood flows. Due to the vegetation establishment, the active channel width decreases by up to 700 per cent. The deposition of fine sediment and organic material changed the sediment texture within the stable bar surface and the whole bed surface texture in Nahal Hebron. The recovery of Nahal Og after the halting of the wastewater flow was relatively fast; within two flood seasons the channel almost returned to pre‐wastewater characteristics. The results of the study could be used to indicate what would happen if wastewater flows were introduced along natural desert streams. Also, the results could be used to predict the consequences of vegetation removal as a result of human intervention within the active channel of humid streams. Copyright © 2001 John Wiley & Sons, Ltd.