Subsurface imaging via fully coupled elastic wavefield extrapolation