Modelling the response of seismic wavefields to sharp lateral variations in crustal discontinuities is essential for the application of seismic tomography and the correction of path effects in earthquake source characterization. Researchers generally focus on travel times and amplitudes of a limited number of phases, simplifying the underlying physics and reducing computational times. However, this approach offers insufficient resolution across oceanic basins while leaving a large part of the information contained in seismic records unexploited. Here, we model the sensitivity of wave propagation to crustal structural variations in the Southern Tyrrhenian basin by fitting the data produced by the Accumoli earthquake (Central Italy, 2016) and recorded at the national INGV (Istituto Nazionale di Geofisica e Vulcanologia) network. As a result, we reconstruct the velocity, Moho, and sediment structures in the basin and improve estimates of earthquake-source characteristics. P-, S- and coda-wave are first modelled and fitted to data acquired by a single station in Sicily; then, the procedure is extended to seismic receivers of the INGV network that contour the basin. The optimal fits are obtained using a deep Moho (~ 18 km) in the middle of the basin and a crustal pinch between the basin and the continental crust in Sicily. The deep Moho corresponds to the Issel Bridge, a portion of continental crust trapped between the Vavilov and Marsili volcanic centres. The Accumoli earthquake is optimally described using a boxcar source time function with a rise time of 6 s, a depth of 7.3 km, and the radiation pattern previously identified by INGV. Correlation analyses performed on data recorded across the network demonstrate the improvement in P-, S- and coda waveform fitting relative to previous Moho models obtained inverting single-phase observations. The early S-wave coda comprises trapped and reverberating phases that are particularly sensitive to crustal interfaces. With our approach, full-waveform simulations become a viable strategy for the inversion of Moho depths across the Southern Tyrrhenian using regional earthquakes.