Dips are calculated for a series of sets of adjacent multiple contribution traces, from seismic data representative of subsurface formations, in the vicinity of a boundary of an aperture of a multiple contribution gather, the seismic data acquired by deploying a plurality of seismic sensors proximate an area of the earth's subsurface to be evaluated, the seismic sensors generating at least one of an electrical and optical signal in response to seismic energy. The boundary of the aperture of the multiple contribution gather is recursively extended, based on the calculated dips. Multiple contribution traces in the multiple contribution gather with the extended aperture boundary are stacked to generate predicted multiple traces. The predicted multiple traces are subtracted from the seismic data to generate surface-related multiple attenuated data useful for imaging the subsurface formations.