Model:SedDepEroder

From CSDMS



SedDepEroder


Metadata

Also known as
Model type Single
Model part of larger framework LandLab
Incorporated models or components:
Spatial dimensions 2D
Spatial extent Regional-Scale, Landscape-Scale
Model domain Terrestrial
One-line model description Compute fluvial erosion using using “tools and cover” theory
Extended model description This module implements sediment flux dependent channel incision following::

E = f(Qs, Qc) * ((a stream power-like term) - (an optional threshold)),

where E is the bed erosion rate, Qs is the volumetric sediment flux into a node, and Qc is the volumetric sediment transport capacity at that node.

This component is under active research and development; proceed with its use at your own risk.

Keywords:

stream power, erosion, tools and cover,


First name Daniel
Last name Hobley
Type of contact
Institute / Organization
Postal address 1
Postal address 2
Town / City Boulder
Postal code 80309
State Colorado
Country United States
Email address daniel.hobley@colorado.edu
Phone
Fax


Supported platforms Unix, Linux, Mac OS, Windows
Other platform
Programming language Python
Other program language
Code optimized Single Processor
Multiple processors implemented
Nr of distributed processors
Nr of shared processors
Start year development 2014
Does model development still take place? No
If above answer is no, provide end year model development 2020
Code development status Only maintenance
When did you indicate the 'code development status'? 2020
Model availability As code
Source code availability
(Or provide future intension)
Through web repository
Source web address https://github.com/landlab/landlab/tree/master/landlab
Source csdms web address
Program license type BSD or MIT X11
Program license type other
Memory requirements
Typical run time


Describe input parameters grid : a ModelGrid

A grid.

K_sp : float (time unit must be *years*)
K in the stream power equation; the prefactor on the erosion equation (units vary with other parameters).

g : float (m/s**2)
Acceleration due to gravity.

rock_density : float (Kg m**-3)
Bulk intact rock density.

sediment_density : float (Kg m**-3)
Typical density of loose sediment on the bed.

fluid_density : float (Kg m**-3)
Density of the fluid.

runoff_rate : float, array or field name (m/s)
The rate of excess overland flow production at each node (i.e., rainfall rate less infiltration).

pseudoimplicit_repeats : int
Number of loops to perform with the pseudoimplicit iterator, seeking a stable solution. Convergence is typically rapid.

return_stream_properties : bool
Whether to perform a few additional calculations in order to set the additional optional output fields, 'channel__width', 'channel__depth', and 'channel__discharge' (default False).

sed_dependency_type : {'generalized_humped', 'None', 'linear_decline', 'almost_parabolic'}
The shape of the sediment flux function. For definitions, see Hobley et al., 2011. 'None' gives a constant value of 1. NB: 'parabolic' is currently not supported, due to numerical stability issues at channel heads.

Qc : {'power_law', 'MPM'}
Whether to use simple stream-power-like equations for both sediment transport capacity and erosion rate, or more complex forms based directly on the Meyer-Peter Muller equation and a shear stress based erosion model consistent with MPM (per Hobley et al., 2011).

If ``sed_dependency_type == 'generalized_humped'``...

kappa_hump : float
Shape parameter for sediment flux function. Primarily controls function amplitude (i.e., scales the function to a maximum of 1). Default follows Leh valley values from Hobley et al., 2011.

nu_hump : float
Shape parameter for sediment flux function. Primarily controls rate of rise of the "tools" limb. Default follows Leh valley values from Hobley et al., 2011.

phi_hump : float
Shape parameter for sediment flux function. Primarily controls rate of fall of the "cover" limb. Default follows Leh valley values from Hobley et al., 2011.

c_hump : float
Shape parameter for sediment flux function. Primarily controls degree of function asymmetry. Default follows Leh valley values from Hobley et al., 2011.

If ``Qc == 'power_law'``...

m_sp : float
Power on drainage area in the erosion equation.

n_sp : float
Power on slope in the erosion equation.

K_t : float (time unit must be in *years*)
Prefactor in the transport capacity equation.

m_t : float
Power on drainage area in the transport capacity equation.

n_t : float
Power on slope in the transport capacity equation.

if ``Qc == 'MPM'``...

C_MPM : float
A prefactor on the MPM relation, allowing tuning to known sediment saturation conditions (leave as 1. in most cases).

a_sp : float
Power on shear stress to give erosion rate.

b_sp : float
Power on drainage area to give channel width.

c_sp : float
Power on drainage area to give discharge.

k_w : float (unit variable with b_sp)
Prefactor on A**b_sp to give channel width.

k_Q : float (unit variable with c_sp, but time unit in *seconds*)
Prefactor on A**c_sp to give discharge.

mannings_n : float
Manning's n for the channel.

threshold_shear_stress : None or float (Pa)
The threshold shear stress in the equation for erosion rate. If None, implies that *set_threshold_from_Dchar* is True, and this parameter will get set from the Dchar value and critical Shields number.

Dchar :None, float, array, or field name (m)
The characteristic grain size on the bed, that controls the relationship between critical Shields number and critical shear stress. If None, implies that *set_Dchar_from_threshold* is True, and this parameter will get set from the threshold_shear_stress value and critical Shields number.

set_threshold_from_Dchar : bool
If True (default), threshold_shear_stress will be set based on Dchar and threshold_Shields.

set_Dchar_from_threshold : bool
If True, Dchar will be set based on threshold_shear_stress and threshold_Shields. Default is False.

threshold_Shields : None or float
The threshold Shields number. If None, implies that *slope_sensitive_threshold* is True.

slope_sensitive_threshold : bool
If True, the threshold_Shields will be set according to 0.15 * S ** 0.25, per Lamb et al., 2008 & Hobley et al., 2011.

flooded_depths : array or field name (m)
Depths of flooding at each node, zero where no lake. Note that the component will dynamically update this array as it fills nodes with sediment (...but does NOT update any other related lake fields).

Input format ASCII
Other input format
Describe output parameters "channel__bed_shear_stress":

{ "dtype": float, "intent": "out", "optional": False, "units": "Pa", "mapping": "node", "doc": "Shear exerted on the bed of the channel, assuming all discharge travels along a single, self-formed channel", }

"channel__depth": { "dtype": float, "intent": "out", "optional": True, "units": "m", "mapping": "node", "doc": "Depth of the a single channel carrying all runoff through the node", }

"channel__discharge": { "dtype": float, "intent": "out", "optional": False, "units": "m**3/s", "mapping": "node", "doc": "Volumetric water flux of the a single channel carrying all runoff through the node", }

"channel__width": { "dtype": float, "intent": "out", "optional": True, "units": "m", "mapping": "node", "doc": "Width of the a single channel carrying all runoff through the node", }

"channel_sediment__relative_flux": { "dtype": float, "intent": "out", "optional": False, "units": "-", "mapping": "node", "doc": "The fluvial_sediment_flux_into_node divided by the fluvial_sediment_transport_capacity", }

"channel_sediment__volumetric_flux": { "dtype": float, "intent": "out", "optional": False, "units": "m**3/s", "mapping": "node", "doc": "Total volumetric fluvial sediment flux brought into the node from upstream", }

"channel_sediment__volumetric_transport_capacity": { "dtype": float, "intent": "out", "optional": False, "units": "m**3/s", "mapping": "node", "doc": "Volumetric transport capacity of a channel carrying all runoff through the node, assuming the Meyer-Peter Muller transport equation", }

"drainage_area": { "dtype": float, "intent": "in", "optional": False, "units": "m**2", "mapping": "node", "doc": "Upstream accumulated surface area contributing to the node's discharge", }

"flow__link_to_receiver_node": { "dtype": int, "intent": "in", "optional": False, "units": "-", "mapping": "node", "doc": "ID of link downstream of each node, which carries the discharge", }

"flow__receiver_node": { "dtype": int, "intent": "in", "optional": False, "units": "-", "mapping": "node", "doc": "Node array of receivers (node that receives flow from current node)", }

"flow__upstream_node_order": { "dtype": int, "intent": "in", "optional": False, "units": "-", "mapping": "node", "doc": "Node array containing downstream-to-upstream ordered list of node IDs", }

"topographic__elevation": { "dtype": float, "intent": "inout", "optional": False, "units": "m", "mapping": "node", "doc": "Land surface topographic elevation", }

"topographic__steepest_slope": { "dtype": float, "intent": "in", "optional": False, "units": "-", "mapping": "node", "doc": "The steepest *downhill* slope", }

Output format ASCII
Other output format
Pre-processing software needed? No
Describe pre-processing software
Post-processing software needed? No
Describe post-processing software
Visualization software needed? No
If above answer is yes
Other visualization software


Describe processes represented by the model
Describe key physical parameters and equations
Describe length scale and resolution constraints
Describe time scale and resolution constraints
Describe any numerical limitations and issues


Describe available calibration data sets
Upload calibration data sets if available:
Describe available test data sets
Upload test data sets if available:
Describe ideal data for testing


Do you have current or future plans for collaborating with other researchers?
Is there a manual available? No
Upload manual if available:
Model website if any https://landlab.github.io
Model forum / discussion board https://github.com/landlab/landlab/issues
Comments


This part will be filled out by CSDMS staff

OpenMI compliant Not yet"Not yet" is not in the list (Yes, No but planned, No but possible, No not possible) of allowed values for the "Code openmi compliant or not" property.
BMI compliant Not yet"Not yet" is not in the list (Yes, No but planned, No but possible, No not possible) of allowed values for the "Code IRF or not" property.
WMT component Not yet"Not yet" is not in the list (Yes, In progress, No but possible, No not possible) of allowed values for the "Code CMT compliant or not" property.
PyMT component Not yet"Not yet" is not in the list (Yes, In progress, No but possible, No not possible) of allowed values for the "Code PyMT compliant or not" property.
Is this a data component
Can be coupled with:
Model info
Citation indices SedDepEroder
Nr. of pubs: --
Citations: 0
h-index: --"--" is not a number.
Qrcode SedDepEroder.png
Link to this page



Introduction

History

References




Citation indices SedDepEroder
Nr. of pubs: --
Citations: 0
h-index: --"--" is not a number.


See more publications of SedDepEroder

Issues

Help

Input Files

Output Files