Skip to main content

Snell's Law Seismic Ray Parameter Calculator

In a 1D layered or vertically inhomogeneous Earth model, seismic body waves bend continuously according to Snell's law of refraction.

P-wave velocity at the surface takeoff point.

Initial ray angle measured from the vertical downward normal.

Wave velocity at target depth (e.g. upper mantle).

Calculated Result
0.0989 s/km

Seismic Ray Parameter (p)

Horizontal Phase Velocity

10.11 km/s

Bottom Turning Velocity

10.11 km/s

Angle at Target Depth

54.2°

Calculation Breakdown

  1. p = sin(i₀) / V₀0.0989 s/km
  2. V_phase = 1 / p10.11 km/s
  3. i_z = arcsin(p · V_z)54.2°

What Is the Snell's Law Seismic Ray Parameter Calculator?

The ray parameter p represents the invariant horizontal phase slowness along a seismic ray.

It links the surface epicentral distance and traveltime derivative (dT/dDelta).

How Does the Snell's Law Seismic Ray Parameter Calculator Work?

Evaluates ray parameter p from initial surface velocity and takeoff angle.

Computes horizontal phase velocity and turning point velocity.

Evaluates ray angle of incidence at deeper structural velocity horizons.

Snell's Law Seismic Ray Parameter Calculator Formula & Variables

The core mathematical equation utilized by this calculator is expressed as:

p = \frac{\sin(i_0)}{V_0} = \frac{\sin(i_z)}{V_z}, \quad V_{turn} = \frac{1}{p}, \quad i_z = \arcsin(p \cdot V_z)

Conservation of horizontal slowness along a seismic acoustic ray trajectory.

How to Use the Snell's Law Seismic Ray Parameter Calculator

  1. Enter the surface seismic wave velocity in km/s.
  2. Specify the initial takeoff angle in degrees from vertical.
  3. Enter the seismic velocity at the depth of interest.

Step-by-Step Example Calculation

Crust-to-Mantle P-Wave Refraction

Input Values:

surfaceVelocityKmS:5.8
takeoffAngleDeg:35
targetDepthVelocityKmS:8.2
Worked Steps: Deep seismic ray bending across the Mohorovicic discontinuity.

Understanding Your Result

Ray parameter p stays constant along the ray.

The turning velocity V_turn is the velocity where the ray becomes purely horizontal (i = 90 deg) and turns back toward the surface.

Factors That Affect the Result

  • Velocity gradients: Positive velocity gradients cause rays to curve concavely upwards.
  • Low velocity zones: Regions with decreasing velocity trap rays and cause shadow zones.

When Should You Use This Calculator?

  • Global seismic tomography and traveltime curve modeling.
  • Determining turning depths for teleseismic P and S body wave phases.

Assumptions & Limitations

  • Assumes flat, isotropic, vertically varying 1D Earth model.
  • Does not incorporate spherical Earth Earth-flattening transformations.

Frequently Asked Questions

Calculation Accuracy & Reference Note

Snell refraction invariance is exact in geometric ray theory.

Explore more tools and calculators in Physics Calculators