Seismic Wave Propagation Simulator Back
Earth Science & Earthquake Engineering

Seismic Wave Propagation Simulator

Real-time animation of P-wave and S-wave propagation, reflection, and refraction through Earth's interior layers. Visualize the shadow zone and core-mantle boundary effects.

Parameters
Epicenter Distance
°
Magnitude
Presets
Arrival Times
P-wave
S-wave
Crust
Mantle
Outer Core
Inner Core
Results
P-wave arrival (s)
S-wave arrival (s)
P-S gap (s)
None
Shadow Zone
Seismic
Theory & Key Formulas
At each layer boundary, waves refract according to:
sin(i) / v₁ = sin(r) / v₂
S-waves vanish in the liquid outer core; only P-waves pass through.

What is Seismic Wave Propagation?

🙋
What exactly are P-waves and S-waves? I see them shooting out from the epicenter in the simulator, but they look different.
🎓
Basically, they're the two main types of body waves that travel through the Earth. P-waves (Primary) are compressional waves, like a slinky being pushed and pulled. They travel fastest. S-waves (Secondary) are shear waves, moving side-to-side, like shaking a rope. In the simulator, try setting the "Epicenter Distance" to a large value. You'll see the P-wave arrive first, followed by the S-wave, creating a time gap.
🙋
Wait, really? So why does that big "shadow zone" appear on the opposite side of the Earth? The waves just... stop?
🎓
Great observation! They don't just stop. The Earth isn't uniform; it has layers with different densities and states. The key is the liquid outer core. S-waves can't travel through liquids, so they vanish there. P-waves are refracted (bent) sharply at the core boundary. This combination creates two shadow zones where direct waves don't reach. Try increasing the "Magnitude" slider—the wavefronts get brighter, but the shadow zone pattern remains the same because it's a structural feature.
🙋
So the bending is Snell's Law, like with light? How does that help us figure out what's inside the Earth if we can't dig there?
🎓
Exactly! In practice, seismologists are like doctors doing a CT scan of the planet. By measuring the arrival times and paths of waves from earthquakes all over the globe, they can work backwards. The specific angles of refraction and reflection, which you see animated at each layer in the simulator, tell us the depths and properties of the crust, mantle, and core. Changing the epicenter distance shows how the wave paths through these layers change.

Physical Model & Key Equations

The fundamental behavior governing how seismic waves change direction at a boundary between two layers with different wave speeds is described by Snell's Law.

$$ \frac{\sin(i)}{v_1}= \frac{\sin(r)}{v_2}$$

Here, $i$ is the angle of incidence, $r$ is the angle of refraction, and $v_1$ and $v_2$ are the seismic wave velocities in the first and second layer, respectively. If $v_2 \gt v_1$, the wave bends away from the normal line.

The existence of the shadow zone is a direct consequence of two physical principles combined with Earth's layered structure.

$$ \text{S-wave speed in fluid}= 0 \quad \text{and}\quad \text{P-wave path deflection} \gt 103^\circ $$

S-waves require shear strength to propagate, which a fluid lacks, so they are absorbed at the core-mantle boundary. The P-wave is refracted so strongly at the core boundary that it emerges on the other side at a specific minimum angle, leaving a gap (the shadow zone) between 103° and 143° from the epicenter where no direct P-waves are detected.

Real-World Applications

Earthquake Early Warning Systems: The predictable difference in speed between P and S-waves is exploited for early warnings. Sensitive detectors can feel the fast, less-damaging P-wave and trigger alarms seconds before the destructive S-wave arrives, allowing trains to brake and gas valves to shut.

Oil & Gas Exploration (Seismic Reflection Surveying): The same principles are used on a smaller scale. Trucks send artificial seismic waves into the ground, and the reflected signals are analyzed to map subsurface rock layers and identify potential hydrocarbon reservoirs, essentially creating an ultrasound of the Earth's crust.

Nuclear Test Ban Treaty Verification: A global network of seismic stations monitors for the unique seismic signatures of underground nuclear explosions, which differ from natural earthquakes. Understanding wave propagation is crucial to pinpoint the location and nature of the event.

Planetary Science: By placing seismometers on other celestial bodies (like the InSight lander on Mars), scientists can use the recorded seismic waves to infer the internal structure, composition, and layer thicknesses of these planets, just as we did with Earth.

Common Misconceptions and Points to Note

There are a few key points you should be aware of when starting to use this simulator. First, understand the premise that wave velocity is fixed per layer. In the actual Earth, pressure, temperature, and density change with depth even within the same mantle, causing velocity to increase continuously. The simulator simplifies this into "layers," which makes the abrupt refraction at boundaries appear more pronounced. Real data often shows smoother ray paths.

Next, don't confuse "epicentral distance" with "distance to the observation point". The simulator's "epicentral distance" is the angular distance from the epicenter to a point on the Earth's surface (e.g., 103 degrees). However, with a deep earthquake source, the path the wave takes before reaching the surface changes. For example, for an earthquake at 600 km depth, the shadow zone's range differs from the simulator's (which assumes a surface source). Always keep the source depth in mind.

Finally, avoid overinterpreting amplitude (magnitude) changes. While you can change the amplitude with the slider, this is a simplified representation. The amplitude of real seismic waves varies greatly not only with travel distance but also due to site amplification effects, scattering, and attenuation. Even if you increase the S-wave amplitude in the simulator, the fundamental principle that it cannot propagate through the liquid outer core remains unchanged. The key to learning is to focus on this qualitative difference of "whether it propagates or not."

How to Use

  1. Set epicentral distance (0–180°) using distSlider; typical regional earthquakes display 15–45° ranges
  2. Adjust magnitude on magSlider (4.0–8.0 Mw) to modify wave amplitude and energy release
  3. Observe real-time P-wave and S-wave fronts propagating through layered crust and mantle; P-wave velocity averages 6.0 km/s in upper crust versus S-wave at 3.5 km/s
  4. Monitor P-wave arrival time, S-wave arrival time, and P-S gap (used for distance calculation in seismology)
  5. Note shadow zone formation beyond 103° epicentral distance where S-waves are completely attenuated by the liquid outer core

Worked Example

A magnitude 6.8 earthquake occurs 45° away from a seismic station. P-wave velocity through continental crust is 6.1 km/s; S-wave velocity is 3.8 km/s. Arc distance = 45° × 111.2 km/degree ≈ 5,004 km. P-wave arrival ≈ 818 seconds (13.6 minutes); S-wave arrival ≈ 1,317 seconds (21.9 minutes). P-S gap = 499 seconds. At 103° epicentral distance, the shadow zone begins; S-wave energy reflects off the liquid core boundary while P-waves refract through it, creating the characteristic seismic gap used in real-time distance estimation.

Practical Notes

  1. P-S gap method: subtract P-arrival from S-arrival to estimate epicentral distance without knowing magnitude; 8-second gaps typically indicate ~60 km local earthquakes
  2. Magnitude variation affects waveform amplitude but not velocity; higher magnitudes (7.0+) generate visible surface waves (Love, Rayleigh) propagating at 4.0 km/s in far-field recordings
  3. Shadow zone absence of S-waves beyond 103° confirms outer core is liquid; use reflected phases (ScS, sScS bouncing off core-mantle boundary) for deep hypocenter studies
  4. In subduction zones (Wadati-Benioff geometry), angle-of-incidence changes alter reflection patterns; dip angles of 30–45° produce complex multiplet sequences