3D velocity models from high-resolution local earthquake tomography in Northern Chile

We provide four model files. The file raw_output_Vp.txt contains the uninterpolated three-dimensional tomographic results obtained using SIMUL23 (Eberhart-Phillips et al., 2024), with a strongly damped Vp/Vs parameter to effectively restrict the inversion to Vp. The file raw_output_VpoverVs.txt contains the corresponding uninterpolated results from a simultaneous inversion of  Vp and Vp/Vs using the same algorithm. The files Vp_interpolated.txt and VpoverVs_interpolated.txt contain the interpolated Vp and Vp/Vs models, respectively. All four files include descriptive headers.

The study region spans from 18.5–23.5°S in latitude and 67.0–71.0°W in longitude. Events are chosen with focal depths less than 165 km. This depth range was chosen to investigate the velocity structure of the subducting slab and overlying mantle wedge, particularly in relation to intermediate-depth seismicity. Only events with at least 20 P-wave and 10 S-wave arrival-time picks were retained for further analysis.