# 3-D data of thermal regime, water content, and slab dehydration in Central America

## 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 (88.10W,10.70N), the x+ azimuth is 212 degrees clockwise from the north , the y+ azimuth is 122 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_coast.vtk (global coastlines, GMT, Wessel and Smith, 1998) 2. paraview_eq_cen.vtk (Earthquake Catalog by Centennial, 1900–2000, Engdahl and Villasenor, 2002) 3. paraview_eq_USGS.vtk (Earthquake Catalog by IRIS, 2000-2010, Trabant et al., 2012, http://ds.iris.edu/ieb/) 4. paraview_plb.vtk (global plate boundaries, Bird, 2003) 5. paraview_slab.vtk ((3-D thermal regime, slab water content and slab dehydration, this study, 0km-24km from plate interface) 6. paraview_state.vtk (global coastlines, GMT, Wessel and Smith, 1998) 7. paraview_volcano.vtk (global volcanoes at NCEI, Siebert et al., 2010)

## Files

## Steps to reproduce

The data include the 3-D temperature field (degrees Celsius), water content (wt%), dehydration rate (wt%/km), and subduction velocity field (cm/yr) of the subducting plate, as well as the coastline and volcano distribution in Central America. The data of the model region have dimensions of 500×1000×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 subduction velocities inside a prescribed 3-D constrained volume of the oceanic lithosphere are given based on the kinematic plate subduction modeling method and the MORVEL plate motion data. The observation of surface heat flow and Curie point depths are used to constrain the model thermal regime. The geophysical calculation is ensured after the subduction thermal regime reaches a steady state. 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 and Centennial, the data of the slab thermal state, water content, and slab dehydration distribution in Central America are calculated.