Particle motion sensor signals and the pressure sensor signals data from a towed marine seismic streamer are combined to generate an up-going pressure wavefield component and a down-going particle motion wavefield component. The down-going particle motion wavefield component is extrapolated from the receiver position depth level to the source position depth level. The up-going pressure wavefield component is multiplied by the extrapolated down-going particle motion wavefield component, generating a first product. Then, n.sup.th order surface related multiples in the pressure wavefield are iteratively calculated utilizing a product of (n-1).sup.th surface related multiple free data and the extrapolated down-going particle motion wavefield component. The calculated n.sup.th order surface related multiples are iteratively subtracted from the recorded pressure wavefield, generating the n.sup.th order surface related multiple free data.