On the methods for determining the transverse dispersion coefficient in river mixing

In this study, the strengths and weaknesses of existing methods for determining the dispersion coefficient in the two-dimensional river mixing model were assessed based on hydraulic and tracer data sets acquired from experiments conducted on either laboratory channels or natural rivers. From the results of this study, it can be concluded that, when the longitudinal dispersion coefficient as well as the transverse dispersion coefficients must be determined in the transient concentration situation, the two-dimensional routing procedures, 2D RP and 2D STRP, can be employed to calculate dispersion coefficients among the observation methods. For the steady concentration situation, the STRP can be applied to calculate the transverse dispersion coefficient. When the tracer data are not available, either theoretical or empirical equations by the estimation method can be used to calculate the dispersion coefficient using the geometric and hydraulic data sets. Application of the theoretical and empirical equations to the laboratory channel showed that equations by Baek and Seo [[3], 2011] predicted reasonable values while equations by Fischer [23] and Boxwall and Guymer (2003) overestimated by factors of ten to one hundred. Among existing empirical equations, those by Jeon et al. [28] and Baek and Seo [6] gave the agreeable values of the transverse dispersion coefficient for most cases of natural rivers. Further, the theoretical equation by Baek and Seo [5] has the potential to be broadly applied to both laboratory and natural channels.     

A global river routing network for use in hydrological modeling

In this paper a relatively simple procedure is presented to construct a global river routing network on a 0.5° latitude–longitude grid. In this network all grid cells in a catchment are coupled and have a flow direction, making it a useful tool in the modeling of river flow on a global scale. The flow directions are based on a digital elevation model and on information on the locations of major rivers (‘stream burning’). The presented river routing network is specifically designed for the assessment of fresh water shortages. We tested the validity of the river routing network by comparing the computed drainage areas with published estimates. This comparison revealed a good similarity and it is concluded that the presented river routing network has sufficient quality to be implemented in global climate models. This could mean a considerable improvement of the surface parameterization in these models.     

Assessment of river flow with significant lateral inflow through reverse routing modeling

The discharge hydrograph estimation in rivers based on reverse routing modeling and using only water level data at two gauged sections is here extended to the most general case of significant lateral flow contribution, without needing to deploy rainfall–runoff procedures. The proposed methodology solves the Saint‐Venant equations in diffusive form also involving the lateral contribution using a “head‐driven” modeling approach where lateral inflow is assumed to be function of the water level at the tributary junction. The procedure allows to assess the discharge hydrograph at ends of a selected river reach with significant lateral inflow, starting from the stage recorded there and without needing rainfall data. Specifically, the MAST 1D hydraulic model is applied to solve the diffusive wave equation using the observed stage hydrograph at the upstream section as upstream boundary condition. The other required data are (a) the observed stage hydrograph at the downstream section, as benchmark for the parameter calibration, and (b) the bathymetry of the river reach, from the upstream section to a short distance after the downstream gauged section. The method is validated with different flood events observed in two river reaches with a significant intermediate basin, where reliable rating curves were available, selected along the Tiber River, in central Italy, and the Alzette River, in Luxembourg. Very good performance indices are found for the computed discharge hydrographs at both the channel ends and along the tributaries. The mean Nash‐Sutcliffe value (NS_q) at the channel ends of two rivers is found equal to 0.99 and 0.86 for the upstream and downstream sites, respectively. The procedure is also validated on a longer stretch of the Tiber River including three tributaries for which appreciable results are obtained in terms of NS_q for the computed discharge hydrographs at both the channel ends for three investigated flood events.     

Estimating longitudinal dispersion in rivers using Acoustic Doppler Current Profilers

The longitudinal dispersion coefficient (D) is an important parameter needed to describe the transport of solutes in rivers and streams. The dispersion coefficient is generally estimated from tracer studies but the method can be expensive and time consuming, especially for large rivers. A number of empirical relations are available to estimate the dispersion coefficient; however, these relations are known to produce estimates within an order of magnitude of the tracer value. The focus of this paper is on using the shear-flow dispersion theory to directly estimate the dispersion coefficient from velocity measurements obtained using an Acoustic Doppler Current Profiler (ADCP). Using tracer and hydrodynamic data collected within the same river reaches, we examined conditions under which the ADCP and tracer methods produced similar results. Since dead zones / transient storage (TS) are known to influence the dispersion coefficient, we assessed the relative importance of dead zones in different stream reaches using two tracer-based approaches: (1) TS modeling which explicitly accounts for dead zones and (2) the advection–dispersion equation (ADE) which does not have separate terms for dead zones. Dispersion coefficients based on the ADE tend to be relatively high as they describe some of the effects of dead zones as well. Results based on the ADCP method were found to be in good agreement with the ADE estimates indicating that storage zones play an important role in the estimated dispersion coefficients, especially at high flows. For the river sites examined in this paper, the tracer estimates of dispersion were close to the median values of the ADCP estimates obtained from multiple datasets within a reach. The ADCP method appears to be an excellent alternative to the traditional tracer-based method if care is taken to avoid spurious data and multiple datasets are used to compute a distance-weighted average or other appropriate measure that represents reach-averaged conditions.     

基于连续方程的河道水流双变量耦合演算模型

针对现有的河道水流洪水演算模型只能模拟单一变量(流量或水位)的问题,以水流连续方程和河段蓄水量的两种不同表达形式(蓄水量等于平均过水断面面积与河段长乘积,蓄水量等于河段平均流量与传播时间的乘积)为基础,对马斯京根模型进行了通用性改进,提出了双变量耦合通用演算模型.选取了四大水系(包括内陆河流和入海河流)的16个河段汛期洪水资料进行模型检验,模型验证考虑了地理范围、不同的河段特征和水力特征、洪水量级等因素,全面地检验了模型结构的合理性和模拟实际洪水的有效性.将双变量耦合通用演算模型与传统的马斯京根法进行了效果比较,结果表明双变量耦合通用演算模型的模拟精度高于马斯京根法,模拟效果比马斯京根法稳定一些,而且具有较好的通用性.     

Analytical solutions for transport of decaying solutes in rivers with transient storage

Analytical solutions are presented for solute transport in rivers including the effects of transient storage and first order decay. The solute transport model considers an advection–dispersion equation for transport in the main channel linked to a first order mass exchange between the main channel and the transient storage zones. In case of a conservative tracer, it is shown that different analytical solutions presented in the literature are mathematically identical. For non-conservative solutes, first order decay reactions are considered with different reaction rate coefficients in the main river channel and in the dead zones. New analytical solutions are presented for different boundary conditions, i.e. instantaneous injection in an infinite river reach, and variable concentration time series input in a semi-infinite river reach. The correctness and accuracy of the analytical solutions is verified by comparison with the OTIS numerical model. The results of analytical and numerical approaches compare favourably and small differences can be attributed to the influence of boundary conditions. It is concluded that the presented analytical solutions for solute transport in rivers with transient storage and solute decay are accurate and correct, and can be usefully applied for analyses of tracer experiments and transport characteristics in rivers with mass exchange in dead zones.     

A multilinear discrete Nash-cascade model for stage-hydrograph routing in compound river channels

ABSTRACT

Stream gauge-based information is the foundation for many hydrological applications in a river basin including the aquatic-habitat conservation. A simple two-parameter model for routing streamflow depth (alternatively, stream–stage) hydrographs and estimating corresponding discharge hydrographs in river channels is proposed using the multilinear approach, based on Nash-type discrete-cascade model. The applicability of this model is investigated by extending its framework to the realm of compound cross-section trapezoidal channels for both in-bank and overbank flows by using 20 flood events of the Tiber River in the Umbria region of Central Italy, and subsequently comparing the simulated results with the corresponding simulations of the HEC-RAS (Hydrologic Engineering Center – River Analysis System) hydrodynamic model and observed flow depth hydrographs. The field application, comparative study, and uncertainty and sensitivity analysis of the results demonstrate that the proposed multilinear discrete Nash-cascade stage-hydrograph (MDNS) routing model has the potential for routing floods in real-world rivers and canal irrigation systems, especially in operational mode.     

Predicting longitudinal dispersion coefficient in natural streams by artificial intelligence methods

In this study, three artificial neural network methods, i.e. feed forward back propagation, the radial basis function neural network, and the generalized regression neural network are employed to compute the longitudinal dispersion coefficient in order to evaluate its behaviour in predicting dispersion characteristics in natural streams. These methods, which use hydraulic and geometrical data to predict dispersion coefficients, can easily be applied to natural streams and are proven to be superior in explaining their dispersion characteristics more precisely than existing equations. This method of predicting the longitudinal dispersion coefficient in river flows was tested on 65 data sets, obtained by researchers from 30 rivers in the USA. Results using the models are compared with results obtained in many other studies, and are shown to be more accurate than the other methods considered. Copyright © 2008 John Wiley & Sons, Ltd.     

Algorithms for solving the diffusive wave flood routing equation

Generally, the diffusive wave equation, obtained by neglecting the acceleration terms in the Saint-Venant equations, is used in flood routing in rivers. Methods based on the finite-difference discretization techniques are often used to calculate discharges at each time step. A modified form of the diffusive wave equation has been developed and new resolution algorithms proposed which are better adapted to flood routing along a complex river network. The two parameters of the equation, celerity and diffusivity, can then be taken as functions of the discharge. The resolution algorithm allows the use of any distribution of lateral inflow in space and time. The accuracy of the new algorithms were compared with a traditional algorithm by numerical experimentation. Special attention was given to the instability caused by the inflow signal which constitutes the upstream boundary condition. For the fully diffusive wave flood routing problem, all three algorithms tested gave good results. The results also indicate that the efficiency of the new algorithms could be significantly improved if the position of the x-axis is modified by rotation. The new algorithms were applied to flood routing simulation over the Gardon d'Anduze catchment (542 km²) in southern France.     

太湖北岸湖滨带观测场水生植物群落特征及其影响因素分析

研究了鄱阳湖流域在1955-2002年间的径流系数的变化,重点分析了它与水循环的两个基本要素:降水量和蒸发量的关系,同时对其原因进行了初步的探讨.经分析,在鄱阳湖流域中,径流系数较大的是饶河流域和信江流域,较小的是抚河流域;在年内变化上,4-6月为五河流域径流系数比较大的月份,这与鄱阳湖流域降水集中期相对应.在空间上,4-6月仍然以饶河流域和信江流域相对较大,而抚河流域较小,特别是8月份的径流系数远小于其他四河;年代际变化上,1990s径流系数增加较为显著.尽管鄱阳湖流域的径流系数除了受气候因子的影响外,还受到水土流失和地形等因素的影响,但是降水量的增加,特别是暴雨频率的增加仍然是其主要影响因素,蒸发量的减小对径流系数的增加也有一定程度的影响.径流系数与气温并无明显的线性相关关系.     

Relationships between three dispersion measures used in flood frequency analysis

Three dispersion measures of a random variable, i.e., the standard deviation, the mean deviation (MD) about the mean and the second L-moment, are analyzed in terms of their properties and mutual relationships. Emphasis is placed on the MD, as it is less recognized than two other dispersion measures. The relationships between the dispersion measures are derived for distributions commonly applied in flood frequency analysis (FFA). For distributions that are unbounded, there is a distribution-dependent constant value of the ratio of dispersion measures, or equivalently of respective coefficients of variation. For two-parameter distributions that are lower-bounded, the relationship between the coefficients of variation is also distribution dependent and is not linear. For lower-bounded three-parameter distributions, the dispersion measure ratios, or equivalently the ratios of coefficients of variation, depend on the coefficient of skewness and show a strong distributional dependence. For selected distributions, the three dispersion measures are compared both in terms of the robustness to the largest samples element and the accuracy of upper quantile estimation. The MD statistics may be highly competitive to the two other dispersion measure statistics if applied in FFA for parameters estimation.     

Variation of turbidity during subsurface abstraction of river water: A case study

Rivers are known to carry sediment loads. Many rivers in Asia including Kosi, Brahmaputra, Bagmati rivers in India and Yellow and Yarlung Tsangpo rivers in China are well known for their carriage of excessive sediment loads. In India, many water abstraction units are located along river banks. Their success depends on how the sediment load is trapped as the sediment laden water flows towards the pumping wells located few hundred meters away from river bank. In this work, based on the analysis of a data of one year on suspended load variations in a canal, located at Haridwar, India, trap efficiency variations are studied during the passage of source water to the abstraction well. Considering the subsurface basin, located within the source water and the abstraction well, as a hypothetical settling basin, variation of trap efficiency is related with the influent concentration. It is observed that the trap efficiency depends on the influent concentration. Appropriate functional forms for the variation of trap efficiency with the influent concentrations are also identified.     

Reduced‐complexity flow routing models for sinuous single‐thread channels: intercomparison with a physically‐based shallow‐water equation model

Reduced‐complexity models of fluvial processes use simple rules that neglect much of the underlying governing physics. This approach is justified by the potential to use these models to investigate long‐term and/or fundamental river behaviour. However, little attention has been given to the validity or realism of reduced‐complexity process parameterizations, despite the fact that the assumptions inherent in these approaches may limit the potential for elucidating the behaviour of natural rivers. This study presents two new reduced‐complexity flow routing schemes developed specifically for application in single‐thread rivers. Output from both schemes is compared with that from a more sophisticated model that solves the depth‐averaged shallow water equations. This comparison provides the first demonstration of the potential for deriving realistic predictions of in‐channel flow depth, unit discharge, energy slope and unit stream power using simple flow routing schemes. It also highlights the inadequacy of modelling unit stream power, shear stress or sediment transport capacity as a function of local bed slope, as has been common practice in a number of previous reduced‐complexity models. Copyright © 2009 John Wiley & Sons, Ltd.     

Sediment size variation in a braided reach of the Sunwapta River,Alberta, Canada

The downstream diminution in sediment size in a braided reach of the proglacial Sunwapta River, Alberta, Canada, was examined statistically to identify the sources of the observed variation about an expected exponential relationship between clast size and distance. Major deviations from this hypothetical relationship, such as a relative increase in grainsize, may be attributed to the effects of tributary sediment inputs and downstream changes in channel behaviour, whilst local variation is associated with complex patterns of sediment deposition observed at a bar scale. A comparison of diminution coefficients, calculated for separate lithologies and for subreaches along the river, with those obtained from previous studies, is used as an indicator of river behaviour and sediment transport processes. It is shown that rates of diminution vary within the reach in response to differing rates of aggradation and to the backwater effects created by tributary alluvial fans. The relatively high values for the calculated diminution coefficients indicate that processes of differential transport are the main cause of the grain size decrease.     

Experimental analysis of reservoir release wave routing in upland boulder bed rivers

Reservoir release wave routing during 33 controlled reservoir releases, along 15 upland boulder bed river channel reaches, on five different regulated rivers were monitored to assess the importance of river channel roughness and reservoir release magnitude on reservoir release wave speeds. Wave speeds varied between 0.52 and 3.01 m s^?1. Reservoir release wave translation, steepening, and attenuation occurred. With high channel roughness values reservoir release wave arrival speed is retarded in comparison to peak stage and wave steepening occurs, but with a reduction in channel roughness reservoir release wave front arrival is accelerated producing attentuation. The threshold between reservoir release wave front attenuation and steepening occurs at a pre-release discharge/channel width of approximately 0.1, an index of channel roughness. The paper also demonstrates, via comparison of observed and calculated reservoir release wave speeds on the River Washburn, Yorkshire, the difficulty of accurately predicting flood wave movement in upland boulder bed channels using existing prediction equations. The calculated values, however, revealed systematic error with pre-release discharge and reservoir release magnitude. Apparently the equations fail to account for the effects of high channel roughness together with pressure gradient forces, induced by rapid rates of stage change on the rising limb of reservoir releases. In order to accurately predict reservoir release wave movement in regulated rivers, this paper demonstrates that hydraulic studies need to be undertaken and pre-release discharges prescribed to determine desired reservoir release wave routing behaviour. Manipulation of the reservoir release pattern at the dam alone, cannot dictate reservoir release wave front form downstream or wave speed.     

Influence of frequency-dependent anisotropy on seismic amplitude-versus-offset signatures for fractured poroelastic rocks

Frequency-dependent amplitude variation with offset offers an effective method for hydrocarbon detections and analysis of fluid flow during production of oil and natural gas within a fractured reservoir. An appropriate representation for the frequency dependency of seismic amplitude variation with offset signatures should incorporate influences of dispersive and attenuating properties of a reservoir and the layered structure for either isotropic or anisotropic dispersion analysis. In this study, we use an equivalent medium permeated with aligned fractures that simulates frequency-dependent anisotropy, which is sensitive to the filled fluid of fractures. The model, where pores and fractures are filled with two different fluids, considers velocity dispersion and attenuation due to mesoscopic wave-induced fluid flow. We have introduced an improved scheme seamlessly linking rock physics modelling and calculations for frequency-dependent reflection coefficients based on the propagator matrix technique. The modelling scheme is performed in the frequency-slowness domain and can properly incorporate effects of both bedded structure of the reservoir and velocity dispersion quantified with frequency-dependent stiffness. Therefore, for a dispersive and attenuated layered model, seismic signatures represent a combined contribution of impedance contrast, layer thickness, anisotropic dispersion of the fractured media and tuning and interference of thin layers, which has been avoided by current conventional methods. Frequency-dependent amplitude variation with offset responses was studied via considering the influences of fracture fills, layer thicknesses and fracture weaknesses for three classes amplitude variation with offset reservoirs. Modelling results show the applicability of the introduced procedure for interpretations of frequency-dependent seismic anomalies associated with both layered structure and velocity dispersion of an equivalent anisotropic medium. The implications indicate that anisotropic velocity dispersion should be incorporated accurately to obtain enhanced amplitude variation with offset interpretations. The presented frequency-dependent amplitude variation with offset modelling procedure offers a useful tool for fracture fluid detections in an anisotropic dispersive reservoir with layered structures.     

Impulse response of a linear diffusion analogy model as a flood frequency probability density function

Abstract

The impulse response of a linear convective-diffusion analogy (LD) model used for flow routing in open channels is proposed as a probability distribution for flood frequency analysis. The flood frequency model has two parameters, which are derived using the methods of moments and maximum likelihood. Also derived are errors in quantiles for these parameter estimation methods. The distribution shows that the two methods are equivalent in terms of producing mean values—the important property in case of unknown true distribution function. The flood frequency model is tested using annual peak discharges for the gauging sections of 39 Polish rivers where the average value of the ratio of the coefficient of skewness to the coefficient of variation equals about 2.52, a value closer to the ratio of the LD model than to the gamma or the lognormal model. The likelihood ratio indicates the preference of the LD over the lognormal for 27 out of 39 cases. It is found that the proposed flood frequency model represents flood frequency characteristics well (measured by the moment ratio) when the LD flood routing model is likely to be the best of all linear flow routing models.     

Scaling aspects of river flow routing

Scaling aspects of river flow routing are studied by comparing two flow routing schemes, one designed for use in coupled general circulation models (GCMs) and operated at large spatial scales (～350 km), and the other designed for use in typical hydrological applications at small spatial scales (～25 km). The same runoff data are used as input into the two routing schemes, and comparisons are made between mean annual, mean monthly and daily streamflow simulated at four locations within the Mackenzie River Basin. The results suggest that for the purpose of realistically modelling monthly streamflow at the mouth of the rivers in GCMs, flow routing at large spatial scales gives similar results. However, the amplitude of the annual streamflow cycle is slightly but characteristically larger, when routing is performed at large spatial scales. Flow routing at large spatial scales also results in overestimation of high flows, while low flows are underestimated. Copyright © 2001 John Wiley & Sons, Ltd.     

Spatial analysis of stream power using GIS: SLk anomaly maps

The stream length‐gradient index (SL) shows the variation in stream power along river reaches. This index is very sensitive to changes in channel slope, thus allowing the evaluation of recent tectonic activity and/or rock resistance. Nevertheless, the comparison of SL values from rivers of different length is biased due to the manner in which the index is formulated, thus making correlations of SL anomalies along different rivers difficult. Therefore, when undertaking a comparison of SL values of rivers of different lengths, a normalization factor must be used. The graded river gradient (K) has already been used in some studies to normalize the SL index. In this work, we explore the relationships between the graded river gradient (K), the SL index and the stream power, proposing the use of a re‐named SLk index, which enables the comparison of variable‐length rivers, as well as the drawing of SLk anomaly maps. We present here a GIS‐based procedure to generate SLk maps and to identify SLk anomalies. In order to verify the advantages of this methodology, we compared an SLk map of the NE border of the Granada basin with both simple river profile–knickpoint identification and with an SL map. The results show that the SLk map supplies good results with defined anomalies and suitably reflects the main tectonic and lithological features of the study area. Copyright © 2008 John Wiley & Sons, Ltd.     

Alkalinity and dissolved inorganic carbon exports from tropical and subtropical river catchments discharging to the Great Barrier Reef,Australia

Dissolved inorganic carbon (DIC) transport by rivers is an important control on the pH and carbonate chemistry of the coastal ocean. Here, we combine DIC and total alkalinity (TAlk) concentrations from four tropical rivers of the Great Barrier Reef region in Australia with daily river discharge to quantify annual river loads and export rates. DIC in the four rivers ranged from 284 to 2,639 μmol kg⁻¹ and TAlk ranged from 220 to 2,612 μmol kg⁻¹. DIC:TAlk ratios were mostly greater than one suggesting elevated exports of free [CO₂*]. This was pronounced in the Johnstone and Herbert rivers of the tropical wet north. The largest annual loads were transported in the two large river catchments of the southern Great Barrier Reef region, the Fitzroy and Burdekin rivers. The carbon stable isotopic composition of DIC suggests that carbonate weathering was the dominant source of DIC in the southern rivers, and silicate weathering was likely a source of DIC in the northern Wet Tropics rivers. Annual loads and export rates were strongly driven by precipitation and discharge patterns, the occurrence of tropical cyclones, and associated flooding events, as well as distinct seasonal dry and wet periods. As such, short-lived hydrological events and long-term (seasonal and inter-annual) variation of DIC and TAlk that are pronounced in rivers of the tropical and subtropical wet and dry climate zone should be accounted for when assessing inorganic carbon loads to the coastal ocean and the potential to buffer against or accelerate ocean acidification.