“…There are two groups of numerical models: depth‐resolving models [ Strauss and Glinsky , ; Yeh et al , ] and layer‐averaged models, which include the three‐equation model (TEM) and its variants [ Fukushima et al , ; Parker et al , ; Zeng and Lowe , ; Choi , ; Imran et al , ; Bradford and Katopodes , ; Kostic and Parker , , ; de Luna et al , ; Toniolo , ; Hu and Cao , ; Kostic et al , ; Eke et al , ; Hu et al , ; Lai and Wu , ; Kostic , ; Elfimov and Khakzad , ] and the Four‐Equation‐Model (FEM) and its variants [ Fukushima et al , ; Parker et al , ; Salaheldin et al , ; Pratson et al , ; Das et al , ; Fildani et al , ; Kostic and Parker , ; Yi and Imran , ; Eke et al , ; Kostic , ; Tracer et al , ]. The FEM differs from the TEM in the way in which the bed shear velocity ( u ∗ ) is computed [ Fukushima et al , ; Parker et al , ]: while the TEM computes u ∗ from the drag exerted on the bed, roughly approximated by where C D >0 is the bed drag coefficient and U the layer‐averaged velocity of the sediment‐water mixture, the FEM computes u ∗ from the assumption that the bed shear stress is proportional to the layer‐averaged turbulent kinetic energy ( k ), where α > 0 is the dimensionless proportionality constant.…”