Abstract: Based on exact dynamic-stiffness matrices and Greenrsquo;s functions of a layered saturated soil, the scattering of SV waves by a canyon of arbitrary shape in a saturated layered half-space is modeled using the indirect boundary element method in frequency domain. The free-field responses are calculated to determine the displacement and stress at canyon surface, and fictitious distributed loads are then applied on the canyon surface in the free field to calculate the Greenrsquo;s functions for displacement and stress. The amplitude of the fictitious distributed loads are determined from the boundary conditions, and the displacements arising from the waves in the free field and from the fictitious distributed loads are summed to obtain the solution. The present method can be degenerated and the result agrees well with that of dry soil without fluid. The freefield responses of three different saturated sites are calculated and numerical results for displacements of the canyon due to incident SV waves are presented. The results of the saturated sites are compared with that of corresponding dry layered site. It is shown that the incident angle, the frequency of incident wave and penetrability of soils have a great effect on displacements. The frequency of incident wave and the internal friction parameter or the penetrability of soil influence the velocities of three waves propagated in layered saturated soils, which result in the apparent difference between the saturated layered half-space and the dry site, especially for the incident angles from 30deg;to 60deg;.The incident angle and frequency of incident wave should be comprehensively considered in dynamic analysis of the saturated sites.