Love wave

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 ${\displaystyle {\frac {1}{\sqrt {r}}}}$, where ${\displaystyle 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.^[1]

This section needs additional citations for verification. Please help improve this article by adding citations to reliable sources in this section. Unsourced material may be challenged and removed. (February 2020) (Learn how and when to remove this message)

The conservation of linear momentum of a linear elastic material can be written as:^[2]

${\displaystyle {\boldsymbol {\nabla }}\cdot ({\mathsf {C}}:{\boldsymbol {\nabla }}\mathbf {u} )=\rho ~{\ddot {\mathbf {u} }}}$

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

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

${\displaystyle 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 ${\displaystyle (x,z)}$ plane. The function ${\displaystyle {\hat {v}}(x,z,t)}$ can be expressed as the superposition of harmonic waves with varying wave numbers (${\displaystyle k}$) and frequencies (${\displaystyle \omega }$). Consider a single harmonic wave, i.e.,

${\displaystyle {\hat {v}}(x,z,t)=V(k,z,\omega )\,\exp[i(kx-\omega t)]}$

where ${\displaystyle i}$ is the imaginary unit, i.e. ${\displaystyle i^{2}=-1}$. The stresses caused by these displacements are

${\displaystyle \sigma _{xx}=0~,~~\sigma _{yy}=0~,~~\sigma _{zz}=0~,~~\tau _{zx}=0~,~~\tau _{yz}=\mu (z)\,{\frac {dV}{dz}}\,\exp[i(kx-\omega t)]~,~~\tau _{xy}=ik\mu (z)V(k,z,\omega )\,\exp[i(kx-\omega t)]\,.}$

If we substitute the assumed displacements into the equations for the conservation of momentum, we get a simplified equation

${\displaystyle {\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 ${\displaystyle (z=0)}$ must be zero. Another requirement is that the stress component ${\displaystyle \tau _{yz}}$ in a layer medium must be continuous at the interfaces of the layers. To convert the second order differential equation in ${\displaystyle V}$ into two first order equations, we express this stress component in the form

${\displaystyle \tau _{yz}=T(k,z,\omega )\,\exp[i(kx-\omega t)]}$

to get the first order conservation of momentum equations

${\displaystyle {\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).^{[citation needed]}

• Longitudinal wave
• Antiplane shear