In order to broaden our knowledge about the deep lithosphere using seismic methods, we develop leading-edge imaging workflows integrating different types of data. Here we exploit the complementary information carried by seismic wavefields, which are fundamentally different in terms of acquisition setting. We cast this information into our processing workflow and build a detailed model of the subduction zone, which is subject to further geological interpretation.