The study presents and solves the 2-D elastodynamic model for seismic in-plane wave propagation in laterally inhomogeneous geological profiles imbedded in a vertically inhomogeneous half-space in which an earthquake source is buried. To this end, an efficient hybrid modal summation- boundary integral equation method (MSM-BIEM) is developed and applied. The MSM is used as a tool for simulating wave propagation from the source position to the multilayered laterally inhomogeneous geological profile where the BIEM is applied. The proposed model and the hybrid tool are used to investigate the phenomena of site effects. In fact, such a methodology has the potential to investigate the combined effects of different physical phenomena like surface topography, lateral inhomogeneity and the existence of water saturation in soils on the estimation of site effects. The model and hybrid computational tool developed are applied to contribute to the seismic risk analysis of the Bulgarian capital Sofia.
