This data repository for the Southern Caribbean and NW South America contains a 3D thermal model computed down to 75 km depth, the modelled hypocentral temperatures and geothermal gradients at the locations of crustal earthquakes, and the crustal seismogenic depths calculated from earthquake statistics, as well as the associated modelled temperatures.

We used the uppermost 75 km of the gravity-constrained structural and density model of Gómez-García et al. (2020, 2021) to derive the 3D thermal configuration of the study area (5°-15° N, 63°-82° W). A steady-state approach was followed, in which upper and lower boundary conditions were set to run the thermal calculations using the software GOLEM (Cacace & Jacquey, 2017; Jacquey & Cacace, 2017).

A catalogue of earthquakes occurred within the study area and surroundings was compiled from public sources. In the database archived here, we provide data of the best located crustal earthquakes within the boundaries of this area, from January 1980 to June 2021. Earthquakes below the magnitude of completeness, or with poorly determined depths, were disregarded. Earthquakes were deemed crustal if their hypocentres were located between the topo-bathymetry from the GEBCO relief (Weatherall et al., 2015) and the Moho depth from the GEMMA model (Reguzzoni & Sampietro, 2015).

We computed the crustal seismogenic depth as the 90th and 95th percentiles (D90 and D95), respectively, of the crustal hypocentral depths. These percentiles were mapped on a latitude-longitude grid, using for each grid node at least the 20 closest earthquakes as sample.

The hypocentral temperatures, the geothermal gradient at the earthquake locations, and the temperatures at the D90 and D95 surfaces were calculated from the lithospheric-scale thermal model. For more details about the modelling approach and interpretation of the results, we kindly ask the reader to refer to the main publication: Gomez-Garcia et al. (2023).

Version history:

XXXX release of version 2023. The main changes are:

• New title
• New earthquake dataset compiled from different sources, with more accurate hypocentral determinations.
• Error propagation in the D90 and D95 depths considering the reported errors of each event.
• D10 values are not computed, as they are prone to large uncertainties.
• The geothermal gradient at each hypocentral location is computed and reported.