The Finite-Volume discretization of the bed-slope term is

1) S_{bki} = \frac{w_i}{A^W} \sum _f \kappa _{bkf} \left| q_{bkf} \right| \nabla _{\bot} z_{bf}^W L_f^W

where

w_{i} : weights applied to each subarea [-]

A^{W} : cell wet area [-]

κ_{bkf} : non-dimensional bed-slope coefficient at face f [-]

|q_{bkf}| : bed-load mass transport magnitude at face [M/L/T]

∇_⊥{z_{bf}}^W : face-normal gradient of wet elevations at face f [L]

{L_{f}}^{W} : length of wet portion of face f [L]

The be-slope term is computed within a subcell by applying depth weighting in order to capture the effect in which deeper parts of cell are expected to experience larger bed change compared to less deep portions of the cell. The weights are computed as

2) w_i = \frac{h_i^b}{\sum_k h_k^b}

where

h_{i} : subarea water depth [L]

b : nondimensional empirical coefficient (b ≥ 0)

A coefficient b = 0 represents applying the same slope-induced bed change to each subregion.