1. A computer-implemented method for accounting for temporal dispersion in low-order finite difference seismic wave propagation, comprising: a. transforming a seismic dataset from time domain to frequency domain to obtain a frequency-domain seismic dataset; b. applying a frequency-domain time varying filter to the frequency-domain seismic dataset to obtain a filtered frequency-domain seismic dataset; and c. transforming the filtered frequency-domain seismic dataset from the frequency domain to the time domain to obtain a time-domain filtered seismic dataset.