Critical Shields values (
τc*) suitable for specific applications are back‐calculated from representative bed load samples in mountain streams and a flow competence/critical flow approach. The general increase of
τc50* (for the bed D50 size) as well as
τc16* and
τc84* (for the bed D16 and D84 sizes) with stream gradient Sx and also the stratification of
τc50* by relative flow depth and relative roughness are confirmed. Critical Shields values
τc16* are shown to exceed
τc50* by about threefold, while those for
τc84* are nearly half of
τc50*. However, it remains unclear to what extent physical processes or numerical artifacts contribute to determining critical Shields values. Critical bankfull Shields values (
τcbf*) back‐computed from the average largest particles mobile at bankfull flow DBmax,bf approach
τc16* at steep gradients and
τc84* at low gradients and therefore increase very steeply with Sx. The relation
τcbf*=f(Sx) is stratified by bed stability (D50/DBmax,bf) and predictable if bed stability can be field categorized. Noncritical Shields values (
τbf50*) computed from bankfull flow depth and the D50 size differ from
τc50* and
τcbf*. Only in bankfull mobile streams where D50/DBmax = 1 can τ*cbf,
τc50*, and
τbf50* be used interchangeably. In highly mobile streams, substituting
τcbf* by
τbf50* overpredicts the DBmax,bf size by up to fivefold and underpredicts DBmax,bf by the same amount in highly stable streams. A value of 0.03 is appropriate for
τcbf* only on low stability beds with Sx ≅ 0.01, but overpredicts DBmax,bf by 30‐fold on highly stable beds with Sx ≅ 0.1. Differences in field and computational methods also affect critical Shields values.