Skip to content

Commit

Permalink
Merge pull request #727 from ljhwang/patch-44
Browse files Browse the repository at this point in the history
Update 16_mass_conserving.md
  • Loading branch information
MFraters authored Jun 13, 2024
2 parents 0916344 + 3f12c28 commit 910de0b
Showing 1 changed file with 6 additions and 6 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -2,10 +2,10 @@
Changing to a mass conserving slab temperature
===============================

In [the section on slab temperature](part:user_manual:chap:basic_starter_tutorial:sec:12_subducting_plate_temp), we added the {cite:t}`McKenzie_1970` (that is the plate model) slab temperature structure. Although this is a good first order approximation of a slab temperature, the recently develop {cite:t}`Billen_Fraters_AGU_2023` temperature model has many advantages over it, with the downside that is it a bit more involved to use. In this chapter we will just switch out the plate model with mass conserving, without going too much into the detail of how to actually use it in practice. If you plan to use the mass conserving temperature model, please first read through both the [Simple Subduction Model: 2D Cartesian](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_cartesian)and [Simple Subduction Model: 2D Chunk](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_chunk) cookbooks.
In the section on [slab temperature](part:user_manual:chap:basic_starter_tutorial:sec:12_subducting_plate_temp), we added the {cite:t}`McKenzie_1970` (that is the plate model) slab temperature structure. Although this is a good first order approximation of a slab temperature, the recently develop {cite:t}`Billen_Fraters_AGU_2023` temperature model has many advantages over it, with the downside that is it a bit more involved to use. In this chapter we will just switch out the plate model with mass conserving without going too much into the detail of how to actually use it in practice. If you plan to use the mass conserving temperature model, please first read through both the [Simple Subduction Model: 2D Cartesian](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_cartesian)and [Simple Subduction Model: 2D Chunk](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_chunk) cookbooks.

## Changing the Subducting oceanic plate temperature
One of the advantages of using the mass conserving slab temperature model is that it doesn't assume a linear temperature structure at the trench. This means we can seamlessly connect a half space cooling model or a plate model to the mass conseving slab. But for that to make sense, we need to change the subucting oceanic plate to a half space cooling model (the default for the mass conserving temperature model). In this case we will put the ridge far away. Because the half space model can affect the temperature at much deeper depths, we also need to change the max depth for the feature and models. Note that we will want to keep the max depth of the composition 100km, so we now need to set a max depth in the lowest layer of the composition model.
One of the advantages of using the mass conserving slab temperature model is that it doesn't assume a linear temperature structure at the trench. This means we can seamlessly connect a half space cooling model or a plate model to the mass conserving slab. But for that to make sense, we need to change the subucting oceanic plate to a half space cooling model (the default for the mass conserving temperature model). In this case we will put the ridge far away. Because the half space model can affect the temperature at much deeper depths, we also need to change the max depth for the feature and models. Note that we will want to keep the max depth of the composition at 100km, so we now need to set a max depth in the lowest layer of the composition model.


::::::{tab-set}
Expand Down Expand Up @@ -59,9 +59,9 @@ One of the advantages of using the mass conserving slab temperature model is tha

## Changing to a mass conserving temperature model

One of the main idea's behind the mass consering model is that the the slab not just warms up, but that the surrounding material cools down as well, thus conserving energy/mass. This means that the subducting plate feature needs to be able to change the temperature outside the slab. Given that a feature can only change the temperature inside it's defined area we need to expand that area. This means we need to increase the thickness, which we set to 300 km. But this only extends the slab feature downwards, while we also want to cool down the material above the slab. For this we can use the `top truncation` parameter. The `top truncation` parameter is designed to make the slab thinner from above, but when set to a negative value we can actually make it thicker! In this case we will set it to -100 km.
One of the main ideas behind the mass conserving model is that the the slab not just warms up, but that the surrounding material cools down as well, thus conserving energy/mass. This means that the subducting plate feature needs to be able to change the temperature outside the slab. Given that a feature can only change the temperature inside it's defined area we need to expand that area. Hence, we need to increase the thickness which we originally set to 300 km. But this only extends the slab feature downwards, while we also want to cool down the material above the slab. For this we can use the `top truncation` parameter. The `top truncation` parameter is designed to make the slab thinner from above but when set to a negative value we can actually make it thicker! In this case we will set it to -100 km.

Now that all the preparatory work is done, we can finally add the mass conserving model itself. The two most important new parameter are the `ridge coordinates`, which should be the same as the subducting oceanic plate ridge coordinates in this case and a parameter called `coupling depth`. The `coupling depth` defines the depth at which the slab surface first comes in contact with the hot mantle wedge.
Now that all the preparatory work is done, we can finally add the mass conserving model itself. The two most important new parameter are the `ridge coordinates` which should be the same as the subducting oceanic plate ridge coordinates in this case, and a parameter called `coupling depth`. The `coupling depth` defines the depth at which the slab surface first comes in contact with the hot mantle wedge.

::::::{tab-set}

Expand Down Expand Up @@ -119,7 +119,7 @@ Now that all the preparatory work is done, we can finally add the mass conservin
:alt: Basic Starter Tutorial section 18 highres result.
:align: center
Basic Starter Tutorial section 16 high resolution result, were a mass conserving slab temperature is used. This has 8 times the resolution then the grid file above.
Basic Starter Tutorial section 16 high resolution result, where a mass conserving slab temperature is used. This has 8 times the resolution than the grid file above.
```

These are only the very basics of what the mass conserving temperature model can do. If you are interested in using this temperature model, please read the [Simple Subduction Model: 2D Cartesian](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_cartesian)and [Simple Subduction Model: 2D Chunk](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_chunk) cookbooks.
These are only the very basics of what the mass conserving temperature model can do. If you are interested in using this temperature model, please read the [Simple Subduction Model: 2D Cartesian](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_cartesian)and [Simple Subduction Model: 2D Chunk](part:user_manual:chap:cookbooks:sec:simple_subduction_2d_chunk) cookbooks.

0 comments on commit 910de0b

Please sign in to comment.