This work is an exploratory study of coupling high-order time integrations to a finite-difference (FD) spatial discretization of the 1-D wave equation that combines eigth-order differencing at grid interior, with lateral formulas of order sixth and fourthat boundary neighborhood. This reduction of spatial accuracy at the grid vecinity of free surfaces is a known stability limitation of FD methods, when coupled to the two-step Leap-frog (LF) time stepping, which is widely used on seismic modeling. We first implement LF time integrations with an arbitrary accuracy order, as given from a standard Lax-Wendroff procedure, and compare results from the fourth-, sixth-, and twelfth order schemes, against the popular second-order LF. Our emphirical analyses establish the CFL stability constraints for propagation on an homogeneous medium, as first results, and then consider velocity heterogeneities when assessing dispersion and dissipation anomalies. Finally, we use a rapid expansion method (REM) to approximate the exponential of the semidiscrete FD discretization operator by a truncated Chebyshev matrix expansion. Althougth, REM has been previouslly applied to peudospectral (PS) wave simulations, this REM-FD scheme is the first reported in the technical literature according to our knowledge.


Article metrics loading...

Loading full text...

Full text loading...


  1. Blanch, J.O. and Robertsson, J.O.A.
    [1997] A modified Lax-Wendroff correction for wave propagation in media described by Zener elements. Geophysical Journal International, 131(2), 381–386.
    [Google Scholar]
  2. Pestana, R.C. and Stoffa, P.L.
    [2010] Time evolution of the wave equation using rapid expansion method. Geophysics, 75(4), T121–T131.
    [Google Scholar]
  3. Rojas, O., Dunham, E., Day, S.M., Dalguer, L.A. and Castillo, J.E.
    [2009] Finite difference modeling of rupture propagation with strong velocity-weakening friction. Geophysical Journal International, 179, 1831–1858.
    [Google Scholar]
  4. Tal-ezer, H., Kosloff, D. and Koren, Z.
    [1987] An accurate scheme for seismic forward modelling. Geophysical Prospecting, 35(5), 479–490.
    [Google Scholar]

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error