P. Arroucau, Graça Silveira, Nuno Dias, Susana Custódio, Antonio Villaseñor, Chiara Civiero, Thomas Bodin, Jordi Diaz, EDF (EDF), Dublin Institute for Advanced Studies (DIAS), Instituto Dom Luiz, Universidade de Lisboa (ULISBOA), Institute of Earth Sciences Jaume Almera, Consejo Superior de Investigaciones Científicas [Madrid] (CSIC), Institute of Marine Sciences / Institut de Ciències del Mar [Barcelona] (ICM), Laboratoire de Géologie de Lyon - Terre, Planètes, Environnement [Lyon] (LGL-TPE), Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Université Claude Bernard Lyon 1 (UCBL), Université de Lyon-Université de Lyon-École normale supérieure - Lyon (ENS Lyon), Science Foundation Ireland, Geological Survey of Ireland, Universidade da Beira Interior, Agencia Estatal de Investigación (España), Diaz, J. [0000-0003-1801-0541], Universidade de Lisboa = University of Lisbon (ULISBOA), Laboratoire de Géologie de Lyon - Terre, Planètes, Environnement (LGL-TPE), École normale supérieure de Lyon (ENS de Lyon)-Université Claude Bernard Lyon 1 (UCBL), Université de Lyon-Université de Lyon-Institut national des sciences de l'Univers (INSU - CNRS)-Université Jean Monnet - Saint-Étienne (UJM)-Centre National de la Recherche Scientifique (CNRS), and Diaz, J.
22 pages, 17 figures, 2 tables, supporting information https://doi.org/10.1093/gji/ggab005, We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (, This publication is supported by the FCT project SPIDER - Seismogenic processes in slowly deforming regions (PTDC/GEO-FIQ/2590/2014). The work presented is a contribution to FCT UIDB/50019/2020 IDL. P. Arroucau acknowledges support from Science Foundation Ireland (grant 13/CDA/2192) and from the Geological Survey of Ireland (grant 2016-PD-06). C. Civiero was supported by the Science Foundation Ireland, the Geological Survey of Ireland, and the Marine Institute (grants 13/CDA/2192 and 16/IA/4598). Figures were plotted using the GMT - Generic Mapping Tools software (Wessel & Smith 1998) and Python Matplotlib and Basemap packages. NonLinLoc was used for the forward computation of arrival times (Lomax et al. 2000). We also wish to thank Catarina Matos for the tests she performed on an early version of the model., With the funding support of the ‘Severo Ochoa Centre of Excellence’ accreditation (CEX2019-000928-S), of the Spanish Research Agency (AEI).