Model help:AgDegNormalSub

The CSDMS Help System


This program is used to calculate the evolution of upward-concave bed profiles in rivers carrying uniform sediment in subsiding basins.

Model introduction

The program computes the approach to mobile-bed equilibrium in a river carrying uniform material and flowing into a subsiding basin. It is a descendant of AgDegNormal.

Model parameters

Parameter Description Unit
Input directory path to input files
Site prefix Site prefix for Input/Output files
Case prefix Case prefix for Input/Output files
Parameter Description Unit
Flood discharge m3 / s
Intermittency -
Grain size of bed material mm
I -
Manning-Strickler coefficient, k -
Slope of forest face -
upstream bed material sediment feed rate during floods m2 / s
L -
Time step days
Iterations per each printout
Number of printout m
Number of fluvial nodes
Manning-Strickler coefficient, r
Coefficient in total bed material relation
Exponent in load relation
Critical Shield stress
Submerged specific gravity of sediment
initial length of fluvial zone m
B -
O -
Y -
Parameter Description Unit
Model name name of the model -
Author name name of the model author -

Uses ports

This will be something that the CSDMS facility will add

Provides ports

This will be something that the CSDMS facility will add

Main equations

  • Exner equation for uniform sediment from a river
[math]\displaystyle{ \left ( 1 - \lambda _{p} \right ) \left ( {\frac{\partial \eta }{\partial t}} + \delta \right ) = - {\frac{I_{f}}{r_{B}}} \left ( 1 + \Lambda \right ) \Omega {\frac{\partial q_{t}}{\partial x}} }[/math] (1)
  • Ratio of depositional to channel width
[math]\displaystyle{ r_{B} = {\frac{B_{d}}{B_{c}}} }[/math] (2)
  • The maximum possible length of the fluvial reach
[math]\displaystyle{ L_{max} = {\frac{I_{f} \left ( 1 + \Lambda \right ) \Omega }{r_{B}}} {\frac{q_{tf}}{\left ( 1 - \lambda _{p} \right ) \delta}} }[/math] (3)


some assumptions: The subsidence rate s is assumed to be constant in time and space. The sediment is assumed to be uniform with size D. A Manning-Strickler formulation is used for bed resistance. A generic relation of the general form of that due to Meyer-Peter and Muller is used for sediment transport. The flow is computed using the normal flow approximation. The river is assumed to have a constant width.

The program is a descendant of AgDegNormal. Three relatively minor changes have been implemented as follows: a) The input parameters have been modified to include the following parameters: subsidence rate σ, ratio of depositional width to channel width rB, ratio of wash load deposited per unit bed material load Λ and channel sinuosity Ω; b) The code has been modified so as to include subsidence in the calculation of mass balance; c) The output has been modified to show the time evolution of not only the profile of bed elevation η, but also the profiles of bed slope S and the ratio qt/qtf, where qtf denotes the volume feed rate of bed material load per unit width.

The ratio of depositional to channel width, rB, has been introduced to model the fact that in an aggrading river sediment deposits not only in the channel itself, but also in a much wider belt (e.g. the floodplain or basin width, due to overbank deposition, channel migration and avulsion). Here channel width is denoted as Bc (which can be taken to be synonymous with bankfull width) and effective depositional width is denoted as Bd.

The parameter Λ that represents the units of wash load deposited per unit of bed material is introduced to consider that in the 1D formulation implemented in this model it is assumed that deposition occurs not only in the channel but on a much wider area (e.g. the floodplain). Sediment deposited in the channel is mostly made of bed material but sediment deposited around the channel contains a significant amount of wash load. A precise mass balance for wash load is beyond the scope of this model. For simplicity it is assumed that for every unit of sand deposited in the system, Λ units of wash load are deposited. It is also assumed that the supply of wash load from upstream is always sufficient for deposition at such a rate. This is not likely to be strictly true, but should serve as a useful starting assumption.

The parameter Ω has been introduced to consider that channels may be sinuous. Here it is assumed that the channel has a sinuosity, Ω, but that the depositional surface across which it wanders is rectangular. In the present formulation the sinuosity is defined as the ratio of downchannel distance per unit of downvalley distance.

Boundary and initial conditions are equal to that implemented for the ancestor model AgDegNormal.

  • Note on model running

Flow is calculated assuming normal flow approximation

If the input channel length is longer than the maximum possible length of the fluvial reach, the program cannot perform the calculation. The maximum possible length of the fluvial reach, Lmax, is defined as the maximum length of basin that the sediment supply can fill; at this length the sediment transport rate out of the basin drops precisely to zero.


An example run with input parameters, BLD files, as well as a figure / movie of the output

Follow the next steps to include images / movies of simulations:

See also: Help:Images or Help:Movies


Gary Parker


Key papers