2D visco-elasto-plastic subduction models in the presence of a high density and viscosity continental block

<p>Three main models with different rheological conditions and density are incorporated for an imposed intrusion between the upper and lower continental crust. Output files for a model with a crust without intrusions are also incorporated.</p> <p>The file contains data on the conditions of plate motion velocity, plate age, plate thickness in Excel (Input_Data_Age_Vel.xlsx), a table with the rheology used to reproduce different numerical models incorporating the intrusion (Table_S3.pdf), and a file High_density_materials.m showing the construction of the intrusion within the continental crust.</p> <p>High_density_materials.m m should be imposed on the original I2ELVIS code provided by Taras Gerya - ETH Z&uuml;rich, Institut f&uuml;r Geophysik, Switzerland - email: taras.gerya@erdw.ethz.ch</p> <p>The file has the main time steps as a function of temperature, viscosity, density and rock composition&nbsp;into the folders. If you want to visualize them please run the following Script PLOT_OUTPUT_DATA.m where,</p> <p>str=string(2040), represent the time_step</p> <p>Files are organized in a matrix form, each folder has the</p> <p>grid_x :coordinate matrix,&nbsp;<br>gridy_: coordinate matrix,&nbsp;</p> <p>matrix_temperature_, density and viscosity.</p> <p>To plot the Rock composition.</p> <p>mx and my, tracer coordinates in x-direction and y-direction<br>matrix_markcom is&nbsp; matrix tracers.</p> <p>&nbsp;</p> <p>Model 1</p> <p>(Time step= 10= 14.99 Ma, 120= 14 Ma, 310= 12 Ma, 670 = 9 Ma, 1120 = 6 Ma, 1620 = 3 Ma, 2040 = 0 Ma)</p> <p>&nbsp;</p> <p>The temperature, viscosity, density, and rock composition for the PHS slab are shown in this animation for Model 1.&nbsp;<br>This model is reproduced under the conditions of plate motion velocity and age from Figure 3 with an initial dip angle for the weak zone of 20&deg;. &nbsp;<br>The high density material between the upper and lower crust in the vicinity of the slab is not included in Model 1.&nbsp;<br>The model evolves over a period of the last 15 Myr and shows a shallow subduction.&nbsp;</p> <p><br>Model 2</p> <p>(Time step= 10= 14.99 Ma, 120= 14 Ma, 380= 12 Ma, 850 = 9 Ma, 1490 = 6 Ma, 2090 = 3 Ma, 2540 = 0 Ma)</p> <p>With high density block<br>The temperature, &nbsp;viscosity, density, and rock composition for the PHS slab are shown in this animation for Model 2. This model is reproduced under the conditions of plate motion velocity and age from Figure 3 with an initial dip angle for the weak zone of 20&deg;. &nbsp;The high density material between the upper and lower crust in the vicinity of the slab is included in Model 2. The model evolves over a period of the last 15 Myr and shows a steep subduction.</p> <p><br>Model 3</p> <p>(Time step= 10= 14.99 Ma, 120= 14 Ma, 310= 12 Ma, 670 = 9 Ma, 1120 = 6 Ma, 1690 = 3 Ma, 2040 = 0 Ma)</p> <p>Without high density block and high viscosity</p> <p>The temperature, viscosity, density, and rock composition for the PHS slab are shown in this animation for Model 3. This model is reproduced under the conditions of plate motion velocity and age from Figure 3 with an initial dip angle for the weak zone of 20&deg;. &nbsp;The initial high viscosity material and without high density between the upper and lower crust in the vicinity of the slab is included in Model 3. The model evolves over a period of the last 15 Myr and shows subduction with a high dip angle.&nbsp;</p> <p><br>Model4</p> <p>(Time step= 10= 14.99 Ma, 120= 14 Ma, 330= 12 Ma, 770 = 9 Ma, 1190 = 6 Ma, 1640 = 3 Ma, 2040 = 0 Ma)<br>&nbsp;&nbsp;<br>With high density block and high viscosity<br>The temperature, viscosity, density, and rock composition for the PHS slab are shown in this animation for Model 4. This model is reproduced under the conditions of plate motion velocity and age from Figure 3 with an initial dip angle for the weak zone of 20&deg;. &nbsp;The high density and high viscosity material between the upper and lower crust in the vicinity of the slab is included in Model 4. The model evolves over a period of the last 15 Myr and shows a steep subduction</p>