Dip change in segment

Dip change in segment#

What we have achieved in the last section was already great, but slabs in the Earth are usually not thought of as straight lines. For example, near the surface they start with a dip angle of zero and then increase their dip with depth to maybe 60 degrees. So, how do we do that in the world builder?

This is where the array for the angle comes in. Last time we provided just one value, but we are also allowed to provide two values. If you provide two values, the first value is the dip angle at the top of the segment and the second value is the dip angle at the bottom of the segment. So to achieve this we only need to add a zero before the 60 in the angle array.

43    { 
44      "model":"subducting plate", "name":"Slab", "dip point":[0,0],
45      "coordinates":[[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
46      "segments":[{"length":300e3, "thickness":[100e3], "angle":[0,60]}],
47      "temperature models":[{"model":"uniform", "temperature":400 }],
48      "composition models":[{"model":"uniform", "compositions":[2]}]
49    }
BST_11_dip_change_in_segment.wb
BST_11_dip_change_in_segment.grid
Paraview v3 state file
 1{
 2  "version": "1.0",
 3  "coordinate system":{"model":"cartesian"},
 4  "features":
 5  [
 6    {
 7       "model":"mantle layer", "name":"upper mantle", "max depth":660e3,
 8       "coordinates":[[-1000e3,0],[-1000e3,1000e3],[3000e3,1000e3],[3000e3,0]],
 9       "composition models":[{"model":"uniform", "compositions":[4]}]
10    },
11    {
12       "model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3, 
13       "coordinates":[[0,0],[0,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
14       "temperature models":
15         [
16           {"model":"half space model", "max depth":100e3, "spreading velocity":0.04, 
17            "ridge coordinates":[[[400e3,-1],[-100e3,2000e3]]]}
18         ],
19       "composition models":[{"model":"uniform", "compositions":[0], "max depth":50e3}]
20    },
21    {
22       "model":"continental plate", "name":"Passive margin", "max depth":[[200e3]], 
23       "coordinates":[[-1000e3,0],[-1000e3,1000e3],[0,1000e3],[0,0]],
24       "temperature models":
25         [
26           {"model":"linear", "max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
27         ],
28       "composition models":
29         [
30           {"model":"uniform", "compositions":[3], 
31            "max depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]},
32           {"model":"uniform", "compositions":[1], 
33            "min depth":[[100e3],[200e3,[[-250e3,0],[-750e3,1000e3]]]]}
34         ]
35    },
36    {
37      "model":"oceanic plate", "name":"Subducting Oceanic plate", "max depth":100e3, 
38      "coordinates":[[2000e3,0],[2000e3,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
39      "temperature models":[{"model":"linear", "max depth":100e3}],
40       "composition models":[{"model":"uniform", "compositions":[3], "max depth":50e3},
41                             {"model":"uniform", "compositions":[1], "min depth":50e3}]
42    },
43    { 
44      "model":"subducting plate", "name":"Slab", "dip point":[0,0],
45      "coordinates":[[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
46      "segments":[{"length":300e3, "thickness":[100e3], "angle":[0,60]}],
47      "temperature models":[{"model":"uniform", "temperature":400 }],
48      "composition models":[{"model":"uniform", "compositions":[2]}]
49    }
50  ]
51}
BST_11_dip_change_in_segment.wb
BST_11_dip_change_in_segment.grid
Paraview v3 state file

Note

You can do the same for the thickness. If we want to start with a thickness of 100km and linearly taper it down to 0 km to form a slab tip, we can do the following: "thickness":[100e3,0]

Todo

The explanation can be significantly improved by adding conceptual figures

Basic Starter Tutorial section 11.

Basic Starter Tutorial section 11. The top part of the figure shows where the composition as been assigned as an object. Currently is shows composition 0 as green, composition 1 as yellow, composition 2 as purple and composition 3 as blue. Composition 4 is not shown to be able to see the slab. The front half of the overriding plate (composition 1) has also been removed to be able to better view the slab. The bottom part shows the temperature as seen slightly from below where only temperatures between 300K and 1600K are shown.#