Shields Parameter Calculator: Sediment Transport Threshold & Bedload

simulator intermediate ~11 min
Loading simulation...
θ = 0.076 — active bedload transport of coarse sand

A 2 mm quartz grain with shear velocity 0.05 m/s gives Shields parameter 0.076, exceeding the critical value of ~0.047. Active bedload transport by saltation and rolling is occurring.

Formula

θ = ρ_f × u*² / ((ρ_s - ρ_f) × g × d)
q_b* = 8 × (θ - θ_cr)^1.5 (Meyer-Peter & Mueller)
Re* = u* × d / ν (grain Reynolds number)

The Threshold of Motion

A river flowing over a sand bed does not always move the sand — there is a critical shear stress below which grains remain stationary, locked in place by friction and their own weight. Albert Shields, in his 1936 doctoral thesis at the Technical University of Berlin, discovered that this threshold could be expressed as a dimensionless parameter comparing fluid force to gravitational resistance. The Shields diagram remains the single most important tool in sediment transport engineering nearly a century later.

Shields Parameter Physics

The Shields parameter theta = tau/((rhos-rhof)gd) compares the destabilizing fluid shear stress tau to the stabilizing submerged weight of grains per unit area. When theta exceeds the critical value theta_cr (approximately 0.03-0.06, depending on grain Reynolds number), grains begin to roll and saltate. The critical value is not a single number but a curve in Shields space, reflecting the transition from viscous-dominated to inertia-dominated grain environments.

Bedload Transport Formulae

Once grains are in motion, the transport rate increases rapidly with excess shear stress. The Meyer-Peter and Mueller formula q* = 8(theta - theta_cr)^1.5 relates dimensionless transport rate to the excess above threshold. This 3/2 power law means that doubling the excess stress nearly triples the transport rate — explaining why floods carry disproportionately more sediment than average flows. Modern formulations by van Rijn, Engelund-Hansen, and others refine this for specific grain-size distributions and flow conditions.

Applications in River Engineering

The Shields criterion guides the design of stable river channels, bridge foundations, and coastal protection. Engineers select armor stone sizes so that the Shields parameter remains below critical under design flood conditions. Conversely, understanding when sediment moves helps predict channel migration, reservoir siltation, and habitat change. In environmental restoration, managed floods are designed to mobilize fine sediment and restore natural channel processes.

FAQ

What is the Shields parameter?

The Shields parameter θ = τ/((ρs-ρf)gd) is a dimensionless ratio of fluid shear stress to submerged grain weight per unit area. When θ exceeds a critical value (θ_cr ≈ 0.045-0.06), grains begin to move. It was introduced by Albert Shields in 1936 and remains the fundamental criterion for sediment entrainment in rivers, estuaries, and the sea.

What is the critical Shields value for sand?

For medium to coarse sand (0.25-2 mm), the critical Shields parameter is approximately 0.03-0.06, depending on grain Reynolds number. Fine sand requires slightly higher θ_cr due to cohesive effects, while well-sorted gravel has θ_cr around 0.045. The Shields curve plots θ_cr versus Re* to capture this dependency.

What is the difference between bedload and suspended load?

Bedload consists of grains rolling, sliding, and saltating (bouncing) along the bed — they maintain contact with the bottom. Suspended load consists of grains carried aloft by turbulent eddies when vertical turbulent velocity exceeds settling velocity. The transition occurs when u*/w_s > 1, where w_s is settling velocity.

How is bedload transport rate calculated?

The Meyer-Peter and Mueller (1948) formula gives dimensionless bedload transport as q* = 8(θ-θ_cr)^1.5. Dimensional transport rate is q_b = q* × sqrt((s-1)gd³), where s=ρs/ρf. More advanced formulas by Einstein, Engelund-Hansen, and van Rijn refine predictions for specific conditions.

Sources

Embed

<iframe src="https://homo-deus.com/lab/sedimentology/sediment-transport/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub