Abstract: Seismic modeling is one of the key problems in reverse time migration and full waveform inversion, because both of them require efficient and accurate simulation of forward and backward wave propagation. In order to improve the accuracy of numerical simulation, the high order finite difference method is widely used, but most of the methods have high-order accuracy in space while only second-order accuracy in time. In this paper, we first introduced basic principle of the finite difference method with high-order accuracy not only in space, but also in time. Then the model is utilized to verify the effectiveness of the method. The results show that the finite difference methods with high-order accuracy in time-space domains have lower dispersion than the conventional staggered-grid finite difference method.