Shields Parameter Sediment Transport Calculator

Understanding the Shields parameter for sediment motion

A Shields calculation is most useful when you want a quick answer to a practical river question: given this depth, slope, water density, and grain size, is the bed likely to stay in place or begin to move? The ratio is built for exactly that kind of screening, because it compares the driving stress from flowing water with the resisting submerged weight of a representative sediment grain. Instead of relying on depth or slope alone, the Shields approach folds the main hydraulic controls into one dimensionless number that is easier to compare across different channels and sediment textures.

This calculator turns your inputs into a compact mobility check. It computes bed shear stress, the Shields parameter, the threshold ratio against your chosen critical value, and the slope that would be needed to reach incipient motion at the same depth. That makes the output useful when you are comparing a riffle to a pool, testing a restoration concept, checking a culvert approach, or making a quick field estimate before you move on to a fuller sediment transport method. The result is a first-pass judgment about mobility, not a complete design answer.

Because of that, the calculator is best used to identify where a reach is comfortably stable, where it is likely mobile, and where the site sits in the uncertain transition zone around the motion threshold. A small change in depth, slope, or grain size can matter a lot, especially in gravel-bed channels where the selected d50 may be close to the size that starts moving during moderate flows. If the answer is near the threshold, that is often a cue to inspect the bed surface, the grain mix, and the local hydraulic setting more carefully.

Estimate whether flow can move sediment

In Shields screening, the key comparison is between the flow's bed shear stress and the submerged weight of the grains you selected. If the ratio stays below the critical value, the bed is more likely to remain stable; if it rises above the critical value, incipient motion becomes more likely. The result is not a promise that every grain will move everywhere at once, but it does tell you whether the hydraulic forcing is strong enough to start mobilizing the chosen representative size.

The calculator uses 0.045 as a practical reference critical Shields parameter for noncohesive sediment. That value is widely used for quick checks, but it is still only a screening reference. Real thresholds vary with grain Reynolds number, hiding and exposure, sorting, bedforms, cohesion, turbulence, and how armored the bed surface has become. When your result lands very close to the threshold, the main lesson is to treat the answer as approximate rather than exact.

Input guidance for Shields sediment checks

Fluid density ρ is usually close to 1000 kg/m³ for freshwater, but that is not the only case worth checking. Slightly saltier water, sediment-laden water, or a laboratory mixture can change the density enough to matter when you are near the mobility boundary. In the Shields calculation, a higher fluid density increases the applied shear stress and also changes the submerged weight term, so it affects both the numerator and the denominator of the ratio.

Sediment density ρ_s describes the material making up the grains. Quartz-rich sand and gravel are often near 2650 kg/m³, which is why that value appears as the default. If the bed includes lighter organic particles, heavier minerals, or engineered material, update this input before interpreting the result. The density contrast between the grain and the fluid controls how much resistance the submerged weight provides.

Flow depth h is the representative hydraulic depth for the reach you want to test. In screening work, that might be a normal-depth estimate, a measured stage, a design depth, or a flood depth. Because the simplified wide-channel form multiplies depth directly into bed shear stress, this input often has a strong effect on the outcome. When you are unsure which depth matters most, it is worth trying several plausible values that span low flow, bankfull-like conditions, and a larger storm event.

Energy slope S should be entered as a decimal, not a percent. A slope of one-tenth of one percent is 0.001. In a near-uniform reach, this may be approximated by the water-surface slope, the energy grade line slope, or the bed slope. Since the formula uses depth times slope as the main driving term, even a small change in slope can shift the Shields number enough to move a coarse surface layer across the critical boundary.

Median grain size d is entered in millimetres, and the calculator converts it to meters internally. Using d50 is a practical shortcut because it gives one representative size for a mixed bed. Even so, real beds are seldom uniform. Hiding, armoring, and a wide size distribution can make a single d50 look more stable or more mobile than the full surface actually is. If the bed is mixed, try a few sizes to see how sensitive the result is to the chosen representative grain.

Critical Shields threshold θ_c is the reference level used to judge whether the computed θ is below, near, or above incipient motion. The default 0.045 is a common screening value for noncohesive material, but it is not universal. Bed structure, grain exposure, angularity, and hydraulic roughness can all shift the effective threshold. If you have a published value from a local study or a site-specific calibration, enter it here and compare the result to that local criterion instead of relying only on the default.

Formula for Shields bed shear stress and incipient motion

For a broad, approximately uniform open channel, the calculator estimates bed shear stress 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 form, the ratio is easiest to read as a competition between the hydraulic product h × S and the grain-size term in the denominator. Deeper water and steeper slope raise θ; larger grains lower it. That is why coarse gravel can remain stable under conditions that would move sand without much difficulty.

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 read the Shields sediment result

The result table is meant to help you interpret the Shields check, not to replace judgment from the field. Bed shear stress τ is shown in pascals, which gives you a dimensional sense of the forcing. Critical shear τ_c is the stress associated with the critical Shields value you entered. The dimensionless Shields parameter θ and the ratio θ/θ_c then show how close your case is to the mobility threshold. In this tool, the category bands are defined as 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 reverse check is useful in planning discussions because it answers a practical question: if the channel were this deep, how steep would it need to be before the selected sediment began to move? For restoration concepts or channel stability checks, that can help separate layouts that are comfortably stable from those that may need roughness, grade control, or larger material.

Field interpretation checklist for Shields screening

Use the checklist below when the Shields result and the reach you are studying do not quite match. Field evidence often shows why a simple screening calculation feels right in one location and wrong in another. A mixed gravel bed, an armored surface, or cohesive material can shift the effective threshold enough that the calculator should be read as a starting point rather than a final verdict.

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 often tells you more than a single deterministic run because channel geometry, roughness, and surface texture rarely stay fixed in the field.

Worked example: gravel-bed flow at the Shields threshold

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 everywhere in the reach. Local turbulence, protruding grains, and finer fractions can still move, especially if the surface is uneven or partially sorted. What the result does say is that the chosen 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 for Shields sediment screening

This calculator is a screening tool. 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 limitation is time. The calculator estimates whether 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 remain near threshold for a long period without obvious geomorphic change, or it can exceed threshold briefly during a sharp flood peak and still produce a visible response. If duration matters, combine this screening with a hydrograph, field observations, and a transport relation suited to the sediment and channel form you are studying.

Design use for river and channel projects

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?