Threshold criterion for debris flow initiation in seasonal gullies

A series of flume experiments were done to investigate the effect of grain composition on the critical gradient and discharge of debris flows initiated in seasonal gullies. The results indicated that the critical gradient and discharge for debris flow initiation decrease initially, and then increase as the mass content of fine particles (<2 mm) increases. As the mass content of fine particles increases, the angle of repose, permeability of widely graded gravel soils, and the incipient motion conditions of the coarse grains in non-uniform sediments decrease at first, and then increase. The mass content of fine particles of all inflection points is the same. The theoretical model based on the combination of hydrodynamic force and shear stress is more applicable to the prediction of the critical gradient for debris flow initiation. The critical discharge model considering the effect of non-homogeneity of the soil and the size of coarse and fine grains provides a more accurate prediction of debris flow initiation than other models based on the mean diameter.     

Classification of natural flow regimes in Iran to support environmental flow management

Abstract

Development of environmental flow standards at the regional scale has been proposed as a means to manage the influence of hydrological alterations on riverine ecosystems in view of the rapid pace of global water resources management. Flow regime classification forms a critical part in such environmental flow assessments. We present a national-scale classification of hydrological regimes for Iran based on a set of hydrological metrics. It describes ecologically relevant characteristics of the natural hydrological regime derived from 15- to 47-year-long records of daily mean discharge data for 539 streamgauges within a 47-year period. The classification was undertaken using a fuzzy partitional method within Bayesian mixture modelling. The analysis resulted in 12 classes of distinctive flow regime types that differ in various hydrological aspects. This classification is being used for further research in regional-scale environmental flow studies in Iran.

Editor D. Koutsoyiannis     

Numerical simulation of groundwater flow patterns using flux as upper boundary

Since the 1960s, most of the studies on groundwater flow systems by analytical and numerical modelling have been based on given‐head upper boundaries. The disadvantage of the given‐head approach is that the recharge into and discharge from a basin vary with changes in hydraulic conductivity and/or basin geometry. Consequently, flow patterns simulated with given‐head boundaries but with different hydraulic conductivities and/or basin geometry may not reflect the effects of these variables. We conducted, therefore, numerical simulations of groundwater flow in theoretical drainage basins using flux as the upper boundary and realistically positioned fluid‐potential sinks while changing the infiltration intensity, hydraulic conductivities, and geometric configuration of the basin. The simulated results demonstrate that these variables are dominant factors controlling the flow pattern in a laterally closed drainage basin. The ratio of infiltration intensity to hydraulic conductivity (R_ic) has been shown to be an integrated pattern‐parameter in a basin with a given geometric configuration and possible fluid‐potential‐sink distribution. Successively, the changes in flow patterns induced by stepwise reductions in R_ic are identical, regardless of whether the reductions are due to a decrease in infiltration intensity or an increase in hydraulic conductivity. The calculated examples show five sequential flow patterns containing (i) only local, (ii) local–intermediate, (iii) local–intermediate–regional, (iv) local–regional, and (v) just regional flow systems. The R_ic was found to determine also whether a particular sink is active or not as a site of discharge. Flux upper boundary is preferable for numerical simulation when discussing the flow patterns affected by a change of infiltration, the hydraulic conductivity, or the geometry of a basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Bed scour by debris flows: experimental investigation of effects of debris‐flow composition

Debris flows can grow greatly in size by entrainment of bed material, enhancing their runout and hazardous impact. Here, we experimentally investigate the effects of debris‐flow composition on the amount and spatial patterns of bed scour and erosion downstream of a fixed to erodible bed transition. The experimental debris flows were observed to entrain bed particles both grain by grain and en masse, and the majority of entrainment was observed to occur during passage of the flow front. The spatial bed scour patterns are highly variable, but large‐scale patterns are largely similar over 22.5–35° channel slopes for debris flows of similar composition. Scour depth is generally largest slightly downstream of the fixed to erodible bed transition, except for clay‐rich debris flows, which cause a relatively uniform scour pattern. The spatial variability in the scour depth decreases with increasing water, gravel (= grain size) and clay fraction. Basal scour depth increases with channel slope, flow velocity, flow depth, discharge and shear stress in our experiments, whereas there is no correlation with grain collisional stress. The strongest correlation is between basal scour and shear stress and discharge. There are substantial differences in the scour caused by different types of debris flows. In general, mean and maximum scour depths become larger with increasing water fraction and grain size, and decrease with increasing clay content. However, the erodibility of coarse‐grained experimental debris flows (gravel fraction = 0.64) is similar on a wide range of channel slopes, flow depths, flow velocities, discharges and shear stresses. This probably relates to the relatively large influence of grain‐collisional stress to the total bed stress in these flows (30–50%). The relative effect of grain‐collisional stress is low in the other experimental debris flows (<5%), causing erosion to be largely controlled by basal shear stress. Copyright © 2016 John Wiley & Sons, Ltd.     

Classifying the flow regimes of Mediterranean streams using multivariate analysis

Rivers in the Mediterranean region often exhibit an intermittent character. An understanding and classification of the flow regimes of these rivers is needed, as flow patterns control both physicochemical and biological processes. This paper reports an attempt to classify flow regimes in Mediterranean rivers based on hydrological variables extracted from discharge time series. Long‐term discharge records from 60 rivers within the Mediterranean region were analysed in order to classify the streams into different flow regime groups. Hydrological indices (HIs) were derived for each stream and principal component analysis (PCA) and then applied to these indices to identify subsets of HIs describing the major sources of variations, while simultaneously minimizing redundancy. PCA was performed for two groups of streams (perennial and temporary) and for all streams combined. The results show that whereas perennial streams are mainly described by high‐flow indices, temporary streams are described by duration, variability and predictability indices. Agglomerative cluster analysis based on HIs identified six groups of rivers classified according to differences in intermittency and variability. A methodology allowing such a classification for ungauged catchments was also tested. Broad‐scale catchment characteristics based on digital elevation, climate, soil and land use data were derived for each long‐term station where these data were available. By using stepwise multiple regression analysis, statistically significant relationships were fitted, linking the three selected hydrological variables (mean annual number of zero‐flow days, predictability and flashiness) to the catchment characteristics. The method provides a means of simplifying the complexity of river systems and is thus useful for river basin management. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrodynamic characteristics of rill flow on steep slopes

Rill erosion is a dominant sediment source on sloping lands. However, the amount of soil loss from rills on steep slopes is vastly more than that on gentle slopes because of differences in rill shape and hydraulic patterns. The aims of this paper are to determine the hydrodynamic characteristics of rills and the friction coefficients in steep slope conditions and to propose modifications of some hydraulic parameters used in soil loss prediction models. A series of inflow experiments was conducted on loess slopes. The results show that the geometric and hydraulic properties of rill on the steep loess slopes, which are characterized by the mean width of cross sections, mean velocity and mean depth of flow, are related to discharge and slope gradient in power functions. However, the related exponents to discharge are 0.26, 0.48 and 0.26, respectively, which are different from the exponents derived in previous studies, which were conducted on gentle slopes. The Manning roughness coefficient ranged from 0.035 to 0.071, with an average of 0.0536, and the Darcy–Weisbach friction coefficients varied from 0.4 to 1.9. The roughness coefficients are closely related to the Reynolds numbers and flow volumes; however, the correlations vary with slope gradient. The roughness coefficients are directly proportional to the Reynolds number and the flow volume on steep slopes, in contrast with the roughness coefficients found on gentle slopes, which decrease as the Reynolds number and flow volume increase. This difference is caused by the interactions among the hydraulics of the flow, the shape of the rills and the sediment concentrations on steep slopes. The results indicate that parameters used in models to predict rill erosion have to be modified according to slope gradient. Copyright © 2015 John Wiley & Sons, Ltd.     

Simplified nonlinear modeling of river flow routing

Management of open-channel flow systems requires accurate models of flow transfer. This article presents a simple nonlinear model representative of the flow transfer in a river reach. The model is obtained through linearization of a physical model, simplification using the cumulant matching method and analytic identification of a nonlinear model coinciding with the linear model around equilibrium points, corresponding to the hydraulic permanent regimes. The methodology is illustrated on the diffusive wave equation and the Saint-Venant equations. The obtained nonlinear models are compared in simulation to the initial models. The nonlinear model is shown to ensure mass conservation, despite the variable delay element of the model. The proposed model can reproduce the nonlinear behavior of the time-delay with discharge variations. It is well-suited for fast simulations, flow forecasting, and for controller design.     

Stormflow generation involving pipe flow in a zero‐order basin of Peninsular Malaysia

Hydrological responses in a zero‐order basin (ZOB), a portion of whose discharge emerged via preferential flow through soil pipes, were examined over a 2‐year period in Peninsular Malaysia to elucidate primary stormflow generation processes. Silicon (Si) and specific conductance (EC) in various runoff components were also measured to identify their sources. ZOB flow response was dependent on antecedent precipitation amount; runoff increased linearly with precipitation during events >20 mm in relatively wet antecedent moisture conditions. Runoff derived from direct precipitation falling onto saturated areas accounted for <0·2% of total ZOB flow volume during the study period, indicating the predominance of subsurface pathways in ZOB flow. ZOB flow (high EC and low Si) was distinct from perennial baseflow via bedrock seepage (low EC and high Si) 5 m downstream of the ZOB outlet. Pipe flow responded quickly to ZOB flow rate and was characterized by a threshold flow capacity unique to each pipe. Piezometric data and pipe flow records demonstrated that pipes located deeper in the soil initiated first, followed by those at shallower depths; initiation of pipe flow corresponded to shallow groundwater rise above the saprolite‐soil interface. Chemical signatures of pipe flow were similar to each other and to the ZOB flow, suggesting that the sources were well‐mixed soil‐derived shallow groundwater. Based upon the volume of pipe flow during storms, the combined contribution of the pipes monitored accounted for 48% of total ZOB flow during the study period. Our results suggest that shallow groundwater, possibly facilitated by preferential flow accreted above the saprolite–soil interface, provides dominant stormflow, and that soil pipes play an important role in the rapid delivery of solute‐rich water to the stream system. Copyright © 2006 John Wiley & Sons, Ltd.     

Analytical findings for power law cross-sections: Uniform flow depth

Analytical results concerning open channel flows are presented, assuming that the cross-section is defined by a power law relationship between the channel width and the channel depth. Explicit equations to compute the normal flow depth are derived by considering the liquid discharge, the channel roughness height and the cross-section geometry (based on knowledge of the power law exponent, the reference width, and the reference depth) as known quantities.Such equations are deduced by writing the physical quantities as a power expansion in the power law exponent and expressing the wetted perimeter using a Gauss hypergeometric function. With the designed procedure, an accurate estimations of the integrals required to invert the uniform flow formula are obtained, at least for cross-sections characterized by aspect ratios of technical interest.Two relationships are proposed between the normal depth and the flow discharge. The first relationship is shown to work well for any discharge, provided that the width to depth ratio is sufficiently large. If this is not the case, the second procedure must be used for non-dimensional discharge larger than a given threshold, while the former procedure remains valid under the threshold.     

Fine‐scale characteristics of groundwater flow in a peatland

Fine‐scale dynamics of groundwater flow were studied in a 1·5 ha peatland in central New York. Measurements of the hydraulic head throughout a detailed network of piezometer clusters revealed spatial and temporal variability in the direction of groundwater flow at a very fine (within a few metres) scale of analysis. Within the small wetland, there were areas of groundwater recharge, discharge and lateral flow. Such patterns of groundwater flow frequently reversed or changed due to fluctuations of only a few centimetres in hydraulic head. Specific conductance, deuterium signatures and calcium concentrations of groundwater corroborated the groundwater flow patterns determined with hydraulic head measurements and illustrated the influence of source water chemistry and evaporation on different layers in the peat column. The control of peat chemistry by such fine‐scale groundwater flow may have important implications for plant community composition and diversity in groundwater‐fed peatlands. Copyright © 1999 John Wiley & Sons, Ltd.     

River flow regimes in a changing climate

Abstract

A river flow regime describes an average seasonal behaviour of flow and reflects the climatic and physiographic conditions in a basin. Differences in the regularity (stability) of the seasonal patterns reflect different dimensionality of the flow regimes, which can change subject to changes in climate conditions. The empirical orthogonal functions (EOF) approach can be used to describe the intrinsic dimension of river flow regimes and is also an adopted method for reducing the phase space in connection to climate change studies, especially in studies of nonlinear dynamic systems with preferred states. A large data set of monthly river flow for the Nordic countries has been investigated in the phase space reduced to the first few amplitude functions to trace a possible signature of climate change on the seasonal flow patterns. The probability density functions (PDF) of the weight coefficients and their possible change over time were used as an indicator of climate change. Two preferred states were identified connected to stable snowmelt-fed and rainfed flow regimes. The results indicate changes in the PDF patterns with time towards higher frequencies of rainfed regime types. The dynamics of seasonal patterns studied in terms of PDF renders it an adequate and convenient characterization, helping to avoid bias connected to flow regime classifications as well as uncertainties inferred by a modelling approach.     

An experimental study of discharge partitioning and flow structure at symmetrical bifurcations

Recent research has examined the factors controlling the geometrical configuration of bifurcations, determined the range of stability conditions for a number of bifurcation types and assessed the impact of perturbations on bifurcation evolution. However, the flow division process and the parameters that influence flow and sediment partitioning are still poorly characterized. To identify and isolate these parameters, three‐dimensional velocities were measured at 11 cross‐sections in a fixed‐walled experimental bifurcation. Water surface gradients were controlled, and systematically varied, using a weir in each distributary. As may be expected, the steepest distributary conveyed the most discharge (was dominant) while the mildest distributary conveyed the least discharge (was subordinate). A zone of water surface super‐elevation was co‐located with the bifurcation in symmetric cases or displaced into the subordinate branch in asymmetric cases. Downstream of a relatively acute‐angled bifurcation, primary velocity cores were near to the water surface and against the inner banks, with near‐bed zones of lower primary velocity at the outer banks. Downstream of an obtuse‐angled bifurcation, velocity cores were initially at the outer banks, with near‐bed zones of lower velocities at the inner banks, but patterns soon reverted to match the acute‐angled case. A single secondary flow cell was generated in each distributary, with water flowing inwards at the water surface and outwards at the bed. Circulation was relatively enhanced within the subordinate branch, which may help explain why subordinate distributaries remain open, may play a role in determining the size of commonly‐observed topographic features, and may thus exert some control on the stability of asymmetric bifurcations. Further, because larger values of circulation result from larger gradient disadvantages, the length of confluence–diffluence units in braided rivers or between diffluences within delta distributary networks may vary depending upon flow structures inherited from upstream and whether, and how, they are fed by dominant or subordinate distributaries. Copyright © 2011 John Wiley & Sons, Ltd.     

Back‐barrier flooding by storm surges and overland flow

Low‐elevation areas within a sandy barrier island are subject to flooding via saturation overland flow following moderate storm surges and rainfall events. Using a high resolution topographic survey and simple hydrology models, we estimate the discharge and velocities from storm surge return flow and saturation overland flow. Results show that return flow velocities are of the same magnitude as the critical velocity necessary to mobilize sand when a hydraulic connection between the watershed and back‐barrier bay is present. Storms of moderate strength and rainfall intensity may be sufficient to keep the return channels open within the back‐barrier, thus providing natural conduits for water exchange from overwash events during extreme storm surges triggered by hurricanes. Copyright © 2011 John Wiley & Sons, Ltd.     

Persistency of flow patterns in a water repellent sandy soil – Conclusions of TDR readings and a time-delayed double tracer experiment

On a former waste water disposal field with water repellent sandy soil under grass vegetation we analyzed the persistency of flow patterns on a 150 m × 25 m plot by (i) continuous TDR-measurements on a 2 m × 1 m transect combined with seasonal soil moisture sampling campaigns, and (ii) a time-delayed double tracer experiment on a second 3 m × 1 m transect. Here, we applied bromide under wettable soil conditions in spring and chloride under water repellent soil conditions in autumn. At the end of the tracer experiment, after a travel time of 328 days for Br and 87 days for Cl, respectively, the transect was excavated and sampled in high spatial resolution. Tracer concentration, water content, water drop penetration times (WDPT), and soil organic matter content (SOM) of each sample were analyzed in order to characterize flow patterns. The TDR readings were used to predict the effective cross section (ECS) of subsurface flow and flow shifts over the season.During summer, when ECS is low and consecutive precipitation events occur, flow paths – once created – persist over time. However, over longer times (from autumn to autumn), the spatial arrangements of the flow paths can change completely. The Cl distribution showed typical fingering structures with high concentrations in the less water repellent flow paths. In contrast, Br was found mostly in the dry, hydrophobic areas indicating that it was transported before the soil became water repellent. Consequently, the flow patterns generated in spring and early summer differ completely from those in autumn and winter because of water repellent structures established during the vegetative period. These structures could be identified using a critical water content (θ_crit) concept, considering both soil water content and SOM.As not all soil parts being active during to season, four flow categories could be identified: about 10% permanent (=stable flow paths), 45% periodic (i.e. water repellent in summer), 40% occasional (water repellent in summer and autumn), and 5% permanent water repellent.     

Critical discharge for entrainment of poorly sorted gravel

Entrainment criteria for gravel-sized particles on river beds are usually defined in terms of shear stress, but some workers have proposed that critical-discharge relationships are more reliable and more convenient. The critical shear stress in poorly sorted gravels is known from recent work to depend much more on relative than absolute grain size. It is shown that this effect can be included in a semi-theoretical equation for the critical unit discharge to move individual grain sizes in a poorly sorted bed. Comparisons with published data on critical discharges in a variety of rivers show that both the sensitivity of critical discharge to grain size, and the absolute levels of critical discharge, are predicted quite well. However, any prediction of critical discharge is sensitive to the assumptions made, first about flow resistance and secondly about the critical stress to move average-sized particles.     

Evolution of overland flow connectivity in bare agricultural plots

Soil surface roughness not only delays overland flow generation but also strongly affects the spatial distribution and concentration of overland flow. Previous studies generally aimed at predicting the delay in overland flow generation by means of a single parameter characterizing soil roughness. However, little work has been done to find a link between soil roughness and overland flow dynamics. This is made difficult because soil roughness and hence overland flow characteristics evolve differently depending on whether diffuse or concentrated erosion dominates. The present study examined whether the concept of connectivity can be used to link roughness characteristics to overland flow dynamics. For this purpose, soil roughness of three 30‐m² tilled plots exposed to natural rainfall was monitored for two years. Soil micro‐topography was characterized by means of photogrammetry on a monthly basis. Soil roughness was characterized by the variogram, the surface stream network was characterized by network‐based indices and overland flow connectivity was characterized by Relative Surface Connection function (RSCf) functional connectivity indicator. Overland flow hydrographs were generated by means of a physically‐based overland flow model based on 1‐cm resolution digital elevation models. The development of eroded flow paths at the soil surface not only reduced the delay in overland flow generation but also resulted in a higher continuity of high flow velocity paths, an increase in erosive energy and a higher rate of increase of the overland flow hydrograph. Overland flow dynamics were found to be highly correlated to the RSCf characteristic points. By providing information regarding overland flow dynamics, the RSCf may thus serve as a quantitative link between soil roughness and overland flow generation in order to improve the overland flow hydrograph prediction. Copyright © 2016 John Wiley & Sons, Ltd.     

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.     

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.     

Experimental investigation of the morphodynamic response of riffles and pools to unsteady flow and increased sediment supply

Gravel-bed rivers characteristically exhibit shallow riffles in wide sections and deeper pools where the channel becomes constricted and narrow. While rivers can adjust to changing flow and sediment supply through some combination of adjustments of channel slope, bed-surface sorting, and channel shape, the degree to which riffle-pools may adopt these changes in response to changing flows and sediment supplies remains unclear. This article presents results from a flume experiment investigating how constant- and variable-width channels adjust their morphology in response to changing flow and increased sediment supply. Two flume geometries were used: (1) constant-width and (2) variable-width, characterized by a sinusoidal pattern with a mean width equal to that of the first channel. The variable-width channel developed bed undulations in phase with the width, representing riffle-pools. The experiment consisted of three phases for each flume geometry: (1) steady flow, constant sediment supply; (2) unsteady flow, constant sediment supply; and (3) unsteady flow, doubled sediment supply. Unsteady flow was implemented in the form of repeated symmetrical stepped hydrographs, with a mean discharge equal to that in the steady flow phase. In all phases the bed and sediment supply were composed of a sand/gravel mixture ranging from 1 to 8 mm. In both the straight and variable-width channels, transitioning from steady flow to repeated hydrographs did not result in significant changes in bed morphology. The two channel geometries had different responses to increased sediment supply: the slope of the straight channel increased nearly 40%, while the variable-width channel reduced the relief between bars and pools and decreased the variability in cross-sectional elevation with a slight slope increase. Bar-pool relief varied with repeat discharge hydrographs. Pool elevation changed twice the distance of bar elevations, emphasizing the relevance of pool scour for riffle-pool self-maintenance in channels with width variations.