LAUSR

High-order staggered-grid finite-difference (SFD) schemes are preferred in elastic wave simulation for geophysical problems because they decrease the accumulation of error from grid dispersion. However, most SFD approaches reach high-order spatial but limited temporal accuracy. To tackle the issue, we develop a novel temporal and spatial high-accuracy elastic SFD scheme by selectively modifying the spatial operators of the original SFD stencil. This modification has three main advantages. First, it facilitates the design of a new SFD stencil with temporal and spatial accuracies to arbitrary even-order by a Taylor-series expansion method. Second, it helps boost the accuracy further by implementing a linear optimization method. Third, the new selectively modified SFD (SMSFD) stencil needs fewer float-point operations (FPOs) than the existing temporal high-order SFD stencil. We compare our new SMSFD scheme with spatial high-order and temporal–spatial high-order SFD schemes and show that our new elastic SMSFD scheme possesses better accuracy and stability and requires fewer FPOs than these methods.