Heterogeneous slab thermal dehydration driving warm subduction zone earthquakes
Description
The data includes the three-dimensional temperature field (degrees Celsius), water content (wt%), dehydration rate (wt%/km), rock phase change distribution, and subduction velocity field (cm/yr) of the subducting plate, as well as the coastline and volcano distribution. The data coordinates are an orthogonal coordinate system, where the origin of the coordinates is (166.3E, 16.4S).The data of the model region have dimensions of 1200×500×400 km (length×width×depth). The geometry of the subducted plate is well constrained by Slab2.0, and the plate ages are provided by EarthByte. The temperature boundary condition agrees with the plate cooling model. The bottom of the slab and the perpendicular plane are prescribed as adiabatic and permeable, and the top surface is set to a fixed temperature (0°C) and is rigid. The data was acquired through 3-D numerical thermal modeling which was developed from originally code stag3d. An anelastic liquid approximation and the equations of conservation of mass, momentum, and energy are used in this model. Based on a three-dimensional thermomechanical model and the collect earthquake catalog from IRIS, the data of the slab thermal state, water content, and slab dehydration distribution in Vanuatu are calculated. The geophysical calculation is ensured after the subduction thermal regime reaches a steady state. The temperature field error range is ±50 degrees Celsius, and the velocity field error range is ±0.1cm/ yr. This data can be further used analyzing the geophysical field of subduction zone.
Files
Steps to reproduce
The data includes the three-dimensional temperature field (degrees Celsius), water content (wt%), dehydration rate (wt%/km), rock phase change distribution, and subduction velocity field (cm/yr) of the subducting plate, as well as the coastline and volcano distribution. The data coordinates are an orthogonal coordinate system, where the origin of the coordinates is (166.3E, 16.4S).The data of the model region have dimensions of 1200×500×400 km (length×width×depth). The geometry of the subducted plate is well constrained by Slab2.0, and the plate ages are provided by EarthByte. The temperature boundary condition agrees with the plate cooling model. The bottom of the slab and the perpendicular plane are prescribed as adiabatic and permeable, and the top surface is set to a fixed temperature (0°C) and is rigid. The data was acquired through 3-D numerical thermal modeling which was developed from originally code stag3d. An anelastic liquid approximation and the equations of conservation of mass, momentum, and energy are used in this model. Based on a three-dimensional thermomechanical model and the collect earthquake catalog from IRIS, the data of the slab thermal state, water content, and slab dehydration distribution in Vanuatu are calculated. The geophysical calculation is ensured after the subduction thermal regime reaches a steady state. The temperature field error range is ±50 degrees Celsius, and the velocity field error range is ±0.1cm/ yr. This data can be further used analyzing the geophysical field of subduction zone. File naming: the data is stored in text paraview format. Data reading method: All ASCII code files in the data set can be directly opened with paraview software. Data list: 1. eqs_USGS.vtk (Earthquake Catalog by IRIS, 2000-2010, Trabant et al., 2012) 2. volcano.vtk (global volcanoes at NCEI, Siebert et al., 2010) 3. map.vtk (map and state boundary, GMT, Wessel and Smith, 1998) 4. plate_boundary.vtk (global plate boundaries, Bird, 2003) 5. mantle.vtk (3-D velocity field in this study). 6. slab.vtk (3-D thermal regime, slab water content and slab dehydration in this study). 7. five pvsm files (state files of paraview, corresponding to the figures in article) 8. readme.txt (parameters and explanation)