The scattering of elastic waves is studied in the vicinity of a vacuum-medium boundary. The Green's function for a half-space is re-derived within the mixed 2D-Fourier representation, which is convenient for studying layered media. Monte-Carlo simulations of elastic wave scattering from random inhomogeneities within a simplified scalar model are performed, accounting for a boundary-induced term in the Green's function. The multiply scattered elastic waves' radiation is shown to decay with distance from the source much slower in vicinity of boundary than in an infinite medium, due to the boundary condition requirements.