Unstructured grid methods are well suited for earthquake problems in complex geometries, such as non-planar and branching faults. Unfortunately they are indefficient in comparsion wiht high-order finite differences. With the use of summation-by-parts (SBP) operators and the SAT penalty method (simultaneous approximation term) it is possible to couple unstructured finite volume methods with high-order finite difference methods in an accurate and stable way. The couple method is more efficient than the unstructured method alone. Another advantage of the SBP and SAT method is that it is possible to prove strict stability, meaning that the semi-discrete solution dissipates energy at a slightly faster rate than the continuous solution so that the error remains bounded in time, which is particulary useful for long time computations.