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 60°
Magnitude M 6.5
Presets
Arrival Times
P-wave arrival (s)
S-wave arrival (s)
P-S gap (s)
None
Shadow Zone
P-wave
S-wave
Crust
Mantle
Outer Core
Inner Core

Snell's Law for Seismic Waves

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 > 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} > 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."

Related Engineering Fields

The concepts of wave propagation behind this simulator form the foundation for various fields that any CAE engineer should know. The first that comes to mind is Non-Destructive Testing (NDT). For instance, in ultrasonic flaw detection used to inspect internal defects in structures, longitudinal waves (similar to P-waves) and transverse waves (similar to S-waves) are used selectively to identify the type and location of flaws. Propagation within metals, including refraction and reflection at boundaries of different materials, follows the exact same physics as seismic waves.

Another is Geotechnical Engineering and Seismic Design of Civil Structures. Here, the "layered structure of the ground" you see in the simulator is directly applied. Calculations are made for how waves traveling through shallow, soft soil layers (analogous to the crust) refract at boundaries with hard bedrock (analogous to the mantle) and how they amplify at the surface. This is the starting point for ground response analysis, which evaluates whether these motions resonate with the natural period of buildings or bridges.

Broadening the perspective further, the principles also connect to Acoustical Engineering and Underwater Sonar. In water, S-waves do not propagate (because it's a liquid), so only P-waves (sound waves) travel. When a submarine's sonar detects objects in the ocean, it must account for layered sound velocity structures (like the SOFAR channel) due to differences in water temperature and salinity, which are also described by Snell's Law. The wave propagation expertise developed in geophysics is truly a common language for multiphysics simulation.

For Further Learning

Once you're comfortable with this simulator, try challenging yourself with the concept of the "Inverse Problem". The simulator calculates wave paths (the forward problem) assuming the Earth's structure (velocity model) is known. However, in practice, you often need to invert the source location or subsurface structure from the P-wave and S-wave arrival time data obtained at observation points. The basic method for solving this inverse problem is analyzing travel-time curves. The shape of the curve plotting arrival time against epicentral distance is used to infer the velocity structure. Your understanding will deepen significantly if you try sketching simple travel-time curves yourself while varying different parameters in the simulator.

If you want to delve a bit deeper into the mathematical background, explore the partial differential equations that describe waves. Seismic wave propagation is expressed by the wave equation derived from the equations of motion for an elastic medium. For a simple one-dimensional case, it takes this form: $$ \frac{\partial^2 u}{\partial t^2} = v^2 \frac{\partial^2 u}{\partial x^2} $$ Here, $u$ is displacement and $v$ is wave velocity. The next step is to investigate how the solutions to this equation take the form of traveling waves, and how the wave velocity $v$ connects to the formulas for $v_P$ and $v_S$ mentioned earlier. Understanding this will allow you to start reading papers on FEM (Finite Element Method) simulations of seismic waves.

Finally, as a concrete next topic, I recommend "Surface Waves (Rayleigh and Love Waves)". The simulator deals only with "body waves" (P and S waves) that travel through the Earth's interior. However, in actual earthquakes, the large shaking and long-period vibrations are caused by surface waves propagating near the surface. Analyzing the dispersive nature of these waves (where velocity changes with period) allows for detailed estimation of crustal thickness and shallow structure. Broadening your perspective from internal structure to surface phenomena will make the applications to earthquake engineering feel much more real.