Zoeppritz Equations¶

The exact plane-wave interface relations that use elastic boundary continuity to partition an incident P or S wave among reflected and transmitted P and S modes as a function of angle and material properties.

Version: v1 · 2026-08-30 · History
Domain-specific #: 3144
Origin domain: geophysics
Subdomain: reflection seismology
Aliases: Zoeppritz relations, Zoeppritz system

Core Idea¶

The Zoeppritz Equations are the exact linear relations governing reflection, transmission, and mode conversion of plane elastic waves at a planar boundary between two homogeneous isotropic elastic media. Given the incident wave type and angle plus each medium's density, compressional-wave velocity, and shear-wave velocity, the equations determine amplitude coefficients for the permitted outgoing P and S waves. An incident P wave, for example, can generate reflected P, reflected SV, transmitted P, and transmitted SV components.

The equations arise by requiring continuity of displacement and traction across a welded interface while all participating waves share compatible horizontal slowness under Snell's law. These boundary conditions provide four linear equations for four unknown amplitude coefficients in the two-dimensional P–SV case. Solving the system partitions the incident disturbance among modes; properly converting amplitude coefficients to energy flux provides the corresponding energy partition.^[1]

The locked identity is two elastic half-spaces + plane interface + specified densities and P/S velocities + one incident plane P or SV wave + common ray parameter + displacement and stress continuity -> angle-dependent reflected and transmitted P/S coefficients. The system is not simply “a reflection equation.” Its distinct contribution is the coupled exact solution including mode conversion and oblique incidence.

Structural Signature¶

upper and lower media — each is represented as homogeneous, isotropic, linearly elastic, and characterized by density, P-wave velocity, and S-wave velocity;
a planar welded interface — the media meet across a boundary at which the relevant displacement and traction components are continuous;
incident mode — a P or vertically polarized S wave approaches from a declared medium;
incidence angle — measured relative to the interface normal under an explicit convention;
common horizontal slowness — Snell's law relates incident, reflected, and transmitted angles;
outgoing modes — reflected and transmitted P and SV waves, subject to propagation or evanescence conditions;
four boundary equations — tangential and normal displacement plus tangential and normal stress conditions;
four amplitude unknowns — coefficients relative to a defined incident amplitude and wavefield convention;
matrix solution — the coefficient vector is obtained by solving the coupled linear system;
critical-angle behavior — transmitted modes can become evanescent when the angle exceeds the relevant condition;
amplitude-versus-angle relation — coefficients change with incidence geometry and elastic contrast;
energy check — lossless-media solutions should conserve normal energy flux after impedance and angle factors are applied.

Conventions matter. Potentials, displacement amplitudes, particle velocity, stress, and energy-normalized coefficients are not interchangeable. A correct matrix copied from one source can appear different from another because of polarity, angle, or normalization conventions. The physics must be compared after those conventions are reconciled.

What It Is Not¶

Not Snell's law alone. Snell's law gives ray-angle compatibility, not amplitudes.
Not the normal-incidence acoustic coefficient. The simple impedance-contrast formula suppresses obliquity, shear waves, and conversion.
Not a generic wave-interface condition. Zoeppritz fixes isotropic elastic P–SV propagation and a particular boundary-value problem.
Not an approximation. Aki–Richards, Bortfeld, and Shuey forms trade exactness for interpretability under stated assumptions.
Not amplitude-versus-offset by itself. AVO is an observational and inversion practice that uses angle-dependent reflection behavior.
Not energy coefficients without conversion. Squaring a displacement-amplitude coefficient alone generally does not yield energy fraction.
Not directly valid for arbitrary anisotropic, attenuating, rough, curved, or layered interfaces. Those cases require extensions or different models.
Not restricted to incident P waves. Reciprocal S-wave incidence has its own coefficient set.
Not proof of reservoir fluid content. Elastic contrasts inferred from seismic amplitudes remain nonunique and processing-dependent.

Scope of Application¶

Reflection seismology uses the equations to model how seismic amplitudes vary with incidence angle and subsurface elastic properties. Acquisition offset serves as a practical proxy for angle after velocity modeling, making the relations foundational to amplitude-versus-offset or amplitude-versus-angle analysis. Changes in P velocity, S velocity, and density influence intercept, gradient, and curvature, which can support lithology and fluid interpretation when combined with geological constraints.

The exact equations also provide benchmarks for numerical wave solvers, synthetic seismograms, interface modeling, and approximation error. At modest contrasts and angles, simplified expressions can expose interpretable combinations of elastic parameters. Shuey's widely used approximation organizes the P–P reflection coefficient into normal-incidence, gradient, and far-angle terms, with a two-term form often used at moderate angles.^[2] Approximation choice must be matched to angle, contrast, anisotropy, and target accuracy.

The same boundary-value framework is relevant to earthquake seismology and elastic ultrasonics, but the node remains an elastodynamic equation system rather than a general prime. Application to fluids requires recognizing that shear velocity vanishes and changes the permitted modes.

Clarity¶

For incident P energy, the amplitude coefficient RPP describes reflected P amplitude relative to incident amplitude under the chosen convention; RPS describes converted reflected S amplitude; TPP and TPS describe the transmitted modes. Because the modes travel at different speeds and angles in media of different density, energy fractions include impedance and direction-cosine factors. A coefficient magnitude above one can occur under some amplitude normalizations without violating energy conservation.

The full matrix is exact only within its modeling assumptions. “Exact Zoeppritz” does not mean exact Earth. Thin beds, gradients, roughness, anisotropy, attenuation, poroelasticity, and finite-frequency effects can make the half-space interface model incomplete. The word exact distinguishes solution of the ideal boundary conditions from approximations to that solution.

Manages Complexity¶

The equations compress a coupled boundary problem into a reusable solver. Instead of reasoning separately about every reflected and refracted mode, the matrix enforces all interface obligations simultaneously. This prevents physically inconsistent coefficient choices and makes sensitivity to angle and material contrast calculable.

They also establish a reference hierarchy. The exact system anchors approximations; approximations expose parameter effects; AVO workflows connect those effects to data. Without this hierarchy, a convenient linear formula can be used outside its range and mistaken for the underlying physics.

Abstract Reasoning¶

If the two media have identical elastic properties, the interface produces no reflected energy in the ideal model.
At normal P incidence, converted SV amplitude vanishes by symmetry and the result reduces toward the acoustic-impedance contrast relation.
At oblique incidence, P–S conversion is generally permitted because boundary conditions couple normal and tangential motion.
If a transmitted angle becomes complex beyond a critical angle, that mode is evanescent rather than a propagating far-field ray.
If a simplified equation neglects its far-angle term, errors should grow as incidence approaches critical regimes.
If amplitudes are compared across processing steps, polarity and normalization conventions must remain consistent.
If coefficient squares do not sum to one, first convert amplitudes to normal energy flux before diagnosing nonconservation.
If anisotropy is material, isotropic inversion can assign directional effects incorrectly to P/S velocity or density contrast.
If a thin bed contains two close interfaces, interference prevents treating the observed amplitude as one isolated Zoeppritz boundary.
If density and velocity contrasts trade off in the data, a unique geological interpretation requires external constraints.

Knowledge Transfer¶

The portable skeleton is incident mode + rule-governed interface + continuity constraints + coupled outgoing modes -> partition coefficients. It informs acoustics, optics, and transmission-line reasoning, but those fields have their own boundary variables and coefficient equations. The Zoeppritz name should remain attached to elastic seismic P/S relations.

The broader modeling lesson is that boundary constraints determine globally coupled outputs. Solving one reflected mode independently can violate conditions that only the complete set satisfies.

Examples¶

P-wave on a sediment boundary: the model returns P–P reflection, P–S conversion, and two transmitted modes;
normal incidence: the system collapses to an impedance-controlled reflection case with no P–S conversion;
AVO synthetic: coefficients are evaluated over angle to predict gather amplitude trends;
Shuey comparison: a two- or three-term approximation is checked against the exact coefficients over the survey's angle range;
critical angle: one transmitted mode becomes evanescent and coefficient phase behavior changes;
fluid boundary: zero shear velocity removes propagating S behavior in the fluid and requires careful limiting treatment;
non-example—ray tracing: travel paths alone do not determine mode amplitudes;
failure—thin-bed use: one-interface coefficients are assigned to a tuning response produced by multiple reflections.

Structural Tensions¶

exact ideal solution vs. approximate Earth — mathematical exactness coexists with restrictive material assumptions;
completeness vs. intuition — the coupled matrix preserves physics while obscuring parameter effects;
amplitude vs. energy — simple coefficients are convenient but require flux normalization for conservation;
large-angle information vs. instability — far angles can add sensitivity and approach critical, anisotropic, or noise-dominated regimes;
mode conversion vs. interpretation — additional observables offer constraints and additional convention risk;
approximation speed vs. validity — simplified formulas support inversion while imposing angle and contrast limits.

Structural–Framed Character¶

Zoeppritz Equations are structural inside geophysics. The media, wave modes, geometry, and boundary conditions determine the system. Sign and normalization conventions are framed representational choices; equivalent formulations preserve the physical solution.

Structural Core vs. Domain Accent¶

The structural core is coupled partition at an interface under continuity constraints. The domain accent is isotropic elastic P/S waves, density and seismic velocities, ray angles, tractions, and seismological amplitude interpretation. Removing it yields Wave or Interface.

Wave — incident and outgoing disturbances propagate through elastic media.
Interface — a boundary couples two media under continuity rules.
Conservation — lossless solutions preserve energy flux across modes.
Constraint Satisfaction — four boundary equations determine four coefficients.
Approximation — operational forms simplify the exact system within controlled regimes.

The minimal prospective DAG uses a composition edge to prime:wave. Wave behavior is load-bearing, while the candidate adds the elastic interface problem and exact coefficient system.

Relationships to Other Abstractions¶

Current abstraction Zoeppritz Equations Domain-specific

Parents (1) — more general patterns this builds on

Zoeppritz Equations is part of Wave Prime

incident and outgoing disturbances propagate through elastic media.

Hierarchy path (1) — routes to 1 parentless root

Zoeppritz Equations → Wave

Neighborhood in Abstraction Space¶

Zoeppritz Equations sits in a sparse region of the domain-specific corpus (92^nd percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.

Family — Unclustered & Miscellaneous (1565 abstractions)

Nearest neighbors

Computed from structural-signature embeddings · 2026-09-08

Not to Be Confused With¶

Snell's law;
Fresnel equations for electromagnetic waves;
normal-incidence acoustic reflection;
Knott's alternative potential formulation;
Shuey, Aki–Richards, or Bortfeld approximations;
amplitude-versus-offset data as such;
ray-path computation without amplitudes;
arbitrary anisotropic interface equations;
reflection coefficient squared interpreted automatically as energy.

References¶

[1] Robert E. Sheriff and Lloyd P. Geldart, Exploration Seismology, 2^nd ed. (Cambridge University Press, 1995), chapters on elastic-wave reflection and transmission. registry ↩

[2] R. T. Shuey, “A Simplification of the Zoeppritz Equations,” Geophysics 50(4) (1985), 609–614, https://doi.org/10.1190/1.1441936. registry ↩

[3] Karl Zoeppritz, “Erdbebenwellen VII B: Über Reflexion und Durchgang seismischer Wellen durch Unstetigkeitsflächen,” Nachrichten von der Königlichen Gesellschaft der Wissenschaften zu Göttingen (1919), historical scan at https://www.digizeitschriften.de/dms/img/?PID=GDZPPN002505290. registry

[4] “Zoeppritz equations,” Wikipedia, frozen revision 1346296855, https://en.wikipedia.org/wiki/Zoeppritz_equations. registry