LAUSR

ABSTRACT We present a new modified wave equation and apply it to develop a smoothing scheme for seismic wave propagation simulations. With mathematical rigour we show that the solution of the new equation, which is derived as an analog of the advection–diffusion equation, can be obtained by the spatial convolution between a solution of the wave equation and the heat kernel and has a finite propagation speed and a diffusion effect. Using numerical experiments we show that the smoothing scheme based on the modified wave equation has the following advantages. Firstly, it preserves the characteristics of the wave equation such as wave propagation speed. Secondly, it selectively removes the short-wavelength components of the solution. Lastly, the energy decreases slowly after the short-wavelength components have been removed. Since our smoothing scheme can be implemented by adding simple correction terms to usual schemes, it can easily be applied to the seismic wave equation.