Identify 2D sediment transport equations

Step-by-step guide

In InfoWorks ICM, empirical equations are used to simulate 2D sediment transport in rivers.

The QM Parameters and the Water quality and sediment parameters are used to apply different equations that model suspended sediment and bed load transport in 2D.

To set the simulation parameters:

  1. In the Run dialog, select QM parameters.

The Run dialog with QM parameters highlighted in red.

  1. In the QM Parameters dialog, select either:
  • Model 2D erosion-deposition to simulate both bed and suspended sediment loads.
  • Model 2D bed load to simulate bed load transport only.

The QM Parameters dialog, with the options for Model 2D erosion-deposition and Model 2D bed load highlighted in red.

To set the water quality parameters:

  1. Select Model > Model parameters > Water quality and sediment parameters.
  2. In the Properties panel, select a Model type.
  • Suspended load models only the suspended load.
  • Total load models the total sediment load, including the suspended load and the bed load.
  • Bed load models only the bed load.

It is important to keep the following in mind:

  • If the Model type is set to Total load, verify that only one option is selected in QM Parameters. Selecting both will cause double counting of bed load transport and will result in simulation failure.
  • If Model 2D bed load is selected in QM Parameters, the only available Model type will be Bed load transport. Simulating suspended sediment is not possible with this option.

The table below shows the various combinations of parameters and the result of each:

Water quality parametersSimulation parametersResult
Suspended loadModel 2D erosion-depositionSuspended material only
Total loadModel 2D erosion-depositionBoth suspended and bed load
Bed loadModel 2D bed loadBed load
Suspended load (N/A)Model 2D bed loadNot possible

Suspended load model

Remember that Suspended load is the suspension of fine particles, such as silt and clay, within the mass of fluid.

A cross-sectional diagram of a river bed illustrating suspended sediment modeling.

For the 2D erosion/deposition model, if the Model type is set to Suspended load, three additional options become available for selecting the corresponding equations.

In Properties, the Model type set to Suspended load and the three corresponding parameters now available, all highlighted in red.

  • Near-bed concentration is the concentration of particles in the fluid closest to the bed of a river.
  • Skin friction is the resistance caused by fluid viscosity.
  • Equilibrium near-bed concentration is the state where particle concentration near the bed remains constant over time. This occurs when the rate of particle deposition onto the bed is equal to the rate of resuspension into the water column.

For Near-bed concentration, two equations are available:

  • The Rouse profile approximation describes the vertical distribution of suspended sediment concentration and considers the balance between upward turbulent diffusion and downward settling of particles.
  • The Lin formulation models the vertical concentration profile and considers specific flow and sediment characteristics for a more tailored concentration distribution.

For Equilibrium near-bed concentration, the available equations are Van Rijn, Zyserman and Fredsoe, and Smith and McLean. The Van Rijn equation, for example, estimates the reference concentration based on sediment size, flow velocity, and bed shear stress.

For Skin friction, available equations are based on experimental tests.

Bed load model

Several equations are available to calculate the carrying capacity of the flow as a volumetric transport rate.

  1. In Properties, under 2D bed load model, expand the Bed load formula drop-down.
  2. Select an equation. For example, Meyer-Peter-Müller is typically selected for sand and fine gravel bed material.

In the 2D bed load model parameters, the Bed load formula drop-down expanded with the list of available equations highlighted in red.

Total load model

To accurately simulate sediment transport dynamics, it is crucial to select the appropriate equation.

  1. In Properties, under 2D erosion/deposition model, expand the Total load model drop-down to select from the available options: Ackers-White, Engelund-Hansen, Westrich-Jurashek, Van Rijn, and Velikanov.

On the left, in Properties, the Model type being set to Total load, and on the right, the Total load model drop-down expanded to show the available options.

Each equation has specific applications based on sediment type, flow conditions, and environmental settings. For example, Van Rijn is a versatile equation that predicts sediment transport rates, including bed load and suspended load in rivers, accounting for sediment size, flow depth, and current velocity.

Settling velocity

This is the terminal velocity at which a particle in still fluid ceases to accelerate as the forces acting on it are balanced.

  1. In Properties, under Sediment fraction 1, expand the drop-down for 2D sf1 settling velocity (m/hr).
  2. Choose from the available formulas: Van Rijn, Soulsby, Wu and Wang, Hallermeier, and Zanke.

Key parameters required for these calculations include particle size, D35, D50, specific gravity, density, porosity, skin roughness, and the angle of repose for sediment.

In Properties, the Sediment fraction 1 parameters, with the drop-down for 2D sf1 settling velocity expanded to show the available options.

Next, in ICM, perform a suspended load and total bed load simulation:

  1. In the Scenarios drop-down, select Suspended.
  2. On the GeoPlan, select the bed channel.
  3. In the Multiple Selection pop-up, select Sediment bed.
  4. Click OK.

This is a polygon, where different types of bed material and the concentration of each layer can be defined.

In the GeoPlan, the bed channel selected, and in Properties, the available Polygon definitions.

  1. In the Explorer, right-click Initial conditions 2D and select Open.
  2. In the Variable column, set the value to Bed sediment.
  3. Define the type of bed layer and the depth, as well as the concentration layer.
  4. Select Model > Model parameters > Water quality and sediment parameters.
  5. In Properties, under 2D erosion/deposition model, set the Model type to Suspended load.

For this example, based on the characteristics of the particles in the bed, set the following parameters:

  1. Set the Equilibrium near-bed concentration method to Van Rijn.
  2. Set Near-bed concentration to Rouse profile.
  3. Set the Skin friction method to Manning.
  4. Under Sediment fraction 1, specify the size of the particle, Specific gravity, Settling velocity, and other relevant properties.

In Properties, the 2D erosion/deposition model options and the Sediment fraction 1 parameters configured for this example and highlighted in red.

  1. From the Explorer, open a pollutograph. In this case, a previously created pollutograph named Q_concentration_suspended is opened.

In the Explorer, the expanded context menu for the pollutograph, with Open called out and being selected.

  1. On the SF1 tab, define the variation of particle concentration over time.
  2. Right-click one of the profiles and select Profile properties.
  3. In the Profile Properties dialog, for Profile type, select Suspended.
  4. Open the Run dialog for the Suspended load run and drag in the pollutograph.
  5. Select Use QM.
  6. Click QM parameters.

The Run dialog, with the Pollutograph added and highlighted in red, and with Use QM and QM parameters highlighted in red.

  1. In the QM Parameters dialog, select SF1.
  2. Verify that Model 2D erosion-deposition is selected.

QM Parameters with SF1 and Model 2D erosion-deposition selected.

This will perform a suspended load and total bed load simulation. To simulate only bed load, select Model 2D bed load instead.

  1. Click OK.
  2. Back in the Run dialog, click Run simulations.
  3. Confirm the selections in the Schedule Jobs dialog and click OK.
  4. In the Explorer, right-click the Suspended load result and select Open to display the plan view and two graphs.
  • The plan view shows a thematic map of the distribution of the carrying capacity values along the river reach.
  • The first graph shows that the values of carrying capacity are increased in one point close to the meander section.
  • The second graph shows suspended concentration at the same point where the carrying capacity was plotted.

The Suspended load result called out in the Explorer, with the resulting plan view and two graphs displayed in the GeoPlan.