Comparing slab viscous dissipation under slab distributed weakening and slab segmentation scenarios

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==================== ## ========================