“…Some of the best-known numerical approaches to assess the seismic response of topographic features can be categorized as follows: i) Finite difference method (FDM), which is a method with computational efficiency and simple code structure, however, could present difficulties with the modeling of irregular geometry, (Oprsal 2002;Republic 2007;Blanch 2012). ii) Finite element method (FEM), which can accurately capture the effect of complex topographies, such as mountainous areas, without refining the grid resolution, but has several disadvantages such as the need for huge power of calculation and the need to apply numerical damping (Amorosi et al 2010;Jahromi and Karkhaneh 2019). iii) Boundary element method (BEM), direct (DBEM) (Chaillat et al 2012;Kamalian et al 2003a, b;Mogi and Kawakami 2007;Sohrabibidar et al 2010), and indirect (IBEM) (Griffiths and Bollinger 1979;Sánchez-Sesma et al 1993;Gil-Zepeda et al 2003;Lee 2013): the reduction of problem dimensionality is the most important advantage of BEM, it means that only boundary line and surface domain are discretized for twodimensional plane and 3D program, respectively.…”