LaMEM tutorial as part of the 2024 ALW

This repository contains material and scripts for the mini tutorial on LaMEM as part of the 2024 Ada Lovelace workshop on Modelling Mantle and Lithosphere Dynamics.

a. Required software

Participants will have to install the following software:

• VSCode This is the most widely used coding environment worldwide and provides support for julia. Also add the julia extension (click on the left side to add extensions).
• Julia We will use julia to create model setups and run LaMEM. The easiest way to install julia these days is through juliaup. Once you did that, you can install the LaMEM, GeophysicalModelGenerator and Plots packages. We currently use julia 1.10
• Paraview - This is a very powerful visualisation tool. Any version will do.

b. Topics

We will use the julia interface to LaMEM - called LaMEM.jl in this workshop as it is by far the easiest way to setup LaMEM models.

1. Julia intro

Installing software and intro to julia - In case you need some more help on julia you can work through these exercises. Note that we will not discuss this in detail during the tutorial (for time reasons).

2. Falling sphere example

Falling Sphere - our first example is a 3D model of a dense sphere sinking through a less dense matrix. This involves a linear viscous material and multigrid solvers and gives you a first idea on how to setup and run simulations.

3. Plume lithosphere interaction

Plume Lithosphere example - in the next 2D example we show how you can setup a model with temperature-dependent nonlinear viscosities, a free surface, viscoelastic rheologies, shear and adiabatic heating along with passive tracers.

4. 3D subduction

3D Subduction - in the next example we reproduce results of a paper by Schellart and coworkers that looked at the effect of a slab width on the trench curvature.

5. Fold and thrust belt

Fold and Thrust - fold and thrust belts are a nice example of how brittle layers may slide overa detachment horizon and cause both folding and faulting. In this example, we reproduce this using visco-elasto-viscoplastic rheologies and an (internal) free surface

c. Additional help

You can find additional material by looking at the

• documentation of LaMEM.jl
• tests of LaMEM.jl - these tests are run every time a new commit is made to the LaMEM repository so are up to date.
• main input file of LaMEM which lists all available options. We try to keep LaMEM.jl up to date with this.
• CHEES2 workshop material this wqs a workshop held in spring 2024 and contains some usful material
• Heidelberg LaMEM course - a shortcourse at the university of Heidelberg held by Nicolas Riel (for BSc students). Note that we have sligthly updated the names of some functions in GeophysicalModelGenerator since that time so you may have to make a few modifications to the scripts

Finally, if you find mistakes or errors, please don't be shy to leave in issue, or perhaps even fix it yourself and make a pull request. The code will only get better with everyone's help!