Home / Science / Earth sciences / Structural geology / Paleostress inversion (fault-slip data)Tool #869

Paleostress inversion: stress axes and stress ratio from fault-slip data, Michael method with bootstrap (free)

Find the principal stress directions σ1, σ2, σ3 and the stress ratio from faults with striations by Michael's linear least-squares method, with bootstrap confidence, the misfit of every fault, P and T axes and a stereonet.

Keep it one tap away for everyday use.

Bookmark Paleostress inversion (fault-slip data)?

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

  1. Faults
  2. Inversion
  3. Stereonet
  4. Report
Nothing is uploadedWorks offline once open

Fault-slip data

Options

Stress

Faults

Misfit: the angle between the slip you measured and the slip the stress would drive on that plane. Misfits above about 30° point to a fault that moved under a different stress, a wrong sense, or a typing error.

Data
Download

What it does

Faults with striations (slickenlines) record how rock slid. If many faults moved under one uniform stress, the directions they slid in tell you that stress. This tool finds the directions of the three principal stresses (σ1 the most compressive, σ3 the least) and the stress ratio, which says whether σ2 is closer to σ1 or to σ3. It uses Michael’s (1984) linear least-squares method, gives bootstrap confidence for every answer, the misfit of each fault, the P and T axes, and draws it all on a stereonet.

How to use it

  1. Pick how each line is written:
    • Strike, dip, rake: rake in the Aki and Richards convention, measured in the fault plane from the strike direction: positive for reverse, negative for normal, 0 for left-lateral and ±180 for right-lateral.
    • Strike, dip, striation trend and plunge, sense: sense is N (normal), R (reverse or thrust), D (dextral) or S (sinistral). A striation a few degrees off the plane is projected onto it; more than 20° off, the line is skipped.
  2. Choose strike (right-hand rule) or dip direction for the planes, then paste or type the faults, or open a CSV or Excel file.
  3. Pick the number of bootstrap resamples and the confidence. Show or hide planes, slip lines, P and T axes, and the bootstrap cloud.
  4. Read the axes, the stress ratio and the misfits, and download the report (image, PDF or Word) or the tables as CSV.

How it works

Every fault gives three linear equations: the shear traction of the stress tensor on the plane should point along the slip and have length 1. That “same shear stress on every fault” assumption is what makes the problem linear. The five unknowns of the deviatoric tensor (trace zero) are solved by least squares, and the eigenvectors of the tensor are the principal axes. The stress ratio is shown both ways round: Φ = (σ2 − σ3)/(σ1 − σ3) and R = 1 − Φ = (σ1 − σ2)/(σ1 − σ3). The bootstrap repeats the inversion on resampled sets of faults (seed 1, so the numbers do not change between visits) and reports the angle within which the chosen share of resampled axes falls.

Worked example

The example holds 24 synthetic faults made from a known stress with σ1 plunging 10° towards 040 and Φ = 0.4, plus a few degrees of noise. The inversion gives σ1 11° → 037, σ2 2° → 128, σ3 78° → 227 and Φ = 0.44 (R = 0.56): reverse or thrust faulting, close to the stress that made the data. The mean misfit is 4.3°. With 1000 resamples, 95 % of the bootstrap σ1 axes fall within 5.2° of the best fit and Φ ranges from 0.37 to 0.49.

Tips

  • Use at least 10 faults with a good spread of orientations. Faults that all strike one way leave part of the stress unknown.
  • Large misfits often mean a wrong sense of slip, a fault from another event, or a typing slip. Look at those rows first; if the data hold two events, split them and invert each set.
  • The P and T axes of single faults are a quick check that does not assume one stress for all faults.

Frequently asked questions

Is this the full Michael method? It is his linear inversion with bootstrap resampling. It does not choose between a fault plane and its auxiliary plane: it uses the planes you give.

Why is my answer different from a grid-search method? Methods weigh misfits differently; small differences within the confidence cones are normal.

Is my data uploaded? No. Everything runs in your browser and nothing leaves your device; the page works offline once loaded.

Where the numbers come from

  • Method: Michael, A. J. (1984), Determination of stress from slip data: faults and folds, Journal of Geophysical Research 89, B13, 11517-11526 (the method is described in words and formulas, no data copied).
  • Rake convention: Aki and Richards. Bootstrap: resampling faults with replacement (Efron).
  • Example: 24 synthetic faults made by fibub from a known stress tensor with a few degrees of noise.

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

Stereonet plotter #33 Plot planes and lines on an equal-area (Schmidt) or equal-angle (Wulff) stereonet. Get great circles, poles, the mean orientation with its 95 % cone and a rose diagram, and download SVG or PNG. Apparent dip and true thickness calculator #867 Work out the apparent dip of a bed in any section, the true dip from two apparent dips, true and vertical thickness from outcrop width on sloping ground, the depth to a bed in a borehole, and strike and dip from three points. Cross-section and topographic profile #509 Draw a topographic profile from a DEM (GeoTIFF or ASCII grid) along a line you tap on the map, or from typed distance and elevation pairs, with vertical exaggeration. Add contacts by strike and dip: each is drawn at its apparent dip in the section's direction and the units between them are filled. Download SVG, PNG, PDF or CSV. Phone clinometer #866 Lay your phone on a rock surface to measure strike and dip (or dip direction and dip), or along a lineation for trend and plunge. Averaged over a second, corrected by the magnetic declination you enter, logged offline in this browser with notes and an optional position, and exported as CSV for the stereonet plotter. 1D Schrödinger solver #923 Choose a textbook potential, type a formula or draw your own, and get the lowest energy levels and wavefunctions of the one-dimensional Schrödinger equation by the matrix (finite-difference) method, with the exact levels alongside for the box, harmonic oscillator, finite well, Morse, linear and Pöschl–Teller potentials. Adduct and neutral loss calculator #984 Type a formula or a neutral monoisotopic mass and get the m/z of 25 common electrospray adducts in positive and negative mode, the ions left after common neutral losses (water, ammonia, CO2, sugars, phosphate and your own), and match a list of observed m/z values within a ppm tolerance. With no molecule, it turns observed m/z values into candidate neutral masses.