Stability Does Not Guarantee Accuracy: CFL-Compliant 1D Acoustic FDTD and Its Consequences for Near-Surface Layered Modeling
DOI:
https://doi.org/10.29303/goescienceed.v7i1.1693Keywords:
FDTD, CFL Condition, Numerical Dispersion, Waveform Accuracy, Near Surface ModelingAbstract
Explicit finite-difference time-domain (FDTD) solvers for the 1D acoustic wave equation are routinely configured by enforcing the Courant–Friedrichs–Lewy (CFL) stability limit. In near-surface settings, this practice can create a false sense of model fidelity because stability constrains time stepping but does not control phase and amplitude errors induced by numerical dispersion. This report isolates the gap between stable time marching and accurate waveform synthesis in layered, strongly contrasting near-surface media. Beyond the well-known stability–accuracy distinction, the specific contribution is a reproducible near-surface verification workflow that couples a transfer-matrix stratified benchmark with controlled bandwidth variation at fixed CFL and converts the resulting waveform misfits into an explicit dispersion budget and pass–fail acceptance gate. The intended deliverable is an independent verification note that demonstrates, with quantitative evidence, that CFL compliance is necessary for stability but insufficient for accuracy, and that near-surface forward modeling requires an explicit dispersion budget expressed in points per minimum wavelength and validated by waveform misfit diagnostics.
References
Courant, R., Friedrichs, K., & Lewy, H. (1967). On the partial difference equations of mathematical physics. IBM Journal of Research and Development, 11(2), 215–234. https://doi.org/10.1147/rd.112.0215
Gao, Y., Tilmann, F., & Rietbrock, A. (2023). A review of misfit functions for adjoint full waveform inversion in seismology. Geophysical Journal International, 235(3), 2794–2827. https://doi.org/10.1093/gji/ggad372
Huang, J.-P., Peng, W.-T., Yang, J.-D., & Lou, L.-F. (2024). Overview of computation strategies on the dispersion analysis for explicit finite difference solution of acoustic wave equation. Petroleum Science, 21(4), 2311–2328. https://doi.org/10.1016/j.petsci.2024.02.003
Kelly, K. R., Ward, R. W., Treitel, S., & Alford, R. M. (1976). Synthetic seismograms: A finite-difference approach. Geophysics, 41(1), 2–27. https://doi.org/10.1190/1.1440605
Liu, Y., Li, Z., Wang, J., Sun, M., & Liu, Q. (2021). A numerical dispersion-suppressed method for shallow seismic migration. Near Surface Geophysics, 19(1), 109–121. https://doi.org/10.1002/nsg.12134
Pled, F., & Desceliers, C. (2022). Review and recent developments on the perfectly matched layer (PML) method for the numerical modeling and simulation of elastic wave propagation in unbounded domains. Archives of Computational Methods in Engineering, 29, 471–518. https://doi.org/10.1007/s11831-021-09581-y
Rao, Y., & Wang, Y. (2019). Dispersion and stability condition of seismic wave simulation in TTI media. Pure and Applied Geophysics, 176, 1549–1559. https://doi.org/10.1007/s00024-018-2063-y
Tao, K., Grand, S. P., & Niu, F. (2017). Full-waveform inversion of triplicated data using a normalized-correlation-coefficient-based misfit function. Geophysical Journal International, 210(3), 1517–1524. https://doi.org/10.1093/gji/ggx249
Wu, B., Tan, W., Xu, W., & Li, B. (2022). Trapezoid-grid finite-difference time-domain method for 3D seismic wavefield modeling using CPML absorbing boundary condition. Frontiers in Earth Science, 9, Article 777200. https://doi.org/10.3389/feart.2021.777200
Xu, Q., Jin, C., Liang, K., & Cao, D. (2025). An optimized 13-point finite-difference operator for 2D frequency-domain acoustic wave modeling. Journal of Geophysics and Engineering, 22(4), 971–985. https://doi.org/10.1093/jge/gxaf051
