A worldwide correlation for exponential bed particle size variation in subaerial aqueous flows

The particle size of the bed sediments in or on many natural streams, alluvial fans, laboratory flumes, irrigation canals and mine waste deltas varies exponentially with distance along the stream. A plot of the available worldwide exponential bed particle size diminution coefficient data against stream length is presented which shows that all the data lie within a single narrow band extending over virtually the whole range of stream lengths and bed sediment particle sizes found on Earth. This correlation applies to both natural and artificial flows with both sand and gravel beds, irrespective of either the solids concentration or whether normal or reverse sorting occurs. This strongly suggests that there are common mechanisms underlying the exponential diminution of bed particles in subaerial aqueous flows of all kinds. Thus existing models of sorting and abrasion applicable to some such flows may be applicable to others. A comparison of exponential laboratory abrasion and field diminution coefficients suggests that abrasion is unlikely to be significant in gravel and sand bed streams shorter than about 10 km to 100 km, and about 500 km, respectively. Copyright © 1999 John Wiley & Sons, Ltd.     

Exponential Longitudinal Profiles of Streams

A theoretical model is presented which shows that streams with low solids concentrations and low lateral inflows, whose bed load sediments undergo either comminution or hydraulic sorting under steady or quasi-steady conditions, have exponential profiles. Similar streams in which both communition and sorting are significant have exponential profiles only if they are short Sediment threshold and flow depth estimates based on the theoretical model are consistent with field and laboratory data from the literature. A comparison of the model and communition and sorting data from the literature strongly suggests that hydraulic sorting and comminution dominate in short and long natural streams, respectively. No examples of natural streams of intermediate length with exponential longitudinal profiles were found, suggesting that neither sorting nor comminution is dominant in such streams. © 1997 by John Wiley & Sons, Ltd.     

Patterns in stream longitudinal profiles and implications for hyporheic exchange flow at the H.J. Andrews Experimental Forest,Oregon, USA

There is a need to identify measurable characteristics of stream channel morphology that vary predictably throughout stream networks and that influence patterns of hyporheic exchange flow in mountain streams. In this paper we characterize stream longitudinal profiles according to channel unit spacing and the concavity of the water surface profile. We demonstrate that: (1) the spacing between zones of upwelling and downwelling in the beds of mountain streams is closely related to channel unit spacing; (2) the magnitude of the vertical hydraulic gradients (VHGs) driving hyporheic exchange flow increase with increasing water surface concavity, measured at specific points along the longitudinal profile; (3) channel unit spacing and water surface concavity are useful metrics for predicting how patterns in hyporheic exchange vary amongst headwater and mid‐order streams. We use regression models to describe changes in channel unit spacing and concavity in longitudinal profiles for 12 randomly selected stream reaches spanning 62 km² in the H.J. Andrews Experimental Forest in Oregon. Channel unit spacing increased significantly, whereas average water surface concavity (AWSC) decreased significantly with increasing basin area. Piezometer transects installed longitudinally in a subset of stream reaches were used to measure VHG in the hyporheic zone, and to determine the location of upwelling and downwelling zones. Predictions for median pool length and median distance between steps in piezometer reaches bracketed the median distance separating zones of upwelling in the stream bed. VHG in individual piezometers increased with increasing water surface concavity at individual points in the longitudinal profile along piezometer transects. Absolute values of VHG, averaged throughout piezometer transects, increased with increasing AWSC, indicating increased potential for hyporheic exchange flow. These findings suggest that average hyporheic flow path lengths increase—and the potential for hyporheic exchange flow in stream reaches decreases—along the continuum from headwater to mid‐order mountain streams. Copyright © 2005 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.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

An improved time series model for monthly stream flows

This study aims to develop an improved time series model to overcome difficulties in modeling monthly short term stream flows. The periodic, serial dependent and independent components of the classical time series models are improved separately by information transfer from a surrounding long term gauging station to the considered flow section having short term records. Eventually, an improved model preserving the mathematical model structure of the classical time series model, while improving general and monthly statistics of the monthly stream flows, is derived by using the improved components instead of the short term model components in the time series modeling. The correlative relationships between the current short term and surrounding long term stations are used to improve periodic and serial dependent behaviors of monthly flows. Independent components (residuals) are improved via the parameters defining their theoretical probability distribution. The improved model approach is tested by using 50 year records of Göksu-Himmetli (1801) and Göksu-Gökdere (1805) flow monitoring stations located on the Ceyhan river basin, in south of Turkey. After 50 year records of the station 1801 are separated into five 10 year sub series, their improved and classical time series models are computed and compared with the real long-term (50 year) time series model of this station to reveal efficiencies of the improved models for each subseries (sub terms with 10 year observation). The comparisons are realized based on the model components, model estimates and general/monthly statistics of model estimates. Finally, some evaluations are made on the results compared to the regression method classically applied in the literature.     

Power-law index for velocity profiles in open channel flows

Turbulence theory has demonstrated that the log law is one of the established theoretical results for describing velocity profiles, which is in principle applicable for the near-bed overlap region, being less than about 20% of the flow depth. In comparison, the power law that is often presented in an empirical fashion could apply to larger fraction of the flow domain. However, limited information is available for evaluating the power-law exponent or index. This paper attempts to show that the power law can be derived as a first-order approximation to the log law, and its power-law index is computed as a function of the Reynolds number as well as the relative roughness height. The result obtained also coincides with the fact that the one-sixth power included in the Manning equation is of prevalent acceptance, while higher indexes would be required for flows over very rough boundaries.     

A depth‐averaged two‐phase model for debris flows over erodible beds

Mass exchange between debris flow and the bed plays a vital role in debris flow dynamics. Here a depth‐averaged two‐phase model is proposed for debris flows over erodible beds. Compared to previous depth‐averaged two‐phase models, the present model features a physical step forward by explicitly incorporating the mass exchange between the flow and the bed. A widely used closure model in fluvial hydraulics is employed to estimate the mass exchange between the debris flow and the bed, and an existing relationship for bed entrainment rate is introduced for comparison. Also, two distinct closure models for the bed shear stresses are evaluated. One uses the Coulomb friction law and Manning's equation to determine the solid and fluid resistances respectively, while the other employs an analytically derived formula for the solid phase and the mixing length approach for the fluid phase. A well‐balanced numerical algorithm is applied to solve the governing equations of the model. The present model is first shown to reproduce average sediment concentrations in steady and uniform debris flows over saturated bed as compared to an existing formula underpinned by experimental datasets. Then, it is demonstrated to perform rather well as compared to the full set of USGS large‐scale experimental debris flows over erodible beds, in producing debris flow depth, front location and bed deformation. The effects of initial conditions on debris flow mass and momentum gain are resolved by the present model, which explicitly demonstrates the roles of the wetness, porosity and volume of bed sediments in affecting the flow. By virtue of extended modeling cases, the present model produces debris flow efficiency that, as revealed by existing observations and empirical relations, increases with initial volume, which is enhanced by mass gain from the bed. Copyright © 2017 John Wiley & Sons, Ltd.     

A bed‐load transport model for rough turbulent open‐channel flows on plane beds

Data from flume studies are used to develop a model for predicting bed‐load transport rates in rough turbulent two‐dimensional open‐channel flows moving well sorted non‐cohesive sediments over plane mobile beds. The object is not to predict transport rates in natural channel flows but rather to provide a standard against which measured bed‐load transport rates influenced by factors such as bed forms, bed armouring, or limited sediment availability may be compared in order to assess the impact of these factors on bed‐load transport rates. The model is based on a revised version of Bagnold's basic energy equation i_bs_b = e_bω, where i_b is the immersed bed‐load transport rate, ω is flow power per unit area, e_b is the efficiency coefficient, and s_b is the stress coefficient defined as the ratio of the tangential bed shear stress caused by grain collisions and fluid drag to the immersed weight of the bed load. Expressions are developed for s_b and e_b in terms of G, a normalized measure of sediment transport stage, and these expressions are substituted into the revised energy equation to obtain the bed‐load transport equation i_b = ω G ^3·4. This equation applies regardless of the mode of bed‐load transport (i.e. saltation or sheet flow) and reduces to i_b = ω where G approaches 1 in the sheet‐flow regime. That i_b = ω does not mean that all the available power is dissipated in transporting the bed load. Rather, it reflects the fact that i_b is a transport rate that must be multiplied by s_b to become a work rate before it can be compared with ω. It follows that the proportion of ω that is dissipated in the transport of bed load is i_bs_b/ω, which is approximately 0·6 when i_b = ω. It is suggested that this remarkably high transport efficiency is achieved in sheet flow (1) because the ratio of grain‐to‐grain to grain‐to‐bed collisions increases with bed shear stress, and (2) because on average much more momentum is lost in a grain‐to‐bed collision than in a grain‐to‐grain one. Copyright © 2006 John Wiley & Sons, Ltd.     

Models for recession flows in the upper Blue Nile River

Stream‐flow recessions are commonly characterized by the exponential equation or in the alternative power form equation of a single linear reservoir. The most common measure of recession is the recession constant K, which relates to the power function form of the recession equation for a linear reservoir. However, in reality it can be seen that the groundwater dynamics of even the simplest of aquifers may behave in a non‐linear fashion. In this study three different storage–outflow algorithms; single linear, non‐linear and multiple linear reservoir were considered to model the stream‐flow recession of the upper Blue Nile. The recession parameters for the linear and non‐linear models were derived by the use of least‐squares regression procedures. Whereas, for the multiple linear reservoir model, a second‐order autoregressive AR (2) model was applied first in order to determine the parameters by the least‐squares method. The modelling of the upper Blue Nile recession flow performed shortly after the wet season, when interflow and bank storage may be contributing considerably to the river flow, showed that the non‐linear reservoir model simulates well with the observed counterparts. The variation related to preceding flow on a recession parameter of the non‐linear reservoir remains significant, which was obtained by stratification of the recession curves. Although a similar stratification did not show any systematic variation on the recession parameters for the linear and multiple linear reservoir models. Copyright © 2004 John Wiley & Sons, Ltd.     

Scenario analysis for assessing the impact of hydraulic fracturing on stream low flows using the SWAT model

Scientists and water users are concerned about the potential impact on water resources, particularly during low-flow periods, of freshwater withdrawals for hydraulic fracturing (fracking). Therefore, the objective of this paper is to assess the potential impact of hydraulic fracturing on water resources in the Muskingum watershed of Eastern Ohio, USA, especially due to the trend of increased withdrawals for hydraulic fracking during drought years. The Statistical Downscaling Model (SDSM) was used to generate 30 years of plausible future daily weather series in order to capture the possible dry periods. The data generated were incorporated in the Soil and Water Assessment Tool (SWAT) to examine the level of impact due to fracking at various scales. Analyses showed that water withdrawal due to hydraulic fracking had a noticeable impact, especially during low-flow periods. Clear changes in the 7-day minimum flows were detected among baseline, current and future scenarios when the worst-case scenario was implemented. The headwater streams in the sub-watersheds were highly affected, with significant decrease in 7-day low flows. The flow alteration in hydrologically-based (7Q10, i.e. 7-day 10-year low flow) or biologically-based (4B3 and 1B3) design flows due to hydraulic fracking increased with decrease in the drainage area, indicating that the relative impact may not be as great for higher order streams. Nevertheless, change in the annual mean flow was limited to 10%.     

Comparing flow resistance law for fixed and mobile bed rills

Rills caused by run‐off concentration on erodible hillslopes have very irregular profiles and cross‐section shapes. Rill erosion directly depends on the hydraulics of flow in the rills, which may differ greatly from hydraulics of flow in larger and regular channels. In this paper, a recently theoretically deduced rill flow resistance equation, based on a power–velocity profile, was tested experimentally on plots of varying slopes (ranging from 9% to 26%) in which mobile and fixed bed rills were incised. Initially, measurements of flow velocity, water depth, cross‐section area, wetted perimeter, and bed slope, carried out in 320 reaches of mobile bed rills and in 165 reaches of fixed rills, were used for calibrating the theoretical flow resistance equation. Then the relationship between the velocity profile parameter Γ, the channel slope, and the flow Froude number was separately calibrated for the mobile bed rills and for the fixed ones. The measurements carried out in both conditions (fixed and mobile bed rills) confirmed that the Darcy–Weisbach friction factor can be accurately estimated using the proposed theoretical approach. For mobile bed rills, the data were supportive of the slope independence hypothesis of velocity, due to the feedback mechanism, stated by Govers. The feedback mechanism was able to produce quasicritical flow conditions. For fixed bed rills, obtained by fixing the rill channel, by a glue, at the end of the experimental run with a mobile bed rill, the slope independence of the flow velocity measurements was also detected. Therefore, an experimental run carried out by a rill bed fixed after modelling flow action is useful to detect the feedback mechanism. Finally, the analysis showed that, for the investigated conditions, the effect of sediment transport on the flow resistance law can be considered negligible respect to the grain roughness effect.     

The GeoClaw software for depth-averaged flows with adaptive refinement

Many geophysical flow or wave propagation problems can be modeled with two-dimensional depth-averaged equations, of which the shallow water equations are the simplest example. We describe the GeoClaw software that has been designed to solve problems of this nature, consisting of open source Fortran programs together with Python tools for the user interface and flow visualization. This software uses high-resolution shock-capturing finite volume methods on logically rectangular grids, including latitude-longitude grids on the sphere. Dry states are handled automatically to model inundation. The code incorporates adaptive mesh refinement to allow the efficient solution of large-scale geophysical problems. Examples are given illustrating its use for modeling tsunamis and dam-break flooding problems. Documentation and download information is available at www.clawpack.org/geoclaw.     

A model of equilibrium bed topography for meander bends with erodible banks

Channel curvature produces secondary currents and a transverse sloping channel bed, along which the depth increases towards the outer bank. As a result deep pools tend to form adjacent to the outer bank, promoting bank collapse. The interaction of sediment grains with the primary and secondary flow and the transverse sloping bed also causes meanders to move different grain sizes in different proportions and directions, resulting in a consistent sorting pattern. Several models have been developed to describe this process, but they all have the potential to over‐predict pool depth because they cannot account for the influence of erodible banks. In reality, bank collapse might lead to the development of a wider, shallower cross‐section and any resulting flow depth discrepancy can bias associated predictions of flow, sediment transport, and grain‐size sorting. While bed topography, sediment transport and grain sorting in bends will partly be controlled by the sedimentary characteristics of the bank materials, the magnitude of this effect has not previously been explored. This paper reports the development of a model of flow, sediment transport, grain‐size sorting, and bed topography for river bends with erodible banks. The model is tested via intercomparison of predicted and observed bed topography in one low‐energy (5·3 W m^?2 specific stream power) and one high‐energy (43·4 W m^?2) study reach, namely the River South Esk in Scotland and Goodwin Creek in Mississippi, respectively. Model predictions of bed topography are found to be satisfactory, at least close to the apices of bends. Finally, the model is used in sensitivity analyses that provide insight into the influence of bank erodibility on equilibrium meander morphology and associated patterns of grain‐size sorting. The sensitivity of meander response to bank cohesion is found to increase as a function of the available stream power within the two study bends. Copyright © 2002 John Wiley & Sons, Ltd.     

引水结构变化对天津于桥水库磷滞留的影响分析与生态水量估算

2014年南水北调中线一期工程通水后,天津市的水源结构发生变化,由单一滦河水源变为双水源.引水结构变化导致于桥水库入库水量变化,从而影响水库的总磷(TP)滞留.一般认为,入库总磷负荷随入库水量减少而降低,且水力停留时间越长,将越有利于总磷滞留,从而水库TP浓度降低.但根据统计数据发现,南水北调通水后于桥水库入库水量降低,TP负荷反而升高,水库TP浓度升高.何种TP滞留机制造成了这样的结果?本文通过分析2001-2018年间入库水量(W_IN)、入库总磷负荷(TPL_IN)在南水北调通水前后的变化规律,采用Vollenweider模型推求水库TP浓度和磷滞留量(R_P),探讨水力停留时间(T)对库区总磷滞留量的影响机制,并估算双水源新情势下控制于桥水库磷滞留量的生态水量.结果表明:TPL_IN随W_IN呈先降低后升高趋势;引水期τ与R_P之间存在显著的正相关关系,随着τ增加,正相关性下降,但同时考虑非引水期和引水期τ时,R_P受较小的TPL_...     

Concentration profiles for fine and coarse sediments suspended by waves over ripples: An analytical study with the 1-DV gradient diffusion model

Field and laboratory measurements of suspended sediments over wave ripples show, for time-averaged concentration profiles in semi-log plots, a contrast between upward convex profiles for fine sand and upward concave profiles for coarse sand. Careful examination of experimental data for coarse sand shows a near-bed upward convex profile beneath the main upward concave profile. Available models fail to predict these two profiles for coarse sediments. The 1-DV gradient diffusion model predicts the main upward concave profile for coarse sediments thanks to a suitable β(y)$\beta (y)$-function (where β$\beta$ is the inverse of the turbulent Schmidt number and y   is the distance from the bed). In order to predict the near-bed upward convex profile, an additional parameter α$\alpha$ is needed. This parameter could be related to settling velocity (α$\alpha$ equal to inverse of dimensionless settling velocity) or to convective sediment entrainment process. The profiles are interpreted by a relation between second derivative of the logarithm of concentration and derivative of the product between sediment diffusivity and α$\alpha$.     

Accuracy of the kinematic wave and diffusion wave approximations for time-independent flows with infiltration included

Errors in the kinematic wave and diffusion wave approximation for time-independent (or steady-state) cases of channel flow with infiltration were derived for three types of boundary conditions: zero flow at the upstream end, and critical flow depth and zero depth gradient at the downstream end. The diffusion wave approximation was found to be in excellent agreement with the dynamic wave approximation, with errors of less than 1·4% for KF²₀≥7·5, and up to 14% for KF²₀≤0·75 for the upstream boundary condition of zero discharge and finite depth, where K is the kinematic wave number and F₀ is the Froude number. The kinematic wave approximation was reasonably accurate except at the channel boundaries and for small values of KF²₀ (≤1). The accuracy of these approximations was significantly influenced by the downstream boundary, both in terms of the magnitude of the error and the segment of the channel reach for which these approximations would be applicable. © 1998 John Wiley & Sons, Ltd.     

Crossover from capillary fingering to compact invasion for two-phase drainage with stable viscosity ratios

Motivated by a wide range of applications from enhanced oil recovery to carbon dioxide sequestration, we have developed a two-dimensional, pore-level model of immiscible drainage, incorporating viscous, capillary, and gravitational effects. This model has been validated quantitatively, in the very different limits of zero viscosity ratio and zero capillary number; flow patterns from modeling agree well with experiment. For a range of stable viscosity ratios (μ_injected/μ_displaced ? 1), we have increased the capillary number, N_c, and studied the way in which the flows deviate from capillary fingering (the fractal flow of invasion percolation) and become compact for realistic capillary numbers. Results exhibiting this crossover from capillary fingering to compact invasion are presented for the average position of the injected fluid, the fluid–fluid interface, the saturation and fractional flow profiles, and the relative permeabilities. The agreement between our results and earlier theoretical predictions [Blunt M, King MJ, Scher H. Simulation and theory of two-phase flow in porous media. Phys Rev A 1992;46:7680–99; Lenormand R. Flow through porous media: limits of fractal patterns. Proc Roy Soc A 1989;423:159–68; Wilkinson D. Percolation effects in immiscible displacement. Phys Rev A 1986;34:1380–90; Xu B, Yortsos YC, Salin D. Invasion Percolation with viscous forces. Phys Rev E 1998;57:739–51] supports the validity of these general theoretical arguments, which were independent of the details of the porous media in both two and three dimensions.