On the use of horizontal acoustic Doppler profilers for continuous bed shear stress monitoring

Continuous monitoring of bed shear stress in large river systems may serve to better estimate alluvial sediment transport to the coastal ocean.Here we explore the possibility of using a horizontally deployed acoustic Doppler current profiler(ADCP) to monitor bed shear stress,applying a prescribed boundary layer model,previously used for discharge estimation.The model parameters include the local roughness length and a dip correction factor to account for sidewall effects.Both these parameters depend on river stage and on the position in the cross-section, and were estimated from shipborne ADCP data.We applied the calibrated boundary layer model to obtain bed shear stress estimates over the measuring range of the HADCP.To validate the results,co-located coupled ADCPs were used to infer bed shear stress,both from Reynolds stress profiles and from mean velocity profiles. From HADCP data collected over a period of 1.5 years,a time series of width profiles of bed shear stress was obtained for a tidal reach of the Mahakam River,East Kalimantan,Indonesia.A smaller dataset covering 25 hours was used for comparison with results from the coupled ADCPs.The bed shear stress estimates derived from Reynolds stress profiles appeared to be strongly affected by local effects causing upflow and downflow,which are not included in the boundary layer model used to derive bed shear stress with the horizontal ADCP.Bed shear stresses from the coupled ADCP are representative of a much more localized flow,while those derived with the horizontal ADCP resemble the net effect of the flow over larger scales.Bed shear stresses obtained from mean velocity profiles from the coupled ADCPs show a good agreement between the two methods,and highlight the robustness of the method to uncertainty in the estimates of the roughness length.     

Comparing different methods of bed shear stress estimates in simple and complex flow fields

Bed shear stress is a fundamental variable in river studies to link ?ow conditions to sediment transport. It is, however, dif?cult to estimate this variable accurately, particularly in complex ?ow ?elds. This study compares shear stress estimated from the log pro?le, drag, Reynolds and turbulent kinetic energy (TKE) approaches in a laboratory ?ume in a simple boundary layer, over plexiglas and over sand, and in a complex ?ow ?eld around de?ectors. Results show that in a simple boundary layer, the log pro?le estimate is always the highest. Over plexiglas, the TKE estimate was the second largest with a value 30 per cent less than the log estimate. However, over sand, the TKE estimate did not show the expected increase in shear stress. In a simple boundary layer, the Reynolds shear stress seems the most appropriate method, particularly the extrapolated value at the bed obtained from a turbulent pro?le. In a complex ?ow ?eld around de?ectors, the TKE method provided the best estimate of shear stress as it is not affected by local streamline variations and it takes into account the increased streamwise turbulent ?uctuations close to the de?ectors. It is suggested that when single‐point measurements are used to estimate shear stress, the instrument should be positioned close to 0·1 of the ?ow depth, which corresponds to the peak value height in pro?les of Reynolds and TKE shear stress. Copyright © 2004 John Wiley & Sons, Ltd.     

Analytical approach for sheet flow transport in purely acceleration-skewed oscillatory flow

An instantaneous analytical approach is developed to predict sheet flow transport in purely acceleration-skewed oscillatory flow. The approach is derived from exponential approximations of velocity and concentration profiles above a mobile seabed, and it particularly considers factors of phase lead; phase lag (i.e. phase residual and phase shift); acceleration modification; and asymmetries in shear stress, roughness height, and boundary layer development. The approach can predict net boundary layer flow above a mobile seabed, and can revert to the classical bedload model. Instantaneous and net sediment transport rates are studied using the approach. The instantaneous sediment transport rate in an onshore flow stage can be approximated by a power function of velocity in which the exponent is confirmed to range between 1 and 5 with a decrease in the phase residual. The net sediment transport rate predicted using the approach is validated using a considerable amount of measured data, and compared with existing instantaneous and half-period type models that consider the phase lag or acceleration modification. For the net sediment transport rate in purely acceleration-skewed oscillatory flow, the phase residual is less important than the acceleration-skewed boundary layer difference between onshore and offshore acceleration stages.     

Simulation and measurement of surface shear stress over isolated and closely spaced transverse dunes in a wind tunnel

Topographic interactions generate multidirectional and unsteady air?ow that limits the application of velocity pro?le approaches for estimating sediment transport over dunes. Results are presented from a series of wind tunnel simulations using Irwin‐type surface‐mounted pressure sensors to measure shear stress variability directly at the surface over both isolated and closely spaced sharp‐crested model dunes. Findings complement existing theories on secondary air?ow effects on stoss transport dynamics and provide new information on the in?uence of lee‐side air?ow patterns on dune morphodynamics. For all speeds investigated, turbulent unsteadiness at the dune toe indicates a greater, more variable surface shear, despite a signi?cant drop in time‐averaged measurements of streamwise shear stress at this location. This effect is believed suf?cient to inhibit sediment deposition at the toe and may be responsible for documented intermittency in sand transport in the toe region. On the stoss slope, streamline compression and ?ow acceleration cause an increase in ?ow steadiness and shear stress to a maximum at the crest that is double that at the toe of the isolated dune and 60–70 per cent greater than at ?ow reattachment on the lower stoss of closely spaced dunes. Streamwise ?ow accelerations, rather than turbulence, have greater in?uence on stress generation on the stoss and this effect increases with stoss slope distance and with incident wind speed. Reversed ?ow within the separation cell generates signi?cant surface shear (30–40 per cent of maximum values) for both spacings. This supports ?eld studies that suggest reversed ?ow is competent enough to return sediment to the dune directly or in a de?ected direction. High variability in shear at reattachment indicates impact of a turbulent shear layer that, despite low values of time‐averaged streamwise stress in this region, would inhibit sediment accumulation. Downwind of reattachment, shear stress and ?ow steadiness increase within 6 h (h = dune height) of reattachment and approach upwind values by 25 h. A distance of at least 30 h is suggested for full boundary layer recovery, which is comparable to ?uvial estimates. The Irwin sensor used in this study provides a reliable means to measure skin friction force responsible for sand transport and its robust, simple, and cost‐effective design shows promise for validating these ?ndings in natural dune settings. Copyright © 2003 John Wiley & Sons, Ltd.     

Examining the physical components of boundary shear stress for water‐worked gravel deposits

It is argued in this commentary that, in order to understand better the physical mechanisms that generate boundary shear stress over water‐worked gravel beds, flow velocity data should be re‐evaluated by spatial averaging the Reynolds equations to produce time‐ and space‐averaged (double‐averaged) momentum equations. A series of laboratory experiments were conducted in which the flow velocities were measured using a PIV system over two water‐worked gravel deposits. Combined with detailed data on the bed surface topography and vertical porosity, the physical components of shear stress were obtained. This enabled the various momentum transfer mechanisms present above, within and at the interface of a porous, fluvial deposit, to be quantified. This included the examination of the relevant contributions of temporal and spatial fluctuations in velocity and surface drag to the overall momentum transfer. It is demonstrated that double‐averaging represents a logical framework for assessing the fluid forces responsible for sediment entrainment and for investigating intragravel flow and sediment–water interface exchange mechanisms within the roughness layer in water‐worked gravel deposits. By considering the physical components of shear stress and their relative sizes it was possible to provide a physically based explanation for existing observations of enhanced mobility of gravel–sand mixtures and the transfer of solutes into porous, gravel deposits. This analysis reveals the importance of obtaining co‐located, high quality spatial data on the flow field and bed surface topography in order to gain a physical understanding of the mechanisms which generate boundary shear stress. Copyright © 2010 John Wiley & Sons, Ltd.     

Changes in the saltation flux following a step‐change in macro‐roughness

The effect of a step change in macro‐roughness on the saltation process under sediment supply limited conditions was examined in the atmospheric boundary layer. For an array of roughness elements of roughness density λ = 0.045 (λ = total element frontal area/total surface area of the array) the horizontal saltation flux was reduced by 90% (±7%) at a distance of ≈150 roughness element heights into the array. This matches the value predicted using an empirical design model and provides confidence that it can be effectively used to engineer roughness arrays to meet sand flux reduction targets. Measurements of the saltation flux characteristics in the vertical dimension, including: saltation layer decay (e‐folding) height and particle size, revealed that with increasing distance into the array, the rate of mass flux change with increasing height decreased notably, and (geometric) mean particle diameter decreased. The distribution of the saltation mass flux in the vertical remains exponential in form with increasing distance into the roughness array, and the e‐folding height increases as well as increasing at a greater rate as particle diameter diminishes. The increase in e‐folding height suggests the height of saltating particles is increasing along with their mean speed. This apparent increase in mean speed is likely due to the preferential removal, or sequestration, of the slower moving particles across the size spectrum, as they travel through the roughness array. Copyright © 2018 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.     

Near‐bed shear stress,turbulence production and dissipation in a shallow and narrow tidal channel

Near‐bed, highly resolved velocity profiles were measured in the lower 0.03 m of the water column using acoustic Doppler profiling velocimeters in narrow tidal channels in a salt marsh. The bed shear stress was estimated from the velocity profiles using three methods: the log‐law, Reynolds stress, and shear stress derived from the turbulent kinetic energy (TKE). Bed shear stresses were largest during ebbing tide, while near‐bed velocities were larger during flooding tide. The Reynolds stress and TKE method gave similar results, while the log‐law method resulted in smaller bed shear stress values during ebbing tide. Shear stresses and turbulent kinetic energy followed a similar trend with the largest peaks during ebbing tide. The maximum turbulent kinetic energy was on the order of 1 × 10^{? 2} m²/s². The fluid shear stress during flooding tide was approximately 30% of the fluid shear stress during ebbing tide. The maximum TKE‐derived shear stress was 0.7 N/m² and 2.7 N/m² during flooding and ebbing tide, respectively, and occurred around 0.02 m above the bed. Turbulence dissipation was estimated using the frequency spectrum and structure function methods. Turbulence dissipation estimates from both methods were maximum near the bed (~0.01 m). Both the structure function and the frequency spectrum methods resulted in maximum dissipation estimates on the order of 4 × 10^{? 3} m²/s³. Turbulence production exceeded turbulence dissipation at every phase of the tide, suggesting that advection and vertical diffusion are not negligible. However, turbulence production and dissipation were within a factor of 2 for 77% of the estimates. The turbulence production and dissipation decreased quickly away from the bed, suggesting that measurements higher in the water column cannot be translated directly to turbulence production and dissipation estimates near the bed. Copyright © 2015 John Wiley & Sons, Ltd.     

General velocity formula of boundary layer above mobile sediment bed induced by asymmetric waves

Sediment movement in the wave boundary layer above a mobile sediment bed is complex.A velocity formula for the boundary layer is proposed for sheet flow induced by asymmetric waves above a mobile sediment bed.The formula consists of a free stream velocity and a defect function which contains a phase-lead,boundary layer thickness and mobile sediment bed.Phase-lag of sediment movement is considered in the formula for the mobile sediment bed.The formula needs six dependent variables about asymmetric wave and sediment characteristics.Asymmetry effects on parameters(orbital amplitude,roughness height,bed shear stress,and boundary layer thickness)are properly considered such that the formula can yield velocity differences among onshore,offshore,acceleration,and deceleration stages.The formula estimates the net boundary layer velocity resulting from the mobile sediment bed and asymmetric boundary layer thickness.In addition,a non-constant phase-lead also contributes to the net boundary layer velocity in asymmetric oscillatory sheet flow.Results of the formula are as good as that of a two-phase numerical model.Sheet flow transport induced by asymmetric waves,and the offshore net sediment transport rate with a large phase-lag under velocity-skewed waves,can be adequately estimated by the formula with a power sediment concentration function.     

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.     

Experimental study of extreme shear stress for shallow flow under simulated rainfall

The raindrop impact and overland flow are two major factors causing soil detachment and particle transportation. In this study, the turbulent characteristics of the shallow rain‐impacted water flow were investigated using a 2‐D fibre‐optic laser Doppler velocimetry (FLDV) and an artificial rainfall simulator. The fluctuating turbulent shear stress was computed using digital data processing techniques. The experimental data showed that the Reynolds shear stress follows a probability distribution with heavy tails. The tail probability increases with an increase of rainfall intensity or raindrop diameter, and it decreases with an increase of Reynolds number. A modified empirical equation was derived using both the raindrop diameter and rainfall intensity as independent variables to provide a better prediction of the Darcy‐Weisbach friction coefficient f under rainfall conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

A field study of mean and turbulent flow characteristics upwind,over and downwind of barchan dunes

Field‐measured patterns of mean velocity and turbulent airflow are reported for isolated barchan dunes. Turbulence was sampled using a high frequency sonic anemometer, deriving near‐surface Reynolds shear and normal stresses. Measurements upwind of and over a crest‐brink separated barchan indicated that shear stress was sustained despite a velocity reduction at the dune toe. The mapped streamline angles and enhanced turbulent intensities suggest the effects of positive streamline curvature are responsible for this maintenance of shear stress. This field evidence supports an existing model for dune morphodynamics based on wind tunnel turbulence measurements. Downwind, the effect of different dune profiles on flow re‐attachment and recovery was apparent. With transverse incident flow, a re‐attachment length between 2·3 and 5·0h (h is dune brink height) existed for a crest‐brink separated dune and 6·5 to 8·6h for a crest‐brink coincident dune. The lee side shear layer produced elevated turbulent stresses immediately downwind of both dunes, and a decrease in turbulence with distance characterized flow recovery. Recovery of mean velocity for the crest‐brink separated dune occurred over a distance 6·5h shorter than the crest‐brink coincident form. As the application of sonic anemometers in aeolian geomorphology is relatively new, there is debate concerning the suitability of processing their data in relation to dune surface and streamline angle. This paper demonstrates the effect on Reynolds stresses of mathematically correcting data to the local streamline over varying dune slope. Where the streamline angle was closely related to the surface (windward slope), time‐averaged shear stress agreed best with previous wind tunnel findings when data were rotated along streamlines. In the close lee, however, the angle of downwardly projected (separated) flow was not aligned with the flat ground surface. Here, shear stress appeared to be underestimated by streamline correction, and corrected shear stress values were less than half of those uncorrected. Copyright © 2011 John Wiley & Sons, Ltd.     

Retrieval of the boundary layer height from active and passive remote sensors. Comparison with a NWP model

In this study, we used boundary layer heights derived from lidar in Romania to validate the Weather Research Forecast (WRF) model improved by ARIA Technologies SA in the framework of ROMAIR LIFE project. Lidar retrievals were also compared to the retrievals from meteorological data, both modeled (Global Data Assimilation System; GDAS) and measured (microwave radiometry). Both the gradient and the wavelet covariance methods were used to compute the boundary layer height (BLH) from the range corrected lidar signal, and their equivalence was shown. The analysis was performed on 102 datasets, spread over all seasons and 3 years (2009–2011). A good agreement was found for the remote sensors (lidar and microwave radiometer) which are co-located and measure simultaneously. The correlation of the measured boundary layer height and the modelled one was 0.66 for the entire dataset, and 0.73 when considering daytime data, i.e., for a well defined boundary layer. A systematic underestimation of the boundary layer height by the WRF during non-convective periods (nocturne, stable atmosphere) was found.     

EXPERIMENTAL STUDY ON BED SHEAR STRESS AND SEDIMENT TRANSPORT IN A MEANDERING LABORATORY CHANNEL

1 INTRODUCTION River erosion is a complex phenomenon. The rate of bank retreat is determined by flow, bed topography, sediment transport, bank properties, and water quality. Prediction of future river planform changes and the knowledge of river erosion and river meandering are required for land use planning in alluvial river valleys and determining locations for bridges and hydraulic structures. The control of riverbank erosion requires prediction of flow and bed features in a meanderin…     

Turbulent flow structures in alluvial channels with curved cross‐sections under conditions of downward seepage

Experimental investigations have been done to analyze turbulent structures in curved sand bed channels with and without seepage. Measures of turbulent statistics such as time‐averaged near‐bed velocities, Reynolds stresses, thickness of roughness sublayer and shear velocities were found to increase with application of downward seepage. Turbulent kinetic energy and Reynolds normal stresses are increased in the streamwise direction under the action of downward seepage, causing bed particles to move rapidly. Analysis of bursting events shows that the relative contributions of all events (ejections, sweeps and interactions) increase throughout the boundary layer, and the thickness of the zone of dominance of sweep events, which are responsible for the bed material movement, increases in the case of downward seepage. The increased sediment transport rate due to downward seepage deforms the cross‐sectional geometry of the channel made of erodible boundaries, which is caused by an increase in flow turbulence and an associated decrease in turbulent kinetic energy dissipation and turbulent diffusion.     

The causes of bedload pulses in a gravel channel: the implications of bedload grain‐size distributions

Unsteady bedload transport was measured in two c. 5 m wide anabranches of a gravel‐bed braided stream draining the Haut Glacier d'Arolla, Switzerland, during the 1998 and 1999 melt seasons. Bedload was directly sampled using 152 mm square Helley–Smith type samplers deployed from a portable measuring bridge, and independent transport rate estimates for the coarser size fractions were obtained from the dispersion of magnetically tagged tracer pebbles. Bedload transport time series show pulsing behaviour under both marginal (1998) and partial (1999) transport regimes. There are generally weak correlations between transport rates and shear stresses determined from velocity data recorded at the measuring bridge. Characteristic parameters of the bedload grain‐size distributions (D₅₀, D₈₄) are weakly correlated with transport rates. Analysis of full bedload grain‐size distributions reveals greater structure, with a tendency for transport to become less size selective at higher transport rates. The bedload time series show autoregressive behaviour but are dif?cult to distinguish by this method. State–space plots, and associated measures of time‐series separation, reveal the structure of the time series more clearly. The measured pulses have distinctly different time‐series characteristics from those modelled using a one‐dimensional sediment routing model in which bed shear stress and grain size are varied randomly. These results suggest a mechanism of pulse generation based on irregular low‐amplitude bedforms, that may be generated in‐channel or may represent the advection of material supplied by bank erosion events. Copyright © 2003 John Wiley & Sons, Ltd.     

Analysis of bedload transport in laminar flows

Bedload transport generally depends on the bed shear stress and Reynolds number. Many studies conducted for the condition of turbulent flows have revealed the dependence of the transport rate on the bed shear stress, while knowledge of the Reynolds number effect on the transport rate is very limited. As an extreme case to reflect the viscous effect on sediment transport, sediment transport in laminar flows is considered in this paper. A stochastic approach is adopted to explore how the transport rate can be associated with characteristics of laminar flows. First, the probability of erosion in the absence of turbulence is assumed to depend only on the randomness of bed particles. The probability is then applied to formulate the sediment transport rate, of which the derivation is made largely based on Einstein’s bedload theory. The theoretical result indicates that the dimensionless transport rate for laminar flows is dependent on the dimensionless shear stress and dimensionless particle diameter or the shear Reynolds number. Comparisons are finally made between the derived formula and an empirical correlation available in the literature.     

Analytical model of the mean velocity distribution in an open channel with double-layered rigid vegetation

An analytical model for predicting the vertical distribution of mean streamwise velocity in an open channel with double-layered rigid vegetation is proposed. The double-layered model was constructed in a laboratory flume with an array of steel cylinders of two heights. For each vegetation layer (i.e., the short- or tall-vegetation layer), the flow is vertically separated into a lower vegetation zone and an upper vegetation zone, and corresponding momentum equations for each zone are formulated. For the lower vegetation zone, a uniform velocity was adopted since turbulent shear is relatively small and the Reynolds stress is ignored. For the upper vegetation zone, a power series was used to solve the momentum equations. For the free-water zone, a new expression was suggested to obtain a zero velocity gradient at the water surface instead of the traditional logarithmic velocity distribution. Good agreement between the analytical predictions and experimental data demonstrated the validity of the model.     

Comparison of different mass transport calculation methods for wind erosion quantification purposes

Quantitative estimation of the material transported by the wind under field conditions is essential for the study and control of wind erosion. A critical step of this calculation is the integration of the curve that relates the variation of the amount of the material carried by the wind with height. Several mathematical procedures have been proposed for this calculation, but results are scarce and controversial. One objective of this study was to assess the efficiency of three mathematical models (a rational, an exponential, and a simplified Gaussian function) for the calculation of the mass transport, as compared to the linear spline interpolation. Another objective of this study was to compare the mass transport calculated from field measurements obtained from a minimum of three discrete sampling heights with measurements of nine sampling heights. With this purpose, wind erosion was measured under low surface roughness conditions on an Entic Haplustoll during 25 events. The rational function was found to be mathematically limited for the estimation of wind eroded sediment mass flux. The simplified Gaussian model did not fit to the vertical mass flux profile data. Linear spline interpolation generally produced higher mass transport estimates than the exponential equation, and it proved to be a very flexible and robust method. Using different sampling arrangements and different mass flux models can produce differences of more than 45% in mass transport estimates, even under similar field conditions. Under the conditions of this study, at least three points between the soil surface and 1·5 m high, including one point as closest as possible to the surface, should be sampled in order to obtain accurate mass transport estimates. Additionally, the linear spline interpolation and the non‐linear regression using an exponential model, proved to be mathematically reliable methods for calculating the mass transport. Copyright © 2010 John Wiley & Sons, Ltd.     

Flow competence: A criticism of a classic concept

The largest grains found in samples of transported sediment are commonly used to estimate flow competence. With samples from a range of flows, a relationship between the flow and the largest mobile grain can be derived and used to estimate the critical shear stress for incipient motion of the different grain sizes in the bed sediment or, inversely, to estimate the magnitude of the flow from the largest grain found in a transport sample. Because these estimates are based on an extreme value of the transport grain-size distribution, however, they are subject to large errors and are sensitive to the effect of sample size, which tends to vary widely in sediment transport samples from natural flows. Furthermore, estimates of the critical shear stress based on the largest sampled moving grain cannot be scaled in a manner that permits reasonable comparison between fractions. The degree to which sample size and scaling problems make largest-grain estimates of fractional critical shear stress deviate from a true relationship cannot be predicted exactly, although the direction of such a deviation can be demonstrated. The large errors and unknown bias suggest that the largest sampled mobile grain is not a reliable predictor of either critical shear stress or flow magnitude. It is possible to define a single flow competence for the entire mixture, based on a central value of the transport grain-size distribution. Such a measure is relatively stable, does not require between-fraction scaling, and appears to be well supported by observation.