A work in progress movie of the coupling between ASPECT and openLEM. This has a rotating continent and a plume rising below it in ASPECT. the velocity field is than being pass to openLEM, which extracts the horizontal rotation and translation of the continental area and applies that to the gird. The residual horizontal velocity and the vertical velocity is then used to compute the uplift, which is passed back to ASPECT. The sea level is determined by a sinus wave with a amplitude of 50m and a period of 20kyr. It is meant as a benchmark and the parameters used, nor the results, are note necessarily physical.

This movie shows a continental extension model where intrusion locations (white lines) are computed based on where and how much melt (red line) is generated. Each timestep is 20.000 years. The intrusions last one timestep and start at around 15Myr. Note the this is only part of the model run. For a more complete model run (without dilation) see https://tube.mfraters.net/w/4EgbqeurD17XPTcCDUq3WW.

This movie shows a continental extension model where intrusion locations (white lines) are computed based on where and how much melt (red line) is generated. Each timestep is 20.000 years. The intrusions last one timestep and start at around 15Myr.

This is a movie of a Work in Progress model, which will be shown at my AGU presentation. It shows about 1.5Ma evolution of the Cascadia subduction zone and the development of Crystal Preferred Orientation (CPO). The model is not ready for actually usage, but is a good showcase of the progress which has been made.

On the left and bottom right there is a top view of the region. The bottom right panel show the viscosity with colors at the top and velocity fields with the colored arrows. The ridges are clearly viable in the light colors. The top right panel show a cross-section through the model, where the cylinders show the elastic fast direction. The two white lines outline 100km and 200km depth. The left figure show the the development of the CPO and the velocity field bewteen 100km and 200km depth, viewed from above.

This movie shows the development of CPO in a simple 3D subduction zone. The numbers on the figures are in meters. The figures on the left are the top view, while the right figures are a cross section. In each sub figure, the top image shows the velocity field with the arrows, where the size of the of the arrows are determined by the velocity magnitude and the color the depth, and on the right side the cross section colors show the viscosity field. In the top left figure, the top 150km has been removed to better show the velocity field below. The bottom figure of each sub figure shows the direction of the largest elastic eigenvector, colored by the active olivine fabric type and scaled by the integrated strain invariant such that the size of the cylinder does not grow beyond a integrated strain invariant of 0.5, and on the right side the cross section colors show the background shows the water content. The background on the right pictures is located at 550km, and the arrows and cylinders are cut off up to 450km.

This movie is part of Fraters and Billen (submitted). Preprint available at doi.org/10.1002/essoar.10506908.1.

CPO development in the mid ocean ridge/corner flow case at different model times. Each particle contains 1000 grains and has a mixture of 70% olivine and 30% enstatite. The bar orientation is the direction associated with the largest eigenvector of the elastic tensor. The background color shading indicates the water content, while the color of the bars indicates the deformation type for this particle at this time in the simulation. Only the top 150~km of the model is shown. Only CPO directions with strain higher than 0.25 are shown. The numbers on the axes are meters.

This movie is part of Fraters and Billen (submitted). Preprint available at doi.org/10.1002/essoar.10506908.1.