A probabilistic model for sediment entrainment:The role of bed irregularity

A generalized probabilistic model is developed in this study to predict sediment entrainment under the incipient motion, rolling, and pickup modes. A novelty of the proposed model is that it incorporates in its formulation the probability density function of the bed shear stress, instead of the near-bed velocity fluctuations, to account for the effects of both flow turbulence and bed surface irregularity on sediment entrainment. The proposed model incorporates in its formulation the collective effects of three para-meters describing bed surface irregularity, namely the relative roughness, the volumetric fraction and relative position of sediment particles within the active layer. Another key feature of the model is that it provides a criterion for estimating the lift and drag coefficients jointly based on the recognition that lift and drag forces acting on sediment particles are interdependent and vary with particle protrusion and packing density. The model was validated using laboratory data of both fine and coarse sediment and was compared with previously published models. The study results show that all the examined models perform adequately for the fine sediment data, where the sediment particles have more uniform gra-dation and relative roughness is not a factor. The proposed model was particularly suited for the coarse sediment data, where the increased bed irregularity was captured by the new parameters introduced in the model formulation. As a result, the proposed model yielded smaller prediction errors and physically acceptable values for the lift coefficient compared to the other models in case of the coarse sediment data.     

Spatial variations in surface sediment structure in riffle–pool sequences: a preliminary test of the Differential Sediment Entrainment Hypothesis (DSEH)

Riffle–pool sequences are maintained through the preferential entrainment of sediment grains from pools rather than riffles. This preferential entrainment has been attributed to a reversal in the magnitude of velocity and shear stress under high flows; however the Differential Sediment Entrainment Hypothesis (DSEH) postulates that differential entrainment can instead result from spatial sedimentological contrasts. Here we use a novel suite of in situ grain‐scale field measurements from a riffle–pool sequence to parameterize a physically‐based model of grain entrainment. Field measurements include pivoting angles, lift forces and high resolution digital elevation models (DEMs) acquired using terrestrial laser scanning, from which particle exposure, protrusion and surface roughness were derived. The entrainment model results show that grains in pools have a lower critical entrainment shear stress than grains in either pool exits or riffles. This is because pool grains have looser packing, hence greater exposure and lower pivoting angles. Conversely, riffle and pool exit grains have denser packing, lower exposure and higher pivoting angles. A cohesive matrix further stabilizes pool exit grains. The resulting predictions of critical entrainment shear stress for grains in different subunits are compared with spatial patterns of bed shear stress derived from a two‐dimensional computational fluid dynamics (CFD) model of the reach. The CFD model predicts that, under bankfull conditions, pools experience lower shear stresses than riffles and pool exits. However, the difference in sediment entrainment shear stress is sufficiently large that sediment in pools is still more likely to be entrained than sediment in pool exits or riffles, resulting in differential entrainment under bankfull flows. Significantly, this differential entrainment does not require a reversal in flow velocities or shear stress, suggesting that sedimentological contrasts alone may be sufficient for the maintenance of riffle–pool sequences. This finding has implications for the prediction of sediment transport and the morphological evolution of gravel‐bed rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of sand content on initial gravel motion in gravel‐bed rivers

When fine sediments are present in gravel streambeds (gravel‐framework beds), the gravel can be more easily removed from its original position, compared with gravel in a streambed without fine sediment but otherwise under the same hydraulic conditions. In this study, the effect of the presence of sand on the initiation of gravel motion in gravel riverbeds was investigated using flume experiments. The relationship between the critical Shields stress for gravel motion initiation and the fraction of sand in the bed was determined experimentally. The results can be summarized as follows. (1) When the fraction of sand in the bed is smaller than about 0.4, the critical Shields stress for the initiation of gravel motion decreases with increasing fraction of sand. The critical Shields stress increases, however, with increasing fraction of sand when it is larger than about 0.4. (2) The difference between the value of the critical Shields stress predicted by the Egiazaroff equation and the value obtained from the experimental data becomes maximum at about 0.4 of the fraction of sand. Here an empirical relation between the critical Shields stress and the fraction of sand is proposed so as to consider the effects of the ratio of the characteristic gravel size to the mean size of the bed material on the critical Shields stress. (3) Gravel in armored beds can be more easily mobilized by supplying sand as part of a sediment augmentation scheme. The sand fraction in the subsurface layer of the bed appears to reduce the friction angle of exposed particles. Sediment augmentation using sand has been recently demonstrated to be a viable alternative for mobilizing gravel for the restoration of gravel‐bed rivers downstream of dams. The quantitative evaluation obtained through the experiments reported here may be useful for the design of augmentation schemes. Copyright © 2017 John Wiley & Sons, Ltd.     

One-dimensional morphodynamic model for retrogressive erosion based on a sediment entrainment theory at high flow velocity

Retrogressive erosion is a high-speed erosion process that usually occurs during the rapid release of stored water in reservoirs built on sandy rivers.Retrogressive erosion has been utilized in the practice of reservoir sedimentation control,but accurate prediction of the bed deformation process by numerical models has rarely been reported.The current study presents a one-dimensional morphodynamic model for simulating the evolution process of retrogressive erosion induced by high-velocity flows on steep slopes.The governing equations apply a Cartesian coordinate system with a vertically oriented z axis.The bed surface gradient and friction terms in the flow equations include correction factors to take account of the effects of high slope on flow movement.The net vertical sediment flux term in the sediment transport and bed deformation equations is calculated using an equation of erosion velocity.Particularly,this equation is based on an empirical relation between the sediment entrainment rate and the Shields parameter in contrast to the traditional sediment transport capacity,and the critical Shields parameter is modified by taking into account the permeability of the sediment layer and the stability of particles on a slope.The feedback of scoured sediment on the flow movement is considered by additional terms in the governing equations.Flume experiments of retrogressive erosion in literature were simulated to validate the model.The temporal variations of the longitudinal profiles of the free surface and channel bed and the sediment transport rate were well predicted.The algorithm calculating sediment entrainment in the proposed model also was validated for an experiment measuring entrainment rate from the literature.More importantly,it was found that the morphodynamic model using the sediment transport capacity equation predicts the trend of cumulative erosion contrary to the measurements,while results of the proposed model can follow a similar trend with the observed data in the retrogressive erosion process.     

Evaluating overland flow sediment transport capacity

An equation for evaluating the sediment transport capacity of overland flow is a necessary part of a physically based soil erosion model describing sediment detachment and transport as distributed processes. At first, for the hydraulic conditions of small-scale and large-scale roughness, the sediment transport capacity relationship used in the WEPP model is calibrated by Yalin and Govers' equation. The analysis shows that the transport coefficient K_t depends on the Shields parameter, Y, according to a semi-logarithmic (Yalin) or a linear (Govers) equation. The reliability of the semi-logarithmic equation is verified by Smart's, and Aziz and Scott's experimental data. Then the Low's formula, whose applicability is also proved by Smart's, and Aziz and Scott's data, is transformed as a stream power equation in which a stream power coefficient, K_SP, depending on Shields parameter, slope, sediment and water-specific weight, appears. A relationship between transport capacity and effective stream power is also proposed. Finally, the influence of rainfall on sediment transport capacity and the prediction of critical shear stress corresponding to overland flow are examined. © 1998 John Wiley & Sons, Ltd.     

Uniform and graded bed-load sediment transport in a degrading channel with non-equilibrium conditions

Bed-load transport plays a critical role in river morphological change and has an important impact on river ecology.Although there is good understanding of the role of the variation of river bed grain size on transport dynamics in equilibrium conditions,much less is understood for non-equilibrium conditions when the channel is either aggrading or degrading.In particular,the relative role of different grain sizes in the promotion and hindering of the transport of coarse and fine fractions in a degrading channel has yet to be investigated.The current study attempts to provide new understanding through a series of flume experiments done using uniform and graded sediment particles.The experiments revealed coarser grain-size fractions for a poorly-sorted sediment,relative to uniform-sized sediment,reduced the transport of finer grains and finer fractions enhanced the transport of coarse grains.This hinderingpromotion effect,caused by relative hiding and exposure of finer and coarse fractions,increased with bed slope and decreased with relative submergence.In particular,as relative submergence increased,the graded fractions tended towards behaving more like their unifo rm-sized counterparts.Also,the bed-load parameter of the graded fractions increased more with a rise in bed slope than observed for the uniformsized counterparts.These results revealed,for degrading channel conditions,such as downstream of a dam,bed-load equations developed for uniform bed sediment are inappropriate for use in natural river systems,particularly in mountain streams.Furthermore,changes in river bed composition due to activities that enhance the input of hill-slope sediment,such as fire,logging,and agricultural development,are likely to cause significant changes in river morphology.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

Formation and erosion of sediment cover in an experimental bedrock‐alluvial channel

Sediment grains in a bedrock‐alluvial river will be deposited within or adjacent to a sediment patch, or as isolated grains on the bedrock surface. Previous analysis of grain geometry has demonstrated that these arrangements produce significant differences in grain entrainment shear stress. However, this analysis neglected potential interactions between the sediment patches, local hydraulics and grain entrainment. We present a series of flume experiments that measure the influence of sediment patches on grain entrainment. The flume had a planar bed with roughness that was much smaller than the diameters of the mobile grains. In each experiment sediment was added either as individual grains or as a single sediment pulse. Flow was then increased until the sediment was entrained. Analysis of the experiments demonstrates that: (1) for individual grains, coarse grains are entrained at a higher discharge than fine grains; (2) once sediment patches are present, the different in entrainment discharge between coarse and fine grains is greatly reduced; (3) the sheltering effect of patches also increases the entrainment discharge of isolated grains; (4) entire sediment patches break‐up and are eroded quickly, rather than through progressive grain‐by‐grain erosion; (5) as discharge increases there is some tendency for patches to become more elongate and flow‐aligned, and more randomly distributed across the bed. One implication of this research is that the critical shear stress in bedrock‐alluvial channels will be a function of the extent of the sediment cover. Another is that the influence of sediment patches equalizes critical shear stresses between different grain sizes and grain locations, meaning that these factors may not need to be accounted for. Further research is needed to quantify interactions between sediment patches, grain entrainment and local hydraulics on rougher bedrock surfaces, and under different types of sediment supply. Copyright © 2016 John Wiley & Sons, Ltd.     

Evaluating spatial variation of suspended sediment rating curves in the middle Yellow River basin,China

Abstract The suspended sediment load in the middle Yellow River basin (YRB) cannot be well predicted by capacity‐based transport formulas because a large fraction of suspended sediment load is composed of wash load. This study evaluated the spatial variations of sediment rating curves (SRCs) in the middle YRB. Both power and linear SRCs were used to fit daily flow and suspended sediment concentration (SSC) historical data at 49 gauging stations throughout the middle YRB. The spatial variation in regression coefficients was investigated, and the relationship between regression coefficients and the physical characteristics of watersheds was discussed. The results indicate that SRC regression coefficients vary with drainage area and basin slope, but their responses to these parameters are remarkably different in watersheds with different underlying surfaces, which indicates the significance of sediment availability, erodibility, and grain size distribution. For power SRCs representing sediment transport in unsaturated flows, the regression coefficients are more closely correlated with the drainage area in loess regions and with the basin slope in rock mountain regions. For linear SRCs representing sediment transport in saturated flows, saturated SSCs vary with coarse (particle size > 0.05 mm) and fine (particle size < 0.01 mm) fractions in suspended sediment. The maximum saturated SSC among the different gauging stations is associated with the optimal grain size composition of suspended sediment, which has been proposed for loess regions in previous studies. This study provides theoretical support for estimating the regression parameters for sediment transport modelling, especially in ungauged basins.     

Stream geomorphology in a mountain lake district: hydraulic geometry,sediment sources and sinks,and downstream lake effects

Lakes are common in glaciated mountain regions and geomorphic principles suggest that lake modifications to water and sediment fluxes should affect downstream channels. Lakes in the Sawtooth Mountains, Idaho, USA, were created during glaciation and we sought to understand how and to what extent glacial morphology and lake disruption of fluxes control stream physical form and functions. First, we described downstream patterns in channel form including analyses of sediment entrainment and hydraulic geometry in one catchment with a lake. To expand on these observations and understand the role of glacial legacy, we collected data from 33 stream reaches throughout the region to compare channel form and functions among catchments with lakes, meadows (filled lakes), and no past or present lakes. Downstream hydraulic geometry relationships were weak for both the single catchment and regionally. Our data show that downstream patterns in sediment size, channel shape, sediment entrainment and channel hydraulic adjustment are explained by locations of sediment sources (hillslopes and tributaries) and sediment sinks (lakes). Stream reaches throughout the region are best differentiated by landscape position relative to lakes and meadows according to channel shape and sediment size, where outlets are wide and shallow with coarse sediment, and inlets are narrow and deep with finer sediment. Meadow outlets and lake outlets show similarities in the coarse‐sediment fraction and channel capacity, but meadow outlets have a smaller fine‐sediment fraction and nearly mobile sediment. Estimates of downstream recovery from lake effects on streams suggest 50 per cent recovery within 2–4 km downstream, but full recovery may not be reached within 20 km downstream. These results suggest that sediment sinks, such as lakes, in addition to sources, such as tributaries, are important local controls on mountain drainage networks. Copyright © 2006 John Wiley & Sons, Ltd.     

EXPERIMENTAL INVESTIGATION AND ANALYSIS OF HEC-6 RIVER MORPHOLOGICAL MODEL

Only comparatively few experimental studies have been carried out to investigate the performance of the HEC-6 river morphological model. The model was developed by the Hydrologic Engineering Center of the US Army Corps of Engineers. In this study, experiments were carried out in a 20 m long concrete flume 0.6 m wide with varying rectangular cross-sections. The channel bed is paved with uniform sand of D₅₀ = 0.9 mm and D₉₀ = 1.2 mm within the test reach of 12 m. Two types of experiments were carried out with sediment transport, one under steady uniform flow and another under steady non-uniform flow conditions. Nine steady uniform flow experiments were carried out to compare the measured equilibrium relationship of flow and sediment transport rate with two bedload formulae, namely, Du Boys and Meyer–Peter and Muller, and with three total load formulae, namely, Toffaleti, Laursen and Yang. It was found that even though the sediment transport consists of a certain portion of bedload, the total load formulae give satisfactory results and better agreement than the two bedload formulae. Five steady non-uniform flow experiments were carried out under various conditions of varying bed profile and channel width and also with sediment addition and withdrawal. The measured transient water surface and bed profiles are compared with the computed results from the HEC-6 model. It was found that the Toffaleti and Yang total load formulae used in the HEC-6 model give the most satisfactory prediction of actual bed profiles under various conditions of non-uniform flow and sediment transport. The effects of Manning's n, variations of sediment inflow, various sediment transport formulae, sediment grain size and the model numerical parameters, i.e. distance interval Δx and numerical weighting factor, on the computed water surface and bed profiles were determined. It was found that the selection of the sediment transport formulae has the most significant effect on the computed results. It can be concluded that the HEC-6 model can predict satisfactorily a long-term average pattern of local scour and deposition along a channel with either a small abrupt change in geometry or gradually varying cross-sections. However, the accuracy of the model prediction is reduced in the regions where highly non-uniform flow occurs.     

Numerical modelling of alternate bar formation,development and sediment sorting in straight channels

A two‐dimensional shallow water hydro‐sediment‐morphodynamic model is applied to investigate alternate bar formation, development and sediment sorting in straight channels. The model is coupled, explicitly incorporating the flow–sediment–bed interactions by using the full mass and momentum conservation equations, which are numerically solved by a well‐balanced version of the finite volume Slope Limiter Centred (SLIC) scheme. The model is first tested against a flume experiment on alternate bars formed over a uniform sediment bed, which clearly exhibits processes of bar formation, migrating and finally approaching an equilibrium state. Then it is applied to another flume experiment on alternate bars due to non‐uniform sediment transport. The computational results are evaluated, with a focus on the longitudinal and vertical sediment sorting. It is argued for the first time that the inconsistent sediment sorting patterns observed in previous studies are determined by different sediment transport conditions, i.e. full versus partial transport. When a condition of full transport is achieved, under which all size fractions are fully mobilized and transported, the longitudinal surface sediment shows a sorting pattern of coarse‐on‐head and fine‐in‐pool, and the vertical substrate sediment exhibits an immobile‐fine‐coarse structure upwards. In contrast, for a partial transport condition, under which only finer fraction participates in the transport process, an opposite longitudinal pattern (i.e. fine‐on‐head and coarse‐in‐pool) and a different vertical structure (i.e. immobile‐coarse‐fine) are observed. Concurrently, numerical experiments with specified conditions show that the critical aspect ratio for the formation of migrating alternate bars is approximately equal to 12. With the increase of the aspect ratio, the bar length grows gradually, while the bar height increases rapidly for moderate values of the aspect ratio and then keeps nearly stable. The bar celerity, however, is weakly sensitive to the variation of this ratio. Copyright © 2016 John Wiley & Sons, Ltd.     

Evolution of Heterogeneous Mixed Siliciclastic/Carbonate Aquifers Containing Metastable Sediments

Mixed carbonate and siliciclastic marine sediments commonly become freshwater aquifers in eastern coastal regions of the United States and many other global locations. As these deposits age, the carbonate fraction of the sediment is commonly removed by dissolution and the aquifer can become a solely siliciclastic system or contain zones or beds of pure quartz sand. During aquifer evolution, the sediment grain size characteristics, hydraulic conductivity, and porosity change. An investigation of these changes using mixed carbonate/siliciclastic sediment samples collected from a modern barrier island beach in southern Florida showed that the average mean grain diameter decreased with removal of the carbonate fraction, but the average hydraulic conductivity and porosity increased slightly, but not to statistical significance. This counterintuitive result occurs because of the change in the pore types from a combined shelter and intergranular pore system producing a dual porosity system in the mixed sediments to a single intergranular pore system in the siliciclastic sediment fraction. Within the mixed carbonate/siliciclastic sediment, in the pure carbonate fraction, large shell fractions form grain‐supported large pores, which become filled with sand‐sized quartz as the shell fragments decrease in size or as the sediment becomes compacted. The hydraulic conductivity increases because the shell fragments that were oriented perpendicular to flow caused an increase in the length of the flow path, or a larger scale tortuosity, compared with the flow through pure quartz sand.     

Effects of sediment grain size and channel slope on the stability of river bifurcations

Channel bifurcations can be found in river network systems from high gradient gravel-bed rivers to fine-grained low gradient deltas. In these systems, bifurcations often evolve asymmetrically such that one downstream channel silts up and the other deepens and, in most cases, they eventually avulse. Past analytical and numerical studies showed that symmetric bifurcations are unstable in high and low Shields stress conditions resulting in asymmetric bifurcations and avulsion, while they can be stable in the mid-Shields range, but this range is smaller for larger width-to-depth ratio. Here, using a one-dimensional (1D) numerical model, we show that effects of sediment grain size and of channel slope are much larger than expected for low-gradient systems when a sediment transport relation is used that separates between bedload and suspended load transport. We found that the range of Shields stress conditions with unstable symmetric bifurcations expanded for lower channel slopes and for finer sediment. In high sediment mobility, suspended load increasingly dominates the sediment transport, which increases the sediment transport nonlinearity and lowers the relative influence of the stabilizing transverse bedslope-driven flux. Contrary to previous works, we found another stable symmetric solution in high Shields stress, but this only occurs in the systems with small width-to-depth ratio. This indicates that suspended load-dominated bifurcations of lowland rivers are more likely to develop into highly asymmetric channels than previously thought. This explains the tendency of channel avulsion observed in many systems.     

Investigating the extrapolation performance of neural network models in suspended sediment data

Suspended sediment modelling is a quite significant issue in hydrology. The prediction of suspended sediment has taken the attention of several scientists in water resources. With extrapolation, the forecasting ability of the employed forecasting model beyond the calibration range is investigated. In the present study, different smoothing parameters are used to differentiate the kurtosis of the local critical points (local minima and maxima). The two models used are an artificial neural network (ANN) model and a multiple linear regression (MLR) model for prediction in order to examine the model extrapolation ability. The ANN model provides closer estimations to the observed peaks, being higher than the corresponding MLR ones. For the local minima, the ANN predictions are higher than the MLR predictions. As there are limited local points, all the remaining ANN predictions are lower than the MLR ones except for one point.     

Implications of the saltation–abrasion bedrock incision model for steady‐state river longitudinal profile relief and concavity

The saltation–abrasion model predicts rates of river incision into bedrock as an explicit function of sediment supply, grain size, boundary shear stress and rock strength. Here we use this experimentally calibrated model to explore the controls on river longitudinal profile concavity and relief for the simple but illustrative case of steady‐state topography. Over a wide range of rock uplift rates we find a characteristic downstream trend, in which upstream reaches are close to the threshold of sediment motion with large extents of bedrock exposure in the channel bed, while downstream reaches have higher excess shear stresses and lesser extents of bedrock exposure. Profile concavity is most sensitive to spatial gradients in runoff and the rate of downstream sediment fining. Concavity is also sensitive to the supply rate of coarse sediment, which varies with rock uplift rate and with the fraction of the total sediment load in the bedload size class. Variations in rock strength have little influence on profile concavity. Profile relief is most sensitive to grain size and amount of runoff. Rock uplift rate and rock strength influence relief most strongly for high rates of rock uplift. Analysis of potential covariation of grain size with rock uplift rate and rock strength suggests that the influence of these variables on profile form could occur in large part through their influence on grain size. Similarly, covariation between grain size and the fraction of sediment load in the bedload size class provides another indirect avenue for rock uplift and strength to influence profile form. Copyright © 2008 John Wiley & Sons, Ltd.     

Entrainment,displacement and transport of tracer gravels

The mass and size distribution of grain entrainment per unit bed area may be measured by replacing a volume of the bed with tracer gravels and observing the mass difference before and after a transport event. This measure of spatial entrainment is relevant to any process involving size-selective exchange of sediment between transport and bed and may be directly used in calculations of sediment transport rate using an elementary relation for fractional transport components presented here. This relation provides a basis for evaluating tracer data collected by different methods and may be used to provide physical insight regarding the expected behaviour of tracer grains. The variation with grain size of total displacement length L_ti depends on the degree of mobilization of the individual fractions on the bed surface: L_ti is independent of D_i for smaller, fully mobile sizes and decreases rapidly with D_i for larger fractions in a state of partial transport (in which a portion of the surface grains remain immobile through the flow event). The boundary between fully and partially mobile grain sizes increases with flow strength. These inferences are supported by values of L_ti calculated from flume experiments and provide a physical explanation for a summary relation between L_ti and D_i based on field data. © 1997 John Wiley & Sons, Ltd.     

Patterns of grain‐size dependent sediment transport in low‐ordered,ephemeral channels

Sediment data were analyzed to determine grain‐size dependant factors affecting sediment transport in a low‐ordered, ephemeral watershed. Sediment and flow samples were collected during 22 flow events at the outlet of a 4·53 ha sub‐watershed within the Walnut Gulch Experimental Watershed in south‐eastern Arizona. Measured concentrations ranged from 4191 to 115 045 mg l^?1 and included grain sizes up to 8·0 mm in diameter. Two grain‐size dependent transport patterns were observed, that of the finer grain‐size fraction (approximately < 0·25 mm) and that of a coarser grain‐size fraction (approximately ≥ 0·25 mm). The concentration of the fine fraction decreased with flow duration, peaking near the beginning of a flow event and declining thereafter. The concentration of the fine fraction showed slight trends with season and recovery period. The concentration of the coarse fraction displayed a slight negative trend with instantaneous discharge and was not correlated with event duration. These patterns typically produced a condition where the majority of the fine fraction of the sediment yield was evacuated out of the watershed before the hydrograph peak while the majority of the coarser sediment was evacuated during the falling limb of the hydrograph. Each grain‐size dependent transport pattern was likely influenced by the source of the associated sediment. At the flow event time scale, the fines were primarily wash load, supplied from the hillslopes and the coarser grains were entrained from the channel bed. Because transport patterns differ based on grain size, attempts to define the total sediment concentration and sediment yield by the behavior of a single grain‐size fraction may lead to erroneous results, especially when a large range of sediment grain sizes are present. Copyright © 2010 John Wiley & Sons, Ltd.