“…In Section 2, we briefly state some relevant results in [21] that we are going to use later. In Section 3, based on the analysis in [21], one can see the Rayleigh waves are corresponding to the solution to Λu = l on the boundary, where l is a source microlocally supported in the elliptic region. Next, inspired by the diagonalization of the Neumann operator for the case of constant coefficients in [20], we diagonalize Λ microlocally up to smoothing operators by a symbol construction in [25].…”