Model help:MARSSIM: Difference between revisions
No edit summary |
No edit summary |
||
Line 1,163: | Line 1,163: | ||
|width=50px align="right"|(28) | |width=50px align="right"|(28) | ||
|} | |} | ||
(C) Total bed sediment discharge | |||
::::{| | |||
|width=500px|<math> Q_{sb} = K_{q} A^ \left (eb \right ) [K_{v} A^ \left (0.6 e \left (1 - b \right ) \right ) S^ \left (0.7\right ) - 1/ \Psi_{c}]^ p </math> | |||
|width=50px align="right"|(29) | |||
|} | |||
::::{| | |||
|width=500px|<math> K_{q} = {\frac{K_{e} \omega \left (1 - \mu \right ) K_{w} K_{a}^b}{\left ( \gamma_{s} - \gamma \right )^p d^ \left ( p - 1 \right )}} </math> | |||
|width=50px align="right"|(30) | |||
|} | |||
::::{| | |||
|width=500px|<math> K_{v} = \gamma \{{\frac{N K_{a}^ \left (1-b\right) }{K_{w} K_{n} K_{p}}}\} ^ \left (0.6\right ) </math> | |||
|width=50px align="right"|(31) | |||
|} | |||
(D) The net erosion of channel | |||
::::{| | |||
|width=500px|<math> {\frac{\partial z}{\partial t}}|_{c} = {\frac{\left (Q_{sbi} - Q_{sbo}\right )}{LW}} </math> | |||
|width=50px align="right"|(32) | |||
|} | |||
<div class="NavFrame collapsed" style="text-align:left"> | <div class="NavFrame collapsed" style="text-align:left"> |
Revision as of 11:17, 23 August 2011
MARSSIM
The MARSSIM model is a landform evolution model primarily focuses on relatively long temporal scales (relative to the timescale for noticeable landform change) through fluvial and mass wasting processes.
Model introduction
The program is designed be computationally efficient such that individual runs can be done on a modern microcomputer in no more than a few tens of hours. The more recent additions to the model have focused on processes relevant to planetary landscapes, including lava flows, groundwater seepage and sapping, impact cratering, surface-normal accretion and ablation, and volatile redistribution by radiation-induced sublimation and recondensation. Individual process formulations vary from completely heuristic to modestly mechanistic. Important limitations for some potential applications are the assumption of a single representative bed material grain size in the fluvial system and no tracking of internal stratigraphy of sedimentary deposits.
Model parameters
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
- Local vertical rate of land surface elevation change
[math]\displaystyle{ {\frac{\partial y}{\partial t}} - U = {\frac{\partial z}{\partial t}} = - \nabla \ast q }[/math] (1)
- Weathering, Rain Slash, and Mass Movement
1) Potential erosion or deposition due to rain splash and regolith mass movement
[math]\displaystyle{ {\frac{\partial z}{\partial t}}|_{m} = -\nabla \ast q_{m} }[/math] (2)
2) The rate of movement
[math]\displaystyle{ q_{m} = [K_{s} \wp \left (S\right ) + K_{f}{\frac{1}{\left (1 - K_{x}|S|^a\right )}}] s }[/math] (3)
[math]\displaystyle{ \wp \left (S\right ) = sin \theta }[/math] (4)
A linear dependency
[math]\displaystyle{ \wp \left (S\right ) = tan \theta = S }[/math] (5)
- Fluvial Processes
1) Potential channel deposition or erosion
[math]\displaystyle{ {\frac{\partial z}{\partial t}}|_{c} = - \nabla \ast q_{s} = - {\frac{\partial q_{s}}{\partial x}} }[/math] (6)
2) Nonalluvial channels (A) Detachment capacity
[math]\displaystyle{ -{\frac{\partial z}{\partial t}}|_{c} = {\frac{\partial q_{s}}{\partial x}} = C_{d} = K_{t} \left (\tau - \tau_{c}\right ) }[/math] (7)
(B) equations for steady, uniform flow
[math]\displaystyle{ \tau = \gamma RS }[/math] (8)
[math]\displaystyle{ V = K_{n}R^\left ({\frac{2}{3}}\right ) S^\left ({\frac{1}{2}}\right ) / N }[/math] (9)
[math]\displaystyle{ Q_{w} = K_{p} RWV }[/math] (10)
(C) Simple power law equations
[math]\displaystyle{ Q_{w} = K_{a} A^e }[/math] (11)
[math]\displaystyle{ W = K_{w}Q_{w}^b = K_{w} K_{a}^b A^ \left (be\right ) }[/math] (12)
(D)
[math]\displaystyle{ {\frac{\partial z}{\partial t}}|_{c} = - K_{t} \left (K_{z} A^g S^h - \tau_{c}\right ) }[/math] (13)
[math]\displaystyle{ g = 0.6 e \left (1 - b \right ) }[/math] (14)
[math]\displaystyle{ h = 0.7 }[/math] (15)
[math]\displaystyle{ K_{z} = \gamma \{{\frac{N K_{a}^\left (1 - b\right )}{K_{p} K_{w} K_{n}}}^ \left ({\frac{3}{5}}\right ) }[/math] (16)
(E) Channel elevation change (a) The net elevation change of the channel when sufficiently large amount of regolith is delivered to the cell (the channel never erodes bedrock)
[math]\displaystyle{ {\frac{\partial z}{\partial t}} = {\frac{\partial z}{\partial t}}|_{m} + F \xi {\frac{\partial z}{\partial t}}|_{c} }[/math] (17)
[math]\displaystyle{ \xi = W / \delta }[/math] (18)
(b) The channel erosion when the channel is capable of eroding all regolith delivered to it
[math]\displaystyle{ \epsilon = {\frac{\partial z}{\partial t}}|_{m} \delta^2 - {\frac{\partial z}{\partial t}}\delta \left ( \delta - W \right ) }[/math] (19)
The fraction of time the channel spends eroding regolith
[math]\displaystyle{ \eta = {\frac{\epsilon}{F \xi \delta^2 {\frac{\partial z}{\partial t}}|_{c}}} }[/math] (20)
The net erosion into bedrock for the channel
[math]\displaystyle{ {\frac{\partial z}{\partial t}} = \left ( 1 - \eta \right ) {\frac{\partial z}{\partial c}}|_{c} = \beta \{{\frac{\partial z}{\partial t}}|_{m} + F \xi {\frac{\partial z}{\partial t}}|_{c} \} }[/math] (21)
[math]\displaystyle{ \beta = {\frac{1}{F \xi + 1 - \xi}} }[/math] (22)
(c) Channel erosion when the volume of regolith delivered to the cell per unit time is negative (the removal of regolith by channel erosion pertains only to the volume surrounding the stream and below the level of removal by slope erosion)
[math]\displaystyle{ \xi = \{{\frac{\partial z}{\partial t}}|_{m} - {\frac{\partial z}{\partial t}}\} \delta \left (\delta - W \right ) }[/math] (23)
The net erosion into bedrock
[math]\displaystyle{ \{{\frac{\partial z}{\partial t}}|_{m} - {\frac{\partial z}{\partial t}}\} = \left ( 1 - \eta \right ) {\frac{\partial z}{\partial t}}|_{c} }[/math] (24)
3) Alluvial channels (A) Potential rate of fluvial erosion
[math]\displaystyle{ {\frac{\partial z}{\partial t}}|_{c} = - \nabla \ast q_{sb} }[/math] (25)
[math]\displaystyle{ \phi = \Im \left ({\frac{1}{\Psi}}\right ) = {\frac{q_{sb}}{\omega d \left ( 1 - \mu \right )}} }[/math] (26)
[math]\displaystyle{ {\frac{1}{\Psi}} = {\frac{\tau}{\left ( \gamma_{s} - \gamma \right ) d }} }[/math] (27)
(B) Bed load or total bedload formulas
[math]\displaystyle{ \phi = K_{e} \{{\frac{1}{\Psi}} - {\frac{1}{\Psi_{c}}}\} ^p }[/math] (28)
(C) Total bed sediment discharge
[math]\displaystyle{ Q_{sb} = K_{q} A^ \left (eb \right ) [K_{v} A^ \left (0.6 e \left (1 - b \right ) \right ) S^ \left (0.7\right ) - 1/ \Psi_{c}]^ p }[/math] (29)
[math]\displaystyle{ K_{q} = {\frac{K_{e} \omega \left (1 - \mu \right ) K_{w} K_{a}^b}{\left ( \gamma_{s} - \gamma \right )^p d^ \left ( p - 1 \right )}} }[/math] (30)
[math]\displaystyle{ K_{v} = \gamma \{{\frac{N K_{a}^ \left (1-b\right) }{K_{w} K_{n} K_{p}}}\} ^ \left (0.6\right ) }[/math] (31)
(D) The net erosion of channel
[math]\displaystyle{ {\frac{\partial z}{\partial t}}|_{c} = {\frac{\left (Q_{sbi} - Q_{sbo}\right )}{LW}} }[/math] (32)
Symbol | Description | Unit |
---|---|---|
U | local tectonic uplift rate | - |
q | spatial divergence of the vector of eroded material flux | - |
\partial y / \partial t | local vertical rate of land surface elevation change | - |
\partial z / \partial t | erosion rate relative to a bedrock-fixed reference frame | - |
qm | spatial divergence of the vector rate of movement | - |
\wp (S) | an increasing function of slope gradient | - |
s | the unit vector in the direction of S | - |
S | value of local slope gradient | - |
Ks, Kx, Kf | constants assumed to be spatially and temporally invariant | - |
a | exponent assumed to be spatially and temporally invariant | - |
θ | sine of the slope angle | - |
c | potential channel deposition or erosion | - |
qs | volumetric sediment transport rate per unit channel width | - |
qsw | wash load (a part of sediment discharge, which is assumed to never be redeposited except in depressions) | - |
qsb | bed load (including bed load and suspended load), which is carried in capacity amounts if an alluvial bed is present | - |
τ | shear stress | - |
τc | critical shear stress | - |
C | actual detachment rate | - |
Cd | intrinsic detachment capacity | - |
γ | unit weight of water | - |
R | the hydraulic radius | - |
N | Manning's resistance coefficient | - |
Kp | a form factor close to unity | - |
Qw | hydraulic geometry for dominant discharge (Assumed to be equal for both bed erosion and sediment transport) | - |
W | channel width | - |
Kt | this constant includes both effects of substrate erodibility as well as magnitude of the dominant discharge | - |
τc | critical Shields number for the onset of sediment motion | - |
ψs | the fraction of bed shear stress | - |
qt * | Einstein number | - |
qt | volume sediment transport rate per unit width | L2 / T |
If | flood intermittency | - |
tf | cumulative time the river has been in flood | T |
Gt | the annual sediment yield | M / T |
ta | the number of seconds in a year | - |
Qf | sediment transport rate during flood discharge | L2 / T |
αt | dimensionless coefficient in the sediment transport equation, equals to 8 | - |
nt | exponent in sediment transport relation, equals to 1.5 | - |
τc * | reference Shields number in sediment transport relation, equals to 0.047 | |
CZ | dimensionless Chezy resistance coefficient. | |
Sl | initial bed slope of the river | - |
ηi | initial bed elevation | L |
τ | shear stress on bed surface | - |
qb | bed material load | M / T |
Δx | spatial step length, equals to L / M | L |
Qw | flood discharge | L3 / T |
Δt | time step | T |
Ntoprint | number of time steps to printout | - |
Nprint | number of printouts | - |
aU | upwinding coefficient (1=full upwind, 0.5=central difference) | - |
αs | coefficient in sediment transport relation | - |
u* | shear velocity | L / T |
αr | coefficient in Manning-Strickler resistance relation | - |
τb * | non-dimensional total shear stress | - |
Notes
Any notes, comments, you want to share with the user
Numerical scheme
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)
References
Key papers