The system under consideration is governed by the equation u_tt = ∇ div u – rot µ rot u in Ω × (0, T); its response operator ("input output" map) R^T plays the role of the inverse data. As in the case of the proper Lamé system, such a model describes a dynamical system with two wave modes (p-waves and s-waves) propagating with different velocities c_p = and c_s = correspondingly. We show that R^2T determines and , where and are the subdomains of Ω filled (at the moment T) with p- and s-waves propagating from ∂Ω. Due to the wave splitting u = ∇p + rot s the problem is reduced to the inverse problems for the acoustical and Maxwell subsystems governed by the equations p_tt = Δp and s_tt = −µ rot rot s with the response operators and determined by R^2T . The first problem can be solved by the BC method (Belishev, 1986), the second one is solved by a version of the method based on a blow up effect. This version is the main subject of the paper. In addition, we derive the inequality, which can be used for approximate determination of the shape of .