Oceanic plate temperature

Oceanic plate temperature#

Using uniform temperatures will only allow you to build very simple models. So let’s make the model a bit more interesting by adding a half space cooling model to the overriding plate. This allows us to add a temperature structure based on the distance along the plate from a ridge.

To be able to do this, we need to provide the following information to the model: model, ridge coordinates, spreading velocity and max depth. Please see /features/items/oneOf/4/temperature models/items/oneOf/2 for more info.

We are already familiar with max depth, but spreading velocity and ridge coordinates are new. The spreading velocity is simply the velocity of how fast one part of the plate is moving away from the ridge. The ridge coordinates are a list of spreading ridges, where each ridge is a list of 2D points. The half space is computed with the distance to the closest point on any ridge and the spreading velocity.

If you are familiar with the half space cooling model, you would also think some kind of top and bottom/mantle temperature is required. These values can be provided, but by default they are set to reasonable values. The top temperature is set to 293.15K and the bottom temperature is set to -1K. Yes, that is negative 1 Kelvin! For this model, if the temperature is set to a negative value, it automatically uses the computed adiabatic temperature at that depth. This will be good enough for our tutorial, so we won’t be changing those values here, but you are free to do so!

 6    {
 7       "model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3, 
 8       "coordinates":[[0,0],[0,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
 9       "temperature models":
10         [
11           {"model":"half space model", "max depth":100e3, "spreading velocity":0.04, 
12            "ridge coordinates":[[[400e3,-1],[-100e3,2000e3]]]}
13         ],
14       "composition models":
BST_06_oceanic_plate_temperature.wb
BST_06_oceanic_plate_temperature.grid
Paraview v2 state file
 1{
 2  "version": "1.0",
 3  "coordinate system":{"model":"cartesian"},
 4  "features":
 5  [
 6    {
 7       "model":"oceanic plate", "name":"Overriding Plate", "max depth":100e3, 
 8       "coordinates":[[0,0],[0,1000e3],[1500e3,1000e3],[1600e3,350e3],[1500e3,0]],
 9       "temperature models":
10         [
11           {"model":"half space model", "max depth":100e3, "spreading velocity":0.04, 
12            "ridge coordinates":[[[400e3,-1],[-100e3,2000e3]]]}
13         ],
14       "composition models":
15         [
16           {"model":"uniform", "compositions":[0], "max depth":50e3}
17         ]
18    }
19  ]
20}
BST_06_oceanic_plate_temperature.wb
BST_06_oceanic_plate_temperature.grid
Paraview v2 state file
Basic Starter Tutorial section 6.

Basic Starter Tutorial section 6. The top part of the figure shows where the composition as been assigned as an object. Currently it only shows composition 0 as green and is now limited to 100km depth. The bottom part shows the temperature as seen slightly from below where only temperatures between 300K and 1600K are shown. This allows to better view the ridge.#