Model help:GravelSandTransition

GravelSandTransition

This model is used to calculate evolution of long profile of river with a migrating gravel-sand transition and subject to subsidence or base level rise.

Model introduction

This program calculates the bed surface evolution at predefined nodes for a river with a gravel to sand transition, as well as calculating the relative location of the transition point and the slopes, bedload transport rate, shear stress, and water depth for plotting.

Model parameters

First tab header

Second tab header

Etc. tab header

Parameter	Description	Unit
First parameter	Description parameter	[Units]

Parameter	Description	Unit
First parameter	Description parameter	[Units]

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

A list of the key equations. HTML format is supported; latex format will be supported in the future

Nomenclature

Symbol	Description	Unit
Qbf	bankfull water discharge at flood	m3 / s
Ifg	flood intermittency, gravel-bed reach	-
Ifs	flood intermittency, sand-bed reach	-
Qgravel,feed	feed rate of gravel	-
Qsand,feed	feed rate of sand	-
Bc,grav	gravel-bed channel width	m
Bc,sand	sand-bed channel width	m
Bd,grav	depositional or floodplain width, gravel-bed reach	m
Bd,sand	depositional or floodplain width, sand-bed reach	m
Ωg	sinuosity of gravel-bed reach	-
Ωs	sinuosity of sand-bed reach	-
Λsg	fraction of sand deposited per unit gravel in depositional zone of gravel-bed reach	-
Λms	fraction of mud deposited per unit sand in depositional zone of sand-bed reach	-
Dg	grain size of gravel	mm
Ds	grain size of sand	mm
Czg	dimensionless Chezy resistance coefficient gravel-bed reach	-
Czs	dimensionless Chezy resistance coefficient sand-bed reach	-
L	reach length	m
Sgsl	initial position of gravel-sand transition (must be < L)	m
Sgl	initial slope of gravel-bed reach	-
Ssl	initial slope of sand-bed reach	-
σ	subsidence rate	mm / year
ξdo	initial water surface base level	m
ξd	rate of base level rise	mm / year
Δt	time step	years
λpg	bed porosity, gravel-bed reach	-
λps	bed porosity, sand-bed reach	-
R	submerged specific gravity of sediment	-
k	initial water surface base level	m
y	year to start sea level rise	yr
Y	year to end sea level rise	yr
p	number of prints	-
i	number of iterations per print	-
d	depositional flood plain width	m
M	number of spatial steps	-
x	downstream coordinate	m
η	bed elevation	m
Sl	bed slope	-
qb	volume bedload transport per unit width	m2 / s
H	water depth	m
τ	shear stress	N / m2
sgs	location of the gravel-sand transition	m
Lmax	maximum reach length	m
εmin	minimum water elevation	m

Output

Symbol	Description	Unit
Ss,feed,n	normal slope associated with sand feed rate	-
Hs,feed,n	normal depth associated with sand feed rate	m
Ss,feed,n	normal slope associated with gravel feed rate
Hs,feed,n	normal depth associated with gravel feed rate	m

Notes

This model calculates the long profile of a river with a gravel-sand transition. The model uses two grain sizes: size D_g for gravel and size D_s for sand. The river is assumed to be in flood for the fraction of time Ifg for the gravel-bed reach and fraction Ifs for the sand-bed reach. All sediment transport is assumed to take place when the river is in flood.

Gravel transport is computed using the Parker (1979) approximation of the Einstein (1950) bedload transport relation. Sand transport is computed using the total bed material transport relation of Engelund and Hansen (1967).

In this simple model the gravel is not allowed to abrade. Both the gravel-bed and sand-bed reaches carry the same flood discharge Q_bf.

Gravel is transported as bed material in, and deposits only in the gravel-bed reach. A small residual of gravel load is incorporated into the sand at the gravel-sand transition. Sand is transported as washload in the gravel-bed reach, and as bed material load in the sand-bed reach.

The model allows for depositional widths Bdgrav and Bdsand that are wider than the corresponding bankfull channel widths Bgrav and Bsand of the gravel-bed and sand-bed channels. As the channel aggrades, it is assumed to migrate and avulse to deposit sediment across the entire depositional width. For each unit of gravel deposited in the gravel-bed reach, it is assumed that Lamsg units of sand are deposited. For each unit of sand deposited on the sand-bed reach, it is assumed that Lamms units of mud are deposited.

The gravel-bed reach has sinuosity Omegag and the sand-bed reach has sinuosity Omegas.

Bed resistance is computed through the use of two specified constant Chezy resistance coefficients; C_zg for the gravel-bed reach and C_zs for the sand-bed reach.

• Note on model running

The values for the spatial step of both gravel and sand max out at 100, allowing 200 spatial steps total (the number of sand spatial nodes, Ms and gravel spatial nodes, Mg must be ≤ 100).

When the slopes are unequal at gravel sand transition the GetData function and the plot data output a slope of 0 at that point.

The reach length, L, must be less than the maximum reach length, L_max; if it is not the program will quit and alert the user what the value was for L_max.

The sand and gravel Froude numbers must be less than 1 (i.e. the flow is assumed subcritical in the sand and gravel portions of the bed); if they are not, the program will quit and alert the user what their values were.

The user inputted downstream water elevation, ξ_d, must be greater than the program calculated minimum water elevation, ξ_min; if this is not the case, the program will quit and the user will be alerted what the value is for ξ_min.

Examples

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:

• Upload file: http://csdms.colorado.edu/wiki/Special:Upload
• Create link to the file on your page: [[Image:<file name>]].

See also: Help:Images or Help:Movies

Developer(s)

Name of the module developer(s)

References

Key papers

Links

[Model:GravelSandTransitio]