We apply two hybrid methods for solving scattering problems affected by resonances, to a four-dimensional reactive surface scattering system. In each method the solution of the problem is divided into two parts: a wave packet propagation, and a resonance calculation; results of the resonance calculation are used to extrapolate the long-time behavior of the system. In the first hybrid method, the propagation is by the multistep Chebyshev method, with calculation of resonances performed by the Lanczos method. In the second, the propagation is done using an implementation of the absorbing boundary condition (ABC) evolution operator, and the resonance calculation by filter diagonalization (FDG). Each method produces accurate scattering results in much less computation time than standard long-time wave packet propagation. The Chebyshev-Lanczos approach proves most capable for the calculation of resonances, but is computationally expensive. The ABC-FDG method is much cheaper to implement, but could not be made to extract accurate data for certain broad, overlapping resonances. This was overcome by propagating longer (still much shorter than for long-time propagation) to allow the elusive resonances time to decay.