Material Collection: `lithospheric_deformation_eb1`

Material List

StickyAir EB1
StickyWater EB1
ClasticSediment EB1
Salt EB1
FelsicContinentalCrust EB1
FelsicContinentalCrustWeak EB1
FelsicContinentalCrustStrongZone EB1
MaficContinentalCrust EB1
MaficContinentalCrustStrongZone EB1
UltramaficContinentalLithosphere EB1
UltramaficContinentalLithosphereWeak EB1
UltramaficContinentalLithosphereStrongZone EB1
UltramaficAsthenosphereDryFertile EB1
UltramaficAsthenosphereWetFertile EB1
OceanicGabbroicCrust EB1
GabbroicMagma EB1
LayeredGabbroicCrust EB1
LayeredGabbroicMagma EB1
SerpentinizedPeridotite EB1

Material Summary

Collection File Path

src/material_library/registries/lithospheric_deformation/lithospheric_deformation_eb1.yaml

Mantle flow laws from HK03.
Base case crustal flow strength from GT95 and R06 (no scaling factor)
Moderate plastic yield strength (20/10; 0/0.1)
Moderate radiogenic heat production (1.8e-6/0.5e-6)

Upper Crustal Rheology (Moderate Strength)

Wet quartzite dislocation creep flow law parameters from Gleason and Tullis (1995). No scaling factor is applied to pre-exponential factor.

Wet anorthite dislocation creep flow law parameters from Rybacki et al. (2006) using a water fugacity of 500 MPa based on fit to Figure 8b at 30 km. No scaling factor is applied to pre-exponential factor.

Mantle Rheology

Dislocation and diffusion creep flow law parameters from Hirth and Kohlstedt (2003) for dry and wet olivine aggregates. A grain size of 6 mm is assumed for the olivine diffusion creep flow law. A water content of 500 ppm H/Si is assumed for the wet olivine aggregates.

Peierls creep law parameters are inferred from Katayama and Karato (2008). The current pre-exponential factor for Peierls creep was adjusted from 6.31e+7 (wet case) to 2.5e+7 approximate dry case and to ensure consistency with the dry dislocation creep law. A visual comparison to Figure 9b of Katayama and Karato (2008) was used to test the adjusted pre-exponential factor at strain rate 1e-5 (~1000 MPa) and 1e-10 (~100 MPa). Note that the activation energy and volume for Peierls creep comes from the current dislocation creep law, which requires adjustment to the Peierls pre-exponential term if different from Katayama and karato (2008).

Frictional Plastic Model

Cohesion: 20 to 10 MPa from plastic strain of 0.0 to 1.0.

Friction angles: 30 to 7 degrees from plastic strain of 0.0 to 1.0. The 7 degree minimum friction angle is from the model presented by Bos and Spiers (2002) based on laboratory experiments (see Figure 6 for the case with a sliding velocity of 1e-3 micrometers/second).

Radiogenic Heat Production

upper crust: 1.8 microW/m^3
lower crust: 0.5 microW/m^3

Material Tables

StickyAir EB1

Property Name	Value	Units	Description
`standard_density`	`1.0`	`kg/m^3`	`Density at standard conditions.`
`thermal_expansion`	`0.0`	`1/K`	`Thermal expansion.`
`compressibility`	`0`	`1/Pa`	`Compressibility.`
`flow_stype`	`Isoviscous`	`None`	`Flow law option.`
`viscosity_iso`	`1.0e18`	`Pa.s`	`Viscosity for isoviscous models.`
`shear_modulus`	`1.0e20`	`Pa`	`Shear modulus.`
`heat_capacity`	`100.0`	`J/kg/K`	`Heat capacity.`
`thermal_conductivity_ref`	`200.0`	`W/m/K`	`Reference thermal conductivity parameter k0.`
`thermal_conductivity_a`	`0.0`	`W/m`	`Reference thermal conductivity parameter a.`
`radiogenic_heat_production`	`0.0`	`W/m^3`	`Radiogenic heat production.`

StickyWater EB1

Property Name	Value	Units	Description
`standard_density`	`1000.0`	`kg/m^3`	`Density at standard conditions.`
`thermal_expansion`	`0.0`	`1/K`	`Thermal expansion.`
`compressibility`	`0.0`	`1/Pa`	`Compressibility.`
`melt_density`	`0.0`	`kg/m^3`	`Melt density.`
`flow_stype`	`Isoviscous`	`None`	`Flow law option.`
`viscosity_iso`	`1.0e18`	`Pa.s`	`Viscosity for isoviscous models.`
`shear_modulus`	`1.0e20`	`Pa`	`Shear modulus.`
`heat_capacity`	`100.0`	`J/kg/K`	`Heat capacity.`
`thermal_conductivity_ref`	`200.0`	`W/m/K`	`Reference thermal conductivity parameter k0.`
`thermal_conductivity_a`	`0.0`	`W/m`	`Reference thermal conductivity parameter a.`
`radiogenic_heat_production`	`0.0`	`W/m^3`	`Radiogenic heat production.`

ClasticSediment EB1

Property Name	Value	Units	Description
`standard_density`	`2710.0`	`kg/m^3`	`Density at standard conditions: Shale50%:Sandstone50%, typical rocks from Hantschel and Kauerauf (2010).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.8e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2400.0`	`kg/m^3`	`Melt density at 1 bar used in linearized Birch-Murnagham density model (Dacite from Lesher and Spera, 2015, Fig5.1).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`0.000857`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for wet quartzite (Gleason and Tullis, 1995).`
`stress_exponent_n_dc`	`4.0`	`None`	`Stress exponent for dislocation creep for wet quartzite (Gleason and Tullis, 1995).`
`activation_energy_dc`	`223.0`	`kJ/mol`	`Activation energy for dislocation creep for wet Quartzite (Gleason and Tullis, 1995).`
`activation_volume_dc`	`0.0`	`J/MPa/mol`	`Activation volume for dislocation creep.`
`shear_modulus`	`3.6e10`	`Pa`	`Shear modulus (Brune et al., 2014, Table S1).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`890.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.3`	`W/m/K`	`Reference matrix thermal conductivity @ 20C: Shale50%:Sandstone50%, typical rocks from Hantschel and Kauerauf (2010).`
`radiogenic_heat_production`	`1.4e-6`	`W/m^3`	`Radiogenic heat production: Shale50%:Sandstone50%, typical rocks from Hantschel and Kauerauf (2010).`
`stype_solidus`	`SedimentGerya2010`	`None`	`Solidus model`
`stype_liquidus`	`SedimentGerya2010`	`None`	`Liquidus model`
`latent_heat`	`300000.0`	`J/kg`	`Latent heat (Liao and Gerya, 2014).`
`porosity_initial`	`0.4`	`fraction`	`Initial porosity.`
`porosity_decay_depth`	`2500.0`	`m`	`Porosity decay depth used in exponential decay law.`

Salt EB1

Property Name	Value	Units	Description
`standard_density`	`2200.0`	`kg/m^3`	`Density at standard conditions.`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.8e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2200.0`	`kg/m^3`	`Melt density.`
`flow_stype`	`Isoviscous`	`None`	`Flow law option.`
`viscosity_iso`	`1.0e19`	`Pa.s`	`Viscosity for isoviscous models.`
`shear_modulus`	`3.6e10`	`Pa`	`Shear modulus (Brune et al., 2014, Table S1).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`860`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`6.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`2.0e-8`	`W/m^3`	`Radiogenic heat production.`

FelsicContinentalCrust EB1

Property Name	Value	Units	Description
`standard_density`	`2860.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.8e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2400.0`	`kg/m^3`	`Melt density at 1 bar used in linearized Birch-Murnagham density model (Dacite from Lesher and Spera, 2015, Fig5.1).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`0.000857`	`1/s/MPa^n`	`Pre-exponential factor for dislocation creep corrected for simple shear: Wet Quartzite, (Gleason and Tullis, 1995).`
`stress_exponent_n_dc`	`4.0`	`None`	`Stress exponent for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_energy_dc`	`223.0`	`kJ/mol`	`Activation energy for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_volume_dc`	`0.0`	`J/MPa/mol`	`Activation volume for dislocation creep.`
`shear_modulus`	`3.6e10`	`Pa`	`Shear modulus (Brune et al., 2014).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al., 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`1.8e-6`	`W/m^3`	`Radiogenic heat production (Hantschel and Kauerauf, 2010).`

FelsicContinentalCrustWeak EB1

Property Name	Value	Units	Description
`standard_density`	`2860.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.8e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2400.0`	`kg/m^3`	`Melt density at 1 bar used in linearized Birch-Murnagham density model (Dacite from Lesher and Spera, 2015, Fig5.1).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`0.000857`	`1/s/MPa^n`	`Pre-exponential factor for dislocation creep corrected for simple shear: Wet Quartzite, (Gleason and Tullis, 1995).`
`stress_exponent_n_dc`	`4.0`	`None`	`Stress exponent for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_energy_dc`	`223.0`	`kJ/mol`	`Activation energy for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_volume_dc`	`0.0`	`J/MPa/mol`	`Activation volume for dislocation creep.`
`shear_modulus`	`3.6e10`	`Pa`	`Shear modulus (Brune et al., 2014).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al., 2017).`
`friction_angle_initial`	`7.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`1.8e-6`	`W/m^3`	`Radiogenic heat production (Hantschel and Kauerauf, 2010).`

FelsicContinentalCrustStrongZone EB1

Property Name	Value	Units	Description
`standard_density`	`2860.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.8e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2400.0`	`kg/m^3`	`Melt density at 1 bar used in linearized Birch-Murnagham density model (Dacite from Lesher and Spera, 2015, Fig5.1).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`0.000857`	`1/s/MPa^n`	`Pre-exponential factor for dislocation creep corrected for simple shear: Wet Quartzite, (Gleason and Tullis, 1995).`
`stress_exponent_n_dc`	`4.0`	`None`	`Stress exponent for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_energy_dc`	`223.0`	`kJ/mol`	`Activation energy for dislocation creep (Wet Quartzite, Gleason and Tullis, 1995).`
`activation_volume_dc`	`0.0`	`J/MPa/mol`	`Activation volume for dislocation creep.`
`shear_modulus`	`3.6e10`	`Pa`	`Shear modulus (Brune et al., 2014).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`2.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al., 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`30.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`1.8e-6`	`W/m^3`	`Radiogenic heat production (Hantschel and Kauerauf, 2010).`

MaficContinentalCrust EB1

Property Name	Value	Units	Description
`standard_density`	`2900.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.6e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`3570.0`	`1/s/MPa^n`	`Tensor-invariant pre-exponential factor for dislocation creep for wet anorthite aggregate, water fugacity = 500 MPa (Rybacki et al., 2006).`
`stress_exponent_n_dc`	`3.0`	`None`	`Stress exponent for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_energy_dc`	`345.0`	`kJ/mol`	`Activation energy for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_volume_dc`	`38.0`	`J/MPa/mol`	`Activation volume for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`shear_modulus`	`4.0e10`	`Pa`	`Shear modulus.`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`5.0e-7`	`W/m^3`	`Radiogenic heat production (Hantschel and Kauerauf, 2010).`

MaficContinentalCrustStrongZone EB1

Property Name	Value	Units	Description
`standard_density`	`2900.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.6e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`3570.0`	`1/s/MPa^n`	`Tensor-invariant pre-exponential factor for dislocation creep for wet anorthite aggregate, water fugacity = 500 MPa (Rybacki et al., 2006).`
`stress_exponent_n_dc`	`3.0`	`None`	`Stress exponent for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_energy_dc`	`345.0`	`kJ/mol`	`Activation energy for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_volume_dc`	`38.0`	`J/MPa/mol`	`Activation volume for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`shear_modulus`	`4.0e10`	`Pa`	`Shear modulus.`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`2.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`30.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`radiogenic_heat_production`	`5.0e-7`	`W/m^3`	`Radiogenic heat production (Hantschel and Kauerauf, 2010).`

UltramaficContinentalLithosphere EB1

Property Name	Value	Units	Description
`standard_density`	`3291.0`	`kg/m^3`	`Density at standard conditions (Phanerozoic depleted model calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`651000.0`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`stress_exponent_n_dc`	`3.5`	`None`	`Stress exponent for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_energy_dc`	`530.0`	`kJ/mol`	`Activation energy for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_dc`	`13.0`	`J/MPa/mol`	`Activation volume for dislocation creep: Dry Olivine (Hirth and Kohlstedt, 2003).`
`pre_exponential_difc`	`0.0104`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for dry olivine aggregate using a grain size of 6 mm (Hirth and Kohlstedt, 2003).`
`activation_energy_difc`	`375.0`	`kJ/mol`	`Activation energy for diffusion creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_difc`	`6`	`J/MPa/mol`	`Activation volume for diffusion creep for dry olivine aggregate, (Hirth and Kohlstedt, 2003)`
`pre_exponential_pei`	`2.5e7`	`s^-m1*MPa^-m2`	`Pre-exponential factor for Peierls creep adjusted from 6.31e+7 for dry conditions and consistency with dislocation creep law (See Katayama and Karato 2008 for wet case)`
`stress_exponent_m1_pei`	`1.0`	`None`	`Stress exponent m1 for Peierls creep (Katayama and Karato, 2008).`
`stress_exponent_m2_pei`	`2.0`	`None`	`Stress exponent m2 for Peierls creep (Katayama and Karato, 2008).`
`peierls_stress`	`9100.0`	`MPa`	`Peierls stress under dry conditions (Katayama and Karato, 2008).`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al,, 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`4.0`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Dry peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

UltramaficContinentalLithosphereWeak EB1

Property Name	Value	Units	Description
`standard_density`	`3291.0`	`kg/m^3`	`Density at standard conditions (Phanerozoic depleted model calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`651000.0`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`stress_exponent_n_dc`	`3.5`	`None`	`Stress exponent for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_energy_dc`	`530.0`	`kJ/mol`	`Activation energy for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_dc`	`13.0`	`J/MPa/mol`	`Activation volume for dislocation creep: Dry Olivine (Hirth and Kohlstedt, 2003).`
`pre_exponential_difc`	`0.0104`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for dry olivine aggregate using a grain size of 6 mm (Hirth and Kohlstedt, 2003).`
`activation_energy_difc`	`375.0`	`kJ/mol`	`Activation energy for diffusion creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_difc`	`6`	`J/MPa/mol`	`Activation volume for diffusion creep for dry olivine aggregate, (Hirth and Kohlstedt, 2003)`
`pre_exponential_pei`	`2.5e7`	`s^-m1*MPa^-m2`	`Pre-exponential factor for Peierls creep adjusted from 6.31e+7 for dry conditions and consistency with dislocation creep law (See Katayama and Karato 2008 for wet case)`
`stress_exponent_m1_pei`	`1.0`	`None`	`Stress exponent m1 for Peierls creep (Katayama and Karato, 2008).`
`stress_exponent_m2_pei`	`2.0`	`None`	`Stress exponent m2 for Peierls creep (Katayama and Karato, 2008).`
`peierls_stress`	`9100.0`	`MPa`	`Peierls stress under dry conditions (Katayama and Karato, 2008).`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al,, 2017).`
`friction_angle_initial`	`7.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`4.0`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Dry peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

UltramaficContinentalLithosphereStrongZone EB1

Property Name	Value	Units	Description
`standard_density`	`3291.0`	`kg/m^3`	`Density at standard conditions (Phanerozoic depleted model calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`651000.0`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`stress_exponent_n_dc`	`3.5`	`None`	`Stress exponent for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_energy_dc`	`530.0`	`kJ/mol`	`Activation energy for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_dc`	`13.0`	`J/MPa/mol`	`Activation volume for dislocation creep: Dry Olivine (Hirth and Kohlstedt, 2003).`
`pre_exponential_difc`	`0.0104`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for dry olivine aggregate using a grain size of 6 mm (Hirth and Kohlstedt, 2003).`
`activation_energy_difc`	`375.0`	`kJ/mol`	`Activation energy for diffusion creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_difc`	`6`	`J/MPa/mol`	`Activation volume for diffusion creep for dry olivine aggregate, (Hirth and Kohlstedt, 2003)`
`pre_exponential_pei`	`2.5e7`	`s^-m1*MPa^-m2`	`Pre-exponential factor for Peierls creep adjusted from 6.31e+7 for dry conditions and consistency with dislocation creep law (See Katayama and Karato 2008 for wet case)`
`stress_exponent_m1_pei`	`1.0`	`None`	`Stress exponent m1 for Peierls creep (Katayama and Karato, 2008).`
`stress_exponent_m2_pei`	`2.0`	`None`	`Stress exponent m2 for Peierls creep (Katayama and Karato, 2008).`
`peierls_stress`	`9100.0`	`MPa`	`Peierls stress under dry conditions (Katayama and Karato, 2008).`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`2.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al,, 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`30.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`4.0`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

UltramaficAsthenosphereDryFertile EB1

Property Name	Value	Units	Description
`standard_density`	`3311.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`651000.0`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`stress_exponent_n_dc`	`3.5`	`None`	`Stress exponent for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_energy_dc`	`530.0`	`kJ/mol`	`Activation energy for dislocation creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_dc`	`13.0`	`J/MPa/mol`	`Activation volume for dislocation creep: Dry Olivine (Hirth and Kohlstedt, 2003).`
`pre_exponential_difc`	`0.0104`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for dry olivine aggregate using a grain size of 10 mm (Hirth and Kohlstedt, 2003).`
`activation_energy_difc`	`375.0`	`kJ/mol`	`Activation energy for diffusion creep for dry olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_difc`	`6`	`J/MPa/mol`	`Activation volume for diffusion creep for dry olivine aggregate, (Hirth and Kohlstedt, 2003)`
`pre_exponential_pei`	`2.5e7`	`s^-m1*MPa^-m2`	`Pre-exponential factor for Peierls creep adjusted from 6.31e+7 for dry conditions and consistency with dislocation creep law (See Katayama and Karato 2008 for wet case)`
`stress_exponent_m1_pei`	`1.0`	`None`	`Stress exponent m1 for Peierls creep (Katayama and Karato, 2008).`
`stress_exponent_m2_pei`	`2.0`	`None`	`Stress exponent m2 for Peierls creep (Katayama and Karato, 2008).`
`peierls_stress`	`9100.0`	`MPa`	`Peierls stress under dry conditions (Katayama and Karato, 2008).`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al,, 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`4.0`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Dry peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

UltramaficAsthenosphereWetFertile EB1

Property Name	Value	Units	Description
`standard_density`	`3311.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`924000.0`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for wet (C[OH] = 500 ppm H/Si) olivine aggregate (Hirth and Kohlstedt, 2003) .`
`stress_exponent_n_dc`	`3.5`	`None`	`Stress exponent for dislocation creep for wet (C[OH] = 500 ppm H/Si) olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_energy_dc`	`480`	`kJ/mol`	`Activation energy for dislocation creep for wet (C[OH] = 500 ppm H/Si) olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_dc`	`10`	`J/MPa/mol`	`Activation volume for dislocation creep for wet (C[OH] = 500 ppm H/Si) olivine aggregate (Hirth and Kohlstedt, 2003).`
`pre_exponential_difc`	`0.00347`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for wet (C[OH] = 500 ppm H/Si) olivine aggregate using a grain size of 6 mm (Hirth and Kohlstedt, 2003).`
`activation_energy_difc`	`335`	`kJ/mol`	`Activation energy for diffusion creep for wet (500 ppm H/Si) olivine aggregate (Hirth and Kohlstedt, 2003).`
`activation_volume_difc`	`4.0`	`J/MPa/mol`	`Activation volume for diffusion creep for wet (500 ppm H/Si) olivine aggregate, (Hirth and Kohlstedt, 2003)`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain (Naliboff et al., 2017).`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain (Naliboff et al,, 2017).`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain (Naliboff et al., 2017).`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain (Naliboff et al., 2017).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`4.0`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Dry peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

OceanicGabbroicCrust EB1

Property Name	Value	Units	Description
`standard_density`	`2900.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.6e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`3570.0`	`1/s/MPa^n`	`Tensor-invariant pre-exponential factor for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`stress_exponent_n_dc`	`3.0`	`None`	`Stress exponent for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_energy_dc`	`345.0`	`kJ/mol`	`Activation energy for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_volume_dc`	`38.0`	`J/MPa/mol`	`Activation volume for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`shear_modulus`	`4.0e10`	`Pa`	`Shear modulus.`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`GabbroGerya2010`	`None`	`Solidus model`
`stype_liquidus`	`GabbroGerya2010`	`None`	`Liquidus model`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

GabbroicMagma EB1

Property Name	Value	Units	Description
`standard_density`	`2750.0`	`kg/m^3`	`Density at standard conditions. Since magma is defined as a melt fraction of 1, the melt_density is used instead of standard density if melt properties are activated.`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.6e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Isoviscous`	`None`	`Flow law option.`
`viscosity_iso`	`1.0e18`	`Pa.s`	`Viscosity for isoviscous models.`
`shear_modulus`	`1.0e11`	`Pa`	`Shear modulus.`
`heat_capacity`	`800.0`	`J/kg/K`	`Heat capacity.`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter k0.`
`stype_solidus`	`GabbroGerya2010`	`None`	`Solidus model`
`stype_liquidus`	`GabbroGerya2010`	`None`	`Liquidus model`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

LayeredGabbroicCrust EB1

Property Name	Value	Units	Description
`standard_density`	`2900.0`	`kg/m^3`	`Density at standard conditions (Calibrated by Theunissen et al, 2022).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.6e-11`	`1/Pa`	`Compressibility (Brune et al., 2014).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`3570.0`	`1/s/MPa^n`	`Tensor-invariant pre-exponential factor for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`stress_exponent_n_dc`	`3.0`	`None`	`Stress exponent for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_energy_dc`	`345.0`	`kJ/mol`	`Activation energy for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`activation_volume_dc`	`38.0`	`J/MPa/mol`	`Activation volume for dislocation creep for wet anorthite aggregate (Rybacki et al., 2006).`
`shear_modulus`	`4.0e10`	`Pa`	`Shear modulus.`
`cohesion_initial`	`2.0e7`	`Pa`	`Cohesion at initial strain.`
`cohesion_final`	`1.0e7`	`Pa`	`Cohesion at final strain.`
`friction_angle_initial`	`30.0`	`Degrees`	`Friction angle at initial strain.`
`friction_angle_final`	`7.0`	`Degrees`	`Friction angle at final strain.`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle.`
`heat_capacity`	`800.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.5`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`GabbroGlacier`	`None`	`Solidus model`
`stype_liquidus`	`GabbroGlacier`	`None`	`Liquidus model`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

LayeredGabbroicMagma EB1

Property Name	Value	Units	Description
`standard_density`	`2750.0`	`kg/m^3`	`Density at standard conditions.`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`1.06e-11`	`1/Pa`	`Compressibility.`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Isoviscous`	`None`	`Flow law option.`
`viscosity_iso`	`1.0e18`	`Pa.s`	`Viscosity for isoviscous models.`
`shear_modulus`	`1.0e11`	`Pa`	`Shear modulus.`
`heat_capacity`	`800.0`	`J/kg/K`	`Heat capacity.`
`thermal_conductivity_ref`	`2.6`	`W/m/K`	`Reference thermal conductivity parameter k0.`
`stype_solidus`	`GabbroGlacier`	`None`	`Solidus model`
`stype_liquidus`	`GabbroGlacier`	`None`	`Liquidus model`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`

SerpentinizedPeridotite EB1

Property Name	Value	Units	Description
`standard_density`	`2900.0`	`kg/m^3`	`Density at standard conditions: Assumed similar density to oceanic crust and consistent with an average bulk serpentinization of ~50%. See Schwartz et al. (2001, Table 1) and Hantschel and Kauerauf, 2010, Table E.1).`
`thermal_expansion`	`3.0e-5`	`1/K`	`Thermal expansion.`
`compressibility`	`7.9e-12`	`1/Pa`	`Compressibility (Turcotte and Schubert. 2014, Table B.6, r=6151km).`
`melt_density`	`2750.0`	`kg/m^3`	`Melt density at 1 bar (1e-4 GPa) used in linearized Birch-Murnagham EOS density model (Lesher and Spera, 2015).`
`flow_stype`	`Composite`	`None`	`Flow law option.`
`pre_exponential_dc`	`1.23e9`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for dislocation creep for dry olivine aggregate (Korenaga and Karato, 2008) x 1000 to approximate strength reduction relative to anhydrous mantle (Schwartz et al., 2001).`
`stress_exponent_n_dc`	`4.94`	`None`	`Stress exponent for dislocation creep for dry olivine aggregate (Korenaga and Karato, 2008).`
`activation_energy_dc`	`610.0`	`kJ/mol`	`Activation energy for dislocation creep for dry olivine aggregate (Korenaga and Karato, 2008).`
`activation_volume_dc`	`13.0`	`J/MPa/mol`	`Activation volume for dislocation creep: Dry Olivine (Korenaga and Karato, 2008).`
`pre_exponential_difc`	`2.14e-7`	`1/s/MPa^n`	`Tensor invariant pre-exponential factor for diffusion creep for dry olivine aggregate using a grain size of 10 mm (Korenaga and Karato, 2008).`
`activation_energy_difc`	`261.0`	`kJ/mol`	`Activation energy for diffusion creep for dry olivine aggregate (Korenaga and Karato, 2008).`
`activation_volume_difc`	`6.0`	`J/MPa/mol`	`Activation volume for diffusion creep for dry olivine aggregate, (Korenaga and Karato, 2008)`
`shear_modulus`	`6.6e10`	`Pa`	`Shear modulus (Turcotte and Schubert, 2014, Table B.6 r=6151km).`
`cohesion_initial`	`1.0e7`	`Pa`	`Cohesion at initial strain (Factor of 2 reduction relative to peridotite, Escartin et al., 2001).`
`cohesion_final`	`5.0e6`	`Pa`	`Cohesion at final strain (Factor of 2 reduction relative to peridotite, Escartin et al., 2001).`
`friction_angle_initial`	`10.0`	`Degrees`	`Friction angle at initial strain (Factor of 2 reduction relative to peridotite, Escartin et al., 2001).`
`friction_angle_final`	`5.0`	`Degrees`	`Friction angle at final strain (Factor of 2 reduction relative to peridotite, Escartin et al., 2001).`
`strain_initial`	`0.0`	`None`	`Initial strain used for cohesion and friction angle.`
`strain_final`	`1.0`	`None`	`Final strain used for cohesion and friction angle (Naliboff et al., 2017).`
`heat_capacity`	`785.0`	`J/kg/K`	`Reference heat capacity in Waples model (Hantschel and Kauerauf, 2010).`
`thermal_conductivity_ref`	`2.6`	`W/m/K`	`Reference thermal conductivity parameter for use in Sekiguchi-Waples model (Hantschel and Kauerauf, 2010).`
`stype_solidus`	`PeridotiteKatz2003`	`None`	`Dry peridotite solidus model from Katz et al. (2003).`
`stype_liquidus`	`PeridotiteKatz2003`	`None`	`Peridotite Liquidus model from Katz et al. (2003).`
`latent_heat`	`400000.0`	`J/kg`	`Latent heat (Turcotte and Schubert, 2014).`