3-D data of thermal state, water content, and slab dehydration in the Colombia-Ecuador subduction zone
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 Nazca Plate, as well as the coastline and volcano distribution. The data coordinates are an orthogonal coordinate system, where the origin of the coordinates is (79.86W,2.61N), the x+ azimuth is 304.2 degrees clockwise from the north, the y+ azimuth is 214.2 degrees clockwise from the north, and the z+ direction is vertical downwards, in kilometers. The data source method is a three-dimensional finite difference numerical simulation. 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. paraview_eqkt2ge762_all_mag_cen.vtk (Earthquake Catalog by Centenial, 1900-2000, Engdahl and Villasenor，2002) 2. paraview_eqkt2ge762_all_mag_USGS.vtk (Earthquake Catalog by IRIS, 2000-2010, Trabant et al., 2012) 3. paraview_kz762.vtk (global volcanoes at NCEI, Siebert et al., 2010) 4. paraview_state762.vtk (state boundary, GMT, Wessel and Smith, 1998) 5. paraview_map762.vtk (global coast, GMT, Wessel and Smith, 1998) 6. paraview_plb762.vtk (global plate boundaries, Bird, 2003) 7. paraview_Kt_2box762.vtk (model boundary) 8. paraview_Kt_2sad762 (3-D thermal state, slab water content and slab dehydration, this study, 0km-24km from plate interface).
Files
Steps to reproduce
The data include the 3-D temperature field (°C), water content (wt%), dehydration rate (wt%/km), and subduction velocity field (cm/yr) of the subducting Nazca Plate, as well as the coastline and volcano distribution in the Colombia-Ecuador region. The data of the model region have dimensions of 800×900×400 km (length×width×depth). The geometry of the subducted plate is well constrained by slab 2.0, and the plate ages are provided by EarthByte. The subduction velocities inside a prescribed 3D constrained volume of oceanic lithosphere are determined based on kinematic plate subduction modeling and MORVEL plate motion data. The observations of surface heat flow and Curie point depths are used to constrain the model thermal regime. Geophysical calculation is performed after the subduction thermal regime reaches a steady state. The data were acquired through 3-D numerical thermal modeling developed from the original stag3d code. This model employs an anelastic liquid approximation and the equations of conservation of mass, momentum, and energy. Using a three-dimensional thermomechanical model and earthquake catalogs from IRIS and Centennial, we calculated the thermal state, water content, and slab dehydration distribution in the Colombia-Ecuador subduction zone.