Models | M@TE

Available Models

Fast modeling results (1.8 cm.y-1 extension speed) at an average strain rate of 3e-15 s-1 and 25 % extension. Colors are the same as in Fig 2. Model a. Only partial melting is allowed (garnet-pyroxene isograde and retrogression into amphibolite are removed). Model b. Partial melting, crystallization of garnet-pyroxene assemblages and retrogression are allowed. The temperature for the transformation of the continental crust into garnet-pyroxene rich rocks is 777 °C (see text for explanation). The reference temperatures for the solidus of the continental crust and the garnet-pyroxene rich crust are 650 °C and 790 °C, respectively. Model c. Same as model b but the temperature for the continental crust solidus is increased to 720 °C. Model d. Zoom on model b illustrating the velocity field (black arrows) when the boundary condition mimicking extension is removed (after 2 Myr of gravity forces operating) showing that the partially molten crust flows while garnet-pyroxene rich rocks do not.

Strain and retrogression partitioning explain long-term stability of crustal roots in stable continents

Uploaded by: Bénédicte Cenki-Tok or Cenki

Uploaded: October 28, 2024

Software: Underworld 2

Tags:

crustal roots retrogression stability

Finite Element Python

Mantle flow and upper plate deformation due to slab step

Subduction dynamics and plate coupling along seismically active margin: the role of slab steps

Uploaded by: Thyagarajulu Gollapalli

Uploaded: October 24, 2024

Software: Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Tags:

Plate boundary Subduction

Finite element Particle-in-cell Python

Marker lines prior to (a) and after 750 years of flow of (b) isotropic and (c) anisotropic ice. The axial plane of the resulting shear fold in isotropic ice mimics the bedrock topography, while it is controlled by shearing along a horizontal shear zone in the case of anisotropic ice. Green: bedrock, flow to the right.

ISMIP-HOM benchmark experiments using Underworld

Uploaded by: Dan Sandiford

Uploaded: August 19, 2024

Software: Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Tags:

Stokes anisotropy benchmark ice-sheet mechanical

Finite-Element Particle-in-cell Python

Figure 5. Evolution of the stress at the reference point (0, 0, −10 km) (a), maximum slip rate along the entire fault zone (b)
and the adaptive time step used in simulation (c) for the reference model.

Numerical Modeling of Earthquake Cycles Based On Navier-Stokes Equations With Viscoelastic-Plasticity Rheology

Uploaded by: Haibin Yang

Uploaded: August 16, 2024

Software: Underworld

Tags:

Cratonic Earthquakes Earthquake Cycles Navier-Stokes Viscoelastic-Plasticity

Finite element Particle-in-cell Python

$Downdip component of strain rate tensor and resolved stress difference from the numerical model, focusing on features within the plate/slab. The resolved stress difference is defined as ($\sigma_{s} - \sigma_{z}$), where $\hat{s}$, and $\hat{z}$ are unit vectors in the downdip and slab normal directions. The fields show, for example, shortening/extension in the downdip direction. Stress profiles at four locations are shown. The blue line ($x_0$) is the first zero crossing based on analysis of the flexural component of the topography. The black line is the location of maximum bending moment.$