A method for modeling shallow water multiples modeling using a predicted sea floor reconstruction is provided. A seafloor reflectivity image is reconstructed using acquired seismic data. Surface multiples are determined by evaluating a seismic wave equation using, as an input, the reconstructed seafloor reflectivity image. The determined surface multiples are subtracted from the acquired seismic data.