A continuous simulation model for design-hydrograph estimation in small and ungauged watersheds

Abstract

A detailed investigation of the behaviour of various hydraulic parameters, using data from rivers in Greece, was conducted in order to explore the universality of features that many natural streams are believed to have in common. Analysis of vertical profiles of temporal mean of horizontal velocities (u) in the longitudinal (river flow) direction and of transverse profiles of depth-mean longitudinal velocities (U) estimated from these vertical profiles, measured at 232 cross-sections of several rivers in Greece, provided valuable information: on the distribution of local roughness coefficients (n_i ) along the wetted perimeter of the cross-sections examined; on the shape of u profiles; on the ratio of maximum to mean cross-sectional velocity, V_max/V_m , and its relation to a dimensionless entropy parameter, M; on the shape of U profiles; and on the normalized intensity, r, of the spatial departure of u velocities from V_m . The similarities among the quantities (u, U, n, V_max/V_m , M, r) analysed in this study and in pertinent literature reveal that the rivers examined exhibit many of the basic features, of rather universal character, shown by other rivers (all over the world) having different geometric and/or other characteristics (aspect ratios, bottom roughness, flow kinematics, etc.). Corresponding differences are also described and explained.     

P‐wave velocity distribution in basalt flows of the Enni Formation in the Faroe Islands from refraction seismic analysis

The main objective of this work is to establish the applicability of shallow surface‐seismic traveltime tomography in basalt‐covered areas. A densely sampled ～1300‐m long surface seismic profile, acquired as part of the SeiFaBa project in 2003 ( Japsen et al. 2006 ) at Glyvursnes in the Faroe Islands, served as the basis to evaluate the performance of the tomographic method in basalt‐covered areas. The profile is centred at a ～700‐m deep well. V_P, V_S and density logs, a zero‐offset VSP, downhole‐geophone recordings and geological mapping in the area provided good means of control. The inversion was performed with facilities of the Wide Angle Reflection/Refraction Profiling program package ( Ditmar et al. 1999 ). We tested many inversion sequences while varying the inversion parameters. Modelled traveltimes were verified by full‐waveform modelling. Typically an inversion sequence consists in several iterations that proceed until a satisfactory solution is reached. However, in the present case with high velocity contrasts in the subsurface we obtained the best result with two iterations: first obtaining a smooth starting model with small traveltime residuals by inverting with a high smoothing constraint and then inverting with the lowest possible smoothing constraint to allow the inversion to have the full benefit of the traveltime residuals. The tomogram gives usable velocity information for the near‐surface geology in the area but fails to reproduce the expected velocity distribution of the layered basalt flows. Based on the analysis of the tomogram and geological mapping in the area, a model was defined that correctly models first arrivals from both surface seismic data and downhole‐geophone data.     

Laboratory flume studies of microflow environments of aquatic plants

The microflow environments of aquatic plants with reference to Myriophyllum and Hydrilla are simulated in a laboratory flume. A Nix Streamflow microflow meter was used to measure the mean velocity profiles of flow at different densities of plants, flow ranges and measurement positions. Each mean velocity profile consists of three hydrodynamic regimes (i.e. within‐canopy zone, above‐canopy zone and a transitional zone between them), which indicate the presence of two benthic boundary layers (internal and external ones). Out of 38 measured mean velocity profiles, most do not fit a logarithmic relationship. The following hydrodynamic parameters are used in characterizing the flow regimes: local shear velocity (u_*), roughness length (z_o), canopy roughness Reynolds number (Re_*), bed shear stress (τ_o) and laminar sublayer (σ). Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

Shear velocity estimates in rough‐bed open‐channel flow

Shear velocity u_* is an important parameter in geophysical flows, in particular with respect to sediment transport dynamics. In this study, we investigate the feasibility of applying five standard methods [the logarithmic mean velocity profile, the Reynolds stress profile, the turbulent kinetic energy (TKE) profile, the wall similarity and spectral methods] that were initially developed to estimate shear velocity in smooth bed flow to turbulent flow over a loose bed of coarse gravel (D₅₀ = 1·5 cm) under sub‐threshold conditions. The analysis is based on quasi‐instantaneous three‐dimensional (3D) full depth velocity profiles with high spatial and temporal resolution that were measured with an Acoustic Doppler Velocity Profiler (ADVP) in an open channel. The results of the analysis confirm the importance of detailed velocity profile measurements for the determination of shear velocity in rough‐bed flows. Results from all methods fall into a range of ± 20% variability and no systematic trend between methods was observed. Local and temporal variation in the loose bed roughness may contribute to the variability of the logarithmic profile method results. Estimates obtained from the TKE and Reynolds stress methods reasonably agree. Most results from the wall similarity method are within 10% of those obtained by the TKE and Reynolds stress methods. The spectral method was difficult to use since the spectral energy of the vertical velocity component strongly increased with distance from the bed in the inner layer. This made the choice of the reference level problematic. Mean shear stress for all experiments follows a quadratic relationship with the mean velocity in the flow. The wall similarity method appears to be a promising tool for estimating shear velocity under rough‐bed flow conditions and in field studies where other methods may be difficult to apply. This method allows for the determination of u_* from a single point measurement at one level in the intermediate range (0·3 < h < 0·6). Copyright © 2013 John Wiley & Sons, Ltd.     

Soil pipe flow tracer experiments: 2. Application of a streamflow transient storage zone model

Soil pipes are important subsurface flow pathways in many soil erosion phenomena. However, limited research has been performed on quantifying and characterizing their flow and transport characteristics. The objectives of this research were to determine the applicability of a streamflow model with transient storage in deriving flow and transport characteristics of soil pipes. Tracer data from pulse inputs were collected in four different soil pipes after a fluorescein dye was injected in the upstream end of each soil pipe network in three branches (west, middle, and east) of a main catchment and a back catchment in Goodwin Creek Experimental Watershed in Mississippi. Multiple sampling stations were positioned along each soil pipe network. The transient storage zone model OTIS‐P was executed inversely to estimate transport parameters by soil pipe reach such as the soil pipe cross‐sectional area (A), soil storage zone cross‐sectional area (A_s), and exchange rate between the soil pipe and the soil storage zone (α_s). Model convergence was achieved, and simulated breakthrough curves of the reaches were in good agreement with actual tracer data for eight of the nine reaches of the three branches of the Main Catchment and five of the seven reaches of the Back Catchment soil pipe. Simulation parameters for the soil pipe networks were similar to the range of values reported for flow and transport characteristics commonly observed in streams. Inversely, estimated soil pipe flow velocities were higher with increased tortuosity, which led to a smaller cross‐sectional areas predicted for the soil pipe flowpaths, while other parameters were not sensitive to tortuosity. In general, application of One‐Dimensional Transport with Inflow and Storage‐P to this unique soil pipe condition suggested larger transient storage (A_s and α_s) compared with most stream systems. This was hypothesized to be because of relatively higher ratio of the wetted perimeter to flow area in the soil pipe, the hydraulic roughness of the soil pipe, potential retention in collapsed portions of the pipe, and interaction with smaller preferential flow systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Effects of patterned Artemisia capillaris on overland flow velocity under simulated rainfall

Vegetation cover is an important factor for erosion control. Laboratory‐simulated rainfall experiments were conducted to quantify the effectiveness of patchy distributed Artemisia capillaris in retarding overland flow velocity. Simulated storms (60, 90, 120, and 150 mm h^?1) were applied on a bare plot (CK) and four different plant patterns, a banded pattern perpendicular to the slope direction (BP), a single long strip parallel to slope direction (LP), small patches distributed like a checkerboard (SP1), and small patches distributed like a letter “X” (SP2). All treatments had three replicates. Each plot underwent two sets of experiments, intact plant plots and root plots (the above‐ground parts were removed, only roots were reserved), respectively. Results showed that flow velocity increased with rainfall intensity, and the lower slope velocity (V_l) was higher than the upper slope velocity (V_u). The removal of grass shoots increased flow velocity. Compared with bare soil plot, intact plants reduced mean flow velocity by 14%–60%, whereas the reduction declined to <40% for the root plots. BP and both SP treatments performed more effectively than LP in retarding flow velocity, whereas no significant differences were identified between BP and SP. The contributions of A. capillaris shoots and roots to the reductions in flow velocity under different rainfall intensities were different. The shoots made greater contribution of 53%–97% at 60 and 90 mm h^–1, and the roots contributed more (51%–81%) at 120 and 150 mm h^–1. Runoff and sediment rate had significant (p < 0.05) linear correlations with mean flow velocity. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Hazard mapping of normalized peak strain in soils during earthquakes: microzonation of a metropolitan area

Seismic hazard maps of the Los Angeles metropolitan area are illustrated for normalized peak strain and for 50 years of exposure. The strain estimates are based on scaling in terms of peak ground velocity. The proportionality factor is the phase velocity with which the wave energy is propagating. A simplified seismicity model is used in which all earthquakes occur on faults represented by buried lines and in one zone of diffused seismicity. Poissonian model of earthquake occurrence is assumed. The same model was used in the 1980's to illustrate a method for microzoning of the same area for response spectral amplitudes. Maps of logarithms of normalized peak strain, cε_max, are presented for probabilities of at least one exceedance p = 0·99, 0·9, 0·5, 0·1 and 0·01. These can be used to construct site specific probability distribution functions of the normalized peak strain, cε_max. Such maps are useful for design of new and for retrofit of existing structures, sensitive to strain and differential ground motions (bridges, tunnels, pipelines, etc.).     

Comparing response of SDF systems to near‐fault and far‐fault earthquake motions in the context of spectral regions

In spite of important differences in structural response to near‐fault and far‐fault ground motions, this paper aims at extending well‐known concepts and results, based on elastic and inelastic response spectra for far‐fault motions, to near‐fault motions. Compared are certain aspects of the response of elastic and inelastic SDF systems to the two types of motions in the context of the acceleration‐, velocity‐, and displacement‐sensitive regions of the response spectrum, leading to the following conclusions. (1) The velocity‐sensitive region for near‐fault motions is much narrower, and the acceleration‐sensitive and displacement‐sensitive regions are much wider, compared to far‐fault motions; the narrower velocity‐sensitive region is shifted to longer periods. (2) Although, for the same ductility factor, near‐fault ground motions impose a larger strength demand than far‐fault motions—both demands expressed as a fraction of their respective elastic demands—the strength reduction factors R_y for the two types of motions are similar over corresponding spectral regions. (3) Similarly, the ratio u_m/u₀ of deformations of inelastic and elastic systems are similar for the two types of motions over corresponding spectral regions. (4) Design equations for R_y (and for u_m/u₀) should explicitly recognize spectral regions so that the same equations apply to various classes of ground motions as long as the appropriate values of T_a, T_b and T_c are used. (5) The Veletsos–Newmark design equations with T_a=0.04 s, T_b=0.35 s, and T_c=0.79 s are equally valid for the fault‐normal component of near‐fault ground motions. Copyright © 2001 John Wiley & Sons, Ltd.     

Transpiration from three dominant shrub species in a desert‐oasis ecotone of arid regions of Northwestern China

Conservation management for the water dependent desert‐oasis ecotone in arid northwest China requires information on the water use of the dominant species. However, no studies have quantified their combined water use or linked species composition to ecotone transpiration. Here, the water use of three dominant shelterbelt shrubs (Haloxylon ammodendron, Nitraria tangutorum, and Calligonum mongolicum) within an ecotone was measured throughout the full leaf‐out period for three shrub species from 30 May to 16 October 2014, with sap flow gauges using the stem heat balance approach. Species‐specific transpiration was estimated by scaling up sap flow velocities measured in individual stems, to stand area level, using the frequency distribution of stem diameter and assuming a constant proportionality between sap flow velocity and basal cross‐sectional area for all stems. The mean peak sap flux densities (J_sn) for H. ammodendron, N. tangutorum, and C. mongolicum, were 40.12 g cm^?2 h^?1, 71.33 g cm^?2 h^?1, and 60.34 g cm^?2 h^?1, respectively, and the mean estimated daily area‐averaged transpiration rates (T_daily) for the same species were 0.56 mm day^?1, 0.34 mm day^?1, and 0.11 mm day^?1. The accumulative stand transpiration was approximately 140.8 mm throughout the measurement period, exceeding precipitation by as much as 42.1 mm. Furthermore, T_daily of these shrubs appeared to be much less sensitive to soil moisture as compared to atmospheric drivers, and the relationship between J_sn and atmospheric drivers was likely uninfluenced by soil moisture regimes in the whole profile (to 1‐m depth), especially for H. ammodendron and C. mongolicum. Results indicate that these shrubs may use deep soil water recharged by capillary rise, or may directly access shallow groundwater. This study provides quantitative data offering important implications for ecotone conservation and water and land resource management. Copyright © 2016 John Wiley & Sons, Ltd.     

The limits to growth: how large can subaqueous,flow-transverse bedforms ultimately become?

Based on field and experimental evidence, the average initial spacing (seed wavelength) of flow-transverse bedforms (ripples and dunes) appears to lie between 80 and 130 grain diameters (L = 80–130D_mm). Starting with an average initial spacing of L = 100D_mm, subsequent bedform growth proceeds by amalgamation of two successive bedforms, which results in a doubling of the spacing in each step. Geometric principles dictate that the combined volume of two smaller bedforms lacks about 40% of the volume required for a fully developed amalgamated bedform. The missing volume is gained by excavation of the troughs, i.e., by lowering the base level. Where base level lowering is prevented by the presence of a coarse-grained armor layer or hard ground pavement, the larger amalgamated bedform remains sediment starved. In its simplest form, bedform growth proceeds by continuous doubling of the spacing in response to increases in flow velocity, the process being reversible in response to flow decelerations. Bedform growth terminates when the shear velocity (u_*) at the crest reaches the mean settling velocity (w_s) of the sediment. At this point, 40% of the bed material is in suspension, at which point the missing volume can no longer be compensated by trough excavation. In shallow water, maximum bedform size is dictated by the water depth, whereas in deep water, bedforms can potentially grow to their ultimate size. Evaluation of bedform data from deep water settings suggests that the largest two-dimensional, flow-transverse bedforms in terms of grain size (phi) can be approximated by the equations: lnL_max = 13.72–4.03D_phi and lnH_max = 9.95–3.47D_phi for grain sizes <  ~ 0.2 mm (> ~ 2.32 phi), with L and H representing bedform spacing and height in meters and D the grain size in phi. For grain sizes >  ~ 0.2 mm (< ~ 3.23 phi), the corresponding relationships are lnL_max = 6.215–0.69 D_phi and lnH_max = 3.18–0.56D_phi, with notations as before, or in terms of grain diameters in mm: L_max = 5 × 10⁵D_mm.

     

Real‐time monitoring and estimation of the discharge of flash floods in a steep mountain catchment

Synchronously and accurately estimating the flood discharges and dynamic changes in the fluid density is essential for hydraulic analysis and forecasting of flash floods, as well as for risk assessment. However, such information is rare for steep mountain catchments, especially in regions that are hotspots for earthquakes. Therefore, six hydrological monitoring sites were established in the main stream and tributaries of the 78.3‐km² Longxi River catchment, an affected region of the Wenchuan earthquake region in China. Direct real‐time monitoring equipment was installed to measure the flow depths, velocities, and fluid total pressures of the flood hydrographs. On the basis of field measurements, real‐time mean cross‐sectional velocities during the flood hydrographs could be derived from easily obtainable parameters: cross‐sectional maximum velocities and the calibrated dimensionless parameter K_h . Real‐time discharges were determined on the basis of a noncontact method to establish the effective rating curves of this mountainous stream, ranging from 1.46 to 386.34 m³/s with the root mean square errors of ≤10.22 m³/s. Compared with the traditional point‐velocity method and empirical Manning's formula, the proposed noncontact method was reliable and safe for monitoring whole flood hydrographs. Additionally, the real‐time fluid density during the flood hydrographs was calculated on the basis of the direct monitoring parameters for fluid total pressures and water depths. During the flood hydrograph, transient flow behaviour with higher fluid density generally occurred downstream during the flood peak periods when the flow was in the supercritical flow regime. The observed behaviour greatly increased the threat of damage to infrastructure and human life near the river. Thus, it is important to accurately estimate flood discharge and identify for fluid densities so that people at risk from an impending flash flood are given reliable, advanced warning.     

Prediction of Remediation of a Heterogeneous Aquifer: A Case Study

Contaminant plumes whose characteristic length is smaller than the horizontal integral scale of the hydraulic conductivity, K, are abundant in shallow, phreatic aquifers. In such cases, the aquifer can be regarded as layered, with K being only a function of the vertical coordinate. The heterogeneity of K has a critical role upon the efficiency of remediation of such sites, for example, by Pump and Treat schemes. The expected efficiency is a random variable, with uncertainty. Quantifying this uncertainty can be of great importance to decision making. In this study, we focus on a case study in the coastal aquifer of Israel and compare two different approaches for constructing realizations of K: continuous and indicator. We observe a significant difference between the constructed realizations, which results in a considerable difference in the predicted remediation efficiency and its uncertainty. Furthermore, we study the effect of conditioning the realizations by a rather limited number of K data points. We find that the conditioning results in a major reduction of the uncertainty. In addition, we compare the results of the transport model to a simplified semi‐analytical solution that is based on assuming radial flow. We find a good agreement with the three‐dimensional numerical model. This result illustrates that the simplified solution can be used for prediction of the remediation efficiency when the flow at the plume vicinity can be regarded as radial.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.     

Using a fluctuating tracer to estimate hyporheic exchange in restored and unrestored reaches of the Truckee River,Nevada, USA

The goal of this research was to compare hyporheic activity in recently restored and adjacent un‐restored reaches of the Truckee River downstream from the Reno/Sparks metropolitan area. The installation of rocky riffles and raised channel bed elevations in the restored reaches may have increased the degree of surface–subsurface interaction. A fluctuating chloride concentration signal served as the tracer, induced by the variable influx of higher salinity water several miles upstream from the study reach. The solute transport model, OTIS, was used in conjunction with the hydrodynamic model, DYNHYD5, to estimate transient storage parameters under unsteady flow conditions. The model was calibrated to chloride concentrations measured over a period of three days at six in‐stream locations representing restored and un‐restored reaches. An automated parameter estimation algorithm (SCE‐UA) was used to optimize parameters for multiple reaches simultaneously and generate a distribution of parameter estimates. Results suggest that the transient storage zone cross‐sectional area (A_s) is larger in the restored reaches than in the unrestored reaches, but the exchange coefficient (α) is smaller, leading to increased hyporheic residence time and hydrologic retention in the vicinity of channel reconstructions. Scenarios were used to simulate the potential effects of increased subsurface residence time on denitrification and in‐stream NO₃‐N concentrations. Monte Carlo analysis was performed to assess uncertainty in the simulation results and show the potential for greater nutrient retention in the lower Truckee River as a result of channel restoration. Copyright © 2009 John Wiley & Sons, Ltd.     

Evaluating Groundwater Interbasin Flow Using Multiple Models and Multiple Types of Data

The groundwater interbasin flow, Q_y, from the north of Yucca Flat into Yucca Flat simulated using the Death Valley Regional Flow System (DVRFS) model greatly exceeds assessments obtained using other approaches. This study aimed to understand the reasons for the overestimation and to examine whether the Q_y estimate can be reduced. The two problems were tackled from the angle of model uncertainty by considering six models revised from the DVRFS model with different recharge components and hydrogeological frameworks. The two problems were also tackled from the angle of parametric uncertainty for each model by first conducting Morris sensitivity analysis to identify important parameters and then conducting Monte Carlo simulations for the important parameters. The uncertainty analysis is general and suitable for tackling similar problems; the Morris sensitivity analysis has been utilized to date in only a limited number of regional groundwater modeling. The simulated Q_y values were evaluated by using three kinds of calibration data (i.e., hydraulic head observations, discharge estimates, and constant‐head boundary flow estimates). The evaluation results indicate that, within the current DVRFS modeling framework, the Q_y estimate can only be reduced to about half of the original estimate without severely deteriorating the goodness‐of‐fit to the calibration data. The evaluation results also indicate that it is necessary to develop a new hydrogeological framework to produce new flow patterns in the DVRFS model. The issues of hydrogeology and boundary flow are being addressed in a new version of the DVRFS model planned for release by the U.S. Geological Survey.     

Inbank and overbank velocity conditions in an arid zone anastomosing river

Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q_2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.     

Entropy‐based equation to assess hillslope sediment production

A new equation to assess hillslope sediment production, based on physical and probabilistic approaches, is presented. The equation, which allows the computation of the delivery ratio for every event, considers the physical variables of travel distance, stream power, settling velocity and gross erosion. The probability density function that arises from the new formulation is solved using the principle of maximum entropy. Based on data from five watersheds in both tropical and temperate zones, the new delivery parameter K_v is calibrated and associated with vegetation cover and conservation practice. The proposed equation is rationally based in relation to parameter K_v. The entropy‐based equation was applied to assess sediment yield in two other experimental watersheds, showing good predictability for the set (mean absolute error of 20·8%). No systematic error was found in the analysed data. The entropy‐based equation showed good predictability for long‐term sets of data and for high‐erosivity events, but did not perform well for the low‐erosivity ones. Copyright © 2007 John Wiley & Sons, Ltd.     

Field experiments constraining the probability distribution of particle travel distances during natural rainstorms on different slope gradients

Rain splash erosion is an important soil transport mechanism on steep hillslopes. The rain splash process is highly stochastic; here we seek to constrain the probability distribution of splash transport distances on natural hillslopes as a function of hillslope gradient and total precipitation depth. Field experiments were conducted under natural precipitation events to observe splash travel on varying slope gradients. The downslope fraction of splash transport on 15°, 25° and 33° gradients were 85%, 96% and 96%, respectively. Maximum splash transport (L_max) was related to the rain splash detachment of soil particles and slope gradient. An empirical relationship of L_max to the precipitation depth and gradient was obtained; it is linearly proportional to hillslope gradient and logarithmically related to precipitation depth. Measured splash distances were calibrated to the fully two‐dimensional (2D) model of splash transport of Furbish et al. (Journal of Geophysical Research 112 : F01001, 2007) that is based on the assumption that radial splash distances are exponentially distributed; calibrated values of mean splash transport distances are an order of magnitude greater than those previously determined in a controlled laboratory setting. We also compared measured data with several one‐dimensional (1D) probability distributions to asses if splash transport distances could be better explained by a heavy‐tailed probability distribution rather than an exponential probability distribution. We find that for hillslopes of 15° and 25°, although a log‐normal probability distribution best describes the data, we find its likelihood is nearly indistinguishable from an exponential distribution based on computing maximum likelihood estimators for all 1D distributions (exponential, log‐normal and Weibull). At 33°, however, we find stronger evidence that measured travel distances are heavy‐tailed. Copyright © 2011 John Wiley & Sons, Ltd.