From geodesic extrapolation to a variational BDF2 scheme for Wasserstein gradient flows