Transient Acceleration and Termination of Japan Sea Opening Controlled by Oblique Subduction Induced Slab Window and Mantle Flow

This repository contains the numerical modelling code used in: "Transient Acceleration and Termination of Japan Sea Opening Controlled by Oblique Subduction Induced Slab Window and Mantle Flow" Model description:- A 3D visco-plastic mechanical model of oblique subduction implemented using Underworld2- The model includes subduction slab, overriding plate, and weakzone structures Additional notes on model setup:- All model geometry, physical parameters, and Underworld2 runtime configurations are defined within a single script.- Different scripts correspond to models with different obliquity angles (45°, 30°, 15°, and 0°), while all other parameters are kept identical.- This suite of experiments is designed as a single-variable test, where only the obliquity angle varies between models.- To test additional angles, users can modify the parameter “ang” in the script. How to run:- Requires Underworld2 (version 2.10 or later) with MPI support- Run the model using: sbatch Run_Cluster.sh- Note: Minor modifications to the submission script may be required depending on the HPC system Outputs:- Outputs are automatically generated to the output directory when reproducing the model- Time evolution of Vrms is recorded in FrequentOutput.dat Visualisation-paraview:- Output files are organised in XMF format, where all timesteps are referenced in a single .xmf file (e.g., XDMF.field.00XX.xmf)- Users can directly open the .xmf file in ParaView to visualise the full set of model fields (e.g., velocity, strain rate, viscosity, density)- To visualise the slab geometry and slab window as presented in the manuscript, apply the Iso Volume filter in ParaView to the viscosity field and retain values between 5-7 × 10²¹ and 1 × 10²⁵ Pa s. This isolates the high-viscosity lithospheric material and makes the slab window clearly visible. Alternative visualisation (Underworld built-in):- For users who are not familiar with ParaView, a simple visualisation workflow is provided using Underworld's built-in visualisation module (Visual_UW.py).- The example script demonstrates how to render material distributions directly from the model swarm. Additional data for figure reproduction:- Supplementary data files (e.g., data_age.txt, data_rate.txt, Fast_Line.csv, Middle_Line.csv, Slow_Line.csv) are provided to reproduce the figures presented in the manuscript. Note:Due to the large size of the full simulation outputs, they are not archived in this repository. All results presented in the manuscript can be reproduced using the provided code and model setup.