Summary: This paper is concerned with time-harmonic elastic scattering from a locally perturbed rough surface in two dimensions. We consider a rigid scattering interface given by the graph of a one-dimensional Lipschitz function which coincides with the real axis in the complement of some compact set. Given the incident field and the scattering interface, the direct problem is to determine the field distribution, whereas the inverse problem is to determine the shape of the interface from the measurement of the field on an artificial boundary in the upper half-plane. We propose a symmetric coupling method between finite element and boundary integral equations to show uniqueness and existence of weak solutions. The synthetic data is computed via the finite element method with the Perfectly Matched Layer (PML) technique. To investigate the inverse problem, we derive the domain derivatives of the field with respect to the scattering interface. An iterative continuation method with multi-frequency data is used for recovering the unknown scattering interface.