Zoeppritz equations: from seismology to medical exploration

Abstract: More than a century ago, Karl Bernhard Zoeppritz derived the equations that determine the reflected and transmitted coefficients at a planar interface for an incident seismic wave. The coefficients so obtained are a function of the elastic parameters of the media on each side of the interface and the angle of incidence. Approximations of the equations have been proposed and used in geophysical exploration, however, full use of the equations and their generalization to multiple layers could offer richer information about the properties of the media and be helpful in medical diagnosis via ultrasound. In this work, we investigate how to extract information from the angle-dependent reflection coefficients, including critical angles and the wave distortion at the interface between two and three media. It is shown that it is possible to separate the effect of density from speed of sound mismatch by measuring amplitudes as a function of angle of incidence (AVA). And examining the critical angle and waveform distortion of the reflected waves can reveal the thickness of an intermediate layer, even with subwavelength resolution. These studies could be integrated into medical imaging and also into the training of artificial intelligence systems that assist in diagnosis. In particular, they could help prevent cerebrovascular accidents by early detection of the formation and hardening of plaque in the arteries that irrigate the brain.

• The paper demonstrates the extension of classical Zoeppritz equations to three-layered media for precise subwavelength estimation in biological tissues.
• It applies amplitude versus angle (AVA) analysis to decouple density and velocity contrasts, revealing key insights in tissue composition.
• Waveform distortion and phase shifts at critical angles are exploited to enhance diagnostic imaging and facilitate AI-driven analysis.

Zoeppritz Equations: Bridging Seismology and Medical Ultrasound

Introduction

The Zoeppritz equations, formulated in 1919, provide a rigorous framework for calculating the amplitudes of reflected and transmitted elastic waves at planar interfaces between different media. While their primary application has been in seismology and geophysical exploration, this paper systematically explores their potential in medical diagnostics, particularly in ultrasound-based imaging. The authors argue that, with modern computational resources, the full, non-approximated Zoeppritz equations can be leveraged to extract richer information from angle-dependent reflection data, enabling subwavelength resolution in tissue characterization and layer thickness estimation.

Theoretical Framework and Methodology

The Zoeppritz equations are derived by enforcing continuity of displacement and stress at the interface between two (or more) elastic media, incorporating Snell's law for wave propagation. Unlike electromagnetic waves, elastic waves in solids can be longitudinal (P-waves) or transverse (S-waves, with both SH and SV polarizations). The equations yield a system of four (for two layers) or eight (for three layers) linear equations relating the amplitudes of incident, reflected, and transmitted waves to the elastic parameters (density, P- and S-wave velocities) and the angle of incidence.

The authors extend the classical two-layer problem to three layers, deriving the full set of equations for arbitrary densities and velocities. This generalization is critical for modeling realistic biological structures, such as vessel walls surrounded by different tissue types.

Amplitude Versus Angle (AVA) Analysis

A central result is the demonstration that amplitude versus angle (AVA) analysis enables the decoupling of density and velocity contrasts at interfaces. For normal incidence, the reflection coefficient depends solely on the acoustic impedance mismatch ($Z = \rho a$), making it impossible to distinguish between changes in density and sound speed. However, at oblique incidence, the angular dependence of the reflection coefficient is sensitive to the individual contributions of density and velocity, as well as to S-wave velocity (shear modulus) differences.

The authors show that, even when acoustic impedances are matched, variations in P- or S-wave velocities manifest as polarity reversals and amplitude changes in the reflected wave as a function of angle. This sensitivity is particularly relevant for detecting compositional or morphological changes in biological tissues, such as the hardening of arterial plaque, which increases the local sound speed without necessarily altering density.

Waveform Distortion and Critical Angle Effects

Beyond amplitude, the phase and shape of reflected wavelets are affected by the critical angle for total internal reflection. When the angle of incidence exceeds the critical angle for a given wave type, the transmitted wave becomes evanescent, and the reflected wave acquires a phase shift. For broadband pulses (wavelets), this results in observable waveform distortion, which can be exploited to infer interface properties.

The authors simulate wavelet reflection at interfaces with sharp velocity contrasts (e.g., soft tissue to bone), showing that the onset of distortion correlates with the critical angle. This effect provides an additional diagnostic observable, particularly for identifying boundaries with large impedance mismatches.

Three-Layered Media and Subwavelength Thickness Estimation

The extension to three-layered media is of particular significance for medical applications. The authors derive the full set of equations for arbitrary elastic parameters and demonstrate, via simulation, that the AVA response of the reflected wave is highly sensitive to the thickness of the intermediate layer, even when the thickness is much less than the wavelength of the probing ultrasound.

This subwavelength sensitivity arises from interference effects between multiple reflections within the intermediate layer, which modulate the amplitude and phase of the composite reflected wave as a function of angle. The authors propose that, by fitting the measured AVA curve (or its shape) to the theoretical model, it is possible to estimate layer thickness with high precision, independent of absolute amplitude calibration.

Implications for Medical Imaging and AI Integration

The practical implications are substantial. The ability to resolve subwavelength features, such as the thickness of arterial walls or the presence of thin plaques, could significantly enhance the early detection of pathologies like atherosclerosis or cancer. The authors suggest that the AVA and waveform distortion signatures can be used as input features for machine learning models, enabling automated classification and diagnosis.

Furthermore, the methodology is compatible with lock-in detection and advanced signal processing techniques to improve SNR and robustness to noise. The approach is not limited to vascular imaging but can be generalized to other layered biological structures, including bone and soft tissue interfaces.

Numerical Results and Claims

The simulations presented demonstrate that:

• The AVA response can distinguish between density and velocity contrasts, even when acoustic impedance is constant.
• S-wave velocity differences, undetectable at normal incidence, become apparent at oblique angles.
• Subwavelength changes in intermediate layer thickness (as small as 0.01 mm for a 1 MHz probe) produce measurable changes in the AVA curve.
• Waveform distortion at and beyond the critical angle provides an independent observable for interface characterization.

These results support the claim that full utilization of the Zoeppritz equations, without recourse to common approximations, enables richer and more precise extraction of tissue properties than current ultrasound imaging techniques.

Limitations and Future Directions

The primary limitation is the requirement for accurate measurement of reflected amplitudes and phases as a function of angle, which may be challenging in vivo due to attenuation, scattering, and limited acoustic windows. The approach also assumes knowledge of the elastic parameters of at least some of the layers, or the ability to constrain them via independent measurements.

Future work should focus on:

• Experimental validation in tissue-mimicking phantoms and clinical settings.
• Development of robust inversion algorithms for parameter estimation from noisy AVA data.
• Integration with AI-based diagnostic systems for automated interpretation.
• Extension to anisotropic and viscoelastic media, which are common in biological tissues.

Conclusion

This work rigorously demonstrates that the full Zoeppritz equations, when applied to angle-dependent reflection data, provide a powerful framework for non-invasive characterization of layered biological media. The ability to decouple density and velocity effects, resolve subwavelength features, and exploit waveform distortion opens new avenues for medical ultrasound imaging and diagnosis. The integration of these physical models with AI-driven analysis holds promise for improved early detection of vascular and oncological pathologies, with potential impact across a range of biomedical applications.