Optional adding mantle layer#
Sometimes you want precise control over what is happening in the mantle. In this case we want to have a specific composition for the mantle.
Note
This is not always needed. For example in the visco-plastic material model in the FEM ASPECT, if no composition is defined, it is assumed to be a background composition. If you don’t define a composition in an area in the world builder, the world builder will tell anyone who wants to know the composition at that location that the fraction of every composition is zero there.
We are going to add the mantle layer as the first object in the feature list, so that other features can override it if necessary. The model name for a mantle layer is mantle layer, and we will name it “upper mantle”. Because we already have used compositions 0,1 and 3, and will use 2 later, we are going to set this to composition 4. We could also specify a different temperature model for this. The default is adiabatic, but you can easily set it to something else, or change the parameters of the adiabatic gradient.
Note
The potential mantle temperature, the thermal expansion coefficient and the specific heat here have a default value of -1. This is, like we have seen before, a special value. In this case it means that the global value is used (as you can read in the documentation of the value). It is strongly recommended to only set the global values so that the whole setup remains self-consistent, but if you do want to change it in individual locations, you can!
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 },
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}
Basic Starter Tutorial section 9. The top part of the figure shows where the composition as been assigned as an object. Currently it shows composition 0 as green, composition 1 as yellow, composition 3 as blue and composition 4 is shown as red. The bottom part shows the temperature as seen slightly from below where only temperatures between 300K and 1600K are shown. We will not show composition 4 in future figures anymore to be able to show other features more clearly.#