This paper presents a step-by-step procedure using the three-dimensional boundary element approach to study the behavior of semi-circular canyons under seismic shear waves. The boundary element code TDASC allows utilization for various canyon geometries, evaluation of concurrent seismic waves and calculation of the ground motions on canyons due to an excitation at any arbitrary point of the incident field. Considering the widening ratio of the canyon (including prismatic, semi-prismatic and non-prismatic canyons), wave characteristics (wavelength, dimensionless period, direction) and maximum amplification pattern, the solution was applied to carry out a series of parametric studies. It was shown that canyon form can significantly affect the displacement amplification, especially at the points located on its edges. By increasing the wave dimensionless frequency (η > 1), the amplification pattern becomes more complex. On the basis of the results from a variety of considered cases, a new expression has been presented for the limiting wavelength beyond which the widening of the canyon will not have a major effect on the displacement amplification. To verify the reliability of the proposed approach, the obtained results, expressed in terms of displacement amplitude, were compared with those from the available published literature and a reasonably good agreement was observed.