DOI: 10.5281/zenodo.18428463
title: Comparing slab viscous dissipation under slab distributed weakening and slab segmentation scenarios

description: This is the README file for the dataset of 'Comparing slab viscous dissipation under slab distributed weakening and slab segmentation scenarios'.

The data file 'Segmentation_40 (60, 80).xlsx' contains the data from the segmentation model where the subducting slab age is 40 Myr (60, 80 Myr).

The data file 'Sdw_40 (60, 80).xlsx' contains the data from the slab distributed weakening model where the subducting slab age is 40 Myr (60, 80 Myr).

The file 'natural_subductionzone_para.xlsx' contains the dataset used to calculate natural subduction-slab dissipation.

The script 'multifactor-regression.py' performs the multi-factor linear regression. The files 'Segmentation_40 (60, 80).xlsx' and 'Sdw_40 (60, 80).xlsx' are used as inputs for this regression.

The script 'set_model.py' is the Underworld input script that we used for our simulations.

The archive 'segmentation(Sdw)-40(60,80).zip' contains post-processing data used to describe the evolution of several slab parameters during subduction. These files share the same structure:

(1) timestep.log — Model evolution time (simulation timestep log).==================== ## ========================

(2) EFP.npy — This dataset includes parameters for subducting-slab energy dissipation and release, subduction velocity, and effective viscosity, enabling analysis of the temporal evolution of the dissipation fraction.

# 1--time                               2--energy_gravity_potential                3-energy_visco_dissipation_slab;# 4-energy_visco_dissipation_mantle;    5-energy_visco_dissipation_crust;          6-energy_visco_dissipation_overriding;# 7-energy_visco_dissipation_litho;     8-energy_visco_dissipation;                9-vel_ad;# 10-vel_sink;                          11-energy_gravity_potential_all;           12-energy_gravity_potential_all_delta;# 13-energy_kinetic;                    14-energy_mantle_diff(temp effect);        15-energy_gravity_potential_1300_mantle;# 16-energy_pressure;                   17-energy_gravity_potential_thermal_slab;  18-rho_shallow;# 19-rho_middle;                        20-rho_deep;                               21-energy_gravity_potential_craton;# 22-energy_gravity_potential_thermal_mantle;# 23-energy_gravity_potential_thermal_overriding;# 24-yita_f(strainRate);                25-yita_f_trenchLeft;                      26-energy_visco_dissipation_slab_upper;# 27-Volume_eclogite                    28-yita_f(Volume)                          29-yita_f_nocrust(strainRate)# 30-yita_f_upper_200km(Volume)         31-yita_f_lower_200km(Volume)              32-vel_trench# 33-yita_f_upper(strainrate)           34-yita_f_lower(strainrate)                35-thickness_slab# 36-thickness_95%_dissipation          37-thckness_90%_dissipation                38-yita_f(Volume and StrainRate) # 39-temperature(slab)==================== ## ========================

(3) Bending.npy — This file contains energy dissipation in different compositional fields during slab subduction, partitioned into the crust and lithospheric mantle, as well as dissipation localized along the subduction interface.

# 1--time;                              2--energy_visco_dissipation_slab_bending;  3-energy_visco_dissipation_crust_bending;# 4-energy_visco_dissipation_litho_bending;# 5-energy_visco_dissipation_slab_iv;   6-energy_visco_dissipation_mantle_iv;      7-energy_visco_dissipation_crust_iv;# 8-energy_visco_dissipation_litho_iv;  9-decoupling_de_c;                         10-decoupling_de_l;# 11-decoupling_de_c_nv;                12-decoupling_de_l_nv;                     13-R_average;# 14-R2==================== ## ========================

(4) Integral.npy — It includes the time evolution of the plate’s curvature-gradient integral and the length of the subducted slab (sum((dk/ds)**2 ds)).

# 1--integral_0(sum);                   2--integral_1(trapz)                       3--integral_2 (dkds**2*yita)# 4--length of subducted slab           5--length of subducted slab (2nd method)   6--trench location# 7--length of unsubduct slab(3rd)      8--length of unsubduct slab (4th method)   9--dkds_upper==================== ## ========================