A crustal tomography of seismic wave velocity was performed in the contact zone between the southern Tyrrhenian, Sicilian and Ionian tectonic units, a zone where the lithospheric structure can be expected to furnish evident signatures of dynamics related to the Tyrrhenian subduction process. A dataset of 10241 P and 5597 S readings from 932 local earthquakes recorded between 1978 and 2001 by stations operating in Sicily and Calabria was inverted by the SIMULPS12 algorithm for simultaneous computation of hypocenter parameters and Vp and Vp/Vs three dimensional distributions. The study brought significant improvement in the knowledge of the local velocity structure, furnishing new information useful to better identify the local tectonic units. The results appear to be compatible with the most recent hypotheses regarding the geodynamics of the study region.