“…The demultiple approach of Amundsen (2001) is a powerful processing step that can jointly suppress the effect of any type and order of free-surface multiples on both source-and receiver-side, apply source signature deconvolution (removing any residual phase in data), and deblend simultaneoussource data. Most implementations approximate multi-dimensional deconvolution of the up-going component with the down-going component by applying a deterministic trace-by-trace deconvolution in tau-p (or f-k) domain by assuming the geology to be horizontally layered (Mispel et al, 2007;Traub et al, 2009). We have instead shown that equation 1 can be successfully solved by direct matrix inversion thus accounting for the three-dimensional nature of the seismic wavefield and lateral heterogeneities in the shallow surface.…”