Seismic Wave Arrival Time Calculator Back
Earth Science & Seismic Engineering

Seismic Wave Arrival Time Calculator

Estimate epicentral distance from the speed difference between P-waves and S-waves. Visualize the Omori formula and the mechanism behind earthquake early warning systems.

Parameter Settings

Seismic Wave Propagation (Epicenter → Station)
P arrival t_P [s]
S arrival t_S [s]
S-P interval Δt [s]
Epicenter distance d [km]

The fast P-wave front (blue) and slower S-wave (red) expand from the epicenter; the station records P first, then S. The travel-time graph below plots distance versus time as two lines. Distance is recovered with Omori's relation d=(t_S−t_P)/(1/V_S−1/V_P).
Verification: V_P=8, V_S=4.6 km/s, d=100 km → t_P=12.5 s, t_S=21.74 s, Δt=9.24 s, recovered d=100.0 km.

Results
P-wave arrival time
25.0 s
S-wave arrival time
42.9 s
P-S time difference Δt
17.9 s
Omori distance estimate
133 km
Travel-Time Curves
Seismogram
Vp/Vs Ratio Analysis
Travel

Blue: P-wave travel-time curve. Red: S-wave travel-time curve. The vertical guide marks the selected source distance.

🎓
Because P-waves are faster than S-waves, the S-wave is still some distance away when nearby seismometers first detect the P-wave. An early-warning network can estimate the source location, depth, and magnitude within seconds, then send alerts to phones and broadcasts before the stronger S-wave arrives. Places farther from the source get more warning time, while locations directly above a shallow source may get almost none.
🙋
The Omori formula calculates distance from the initial motion duration, but where does the number 7.42 come from?
Theory & Key Formulas

$t_P = D / V_P$, $\quad t_S = D / V_S$

Initial motion duration (P-S time difference)
$\Delta t = t_S - t_P = D \left(\dfrac{1}{V_S} - \dfrac{1}{V_P}\right)$

Omori formula (approximation)
$D \approx 7.42 \cdot \Delta t$ [km]
🎓
Fusakichi Omori derived the empirical relation from Japanese earthquake records around 1900. Starting with the P-S time difference, $\Delta t = D(1/V_S - 1/V_P)$, solving for $D$ gives $D = \Delta t / (1/V_S - 1/V_P)$. If Vp is about 6 km/s and Vs about 3.5 km/s, the coefficient is close to 7.42. Modern early-warning systems use more detailed crustal velocity models, but Omori's formula is still useful for education and quick estimates.
🙋
Can knowledge of seismic waves be used in structural analysis with CAE?
🎓
Yes. Seismic design compares the input response spectrum with the building's natural period, and that spectrum depends strongly on site amplification controlled by S-wave velocity, often summarized by VS30. In Abaqus or Ansys seismic-response analysis, engineers may use a soil velocity profile and time-history ground motion as input. In exploration geophysics, P-wave reflection data are also used to image subsurface structure, almost like building a CAE mesh of the ground.
🙋
Why don't S-waves travel through liquids?
🎓
S-waves rely on shear deformation. Solids have shear rigidity, so neighboring particles can pull each other sideways and transmit a transverse wave. Fluids have essentially zero shear rigidity, so they cannot carry S-waves. This fact helped reveal Earth's internal structure: the S-wave shadow zone showed that Earth has a liquid outer core made mostly of molten iron and nickel.

Frequently Asked Questions

What does the P-wave to S-wave velocity ratio (Vp/Vs) indicate?

Vp/Vs is directly related to Poisson's ratio: $\nu = (Vp^2 - 2Vs^2) / (2(Vp^2 - Vs^2))$. The ratio changes with rock type and fluid content, such as granite with Vp/Vs around 1.73 and saturated sandstone often above 2. Seismic tomography uses this ratio to detect magma reservoirs and fluid-rich zones.

How is earthquake focal depth estimated?

Focal depth is estimated by inverting arrival-time differences from multiple seismometers in three dimensions. Source distance combines epicentral distance and depth as $D = \sqrt{r^2 + h^2}$. Deep earthquakes, often deeper than 300 km, occur within subducting plates and define Wadati-Benioff zones.

Can magnitude be estimated from the initial P-S interval?

No. The P-S interval mainly estimates distance. Magnitude estimation requires P-wave amplitude, peak displacement, and frequency content near the onset. Early-warning systems combine magnitude estimates from the first few seconds of P-wave data with source-distance estimates to predict shaking intensity.

What is seismic tomography?

Seismic tomography images Earth's three-dimensional velocity structure using travel-time data from many earthquake-station pairs. It is similar in spirit to a medical CT scan. Cold subducting slabs usually appear as high-velocity regions, while hot mantle plumes appear as low-velocity regions.

How accurate is this calculator?

This tool assumes a uniform one-layer velocity structure. The real Earth has crust, upper mantle, transition zone, lower mantle, outer core, and inner core layers, so seismic waves refract and reflect at depth. Precise travel-time calculations use models such as Jeffreys-Bullen tables or the IASP91 velocity model.

What is Seismic Wave Arrival Time?

Seismic Wave Arrival Time is a fundamental topic in engineering and applied physics. This interactive simulator lets you explore the key behaviors and relationships by directly manipulating parameters and observing real-time results.

By combining numerical computation with visual feedback, the simulator bridges the gap between abstract theory and physical intuition — making it an effective learning tool for students and a rapid-verification tool for practicing engineers.

Physical Model & Key Equations

The simulator is based on the governing equations behind Seismic Wave Arrival Time Calculator. Understanding these equations is key to interpreting the results correctly.

Each parameter in the equations corresponds to a slider in the control panel. Moving a slider changes the equation's solution in real time, helping you build a direct connection between mathematical expressions and physical behavior.

Real-World Applications

Engineering Design: The concepts behind Seismic Wave Arrival Time Calculator are applied across mechanical, structural, electrical, and fluid engineering disciplines. This tool provides a quick way to estimate design parameters and sensitivity before committing to full CAE analysis.

Education & Research: Widely used in engineering curricula to connect theory with numerical computation. Also serves as a first-pass validation tool in research settings.

CAE Workflow Integration: Before running finite element (FEM) or computational fluid dynamics (CFD) simulations, engineers use simplified models like this to establish physical scale, identify dominant parameters, and define realistic boundary conditions.

Common Misconceptions and Points of Caution

Model assumptions: The mathematical model used here relies on simplifying assumptions such as linearity, homogeneity, and isotropy. Always verify that your real system satisfies these assumptions before applying results directly to design decisions.

Units and scale: Many calculation errors arise from unit conversion mistakes or order-of-magnitude errors. Pay close attention to the units shown next to each parameter input.

Validating results: Always sanity-check simulator output against physical intuition or hand calculations. If a result seems unexpected, review your input parameters or verify with an independent method.

How to Use

  1. Enter epicenter distance (0–500 km) using distSlider or input distVal directly
  2. Set P-wave velocity (vpVal): typical range 5.0–7.0 km/s for crustal rock; granite averages 6.0 km/s
  3. Set S-wave velocity (vsVal): typically 58–70% of P-wave velocity; adjust vpSlider and vsSlider to observe arrival time separation
  4. Simulator calculates arrival times using tp = distance/vp and ts = distance/vs, then displays lag time (ts – tp) in seconds

Worked Example

Distance = 150 km, vp = 6.0 km/s (granite crust), vs = 3.5 km/s. P-wave arrival: 150/6.0 = 25.0 s. S-wave arrival: 150/3.5 = 42.9 s. Lag (S–P interval) = 17.9 seconds. Seismograph networks use this 17.9 s offset to triangulate the epicenter; a 1 km error in distance estimate produces ~0.17 s timing error.

Practical Notes

  1. S–P lag time scales linearly with distance; doubling distance doubles the separation, enabling rapid field estimation before digital networks process data
  2. Velocity varies by geology: oceanic basalt vp ≈ 6.8 km/s vs. sedimentary layers vp ≈ 4.5 km/s; adjust parameters to match local lithology
  3. Early warning systems exploit this lag; a 200 km epicenter yields ~34 s warning before S-wave damage at the station
  4. Poisson's ratio affects vs/vp ratio; soft sediments (ratio 0.65) versus rigid shield cratons (ratio 0.58) shift arrival curves noticeably