Because the Earth’s crust is quite heterogeneous and consists of various types of rocks, many faults are located near material interfaces. Recent observational studies have revealed that some faults cross the interfaces, but theoretical treatment of their dynamic behaviour is still undeveloped. In this paper, we develop a boundary integral equation method (BIEM) to analyse the dynamic behavio…