Love wave

In elastodynamics, Love waves, named after Augustus Edward Hough Love, are horizontally polarized surface waves. The Love wave is a result of the interference of many shear waves (S-waves) guided by an elastic layer, which is welded to an elastic half space on one side while bordering a vacuum on the other side. In seismology, Love waves (also known as Q waves (Quer: German for lateral)) are surface seismic waves that cause horizontal shifting of the Earth during an earthquake. Augustus Edward Hough Love predicted the existence of Love waves mathematically in 1911. They form a distinct class, different from other types of seismic waves, such as P-waves and S-waves (both body waves), or Rayleigh waves (another type of surface wave). Love waves travel with a lower velocity than P- or S- waves, but faster than Rayleigh waves. These waves are observed only when there is a low velocity layer overlying a high velocity layer/ sub–layers.

Description
The particle motion of a Love wave forms a horizontal line perpendicular to the direction of propagation (i.e. are transverse waves). Moving deeper into the material, motion can decrease to a "node" and then alternately increase and decrease as one examines deeper layers of particles. The amplitude, or maximum particle motion, often decreases rapidly with depth.

Since Love waves travel on the Earth's surface, the strength (or amplitude) of the waves decrease exponentially with the depth of an earthquake. However, given their confinement to the surface, their amplitude decays only as $$\frac{1}{\sqrt{r}}$$, where $$r$$ represents the distance the wave has travelled from the earthquake. Surface waves therefore decay more slowly with distance than do body waves, which travel in three dimensions. Large earthquakes may generate Love waves that travel around the Earth several times before dissipating.

Since they decay so slowly, Love waves are the most destructive outside the immediate area of the focus or epicentre of an earthquake. They are what most people feel directly during an earthquake.

In the past, it was often thought that animals like cats and dogs could predict an earthquake before it happened. However, they are simply more sensitive to ground vibrations than humans and are able to detect the subtler body waves that precede Love waves, like the P-waves and the S-waves.

Basic theory
The conservation of linear momentum of a linear elastic material can be written as
 * $$\boldsymbol{\nabla}\cdot(\mathsf{C}:\boldsymbol{\nabla}\mathbf{u}) = \rho~\ddot{\mathbf{u}} $$

where $$\mathbf{u}$$ is the displacement vector and $$\mathsf{C}$$ is the stiffness tensor. Love waves are a special solution ($$\mathbf{u}$$) that satisfy this system of equations. We typically use a Cartesian coordinate system ($$x,y,z$$) to describe Love waves.

Consider an isotropic linear elastic medium in which the elastic properties are functions of only the $$z$$ coordinate, i.e., the Lamé parameters and the mass density can be expressed as $$\lambda(z), \mu(z), \rho(z)$$. Displacements $$(u,v,w)$$ produced by Love waves as a function of time ($$t$$) have the form

u(x,y,z,t) = 0 ~, v(x,y,z,t) = \hat{v}(x,z,t) ~, w(x,y,z,t) = 0 \,. $$ These are therefore antiplane shear waves perpendicular to the $$(x,z)$$ plane. The function $$\hat{v}(x,z,t)$$ can be expressed as the superposition of harmonic waves with varying wave numbers ($$k$$) and frequencies ($$\omega$$). Consider a single harmonic wave, i.e.,

\hat{v}(x,z,t) = V(k, z, \omega)\,\exp[i(k x - \omega t)] $$ where $$i$$ is the imaginary unit, i.e. $$i^2 = -1$$. The stresses caused by these displacements are

\sigma_{xx} = 0 ~, \sigma_{yy} = 0 ~, \sigma_{zz} = 0 ~, \tau_{zx} = 0 ~, \tau_{yz} = \mu(z)\,\frac{dV}{dz}\,\exp[i(k x - \omega t)] ~, \tau_{xy} = i k \mu(z) V(k, z, \omega) \,\exp[i(k x - \omega t)] \,. $$ If we substitute the assumed displacements into the equations for the conservation of momentum, we get a simplified equation

\frac{d}{dz}\left[\mu(z)\,\frac{dV}{dz}\right] = [k^2\,\mu(z) - \omega^2\,\rho(z)]\,V(k,z,\omega) \,. $$ The boundary conditions for a Love wave are that the surface tractions at the free surface $$(z = 0)$$ must be zero. Another requirement is that the stress component $$\tau_{yz}$$ in a layer medium must be continuous at the interfaces of the layers. To convert the second order differential equation in $$V$$ into two first order equations, we express this stress component in the form

\tau_{yz} = T(k, z, \omega)\,\exp[i(k x - \omega t)] $$ to get the first order conservation of momentum equations

\frac{d}{dz}\begin{bmatrix} V \\ T \end{bmatrix} = \begin{bmatrix} 0 & 1/\mu(z) \\ k^2\,\mu(z) - \omega^2\,\rho(z) & 0 \end{bmatrix} \begin{bmatrix} V \\ T \end{bmatrix} \,. $$ The above equations describe an eigenvalue problem whose solution eigenfunctions can be found by a number of numerical methods. Another common, and powerful, approach is the propagator matrix method (also called the matricant approach).