We consider the numerical approximation of retarded potentials arising in wave scattering problems. For the temporal approximation we consider a generalization of Lubich's Convolution Quadrature which allows for variable steps. The algorithmic realization of the new method relies on contour integral techniques in the complex plane. Numerical experiments are provided to show the potential of our approach.