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