Dominant role of slab geometry in mediating carbon and water fluxes in subduction systems: Implications from 3-D thermal modeling of the Ryukyu trench

Published: 7 March 2024| Version 1 | DOI: 10.17632/src9zx425d.1
Contributor:
Yingfeng Ji

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 of the model region have dimensions of 950×475×380 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 of the model region have dimensions of 950×475×380 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 including: 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. pvsm files (state files of paraview, corresponding to the figures in article) 8. vtk data nanlysis (excel)

Institutions

Institute of Tibetan Plateau Research Chinese Academy of Sciences

Categories

Geophysics, Geodynamics

Licence