# 5.1: Basics of Wave Propagation

- Page ID
- 3540

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

To understand some of the more complex aspects of seismology, we must first start at the beginning and get a handle on the basics of wave propagation. In this section, we will examine three primary concepts:

- The Basics of Waves
- Types of Seismic Waves
- Optics: Reflection, Transmission (Refraction), and Snell's Law

## The Basics of Waves

In the figure above, \(\lambda\) is the wavelength in meters and \(A\) is the amplitude in \(\mu m-cm\). If you were to stand at x_{1} and watch the wave go by, you would see Figure \(\PageIndex{1}\):

where T is the period in s and \(f\) is the frequency in Hz. Frequency and period have the relation \(f=\frac{1}{T}\), and thus \(Hz=\frac{1}{s}\).

In terms of the relation between sound velocity and seismic velocity:

\[v=f\lambda\; \left[\dfrac{m}{s}\right]\]

Light has a similar relationship between frequency and wavelength,

\(c=f\lambda\), where different values for \(f\lambda\) give different colors of light.

Figure \(\PageIndex{2}\): Period

## Types of Seismic Waves

One category of seismic waves are body waves. Body waves are waves you have likely heard of before, P waves and S waves. P waves act like an accordion, and move parallel to the propagation direction.

S waves can have two components of motion, vertical and horizontal. Most S waves usually have both components, S-vert and S-horiz, which can be polarized.

The other category of seismic waves are surface waves, which decay with depth. Thus, as the name suggests, they are the strongest at the surface. There are two main surface waves, Love waves and Rayleigh waves. Love waves have particle motion like that of the S-H component.

Rayleigh waves have the components of a P wave and a S-V wave. The particle motion is retrograde.

How fast do these waves actually move? Wave speed depends on two things, the density of the material the waves are traveling through and the elastic constant or 'stiffness'. Density is key to wave speed because the wave has to move mass. Imagine moving two materials in a wave-like pattern, a string and a metal cable. It takes more energy to move the more massive metal cable than it does the string because the cable is much denser. To understand the elastic constant, imagine moving a yarn and a stiff rope. In the yarn, it is hard to get the wave to move because the material lacks stiffness. With these two properties, we can write a relationship for wave velocity:

\[v\propto\frac{\text{Elastic constant}}{\text{Density}}\]

As v\(\uparrow\), \(\rho\downarrow\), the elastic constant \(\uparrow\).

Additional useful relationships when dealing with stress-strain due to seismic waves are

\[\sigma=E\epsilon\]

where \(E\) is **Young's modulus** that we discussed in Chapter 1.

\[F=-k\Delta x\]

Where \(k\) is the spring constant.

\[\nu=\frac{\epsilon_1}{\epsilon_3}\]

In addition to stress-strain deformation, we also have **shear deformation**.

\[\epsilon_s=\frac{1}{2}\frac{\Delta x}{y}\]

\[\sigma_s=2G\epsilon_s\]

where \(G\) is the shear modulus.

\[G=\frac{\sigma_s}{2\epsilon_s}\left[\frac{N}{m^2}\right]\]

For reference, sometimes the shear modulus is represented by \(\mu\).

We can also have **volumetric strain**.

\[\frac{\Delta V}{V}=\frac{V_f-V_o}{V_o}\]

\[P=\kappa\frac{\Delta V}{V}\]

Where \(\kappa\) is the bulk Modulus

\[\kappa=\frac{P}{(\frac{\Delta V}{V})}[\frac{N}{m^2}]\]

To summarize and review some of the constants that we have learned so far:

- E-Young's Modulus (elongation dilation)
- \(\nu\)-Poisson's Ratio (compensation deformation)
- G or \(\mu\) - shear muodulus
- \(\kappa\)- bulk modulus

In seismic calculations, you can use one set (E and \(\nu\)) or the other (G and \(\kappa\)).

Now let's learn how to actually calculate seismic velocities. For P-waves:

\[v_p=\sqrt{\frac{\kappa+\frac{4}{3}G}{\rho}}\;or\;\sqrt{\frac{E}{\rho}\frac{(1-\nu)}{(1-2\nu)(1+\nu)}}\]

For S-waves, we can calculate velocity as:

\[v_s=\sqrt{\frac{G}{\rho}}\;or\;\sqrt{\frac{E}{\rho}\frac{1}{(2(1+\nu))}}\]

To find the ratio of P and S-wave velocities:

\[\frac{v_p}{v_s}=\sqrt{\frac{2(1-\nu)}{(1-2\nu)}}\]

Thus, we can see that this only depends on \(\nu\). This ratio is a good indicator of gas/liquids in exploration geology

- \(\frac{v_p}{v_s}<2\) indicates gas+sand
- \(\frac{v_p}{v_s}>2\) indicates sand only

In general, we find that v_{p}~1.7v_{s} (~50% faster).

Unfortunately, for surface waves, we cannot easily write down equations for \(v_R\) and \(v_L\) because of depth dependence. Thus, to determine their velocity, we must solve the wave equation with boundary conditions.

## Optics: Reflection, Transmission (Refraction), and Snell's Law

We can use the ray path to think about what happens to a wave when it encounters a boundary.

The most important information you can take from the above figure is that \(\theta_i=\theta_r\), i.e. the angle of incidence=the angle of reflection. We've now considered reflection of waves, but what about refraction?

\[T_{p-p'}=\frac{d}{v_1}+\frac{e}{v_1}=\frac{\sqrt{a^2+x^2}}{v_1}+\frac{\sqrt{b^2+(c-x)^2}}{v_2}\]

Fermat's Principle of least time states that the ray will follow the path that takes the __shortest__ time. We can use this principle to derive one of the most important laws in seismology, Snell's Law.

\[\begin{align*} \frac{dT}{dx}=0 \\[4pt] &=\frac{x}{v_1\sqrt{a^2+x^2}}-\frac{(c-x)}{v_2\sqrt{b^2+(c-x)^2}} \end{align*}\]

where \(\sin \theta_i=\frac{x}{\sqrt{a^2+x^2}}\) and \(\sin \theta_r=\frac{(c-x)}{\sqrt{b^2+(c-x)^2}}\)

Now we can rewrite the equation as:

\[\frac{dT}{dx}=\frac{\sin \theta_i}{v_1}-\frac{\sin \theta_r}{v_2}\]

\[\frac{\sin \theta_i}{v_1}=\frac{\sin \theta_r}{v_2}\]

The above equation is Snell's Law. To better understand Fermat's Principle and why a ray takes the path of shortest time, we can also say that \(\frac{\sin \theta_i}{v_1}\)=p=slowness. This slowness is constant along a ray path.

\[\frac{sin(i)}{v}=constant=p\]

The ray path curves to adjust for changing velocity between different mediums. Thus, the ray takes the shortest __time__ path. Let's now take a better look at Snell's Law, what it tells us, and what makes it so useful.

The figure shows the difference between reflection and refraction when a ray encounters a boundary.

\[\frac{\sin \theta_i}{v_1}=\frac{\sin \theta_{refr}}{v_2}\]

If we keep increasing v_{2}, \(\theta_r\) increases up to 90^{o}. At 90^{o}, we have critical refraction, where \(sin(90^o)=1\). Plugging this into Snell's Law we get:

\[\frac{\sin \theta_i}{v_1}=\frac{1}{v_2}\]

\[\theta_{ic}=sin^{-1}(\frac{v_1}{v_2})\]

If \(\theta_{ir}>\theta_{ic}\), we get total reflection. What does Snell's Law actually tell us? (slowness)

v increases with depth, i increases.

v decreases with depth, i decreases.

Figure \(\PageIndex{14}\): Changes in i with respect to v