Model help:MARSSIM

From CSDMS
Revision as of 15:30, 22 August 2011 by Happyfei (talk | contribs)
The CSDMS Help System

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

Parameter Description Unit
Input directory - -
MARSSIM application title - -
Variable information filename - -
Parameter Description Unit
Random seed an integer random seed -
Switch for new simulation If set to 1, start a new simulation. Otherwise, continue from a saved state -
Switch for fluvial and slope modeling If set to 1, fluvial and slope processes are modeled -
Switch for impact crater modeling If set to 1, impact craters are modeled -
Switch for lava inundation If set to 1, lava inundation is modeled -
Switch for eolian deposition If set to 1, eolian transport is modeled -
Switch for ocean processes If set to 1, ocean processes is modeled -
Switch for accretion and ablation If set to 1, model accretion and ablation -
Switch for lake evaporation If set to 1, model runoff and evaporation with lakes -
Number of columns for DEM - -
Number of rows for DEM - -
Number of vertical cells - -
Grid cell size Grid cell size arbitrary units. See conversion factor below -
Vertical scale units Elevation units (vertical). Must as same as grid cell size units (horizontal) -
Convert-to-meters factor Factor to convert grid cell size and vertical scale units to meters -
Switch for variable erosion rate If set to 1, erosion rate will vary in time in the specified manner. Otherwise, model parameters are constant -
Boundary lowering rate Rate at which the boundary is lowering over time (meter/year). To be used with a nonperiodic vertical boundary -
Channel timestep scale factor channel process timestep scale factor -
Default channel time step channel process timestep years
Maximum timestep Maximum allowed timestep years
Minimum timestep Minimum allowed timestep years
Sediment timestep scale factor - -
Mass wasting timestep scale factor - -
I-coordinate for center of debug window Used for debug mode -
J-coordinate for center of debug window Used for debug mode -
Size of debug window, I-direction Used for debug mode -
Size of debug window, J-direction Used for debug mode -
Switch to print debug info If set to 1, print info in specified debug window, Affects performance -
Switch to calculate morphometric statistics If set to 1, calculate morphometric statistics -
Switch to use horizontal lower boundary If set to 1, southern (lower ) boundary is level. Do not use with a periodic lower boundary -
Switch to use non-eroding lower boundary If set to 1, use a non-eroding lower/southern boundary can be used in conjunction with a horizontal lower boundary -
Switch for periodic Y axis If set to 1, use periodic boundary condition along the Y-axis -
Switch for periodic X axis If set to 1, use periodic boundary condition along the X-axis -
Switch to allow flow across all boundaries All boundaries are potential flow exits and are non-eroding -
Parameter Description Unit
Starting iteration Generally this is zero for new simulatins -
Present time Generally this is zero for new simulatins years
Maximum iteration The number of timesteps or iterations to be used in the simulation -
Maximum simulation time The maximum time that the simulation is allowed to run years
Iterations between prints The number of iterations (or timesteps) between writing out a lot of info -
Simulation progress interval How often to summarize simulation progress. Generally should be the same as elevinterval -
Recalculate gradient interval Number of interations before recalculating gradients and contributing areas (Currently hardwired) -
Write change interval Number of iterations before printing out information on the rate of change in gradients and flow directions (Currently hardwired) -
Morphometric channel threshold This has to do with the definition of where channels start when printing out morphometric statistics -
Matrices to output What matrices to output at Elevation_print_interval times. Each digit represents a flags are (1) elevation; (2) relative elevation; (3) D8 gradient; (4) erode_slope; (5) cfne; (6) cfnw; (7) cfw; (8) cfn; (9) cfne; (10) drainage_area; (11) relative_gradient; (12) regolith and relative resistance. Generally, only first flag is set (elevation) -
Switch for absolute or relative elevations If set to 1, write out absolute vs. relative elevations -
Interval to print out shaded relief images Shaded relief images are typically generated and saved every 50 to 100 iterations -
Divergence Interval - -
Parameter Description Unit
Switch for wet or dry depressions If set to 1, all depressions fill with water and overflow. Otherwise, all depressions are infinite sinks for water. Most Mars simulations have used 0, but some simulatins with non-periodic BCs use 1 -
Switch for rescale discharges IF set to 1, then effective discharges are rescaled as a power function of contributing area (e.g., Q = K*A**E, where E is typically in [0.5,1.0]) -
Discharge constant Used with option to rescale discharges. In that case, Discharge = discharge_constant*Area**discharge_exponent m3 / s
Discharge exponent Used with option to rescale discharges. In that case, Discharge = discharge_constant*Area**discharge_exponent m3 / s
Area factor Generally not used -
Rain depth - -
Rain standard deviation - -
Switch for divergence-dependent runoff If set to 1, runoff depth functionally depends on topographic divergence (VARYIELDUSE; 0, or 1), usually not set -
Mean convergence runoff Only used if 'Divergence dependent runoff' is toggled on -
High convergence runoff Only used if 'Divergence dependent runoff' is toggled on -
Low convergence runoff Only used if 'Divergence dependent runoff' is toggled on -
Iterations between recalculating evaporation rates Only used if 'Lake evaporation' is toggled on -
Iterations between stocahstic changes in evaporation rates Only used if 'Lake evaporation' is toggled on -
Mean yearly evaporation depth - meters
Std. deviation of mean evaporation depth - meters
Switch to model pelagic deposition If set to 1, deposition of suspended sediment in enclosed basins is modeled (IMODEL_PELAGIC_DEPOSITION, 0 or 1) -
Wash load fraction Fraction of sediment delivered to basins that is deposited as pelagic sediment -
Amount of post-deposition A diffusivity that determines how much post-deposition diffusion occurs in recently-deposited pelagic sediment -
Switch to model influent rivers IF set to 1, incoming rivers are modeled -
Parameter Description Unit
Switch to model fluvial detachment This should only set to 0 when slope processes are being modeled and fluvial erosion, transport and deposition are not -
Switch for explicit or implicit channel bed If set to 1, use EXPLICIT method to model regolith depth (including bare rock). Otherwise use IMPLICIT method -
Vegetation upland resistance Only used if 'variable vegetation resistance' flag is set -
Vegetation channel resistance Only used if 'variable vegetation resistance' flag is set -
Vegetation area minimum Only used if 'variable vegetation resistance' flag is set -
Vegetation area maximum Only used if 'variable vegetation resistance' flag is set -
Parameter Description Unit
Switch to write sediment diagnostics If set to 1, detailed sediment routing information is printed (ISEDDEBUGT, 0 or 1) -
Switch for sediment transport If set to 1, alluvial channels, delta and fans are modeled (ISEDIMENT, 0 or 1) -
Switch for sediment routing If set to 1, sediment routing is modeled (ISEDIMENT, 0 or 1) -
Switch for sediment diffusion If set to 1, sediment diffusion is modeled (ISEDIMENT, 0 or 1) -
Switch for no-flux lower boundary If set to 1, a no-flux lower boundary is used (IREFLECT, 0 or 1) -
Switch for alluvial smoothing If set to 1, alluvial depositional surfaces (fans and deltas) are smoothed during deposition (SMOOTHSEDUSE, 0 or 1) -
Alluvial smoothing factor only used if 'alluvial rerouting' switch is used -
Switch to read alluvial thickness at start If set to 1, alluvium thickness is read at beginning of simulation (IREADALLUV, 0 or 1), generally 0 -
Switch to use an ocean If set to 1, an ocean is used, as when modeling delta (WATERUSE, 0 or 1), set to 1 if USEWET is set to 1 -
Ocean elevation If modeling a global ocean, set this to your assumed ocean level. Otherwise, set it to a very low value -
Delta forest gradient This is the forest bed gradient -
Sediment 1 exponent Used in sediment transport relationship -
Sediment 2 exponent Used in sediment transport relationship -
Effective discharge ratio Used in sediment transport relationship -
Sediment gradient exponent Used in sediment transport relationship -
Sediment transport exponent Used in sediment transport relationship -
Transport factor Generally, 40.0 for sand bed, 8.0 for gravel bed -
Flow fraction Set to 1 when modeling erosion by groundwater seepage -
Sediment porosity - -
Sediment specific gravity - -
Grain size Typical value are 0.0002 for sand and 0.02 for gravel meters
Transport critical dim shear Generally not changed -
Switch to use sticky sediment routing If set to 1, use 'sticky' VS. normal sediment. Used in an experiment to model birdfoot delta. Generally not used -
Sticky routing critical value Use in an experiment to model birdfoot delta. Generally, not used -
Bedload fraction Fraction of eroded sediment that is transported as bedload. Used for calculation of alluvial channel gradients -
Parameter Description Unit
Rock weathering rate Rate of bare rock weathering meters/year
Weather decay rate Determines how rapidly weathering rate decreases with increasing regolith thickness, based on terrestrial values -
Initial regolith thickness - meters
Switch for two-term weathering If set to 1, a humped weathering function is used (WEATHER2USE, 0 or 1) -
Weathering term 2 only used if 'two term weathering' flag is set -
Weathering decay 2 only used if 'two term weathering' flag is set -
Switch for seepage weathering If set to 1, seepage-dependent weathering is modeled (ISSEPAGEWEATHER, 0 or 1) -
Seepage weathering scale factor only used if 'Seepage weathering' flag is set -
Seepage weathering exponent only used if 'Seepage weathering' flag is set -
Critical bedrock gradient Maximum steepness of bedrock slope -
Weather multiplier This scales the rate of bedrock mass wasting -
Weather divergence If greater than 0, then rate of bedrock weathering will depend on topographic divergence -
Switch to read initial regolith thickness If set to 1, read initial regolith thickness from file (READREOLITH, 0 or 1) -
Parameter Description Unit
Switch for slope erosion modelling If set to 1, use slope erosion modelling (ISLOPEUSE, 0 or 1) -
Slope diffusity Determines relative rate of downslope creep (m2/yr) try values between 0.0002 and 0.02 -
Proportion of slope constant in alluvial creep Relative diffusivity of alluvium mass wasting to slope mass wasting. Generally set to 1 -
Switch to use critical slope gradient If use to 1, use a critical slope gradient (Critical_gradient_use, 0 or 1) -
Switch to use Roering mass wasting If set to 1, use the Roering formula for mass wasting. Otherwise, use the Howard formula. The Roering formula allows a threshold maximum slope gradient in regolith. The Howard formula uses linear creep -
Critical slope gradient Maximum stable regolith gradient -
Slope gradient exponent Used for Howard two-term relationship -
Parameter Description Unit
Switch to use bistable fluvial erosion If set to 1, use bistable fluvial erosion (HIGHRATEUSE, 0 or 1) used for modelling gully development -
Switch to use bistable bedrock If set to 1, use bistable bedrock (BEDROCKHIGH, 0 or 1) used for modeling gully development -
Low erosion threshold Only used with bistable erosion and gully development -
High erosion threshold Only used with bistable erosion and gully development -
Erosion rate change lag Only used with bistable erosion and gully development -
Bistable critical shear Only used with bistable erosion and gully development -
Bistable runoff factor Only used with bistable erosion and gully development -
Bistable bedrock erodibility Only used with bistable erosion and gully development -
Parameter Description Unit
Switch to use random critical shear stress If set to 1, use random critical shear stress (RANDTHRESHUSE, 0 or 1) -
Switch to use random discharge If set to 1, use random discharge (RANDDISCHUSE, 0 or 1) -
Critical shear variablity Scale factor for random critical shear stress -
Discharge coeff variation Standard error of discharge, assuming a lognormal distribution -
Omega weight Governs inherientance from past value, similar to Hurst parameter set to 1 for no inheritance -
Parameter Description Unit
Switch to use 3D variable resistance weathering If set to 1, use 3D variable resistance weathering (SLOPEVARUSE, 0 or 1) -
Diffusivity variability - -
Switch to use a resistant surface layer If set to 1, use a resistant surface layer (CRUSTUSE, 0 or 1). Experimental, for simulating the effect of development of a "duricrust" -
Surface layer thickness Crust thickness, using the same units as cell size -
Surface layer resistance Crust relative shear stress factor, greater than 1 for a resistant crust -
Vertical resistance scale factor Ratio of horizontal to vertical scale's (1.0 if equal). Only used for 3D varing rock resistance -
Parameter Description Unit
Switch to use rock deformation DEFORMUSE, 0 or 1, not implemented in current version -
Deformation scale factor Not implemented in current version -
Parameter Description Unit
Switch to model groundwater seepage If set to 1, groundwater seepage is modeled (SEEPUSE, 0 or 1) -
Seepage iteration interval Number of iterations between recalculation of groundwater flow -
Switch for permeability rescaling If set to 1, permeability in exponential vertical decay of permeability is scaled to the present land surface elevation. Otherwise it is scaled to the original land surface elevation (at the start of the simulation) -
Switch to show groundwater calculations If set to 1, show groundwater calculations -
Switch for exponential permeability decay If set to 1, permeability is constant through an aquifer of constant thickness. Otherwise, it decays with depth below the surface -
Yearly recharge rate - meters/year
Viscosity of water - -
Permeability, in darcies Permeability darcies
Groundwater depth scale For exponential permeability, this is the depth to the half value of permeability. For constant permeability, it is the acquiter thickness meters
Groundwater flow fraction The fraction of surface water flow that becomes groundwater flow, set to 1 for all groundwater flow -
Initial groundwater depth Initial depth beneath the land surface for calculation of steady state groundwater surface meters
Depth decay of permeability power Set this to 1 for normal exponential decay -
Max. groundwater iterations Max number of iterations, permitted in calculating steady state groundwater table -
Max. groundwater error Max residual error permitted in calculating steady state groundwater table -
Groundwater relaxation coeff Coefficient used for SOR (Successive Over Relaxation) -
Switch to use groundwater flux If set to 0, the groundwater flow term used for seepage calculation of surface flows and weathering rate is the groundwater flux divergence (Seepage rate to surface). Otherwise, it is the groundwater flow rate per unit aquifer width -
Switch to use seepage averaging If set to 1, the groundwater flow term used in further calculations is a 9 point average value of raw calculated value -
Parameter Description Unit
Eolian event probability Probability of an eolian deposition/erosion event per unit time (year or iteration) -
Eolian time increment Scales the overall Eolian Erosion/Deposition rate -
Switch to use total exposure If set to 1, all cells within the calculation window are used to compute the exposure index. Otherwise, only cells visible to the local cell are used -
Switch to use default eolian process If set to 1, eolian erosion and deposition occur normal to the topographic surface. Otherwise, erosion and deposition are modeled as vertical additions or subtractions -
Minimum eolian deposition rate This can be negative (If eolian erosion of exposed locations) -
Maximum eolian deposition rate Generally unity and scaled by Eolian_time_increment determines Eolian_constant_1 and Eolian_constant_2 -
Eolian parameter code Code that determines how next 2 parameters are used: (1) Exposure_10_percent, Exposure_90_percent; (2) Exposure_50_percent, Exposure_90_percent; (3) Zero_percent_exposure, Exposure_90_percent; (4) Rate0, exposure_50_percent -
Exposure parameter 1 Used to determine Eolian_constant_3 that defines the shape of eolian erosion/deposition curve as related to the exposure index -
Exposure parameter 2 Used to determine Eolian_constant_3 that defines the shape of eolian erosion/deposition curve as related to the exposure index -
Distance decay factor In calculation of exposure index, this determines how rapidly the weighting of surrounding cells decays with distance. Also used for accretion/ablation modeling -
Weighting calculation distance The maximum distance (in number of cells) that elevations are used to calculate the exposure index. Also used for accretion/ablation modeling -
Parameter Description Unit
Number of lava sources Number of sources (vents, volcanoes) that are present on the surface (Multiple sources not yet supported) -
Lava source I coordinate I coordinate column of lava source in DEM -
Lava source J coordinate J coordinate column of lava source in DEM -
Lava event probability Probability of a lava flow event during a single iteration -
Minimum lava flow slope Minimum gradient for laval flow at the flow source -
Lava flow thickness Assumed thickness of individual lava flow deposits meters
Minimum lava flow thickness Minimum thickness of a lava flow that can flow into adjoining cells meters
New segment interval Number of interations before starting a new lava flow source -
Source segment interval Number of iterations between changeover to a different lava flow source -
Eruption stop probability Probability per iteration, that the existing flow will solidify and stop being active. If this happens, a new flow starts at the source -
No flow probability The lower limit of probability for a cell to be a source for a new flow segment. If the probability drops below this value then the cell is no longer considered to be a possible flow source -
Lava gradient weight Determines how much the gradient between the edge of a flow and the neighboring point determines the probability of flow in that direction -
Lava duration weight Determine how rapidly a new cell diminishes in probability that it can be the source of a flow into a neighboring cell -
Parameter Description Unit
Impact probability probability of an impact event, per year -
Switch for ejecta wraparound If set to 1 and if the domain is periodic in both X and Y, then ejecta deposition can carry over onto the opposite side (IFOLD, 0 or 1) -
Switch for regolith crater If set to 1, then crater slopes and ejecta are soft (regolith). Otherwise they are considered to be initially bedrock -
Large crater depth scale Based upon scaling of fresh impact craters in Forsberg et al. (2004) and Howard (2007) meters
Large crater depth exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Large crater rim scale Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) meters
Large crater rim exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007)
Transition diameter Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) meters
Small crater depth scale Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) meters
Small crater depth exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Small crater rim scale Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) meters
Small crater rim exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Large crater shape scale Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Large crater shape exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Small crater shape scale Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Small crater shape exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Crater frequency exponent Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Frequency cutoff scaling Based on scaling of fresh impact craters in Forsberg et al (2004) and Howard (2007) -
Smallest possible crater Depends on the scale of the simulation, generally in meters, Generally keep at about (4*cellscale) -
Smallest modelled crater Depends on the scale of the simulation, generally in meters, Generally keep at about (4*cellscale) -
Largest modeled crater Depends on the scale of the simulation, generally in meters. Generally keep at about (cell_scale*ncols/2) -
Ejecta thickness variability - -
Noisesd Standard deviation of random noise used for ejecta thickness -
Inheritance parameter - -
Maximum rim gradient Maximum gradient of crater rim -
Parameter Description Unit
Switch for variable ocean elevation If set to 1, ocean elevation varies in time and the times and corresponding sea levels are read in from the file "oceanlevels.dat" -
Parameter Description Unit
Accretion rate Rate of non-fluvial and non-eolian surface accretion and degradation -
Switch for exposure dependent creep If set to 1, the mass wasting creep rate depends on exposure -
Switch for solar radiation If set to 1, sublimation due to reflected radiation is modeled -
Switch for top exposure - -
Switch for inverse exposure - -
Switch for exposure smoothing - -
Radiation constant This scales the rate of radiation-dependent sublimation/deposition -
Radiation dust factor Determines the relative amount of re-emitted thermal radiation from dust-covered surface relative to bedrock surface -
Radiation threshold convexity Critical value of exposure index for redeposition of sublimated ice -
Radiation deposit rate Rate scaling for redeposition of sublimated ice on less-exposed surfaces -


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}\right )|S|^a}}] 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

[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)

3) 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)

4) 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)

5)

[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)