A method is disclosed that uses (a) source-receiver reciprocity and (b) a method such as ISR DMO (U.S. Patent Application Publication No. 2006/0098529) that allows reconstruction of densely sampled gathers at arbitrary surface locations, to efficiently predict multiple reflections, either surface related or interbed multiples, in seismic data. For each reconstructed gather and each output (field) trace, two traces are extracted corresponding to the field source and receiver positions of the output trace, then they are convolved and the convolution result is added (summed) to that obtained by applying this procedure to previously reconstructed gathers. The efficiency results from the fact that the convolutions are performed by looping over all traces for each “bounce” point, with the outer loop being over bounce points. Once all the reconstructed gathers are processed, multiple predictions are obtained for the whole survey by conventional means.