A method for processing seismic data, the method including receiving seismic data acquired by at least one receiver over a water-covered subsurface formation; generating, based on the seismic data, a down-going wavefield; generating, based on the seismic data, a partial down-going wavefield with attenuated water-wave; estimating a subsurface reflectivity R using multi-dimensional deconvolution, which equates (1) a convolution of the down-going wavefield with the subsurface reflectivity R to (2) the partial down-going wavefield; and generating an image of the water-covered subsurface formation based on the subsurface reflectivity R.