Subbasalt imaging has gained significant interest in the last two decades, driven by the urge to better understand the geologic structures beneath volcanic layers, which can be up to several kilometers thick. This understanding is crucial for the development and risking of hydrocarbon play models in these areas. However, imaging based on the reflection seismic data alone suffers from severe amplitude transmission losses and interbed multiples in the volcanic sequence, as well as from poor definition of the subbasalt velocity structure. We have considered a sequential imaging workflow, in which the resistivity model from joint controlled-source electromagnetic and magnetotelluric data inversion was used to update the velocity model and to improve the structural definition in the migrated seismic image. The quantitative link between resistivity and velocity was derived from well data. The workflow used standard procedures for seismic velocity analysis, electromagnetic data inversion, and well analysis, and thereby allowed detail control and input based on additional geophysical knowledge and experience in each domain. Using real data sets from the Faroe-Shetland Basin, we can demonstrate that the integration of seismic and electromagnetic data significantly improved the imaging of geologic structures covered by up to several-kilometer-thick extended volcanic sequences. The improved results might alter the interpretation compared with the imaging results from seismic data alone.