Two-dimensional coupled seismic waves, satisfying the equations of linear isotropic elasticity, on a rectangular domain with initial conditions and periodic boundary conditions, are considered. A quantity conserved by the solution of the continuous problem is used to check the numerical solution of the problem. Second order spatial derivatives, in the x direction, in the y direction and mixed derivative, are approximated by finite differences on a uniform grid. The ordinary second order in time system obtained is transformed into a first order in time system in the displacement and velocity vectors. For the time integration of this system, second order and fourth order exponential splitting methods, which are geometric integrators, are proposed. These explicit splitting methods are not unconditionally stable and the stability condition for time step and space step ratio is deduced. Numerical experiments displaying the good behavior in the long time integration and the efficiency of the numerical solution are provided.
