LithoStrongZones

The LithoStrongZones geometry is composed of a central upper weak zone in the lithosphere surrounded by strong lithospheric zones.

Initialization

EarthBox.MaterialGeometry.LithoStrongZones.initialize! — Function

Initialize lithosphere strong zones geometry.

This material geometry involves a central lithospheric zone with different frictional plastic properties than the surrounding lithosphere. The left and right boundaries of the central lithospheric zone are controlled by the parameters x_left_strong and x_right_strong, respectively.

Material properties of the central and lateral zones are defined using the material domains for the upper continental crust strong zone, lower continental crust strong zone, and lithospheric mantle strong zone in either a material.yaml file or Material.jl script.

Keyword arguments:

  • iuse_strong_zones::Int:
    • Use strong zones (0=off, 1=on)
  • x_left_strong::Union{Float64, Nothing}:
    • x-coordinate in meters of the left edge of the lithospheric strong zone
  • x_right_strong::Union{Float64, Nothing}:
    • x-coordinate in meters of the right edge of the lithospheric strong zone
source

Initialization Scenarios

Initialization involves defining the left and right boundaries of the central lithospheric zone, x_left_strong and x_right_strong, respectively. Also, the flag iuse_strong_zones argument must be set to 1 to enable strong zones.

Initialization of these parameters depends on whether T-type refinement parameters xo_highres and xf_highres were previously defined during model initialization with EarthBoxState or during staggered grid initialization with StaggeredGrid.initialize!.

If T-type refinement parameters xo_highres and xf_highres were previously defined then x_left_strong and x_right_strong of the central lithospheric zone will be set equal to xo_highres and xf_highres, respectively.

If T-type refinement parameters xo_highres and xf_highres were not previously defined then x_left_strong and x_right_strong will be set to equal to the values provided as arguments to initialize!.

Examples

Example 1: Direct input of strong zone boundaries

eb = EarthBoxState(xnum = 100, ynum = 100, xsize = 500000.0, ysize = 120000.0)
model = eb.model_manager.model # ModelData container
EarthBox.MaterialGeometry.LithoStrongZones.initialize!(
    model, 
    iuse_strong_zones = 1,
    x_left_strong = 200000.0,
    x_right_strong = 300000.0
)

Example 2: Using T-type refinement parameters to set strong zone boundaries

ttype_refinement_parameters = Dict(
    "xo_highres" => 200000.0,
    "xf_highres" => 300000.0,
    "yf_highres" => 50000.0,
    "dx_highres" => 500.0,
    "dy_highres" => 500.0,
    "dx_lowres" => 2000.0,
    "dy_lowres" => 2000.0
)
eb = EarthBoxState(
    xsize = 500000.0, ysize = 120000.0, 
    ttype_refinement_parameters = ttype_refinement_parameters
)
model = eb.model_manager.model
EarthBox.MaterialGeometry.LithoStrongZones.initialize!(model, iuse_strong_zones = 1)