Model:SedDepEroder

Edit Model information

SedDepEroder

Metadata

Summary Contact Technical specs In/Output Process Testing Other Component info Also known as Model type Single Model part of larger framework LandLab Note on status model Date note status model 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, Name Daniel 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 Authors Daniel Hobley Source code Go to external source code site Model citations Nr. of publications: -- Total citations: 0 h-index: --"--" is not a number. m-quotient: 0 QR-code Link to this page Other models by author ChiFinder FastscapeEroder LakeMapperBarnes PerronNLDiffuse PotentialityFlowRouter SedDepEroder SinkFiller SpatialPrecipitationDistribution SteepnessFinder

References

Automatically enter Reference by DOI

Manually enter Reference

Nr. of publications:	--
Total citations:	0
h-index:	--"--" is not a number.
m-quotient:	0

See more publications of SedDepEroder