Default Behavior#
As Cratermaker is designed to be easy to use, all of its component classes are built to be invoked without any arguments. When a Simulation object is created as in the examples above, a set of component objects are created with their default parameters. It is important to understand what the default behavior of these classes is, as this will affect the results of your simulation.
Target: The default target is the Moon. There are a number of known bodies that can be selected as targets, including some jovian satellites, and small bodies, but currently few of them have been tested.
Scaling: The default scaling model is called
"montecarlo". It is similar to the one used by the Cratered Terrain Evolution Model (CTEM) that is a progenitor to Cratermaker [1], but includes some variability in scaling parameters using a monte carlo method. The projectile to transient scaling model is mostly based on Holsapple (1993) [2] with some additional scaling parameters for ice given by Kraus et al. (2011) [3]. An additiona scaling model is called"ctem", which is identical to “montecarlo” but with the variability turned off.Production: There are two production function models available:
"neukum"and"powerlaw. The default is"neukum", and there are three versions available of this model: “Moon” [4], “Mars” [5], and “Projectile” [6]. The version will be selected based on the target body. When the target body is either “Moon” or “Mars”, then one of the two crater-based production functions are selected, and projectile sizes are determined by the Scaling model. For other bodies, “Projectile” is chosen, and the crater size is determined by the Scaling model.Projectile: There are three projectile velocity and density models available:
"asteroids","comets", and"generic". The default projectile model is determined by the target body. If the body is an inner solar system body, then"asteroids"is chosen, otherwise"comets"is chosen. The asteroid velocity distributions for the Moon are from Yue et al. (2013) [7], and for all other inner planets from Minton et al. (2010) [8]. Comet velocities are from Zahnle et al. (2003) [9].Surface: There are four surface models available:
"icosphere","arbitrary_resolution", and"hireslocal", and"datasurface". The default is"icosphere", which builds fast an efficient representation of a sphere. The resolution of the grid (the number of faces of the mesh) is determined by the formula \(20 \times 4^n\), where \(n\) is given by the argument “gridlevel” with a default value of 8.Morphology: Currently one morphology model is available:
"basicmoon". This is a model that similar to that used by CTEM. Most of the parameters are taken from Pike (1977) [10], except for simple crater profiles, which use a model from Fassett and Thomson (2014) [11]. Ejecta blanket scaling is from McGetchin et al. (1973) [12].Counting: Currently one counting model available:
"depthcount". This is a model that is similar to that used by CTEM, but with some modifications. It is based on a depth-to-diameter cutoff method used in Minton et al. (2019) [13] with corrections for complex craters given by Riedel et al. (2020) [14]
Because the Simulation class contains all other components, the defaults for all of the components can be viewed from printing a simulation object to the console, or by inspecting the cratermaker.yaml configuration file of a simulation with no arguments passed to it.
In [1]: import cratermaker as cm
In [2]: sim = cm.Simulation()
Creating a new grid
Generating a mesh with icosphere level 8.
In [3]: print(sim)
<Simulation>
simdir : .
<Current state>
Interval : 0
<Components>
Target: <Target: Moon>
Material: Soft Rock
Diameter: 3475.1 km
Mass: 7.35e+22 kg
Surface density: 2250.0 kg/m³
Transition Type: silicate
Escape Velocity: 2.376 km/s
Gravity: 1.624 m/s²
Surface: <IcosphereSurface>
Target: Moon
Grid File: surface/grid.nc
Number of faces: 655362
Number of nodes: 1310720
Grid Level: 8
Effective pixel size: 7.602 km +/- 319.2 m
Morphology: <BasicMoonMorphology>
Counting : <depthcount>
Do slope collapse: True
Do subpixel degradation: True
Production: <neukum>
Ejecta Truncation: Off
Ejecta Rays: False
Production: <NeukumProduction>
Generator type crater
Diameter range: 7.161 km to 3475.1 km
N(D) units: N>D(km) per 10⁶ km²
Number of quasi-Monte Carlo craters: 0
Version: Moon
Valid Time Range: 0 My - 4.500 Gy
Valid Diameter Range: 10.00 m - 1000.0 km
Projectile: <AsteroidProjectiles>
Sample from distributions: True
Mean Velocity: 22.10 km/s
Density: 2250.0 kg/m³
Target: Moon
Scaling: <MonteCarloScaling>
Target: Moon
Projectile: <asteroids>
Material: Soft Rock
K1: 0.200
mu: 0.550
Ybar: 7.600 MPa
Target density: 2250 kg/m³
Projectile density: 2250 kg/m³
Nominal simple-complex transition diameter: 12.19 km
Monte Carlo Scaling: True
Counting: <DepthCount>
Number of observed craters: 0
Number of emplaced craters: 0
Surface: <icosphere>