Shields Parameter Sediment Transport Calculator

Introduction

One of the most practical questions in river engineering is also one of the most visual: will the bed stay put, or will the flow start rolling grains downstream? The Shields parameter gives a quick screening answer by comparing the force applied by flowing water with the submerged weight that resists grain motion. Instead of asking only how deep the water is or how steep the channel is, the Shields approach combines the main hydraulic controls into a dimensionless ratio that is easier to interpret across different grain sizes and sediment types.

This calculator is designed for that early decision point. You enter fluid density, sediment density, flow depth, energy slope, a representative grain size, and a reference critical Shields threshold. The tool then computes bed shear stress, the resulting Shields parameter, the threshold ratio, and the critical slope required to reach incipient motion at the depth you supplied. In plain language, it tells you whether the selected conditions appear too weak to move the bed, close enough to the threshold that uncertainty matters, or strong enough that sediment motion is likely.

That makes the calculator useful in many settings: a stream restoration concept review, a culvert or grade-control screening, a classroom exercise, a field comparison between baseflow and flood conditions, or a quick check on whether a proposed channel lining is likely to remain stable. It is intentionally simple, so it should be treated as a first-pass mobility estimate rather than a complete sediment transport design model. Even so, it captures the core idea behind many transport assessments: moving bigger grains requires more shear, and mild changes in depth or slope can push a bed from stable to mobile surprisingly quickly.

Estimate whether flow can move sediment

The Shields parameter helps river engineers, geomorphologists, and restoration planners decide whether water flow is strong enough to mobilize bed material. It compares the shear stress applied by the flow to the submerged weight resisting motion for a representative grain size. Because it is dimensionless, it is especially handy when you want to compare cases with different grain diameters without getting lost in raw force units alone.

A value below the critical threshold suggests the bed is likely stable. A value above the threshold suggests grains can begin to move as bed load. The threshold is not a fixed law of nature, so this calculator uses 0.045 as a practical reference value and shows the margin above or below that value. That margin matters. A result barely under the threshold is not the same as a result far below it, and a result well above the threshold can imply a stronger transport tendency than a simple yes-or-no stability label would suggest.

Inputs explained

Fluid density ρ is usually close to 1000 kg/m³ for fresh water. Saline or sediment-laden water can be slightly heavier, so if you are checking an estuary, hyperconcentrated flow, or laboratory mixture, use the density that best matches your fluid. A higher fluid density slightly increases the applied shear stress and also changes the submerged weight term in the denominator of the Shields expression.

Sediment density ρ_s is the density of the grains themselves. Quartz-rich sediment is commonly around 2650 kg/m³, which is why that value appears as the default. If you are evaluating lighter organic particles, denser heavy minerals, or manufactured channel material, changing this input can noticeably affect the mobility estimate because the resisting submerged weight depends on the difference between sediment density and fluid density.

Flow depth h is the representative hydraulic depth in meters for the case you want to test. In a screening calculation, this is often a normal depth, a measured stage, or a design-event depth. Doubling depth doubles the bed shear stress in this simplified wide-channel form, so depth is often one of the most influential inputs. If you are unsure, it is a good idea to test several depths that represent low flow, a bankfull-like event, and a larger storm condition.

Energy slope S should be entered as a decimal rather than a percent. For example, a one-tenth of one percent slope is 0.001. In practice, this slope is often approximated by water-surface slope, energy grade line slope, or bed slope for a near-uniform reach. Because the formula multiplies depth by slope, a small change in slope can shift the result meaningfully, especially in gravel-bed channels where the mobility threshold may be crossed by only a modest increase in driving force.

Median grain size d is entered in millimetres, but the calculator converts it to meters internally before computing the Shields parameter. Using d50 is common for a first pass because it gives a single representative size for a mixed bed. Still, mixed sediment can hide or expose grains in ways that a single d50 cannot fully describe. If the site contains a broad gravel range, try several grain sizes such as a surface d50 and a coarser percentile to see how sensitive the conclusion is.

Critical Shields threshold θ_c is the reference value used to judge whether the computed θ is below, near, or above incipient motion. The default 0.045 is a widely used screening value for noncohesive sediment, but it is not universal. Bed armoring, exposure, fine sediment hiding, grain angularity, and hydraulic roughness can all shift the effective threshold. Entering a different value here is a simple way to run a sensitivity check without changing the rest of the calculation.

Formula

Bed shear stress for a wide, uniform open channel is estimated as:

Plain-text formulas: bedShearStress = fluidDensity * g * flowDepth * energySlope; shieldsTheta = bedShearStress / ((sedimentDensity - fluidDensity) * g * grainDiameter).

The Shields parameter is:

where ρ is fluid density, ρ_s is sediment density, h is flow depth, S is water-surface or energy slope, and d is median grain size converted from millimeters to meters. In this particular form, gravitational acceleration appears in both the shear-stress calculation and the resisting submerged-weight term. That means the screening outcome is often easiest to think of as a competition between the hydraulic product h × S and the grain-size term in the denominator. All else equal, deeper flow and steeper slope increase θ, while larger grains decrease it.

The critical Shields parameter depends on particle Reynolds number, grain shape, sorting, bed structure, hydraulic conditions, and measurement uncertainty. The default value 0.045 is a screening reference, not a site-specific design threshold. If you have a published threshold from a local study, flume calibration, or agency guidance for a similar bed material and hydraulic regime, enter that value and compare the results directly.

How to use the result

The result table is intended to support interpretation rather than replace it. Bed shear stress τ is shown in pascals, which is useful if you want a dimensional force measure. Critical shear τ_c tells you the approximate stress associated with the threshold you selected. The dimensionless Shields parameter θ and the ratio θ/θ_c then tell you how close your case is to the mobility reference point. The category bands in this tool are set to below threshold when θ is less than 90% of the critical value, near threshold when it falls within 90% to 110%, and above threshold when it exceeds 110%.

• θ < 0.045: sediment is probably stable for the selected conditions.
• θ near 0.045: conditions are near incipient motion; field uncertainty matters.
• θ > 0.045: bed material is likely mobile, especially if turbulence and sorting expose the grains.

The result also reports the critical bed shear stress and the approximate slope needed to reach the threshold at the depth you entered. That last value is especially useful in planning discussions because it answers a common reverse question: if the channel were this deep, how steep would it need to be before the selected sediment started moving? For restoration concepts or channel checks, that can help separate stable layouts from ones that may need roughness, grade control, or coarser material.

Field interpretation checklist

Small changes in slope, depth, or grain size can move the ratio across the threshold. Result categories use the threshold you enter: below threshold, near threshold, and above threshold. For design work, run a sensitivity set with low, best-estimate, and high values for depth and d50. That approach often tells you more than a single deterministic run because channel geometry, roughness, and surface texture rarely stay fixed in the field.

Worked example

Suppose the channel carries fresh water with density 1000 kg/m³, the bed material is quartz-rich gravel with density 2650 kg/m³, the representative depth is 0.8 m, the energy slope is 0.001, and the median grain size is 20 mm. The calculator gives a bed shear stress of about 7.85 Pa. Converting the grain size to meters and applying the Shields formula yields θ of about 0.024. Compared with a critical value of 0.045, the threshold ratio is only about 0.54, which places the case below threshold in this screening framework.

That does not prove the bed is motionless in every patch of the reach. Local turbulence, protruding grains, and finer fractions may still move. What it does say is that the selected representative grain size is not strongly driven above the usual incipient-motion reference under the assumed hydraulic conditions. If you increase the slope, increase the depth, or reduce the grain size, the ratio rises. If you test a coarser gravel or a more armored threshold, the ratio falls.

Assumptions and limitations

This is a screening calculation. It does not account for grain Reynolds number, hiding and exposure in mixed beds, cohesive sediment, bedforms, unsteady floods, channel curvature, vegetation, armoring, or local scour around structures. Use measured hydraulic data and a sediment transport method appropriate to your site before making design decisions. The wide-channel shear-stress approximation is most defensible when the section is reasonably uniform and the hydraulic radius can be represented well by depth.

Another important limitation is time. The calculator estimates whether the flow is capable of initiating motion, not how much sediment will be transported over an event or how quickly a bar will reshape. Sediment can be near threshold for a long period without large geomorphic change, or it can exceed threshold briefly during a sharp flood peak. If duration matters, combine this type of screening with a hydrograph, field observations, and a transport relation suited to the sediment and channel form you are studying.

Design use

For restoration or culvert work, pair the Shields result with site evidence. Fresh bars, exposed roots, embedded gravels, scour pools, and depositional wedges can confirm whether the modeled mobility agrees with observed behavior. If the model says stable but the channel shows recent transport, the depth, slope, grain size, or roughness assumptions may be too mild. If the model predicts strong mobility but the bed is clearly armored, your effective threshold may need to be raised for that surface layer.

Report the selected critical value with the result. A threshold of 0.045 is a common screening reference, but published curves vary. Stating the threshold keeps reviews transparent and makes it easier to rerun the same site with a different incipient-motion criterion. For monitoring, rerun the calculation after major floods or channel work using updated depth, slope, and grain-size observations. A channel that was stable before restoration can become mobile if grade control, vegetation, or bed material changes.

Used thoughtfully, this calculator is a compact way to connect field measurements with a physically meaningful mobility indicator. It does not replace surveying, hydraulic modeling, pebble counts, or sediment transport analysis, but it does help frame the right next question: are you comfortably below the motion threshold, operating in an uncertain transition zone, or already in conditions where bed load should be expected?