Carcione and Ursin (2016) proved that Fermat's variational principle and Snell's law are equivalent for ray propagation in homogeneous lossy (viscoelastic) anisotropic media when rays are defined through the stationary complex slowness vector, but they explicitly left open whether this equivalence extends to the heterogeneous case — source and receiver lying in different media, so that the ray is refracted at the interface (or interfaces) between homogeneous lossy anisotropic media — where the three velocity definitions used to construct rays (ray, envelope, and energy velocities) are known to disagree. This problem asks for a rigorous extension — or refutation — of the Fermat–Snell equivalence to such heterogeneous lossy anisotropic media, together with an explanation of the discrepancies the source reports among the competing ray-construction approaches regarding the refraction point and the value of the minimum traveltime. Recent work (2024–2026) has made ray tracing in such media numerically feasible and has quantified interface discrepancies, but the formal equivalence question remains unsettled.

Background

A lossy or anelastic medium is one in which wave energy dissipates during propagation; in geophysics such media are described by viscoelasticity, and the dissipation is commonly summarized by the quality factor QQ (the number of cycles a wave needs for its energy to decay by a factor e2πe^{2\pi}; small QQ means strong attenuation). In anisotropic media the wave speed depends on propagation direction. Combining loss and anisotropy has a subtle consequence: the slowness vector (the vector reciprocal to phase velocity, pointing normal to the wavefront) becomes complex, and its real and imaginary parts need not be parallel. The angle γ\gamma between them is called the inhomogeneity angle; γ=0\gamma=0 defines a homogeneous wave, the direct generalization of an elastic plane wave. Carcione and Ursin (Geophysics 81, T107–T116, 2016) studied ray theory in this setting. Fermat's principle of least action states that a ray between two points makes the traveltime stationary; Snell's law states that the tangential component of the slowness vector is continuous across an interface (in its complex form, for both real and imaginary parts). In lossy anisotropic media, several distinct velocities appear in the literature — ray, envelope, and energy velocities, which replace the group velocity because group velocity loses physical meaning for anelastic waves — and rays can be constructed from any of them, with the rigorous choice based on the stationary complex slowness vector, i.e. the slowness that makes the (complex, frequency-dependent) traveltime stationary. Carcione and Ursin proved that in homogeneous media all these constructions coincide for stationary rays, so Fermat's principle and Snell's law are equivalent there. However, in heterogeneous media the three velocity definitions disagree, producing different refraction points and different minimum traveltimes; the authors found the differences negligible for real earth materials (Q>5Q>5, wavelike regime) but could not adjudicate them rigorously, because direct numerical or full-wave comparison is impractical: the predicted differences are extremely small for wavelike QQ, and for diffusion-like attenuation (Q≲5Q\lesssim 5) traveltime evaluation itself becomes ill-posed. Later work by Liu, Zhou, Wang and collaborators developed a g*-Hamiltonian method (a reformulation of the ray equations using a modified Hamiltonian that keeps rays real in complex media) for real ray tracing, compared real and complex versions of Snell's law at interfaces between viscoelastic anisotropic materials (Liu et al., Geophys. J. Int. 240, 13–30, 2024; Wang et al., Geophysics 90, C137–C150, 2025), extended real ray tracing to heterogeneous viscoelastic VTI models (media with a vertical symmetry axis of transverse isotropy; Liu, Geophysics 90, C67–C80, 2025), and derived theoretical conditions on the quality factors under which quasi-compressional (qP) and quasi-shear (qSV, qSH) waves in such media are homogeneous (Wang et al., Geophys. Prospect., 2026). These advances are numerical and partial: none establishes the variational foundation itself for heterogeneous media.

Problem Statement

Within the stationary-complex-slowness framework of Carcione and Ursin (2016), establish whether the equivalence between Fermat's principle (stationarity of the complex traveltime functional) and Snell's law (continuity of the tangential complex slowness across an interface), proved for homogeneous lossy anisotropic media, extends to heterogeneous lossy anisotropic media in their sense — configurations in which the source and receiver lie in different media, so that the ray is refracted at the interface (or interfaces) between homogeneous lossy anisotropic media. A complete answer should (1) prove the equivalence for such configurations, or exhibit a heterogeneous configuration in which the Fermat-stationary ray (constructed from the stationary complex slowness that makes the complex traveltime stationary) and the Snell-constructed ray (constructed from continuity of the tangential complex slowness across the interface or interfaces) genuinely differ; and (2) account, within the same equivalence question, for the discrepancies reported in the source among the ray, envelope, and energy velocity constructions in the heterogeneous setting, namely the non-unique point at which the ray crosses the interface and the inconsistent values of the minimum traveltime, explaining the origin of these discrepancies.

References

  1. José M. Carcione and Bjorn Ursin, "On Fermat's principle and Snell's law in lossy anisotropic media", Geophysics 81(3), T107–T116 (2016), DOI 10.1190/geo2015-0585.1, https://doi.org/10.1190/geo2015-0585.1
  2. José M. Carcione, Xu Liu, Stewart Greenhalgh, Marco A. B. Botelho, and Jing Ba, "Quality Factor of Inhomogeneous Plane Waves", Acoustical Physics 66(6), 598–603 (2020), DOI 10.1134/S1063771020060111, https://doi.org/10.1134/s1063771020060111
  3. Xu Liu, Bing Zhou, Abdullatif A Al-Shuhail, Wei Zhou, Bo Liu, and Qinzhuo Liao, "Ray and energy-flux velocities at a contact of two viscoelastic anisotropic materials", Geophysical Journal International 240(1), 13–30 (2024), DOI 10.1093/gji/ggae363, https://doi.org/10.1093/gji/ggae363
  4. Xu Liu, "Real ray tracing in viscoelastic anisotropic media with homogeneous ray and energy velocities", Geophysics 90(4), C67–C80 (2025), DOI 10.1190/geo2024-0706.1, https://doi.org/10.1190/geo2024-0706.1
  5. Pengde Wang, Wei Zhou, and Xu Liu, "Continuous reflection and transmission coefficient curves and ray tracing at an interface of viscoelastic anisotropic media based on complex Snell's law", Geophysics 90(6), C137–C150 (2025), DOI 10.1190/geo2024-0816.1, https://doi.org/10.1190/geo2024-0816.1
  6. Pengde Wang, Kejia Pan, Xu Liu, Jianping Liao, Chao Jin, and Hexiu Liu, "Theoretical Verification of Homogeneous Seismic Body Waves in Viscoelastic Isotropic and Transversely Isotropic Media", Geophysical Prospecting 74(7) (2026), DOI 10.1111/1365-2478.70239, https://doi.org/10.1111/1365-2478.70239
  7. José M. Carcione, "Anisotropic anelastic media", chapter in Wave Fields in Real Media (2022), DOI 10.1016/b978-0-32-398343-3.00013-1, https://doi.org/10.1016/b978-0-32-398343-3.00013-1