Home / Science / Earth sciences / Geophysics and seismology / Seismic refraction interpreterTool #882

Seismic refraction: layer speeds and depths from first-arrival picks by intercept-time and crossover methods, dipping layers and travel-time curves

Paste first-arrival picks from a forward shot and, if you have one, a reverse shot. The page splits them into 2 or 3 layers, fits each line, and gives speeds, depths by the intercept-time and crossover methods, a dipping-layer solution from the two shots, reciprocal time check and a forward model to compare.

Keep it one tap away for everyday use.

Bookmark Seismic refraction interpreter?

Press Ctrl + D now. Your browser opens its own bookmark box, where you choose the name and folder.

  1. Picks
  2. Layers
  3. Lines
  4. Depths
Nothing is uploadedWorks offline once open

First-arrival picks

Automatic breaks try every split and keep the best straight-line fit in which speeds increase with depth. Check them on the plot and type your own if a point is on the wrong line. The example picks are made up.

Layers

Download

Forward model (flat layers)

What it does

In a refraction survey you set off a small source (a hammer on a plate, a weight drop or a charge) and record the first arrival at a line of geophones. Near the shot the first wave travels straight through the top layer. Further out, waves that run along the top of a faster layer underneath overtake it. Plotted as time against distance, the picks fall on straight lines, one per layer, and the slopes and intercepts give the speed and depth of each layer. This page does the whole interpretation:

  • Layer breaks. It splits the picks into 2 or 3 straight lines by trying every split and keeping the best fit in which speeds increase with depth. You can type your own breaks instead.
  • Speeds from the slopes (speed = 1 ÷ slope), with the direct-wave line forced through zero if you wish.
  • Thicknesses by the intercept-time method, for flat layers: ti = Σ 2h cos θ / V, solved layer by layer.
  • Thicknesses by the crossover method, which uses only the speeds and the offsets where the lines cross: for two layers h = (xc/2)√((V2 − V1)/(V2 + V1)).
  • Dipping layer. With a reverse shot at the far end, the slower apparent speed is shot down-dip and the faster one up-dip. The page gives the true refractor speed, the dip, and the depth under each shot (vertical and at right angles to the layer).
  • Reciprocal times. The time from one shot to the other should be the same both ways, a good check on picks and on the shot positions.
  • Forward model. Type speeds and thicknesses of flat layers to get intercepts, critical distances, crossover offsets, the travel-time curve on the plot and its RMS difference from your forward picks.

How to use it

  1. Choose 2 or 3 layers and enter the distance between the shots.
  2. Paste or open the forward picks: offset from the shot in metres and time in milliseconds. Add the reverse shot’s picks, measured from the reverse shot, if you have them.
  3. Check the plot: each pick should sit on its own line. If one is on the wrong line, type the break offsets (the distances between layers), for example 12, 55.
  4. Read the speeds, depths and the dipping-layer result. Download the picks with their layers and fitted times as CSV.

Worked example

The example picks were made from a three-layer model (500, 1800 and 4000 m/s) whose interfaces lie 5 m and 20 m below the forward shot and dip 2° and 3° toward the reverse shot 125 m away, with small timing errors added. The forward shot gives 492, 1631 and 3090 m/s and layer tops at 5.0 m and 19.3 m; the reverse shot gives 505, 1997 and 6042 m/s. The different apparent speeds show the layers dip. The dipping-layer solution gives V2 = 1795 m/s and a dip of 1.67°, with the first refractor 5.06 m below the forward shot and 9.19 m below the reverse shot (the model has 5.0 and 9.4 m). The deeper refractor’s speed is about 4089 m/s, and the reciprocal times agree to 0.1 ms.

Good to know

  • The intercept and crossover methods give the same depth when the direct-wave line goes through zero. Untick that box if your trigger time is uncertain; the two methods then differ, which shows how sensitive the depth is to the picks.
  • A layer that is slower than the one above it makes no head wave and cannot be seen (a hidden layer); depths below it come out wrong.
  • A thin layer may never give first arrivals (a blind zone). Surface checks or boreholes help.
  • The dipping solution here is for the first refractor. For a deeper dipping refractor, the average of the two apparent speeds (2·Vf·Vr/(Vf + Vr)) is a fair estimate of its speed for small dips.

FAQ

How many picks do I need per layer? At least two, but four or more give a much steadier line.

What if the ground is not layered? Strongly varying ground needs methods such as the generalised reciprocal method or tomography. This page is for layered ground.

Is my data sent anywhere? No. Everything is worked out on this page and nothing leaves your device.

Where the numbers come from

  • Straight-line travel times of direct and head waves, the intercept-time and crossover-distance depth formulas and the two-shot dipping-layer solution are standard refraction seismology (for example Telford, Geldart and Sheriff, Applied Geophysics; Kearey, Brooks and Hill, An Introduction to Geophysical Exploration), written by fibub from the formulas.
  • Checked against first-arrival times found independently by Fermat's principle (SciPy minimisation over the ray path) for flat and dipping layers.
  • The example picks are made up by fibub from a dipping three-layer model with small errors added.

Standard scientific values, such as constants, atomic weights and melting points, are facts found in every textbook; we name where ours were checked. If you are an author or publisher and want a value, a credit or a table changed or removed, tell us and we will act quickly.

Last updated 2026-11-09.

Related tools

Level book reducer #133 Paste or open a levelling field book (back, intermediate and fore sights). Get reduced levels by both rise-and-fall and height of collimation, every arithmetic check, the closing error against a benchmark or loop, the allowed misclosure from your own k, adjusted levels and a profile. Publication plot maker #27 Paste columns or open a CSV or Excel file and get a clean scientific figure with error bars, sized for a journal column. Download vector SVG or PDF, or PNG at 300 or 600 dpi. Your data stays in your browser. GPR radargram viewer #513 Open ground-penetrating radar profiles in SEG-Y (rev 0 or 1, IBM or IEEE floats, integers, either byte order) or GSSI DZT format in your browser. Set time zero, dewow, remove the background, add linear, exponential or AGC gain, show depth from a wave speed, fit a hyperbola for the speed, view single traces and save the radargram as PNG. Earthquake locator and magnitude calculator #878 Paste P and S arrival times from three or more seismic stations. Get each station's S−P distance, the epicentre and origin time by least squares (Geiger's method, on a flat map or the globe), optional depth, error estimates, a map with distance circles, local magnitude ML from amplitudes, and Mw from seismic moment or fault size. Seismogram viewer #879 Open SAC and miniSEED seismogram files (Steim-1, Steim-2, integers and floats) in your browser. See several traces at once with zoom and pan, remove the trend, taper, apply a Butterworth bandpass (one pass or zero phase), see the amplitude spectrum and a spectrogram, tap to pick P and S, and download CSV. Gravity survey corrections #880 Paste a day of gravity meter readings with base station repeats. Get drift-corrected observed gravity, normal gravity from GRS80 or WGS84, free-air and Bouguer slab corrections at your density, ring terrain corrections, and free-air, simple and complete Bouguer anomalies as a table, a profile plot and a CSV. Resistivity sounding (VES) inversion #883 Paste Schlumberger or Wenner sounding readings (AB/2, MN/2, ΔV and I, or apparent resistivity). Get apparent resistivities, a layered-earth model of 2 to 6 layers by damped least-squares inversion, the log-log sounding curve, a resistivity-depth plot, error factors and a CSV.