We have developed a new methodology for predicting and removing multiples in the migration depth domain based on wavefield extrapolation and attribute based subtraction. The input for the prediction part is a 3D prestack depth migration volume and the corresponding velocity field. The output is the predicted multiple model (surface related or inter-bed) in the migration depth domain. The residual multiple removal technique combines the multiple prediction part with the recently developed attribute-based subtraction. Applications to both marine and land data have proven this methodology to be very effective in further reducing the residual multiples in the final migration images.