Numerical method for the computation of first-arrival time of seismic waves with general anisotropy

We present a numerical method for computing the first arrival time of seismic waves. The scheme can tackle three dimensional media with elastic properties defined by a fully anisotropic Hooke tensor. Our numerical method proceeds in a single pass over the domain, similar to the Fast Marching method, and achieves third order accuracy in smooth test cases, with quasi-linear computation time. We also validate our solver by comparing it with the solution to the elastic wave equation on a medium with general anisotropy.