Fact-checked by Grok 2 weeks ago

Rayleigh wave

A Rayleigh wave is a type of that propagates along the of a homogeneous isotropic elastic solid, arising from the interference of compressional () and vertically polarized () waves that satisfy the boundary conditions at the surface. These waves were theoretically predicted in 1885 by the British physicist Lord Rayleigh in his seminal paper "On Waves Propagated along the Plane Surface of an Elastic Solid," where he derived the showing their existence as a localized with slightly less than the shear wave speed. The particle motion associated with Rayleigh waves is elliptical and retrograde (opposite to the direction of propagation) in the sagittal plane at the surface, transitioning to prograde deeper in the medium, with amplitudes decaying exponentially away from the surface. In , Rayleigh waves play a critical role as one of the two primary types of surface waves (alongside Love waves) generated during earthquakes, traveling along the with velocities typically 90-95% of the local shear wave speed, making them slower than body waves but capable of producing prolonged ground shaking. Their dispersive nature—where varies with —allows them to carry information about subsurface structure, enabling applications in earthquake estimation, crustal , and site characterization for purposes. Due to their concentration of energy near the surface and complex rolling motion, Rayleigh waves often contribute significantly to structural damage in populated areas during seismic events, amplifying ground motions compared to deeper-penetrating body waves. Beyond , Rayleigh waves have broad applications in , , and acoustics, where they are used to assess surface defects, measure elastic properties, and even in modern technologies like (SAW) devices for in . Their sensitivity to near-surface heterogeneities also makes them valuable in environmental and engineering for mapping soil stiffness and groundwater conditions.

Fundamentals

Definition and basic principles

Rayleigh waves are a type of surface wave that propagates along the free boundary of a homogeneous, isotropic half-space, characterized by particle displacements that follow retrograde elliptical paths in the vertical to the of . These waves arise from the interference of compressional () and () waves reflected at the surface, resulting in a motion where the horizontal and vertical components are out of phase by 90 degrees, with the vertical component leading. The amplitude of the particle motion decays exponentially with depth into the solid, typically penetrating to a depth of about 0.4 wavelengths below the surface, confining the wave's energy primarily near the boundary. The existence of Rayleigh waves relies on the principles of theory, which describes the deformation of solids under stress through the stress-strain relations given by , stating that stress is proportional to within the elastic limit: \sigma_{ij} = C_{ijkl} \epsilon_{kl}, where C_{ijkl} is the stiffness tensor. In isotropic media, this simplifies to two independent constants, such as \lambda and \mu, enabling the derivation of wave equations from Newton's second law applied to . These waves require a boundary condition where the normal and stresses vanish, ensuring no external forces act on the , and they propagate only in elastic solids that are homogeneous and isotropic, without or complicating the mode. Unlike body waves, such as P-waves (which involve purely compressional motion and travel through the bulk of the medium) and S-waves (which involve shear motion and are transverse), Rayleigh waves are confined to the near-surface region and do not propagate into the interior of the solid. Their is slower than that of S-waves in the same medium, typically about 90-95% of the shear wave speed, making them the slowest major seismic wave type. This surface localization results in greater potential for damage in near-surface structures compared to body waves, as the energy is not dispersed volumetrically. The particle motion for Rayleigh waves is retrograde elliptical, meaning particles move in an ellipse oriented in the (containing the and the surface ), with the sense of opposite to the of wave advancement at . This ellipticity arises from the 90-degree phase difference between the longitudinal and vertical components, producing a closed orbital path that diminishes in size with depth. The energy of the wave is predominantly carried within the solid near , decaying exponentially away from the , with negligible into an adjacent low-impedance medium like air due to the impedance mismatch at the boundary.

Historical background

The Rayleigh wave, a type of propagating along the boundary of an elastic solid, was first theoretically predicted in 1885 by John William Strutt, the 3rd Baron Rayleigh, in his seminal paper "On Waves Propagated along the Plane Surface of an Elastic Solid," published in the Proceedings of the London Mathematical Society. This work derived the existence of waves confined to the surface, decaying exponentially with depth, arising from the interference of longitudinal and shear waves at a free boundary. Rayleigh's prediction stemmed from solving the for an isotropic, homogeneous elastic half-space, highlighting their potential role in seismic phenomena. Initial experimental confirmation of Rayleigh waves occurred in the early 20th century through observations of seismic records from distant earthquakes, where signatures matched the predicted elliptical particle motion and slower velocity relative to body waves. In 1911, extended the theoretical framework for surface waves by developing the mathematical model for Love waves—horizontally polarized shear waves in layered media—further advancing understanding of interfacial wave propagation in geodynamics. Post-World War II advancements in and ultrasonics significantly popularized Rayleigh wave studies, driven by improved instrumentation for recording long-period waves and applications in non-destructive testing. By the , Rayleigh waves became integral to earthquake monitoring, enabling global analysis of crustal structure through measurements from teleseismic events. Key contributions include A. Viktorov's 1967 monograph Rayleigh and Lamb Waves: Physical Theory and Applications, which synthesized physical principles and ultrasonic applications, establishing a foundation for practical implementations. Jan D. Achenbach's influential 1973 text Wave Propagation in Elastic Solids further developed scattering theory for Rayleigh waves, influencing advancements in wave interaction with defects and boundaries. The evolution from theoretical curiosity to practical tool accelerated in the with the surge in (SAW) devices, leveraging Rayleigh waves for in filters, delay lines, and sensors due to their high-frequency generation via interdigital transducers on piezoelectric substrates. This period marked widespread adoption in telecommunications and electronics, transforming Rayleigh waves into a of acousto-electronic .

Physical Characteristics

Propagation and particle motion

Rayleigh waves exhibit particle motion characterized by retrograde elliptical orbits in the sagittal plane, which is the plane defined by the direction of wave propagation and the normal to the free surface. At the surface, the particles trace ellipses where the vertical displacement amplitude is approximately twice that of the horizontal displacement, resulting from a phase difference between the vertical and horizontal components that causes the motion to be retrograde—meaning the particle moves backward relative to the propagation direction at the crest of the ellipse. This elliptical polarization arises from the superposition of compressional (P) and shear vertical (SV) waves interfering constructively at the surface while satisfying the free surface boundary condition of zero traction. For shallow penetration near the surface, the motion can be approximated as linear in certain simplified models, though the full elliptical nature dominates. The amplitude of particle displacement in Rayleigh waves decays exponentially with depth below the surface, typically penetrating to a depth of about one to two wavelengths before becoming negligible. This depth-dependent attenuation confines the wave energy primarily to the near-surface region, distinguishing Rayleigh waves from body waves that propagate through the full volume. The elliptical orbits transition from near the surface to prograde (forward-directed) at depths greater than approximately 0.2λ (where the horizontal is zero), reflecting the changing relative phases of the contributing wave components. At free surfaces, Rayleigh waves interact through and partial mode conversion upon encountering boundaries, such as edges or interfaces, where a portion of the may convert to other wave types like head waves or body waves. In layered or inhomogeneous , such as those with velocity contrasts, the can lead to , modifications, or generation of interfacial waves, altering the standard half-space behavior. These interactions are particularly pronounced at geometric discontinuities, where coefficients depend on the incident and material properties. In comparison to Love waves, which feature purely transverse horizontal particle motion perpendicular to the propagation direction with no vertical component, Rayleigh waves' combined vertical and longitudinal motion in the produces a rolling-like disturbance more akin to ocean waves, though with the characteristic retrograde ellipticity. This distinction aids in identifying wave types in seismic records, where Rayleigh waves appear on both vertical and radial components, while Love waves are confined to the transverse.

Speed, dispersion, and attenuation

In homogeneous isotropic , Rayleigh waves are non-dispersive, propagating at a constant speed that is approximately 0.87 to 0.95 times the shear wave speed, with a typical value of about 0.92 times the shear wave speed for a of 0.25. This speed arises from the boundary conditions at the and depends on the material's moduli and , remaining independent of in such ideal conditions. In layered or inhomogeneous media, Rayleigh waves exhibit dispersion, where the varies with due to the constructive of waves reflecting between layers of differing properties. This dependence results in wave packets that spread out over distance, with higher frequencies probing shallower depths and lower frequencies penetrating deeper into the structure. Attenuation of Rayleigh waves occurs through several mechanisms, including viscoelastic damping from internal friction in the material, geometric spreading where amplitude decreases as $1/\sqrt{r} with propagation distance r due to energy distribution over a cylindrical wavefront, and scattering losses from material heterogeneities. The quality factor Q quantifies anelastic attenuation, with the amplitude decay expressed as A(x) = A_0 e^{-\pi f x / (Q c_R)}, where f is frequency, x is distance, and c_R is the phase velocity; higher Q values indicate lower intrinsic damping. Material-specific Rayleigh wave speeds vary significantly: in rocks like , speeds range from 2 to 4 km/s depending on composition and confining pressure; in metals such as , the speed is approximately 3 km/s; and in soft tissues, speeds are much lower, typically 10 to 100 m/s, reflecting the materials' lower moduli. Experimental measurement of Rayleigh wave speeds often employs laser interferometry to detect surface displacements with high precision or piezoelectric sensors to capture vibrational responses, enabling velocity profiling across frequencies and depths.

Mathematical Formulation

Derivation of the Rayleigh wave equation

The derivation of the Rayleigh wave equation for propagation along the free surface of a homogeneous, isotropic elastic half-space (z ≥ 0) starts from the equations of linear elastodynamics. The displacement vector \mathbf{u}(\mathbf{x}, t) satisfies Navier's equation: \rho \frac{\partial^2 \mathbf{u}}{\partial t^2} = (\lambda + 2\mu) \nabla (\nabla \cdot \mathbf{u}) - \mu \nabla \times (\nabla \times \mathbf{u}), where \rho is the mass density and \lambda, \mu are the Lamé constants. This equation governs the motion in the absence of body forces. To solve for surface waves, the displacement is decomposed into irrotational (P-wave) and equivoluminal (S-wave) parts using Helmholtz potentials: \mathbf{u} = \nabla \phi + \nabla \times \boldsymbol{\psi}, where \nabla \cdot \boldsymbol{\psi} = 0. Substituting into Navier's equation yields decoupled wave equations for the scalar potential \phi and the vector potential \boldsymbol{\psi}: \nabla^2 \phi = \frac{1}{c_P^2} \frac{\partial^2 \phi}{\partial t^2}, \quad \nabla^2 \boldsymbol{\psi} = \frac{1}{c_S^2} \frac{\partial^2 \boldsymbol{\psi}}{\partial t^2}, with P-wave speed c_P = \sqrt{(\lambda + 2\mu)/\rho} and S-wave speed c_S = \sqrt{\mu/\rho}. For two-dimensional propagation in the x-z plane (no y-dependence), \boldsymbol{\psi} = (0, \psi(x,z,t), 0), and the relevant displacement components are u_x = \frac{\partial \phi}{\partial x} - \frac{\partial \psi}{\partial z}, u_z = \frac{\partial \phi}{\partial z} + \frac{\partial \psi}{\partial x}. Assume plane-wave solutions propagating in the +x direction with wavenumber k and angular frequency ω: \phi(x,z,t) = f(z) e^{i(kx - \omega t)}, \psi(x,z,t) = g(z) e^{i(kx - \omega t)}. The z-dependent functions satisfy ordinary differential equations f'' + (\omega^2/c_P^2 - k^2) f = 0 and g'' + (\omega^2/c_S^2 - k^2) g = 0. For evanescent waves decaying with depth (z > 0), select the bounded solutions f(z) = A e^{-\alpha z}, g(z) = B e^{-\beta z}, where \alpha = \sqrt{k^2 - \frac{\omega^2}{c_P^2}}, \quad \beta = \sqrt{k^2 - \frac{\omega^2}{c_S^2}}, assuming k > ω/c_P > ω/c_S for phase velocity (real α, β > 0). The free-surface conditions at z = 0 require vanishing and stresses: \sigma_{zz} = 0 and \sigma_{zx} = 0, where \sigma_{zz} = \lambda (\nabla \cdot \mathbf{u}) + 2\mu \frac{\partial u_z}{\partial z}, \quad \sigma_{zx} = \mu \left( \frac{\partial u_x}{\partial z} + \frac{\partial u_z}{\partial x} \right). At z=0, \nabla \cdot \mathbf{u} = (\alpha^2 - k^2) \phi = -\frac{\omega^2}{c_P^2} \phi, \frac{\partial u_z}{\partial z} = \alpha^2 \phi - i k \beta \psi, \frac{\partial u_x}{\partial x} = -k^2 \phi + i k \beta \psi. Thus, \sigma_{zz} = \lambda (-\frac{\omega^2}{c_P^2} \phi) + 2\mu (\alpha^2 \phi - i k \beta \psi) = 0. For the shear stress, \frac{\partial u_x}{\partial z} + \frac{\partial u_z}{\partial x} = -2 i k \alpha \phi - (k^2 + \beta^2) \psi, so \sigma_{zx} = \mu [-2 i k \alpha \phi - (k^2 + \beta^2) \psi] = 0. These form a homogeneous system for A and B (with ϕ = A, ψ = B at z=0). For non-trivial solutions, the determinant must vanish, leading to the Rayleigh secular equation: (2k^2 - \kappa^2)^2 - 4 k^2 \alpha \beta = 0, where \kappa^2 = \omega^2 / c_S^2, and k^2 + \beta^2 = 2k^2 - \kappa^2. This is a biquadratic equation in k^2, with the physical root giving real k > 0 and c = ω/k slightly less than c_S (between c_S and 0, but typically ~0.92 c_S for ν ≈ 0.25). The exact solution requires numerical root-finding.

Velocity expressions and approximations

The of Rayleigh waves in an isotropic, homogeneous half-space is obtained by solving the Rayleigh equation, with a common approximation given by \frac{c_R}{c_S} \approx \frac{0.862 + 1.14 \nu}{1 + \nu}, where c_R is the Rayleigh wave phase velocity, c_S is the shear wave velocity, and \nu is Poisson's ratio. This formula arises from fitting the numerical solution of the boundary conditions at the free surface and provides the non-dispersive speed for undamped propagation in such media. For practical computations, this approximation is used, often yielding c_R \approx 0.92 \, c_S for typical Poisson's ratios around 0.25, with errors less than 0.5%. In anisotropic media, perturbation methods address slight deviations from isotropy by expanding the velocity around the isotropic solution, incorporating elastic constants to first order for weakly anisotropic crystals. In non-homogeneous cases, such as layered structures, the v_g = \frac{d\omega}{dk} is derived from the \omega(k), where \omega is and k is , capturing energy transport along varying velocity profiles. For complex media like layered or heterogeneous solids, numerical methods such as finite element analysis solve the wave equations to compute velocity profiles, enabling accurate predictions of c_R and without analytical closed forms.

Applications in Geophysics

Generation during earthquakes

Rayleigh waves are primarily generated during earthquakes through the dynamic process of fault rupture, where the abrupt release of accumulated at the excites coupled compressional and shear modes that evolve into surface-trapped waves as they reach the . This excitation occurs as the rupture propagates along the fault plane, converting the localized slip into propagating disturbances that decay exponentially with depth, characteristic of Rayleigh wave motion. In shallow crustal events, the proximity of the source to the surface enhances this mode conversion, making Rayleigh waves a dominant component of the far-field seismic signal. Seismic source models, particularly the double-couple representation of fault slip, describe the radiation patterns that preferentially excite at larger epicentral distances greater than 100 km, where the azimuthal dependence of the source leads to nodal lines and lobes that align with propagation paths. This pattern arises from the orientation of the fault and slip , with shallow sources producing stronger amplitudes in directions to the fault strike compared to body waves, which attenuate more rapidly over distance. For instance, teleseismic recordings often show dominating the seismograms beyond regional scales due to this favorable . Amplification of Rayleigh wave amplitudes occurs through site effects involving near-surface sediments and , where low-velocity layers trap and resonate the waves, increasing ground motion intensities. Sedimentary basins can amplify Rayleigh waves 2 to 3 times greater than incident waves at frequencies around 0.3–0.7 Hz, as observed in basin models during events like the 2010 El Mayor-Cucapah . Topographic features, such as ridges, further enhance horizontal displacements at crests with amplification ratios reaching 3.8 at 1 Hz for specific structures like Robinwood Ridge, , while valleys may deamplify or alter particle motion. These effects arise from wave and at irregular boundaries, significantly influencing local shaking. Historical earthquakes illustrate these generation mechanisms vividly. The (M_w 7.9) produced prominent phases recorded at 96 global observatories, with high-quality data from 12 sites revealing extensive propagation and aiding early estimates. Similarly, the 2011 Tohoku earthquake (M_w 9.0) excited intense through its megathrust rupture, with back-projection analysis of regional records showing contributing to coastal ground motions that exacerbated impacts. In shallow events like these, play a key role in long-distance hazard propagation.

Analysis in seismology

Rayleigh waves play a crucial role in seismological analysis for inferring Earth's subsurface structure and characterizing seismic events, leveraging their sensitivity to shear-wave velocities and attenuation properties. In surface wave tomography, dispersion curves of Rayleigh waves—obtained from phase or group velocity measurements across a range of frequencies—are inverted to derive shear-wave velocity (Vs) profiles, particularly in the crust and upper mantle. This inversion process typically involves nonlinear least-squares optimization or Bayesian methods to match observed dispersion data with synthetic curves generated from layered Earth models, resolving lateral variations in Vs down to depths of several kilometers. For instance, Vs30 maps, which represent the average shear-wave velocity in the top 30 meters, are routinely derived from Rayleigh wave inversions to assess seismic site amplification for hazard mapping, as demonstrated in studies using active and passive seismic arrays in regions like the Los Angeles Basin. Array processing techniques enhance the measurement of Rayleigh wave phase velocities by exploiting spatial coherence across seismometer arrays. Beamforming, a delay-and-sum method, scans possible slownesses (apparent velocities and back-azimuths) to maximize signal power, effectively isolating plane-wave arrivals and suppressing noise. Slowness analysis complements this by estimating the horizontal slowness vector from array beam patterns, allowing precise tracking of Rayleigh wave directions and velocities, often at periods of 5–50 seconds for crustal studies. These methods are particularly effective for regional arrays spanning 10–100 km, where they resolve multimode Rayleigh waves and provide input for tomographic inversions, as applied in global and continental-scale networks since the . Seismic using ambient cross-correlations recorded at pairs of stations retrieves empirical Green's functions approximating Rayleigh wave propagation between them. This passive method extracts inter-station Rayleigh wave signals from the causal and acausal parts of the , enabling measurements without active sources. The retrieved Green's functions are then inverted for crustal velocity structures, facilitating high-resolution imaging of sedimentary basins and tectonic features, such as in the where multimode waves from noise correlations have mapped Vs variations to depths of 50 km. This approach has revolutionized crustal imaging by providing dense path coverage from continuous noise fields dominated by ocean microseisms. Attenuation tomography employs Rayleigh wave amplitudes to map the quality factor Q, which quantifies anelastic energy loss in the mantle due to temperature, composition, and partial melting. By measuring amplitude decay along propagation paths and correcting for geometric spreading and focusing effects, Q models are constructed through tomographic inversion, revealing low-Q zones associated with asthenospheric weakening. For example, global Qs models derived from long-period Rayleigh waves (60–400 s) highlight azimuthal variations in upper mantle anelasticity, with Q values as low as 80–100 in tectonically active regions indicating enhanced dissipation. These models, often combined with velocity tomography, provide insights into mantle dynamics and thermal structure. In event characterization, the Global Moment Tensor (CMT) catalog incorporates wave phases in moment tensor inversions to determine source parameters, including location, magnitude, and fault orientation. Since the , the CMT method has used long-period (100–300 s) surface waves, including , in a least-squares inversion of transverse and vertical components against synthetic seismograms in a spherically symmetric model, achieving robust solutions for events above Mw 5.5. Case studies from the catalog demonstrate how wave data constrain the centroid time and depth, particularly for shallow events, improving resolution over body-wave-only inversions in heterogeneous media.

Engineering and Technological Applications

Non-destructive testing

Rayleigh waves are widely employed in non-destructive testing (NDT) for inspecting material integrity, particularly for detecting surface and near-surface defects such as cracks and voids in industrial components. These are generated using ultrasonic methods, including transducers that couple longitudinal or waves at critical angles to produce modes, or non-contact excitation techniques that induce thermoelastic expansion for wave propagation. When these waves encounter flaws, they enable flaw detection by analyzing scattered signals, with -generated waves particularly effective for remote inspection without surface preparation. Defect sizing in Rayleigh wave NDT relies on measuring reflection coefficients and mode conversions at discontinuities, where incident partially reflect and convert to other upon interacting with cracks or voids. This approach provides to surface-breaking flaws extending up to several wavelengths deep, allowing of defect through time-of-flight and variations in the backscattered signals. In practice, the is typically limited to 1-2 wavelengths for optimal resolution, but advanced can extend detection to deeper subsurface features. In applications to pipelines and welds, systems utilizing waves facilitate mapping and weld integrity assessment by steering beams to cover large areas efficiently. These systems, often employing electromagnetic acoustic transducers (EMATs) compliant with standards like ASTM E1774, generate and detect waves without couplant, enabling rapid scanning for defects in challenging environments such as coated pipelines. For instance, inspections have been used to map wall loss due to in oil and gas pipelines, providing quantitative thickness profiles. Compared to waves, Rayleigh waves offer superior near-surface resolution due to their confinement to the material surface, typically within one , which enhances detection of shallow flaws that modes might overlook. Additionally, portable EMAT-based setups for Rayleigh wave testing allow in-situ inspections of rails and components without disassembly, reducing downtime in field applications. Post-2010 advancements include the use of nonlinear for early detection, where material degradation induces higher harmonics in the wave spectrum due to microstructural changes like microcracks. These nonlinear effects, generated via quasi-static or cyclic loading, enable sensitive monitoring of early-stage damage in metals such as , with harmonic amplitudes correlating to accumulated plastic strain. Air-coupled detection techniques further enhance this method by eliminating contact requirements, improving applicability in harsh industrial settings.

Acoustic wave devices

Surface acoustic wave (SAW) devices leverage Rayleigh waves for signal processing and sensing by confining acoustic energy to the surface of a substrate, enabling compact and efficient transduction. These devices typically employ interdigital transducers (IDTs) fabricated on piezoelectric materials such as lithium niobate (LiNbO₃) to generate and detect Rayleigh modes through the converse and direct piezoelectric effects. The IDT consists of interleaved metallic electrodes patterned via photolithography, where the electrode pitch p determines the wavelength \lambda = 2p, allowing precise control of the operating frequency f = v_R / \lambda, with v_R being the Rayleigh wave velocity on the substrate. This configuration was pioneered by White and Voltmer in 1965, who demonstrated direct piezoelectric coupling to surface elastic waves using IDTs on quartz, revolutionizing acoustic wave technology. In , Rayleigh SAW devices serve as bandpass filters and resonators, particularly in mobile phones for , , and applications since the 1990s, where they provide sharp selectivity and low , typically below 3 for preselect filters. For instance, SAW delay lines using LiNbO₃ substrates enable signal buffering with minimal , while resonators exploit reflections from structures to achieve high factors ( > 1000) for stabilization. These devices outperform early discrete implementations by integrating multiple functions on a single chip, supporting bandwidths from hundreds of MHz to several GHz. As of 2025, advancements such as laterally leaking surface acoustic waves (LLSAW) on engineered lithium niobate-on-silicon carbide (LN-on-SiC) substrates have improved performance for front-end modules, offering higher frequencies and better power handling. However, compared to bulk acoustic wave () devices, SAW filters exhibit higher at frequencies above 2 GHz due to increased in the surface-confined , though they excel in surface and ease of fabrication for lower frequencies. Rayleigh SAW sensors detect environmental changes via mass loading, where adsorbed molecules alter the wave velocity or attenuation, enabling chemical and biosensing applications. Rayleigh SAW sensors on ST-cut substrates minimize with a first-order of frequency (TCF) near zero at , ensuring stable operation without active compensation. For example, these sensors achieve sub-ng/cm² mass resolution for vapor detection by coating the propagation path with selective polymers, with enhanced by the wave's confining interactions to the surface. Fabrication involves standard for IDT patterning on 0.5-1 mm thick substrates, followed by thin-film deposition for sensing layers, yielding devices with lifetimes exceeding 10 years in controlled environments. Limitations include vulnerability to viscous in media, often addressed by designs, but overall, Rayleigh SAW sensors provide higher than bulk modes for surface-specific analytes.

Biological and Sensory Interactions

Animal detection and reactions

Animals, particularly certain mammals, exhibit remarkable sensitivity to low-frequency Rayleigh waves, which propagate as surface vibrations in the ground. detect these seismic signals through specialized mechanoreceptors in their feet and trunks, allowing them to perceive vibrations in the 10-40 Hz range that are transmitted efficiently through substrates. Similarly, such as blind mole-rats sense seismic waves (typically 250-300 Hz) via mechanoreceptors in their paws, enabling spatial orientation and environmental monitoring in subterranean habitats. Behavioral responses to Rayleigh waves often manifest as agitation or evasion, especially in seismic events like earthquakes or tsunamis. Prior to the Indian Ocean tsunami, zoo animals in , including and other species, evacuated their enclosures and sought higher ground several hours in advance, suggesting detection of precursor ground vibrations. These reactions highlight how Rayleigh waves, traveling at low speeds through soils, can trigger instinctive flight behaviors in mammals attuned to substrate disturbances. Evolutionary adaptations in burrowing species underscore the role of Rayleigh waves in survival strategies. Golden moles, for instance, utilize substrate-borne vibrations to locate prey such as , detecting seismic cues generated by insect movements through their hypertrophied structures specialized for low-frequency conduction. This adaptation allows mammals to exploit substrate-borne vibrations for . Laboratory experiments provide of vibration impacts on behavior. In studies using vibration platforms to simulate seismic conditions, mice displayed avoidance behaviors, such as freezing or fleeing, when exposed to accelerations exceeding 0.01 g (approximately 0.1 m/s²) at frequencies around 70-100 Hz, indicating a for responses in . These controlled setups demonstrate how low-amplitude s elicit measurable ethological changes, reinforcing the sensory acuity of mammalian foot pads. Despite these observations, research gaps persist regarding Rayleigh waves' specific role in earthquake prediction. Studies from the 2010s, including meta-analyses of over 180 reports, conclude that anomalous animal behaviors are inconclusive for forecasting and often attributable to faster-propagating P-wave precursors rather than slower Rayleigh surface waves alone. The U.S. Geological Survey emphasizes that while animals may detect P-waves seconds before stronger shaking, long-term predictive links to Rayleigh waves lack robust verification.

Human sensory perception

Humans detect Rayleigh waves through tactile and vestibular mechanisms, with the playing a primary role in sensing low-frequency ground motions typical of seismic events. The vestibular apparatus in the is highly sensitive to whole-body vibrations at frequencies below 20 Hz, including those associated with Rayleigh waves, with detection thresholds as low as 0.0003 g at around 100 Hz, though sensitivity extends to infrasonic ranges relevant to earthquakes. This sensitivity arises from otolith organs responding to linear accelerations and to angular motions, allowing perception of the elliptical particle trajectories characteristic of Rayleigh waves. Tactile detection via skin mechanoreceptors, such as Pacinian corpuscles, contributes when higher-frequency components (80–450 Hz) are present, with vibrotactile thresholds minimized at 200–300 Hz; however, these receptors are less dominant for the predominant low-frequency seismic Rayleigh waves (0.01–10 Hz). Auditory perception of occurs indirectly through ground-to-air coupling, where surface vibrations generate (<20 Hz) that can be felt as changes or low rumbles during strong earthquakes, though typically below the audible unless amplitudes are sufficient to produce frequencies up to 70 Hz. For instance, the Earth's surface acts as a , converting seismic energy into atmospheric waves that propagate as , sometimes preceding perceptible shaking. Overall perception limits for vertical accelerations from are around 0.001 g, below which sensations are imperceptible; microseisms—ocean-generated at 0.05–0.5 Hz—occasionally exceed this in coastal areas, manifesting as a subtle or felt through floors or structures. The elliptical particle motion of Rayleigh waves can induce psychological effects, including and disorientation, particularly in buildings where amplified swaying creates sensory conflicts between , visual, and proprioceptive inputs. Studies following major earthquakes, such as the 2016 Kumamoto event, link these low-frequency vibrations (0.1–3.5 Hz) to post-event and , resembling symptoms. Human-subject experiments using shake tables to simulate Rayleigh-like motions have, since the 1990s, demonstrated responses, including altered vestibulo-ocular reflexes and balance disruptions at accelerations as low as 0.005 m/s² (≈0.0005 g), informing thresholds for seismic and impacts.

References

  1. [1]
    [PDF] Seismic Waves and Earthquake location
    Rayleigh (1885) demonstrated that the surface boundary condition can be satisfied leading the existence of a 'coupled' and 'trapped' P-SV wave travelling along ...
  2. [2]
    On Waves Propagated along the Plane Surface of an Elastic Solid
    On Waves Propagated along the Plane Surface of an Elastic Solid. Lord Rayleigh DCL, FRS, Lord Rayleigh DCL, FRS. Search for more papers by this author.
  3. [3]
    Lord Rayleigh and the Love Waves - UC Berkeley Seismology Lab
    Jul 27, 2009 · In a Rayleigh wave, the particles make an elliptical movement against the propagation direction. Hence, their motion is retrograde (see Figure 1) ...
  4. [4]
    Seismic Wave Motions—4 waves animated - IRIS
    Rayleigh Waves—surface waves that move in an elliptical motion, producing both a vertical and horizontal component of motion in the direction of wave ...
  5. [5]
    Seismic Waves and Earth's Interior
    Rayleigh waves are similar to water waves in the ocean (before they "break" at the surf line). As a Rayleigh wave passes, a particle moves in an elliptical ...
  6. [6]
    Rayleigh - Purdue University
    Rayleigh waves have elliptical particle motions and would be recorded on the vertical and radial particle motions away from the source. The equation R(p) = 0 ...
  7. [7]
    Surface Waves - Michigan Technological University
    Rayleigh wave shows rolling motion of a particle below the surface. A Rayleigh wave travels through a medium. Particles are represented by cubes in this model.Missing: discovery | Show results with:discovery
  8. [8]
    [PDF] Chapter 6 Earthquakes - SDSU Geological Sciences
    Because of the vertical motion produced by Rayleigh waves they tend to be the most destructive while the back and forth motion of Love waves can knock buildings ...<|control11|><|separator|>
  9. [9]
    [PDF] Seismic Wave Behavior—Effect on Buildings - IRIS
    Two principal types of surface waves, Love and Rayleigh waves, are generated during an earthquakes. Rayleigh waves cause both vertical and horizontal ground ...
  10. [10]
    Wave Motion in Mechanical Medium - Graduate Program in Acoustics
    Rayleigh waves in an elastic solid are different from surface waves in water in a very important way. In a water wave all particles travel in clockwise circles.
  11. [11]
    [PDF] Chapter 1: - MAE Center
    The dispersive nature of Rayleigh waves enables their use in detecting changes in soil stiffness with depth by measuring the wavelength and velocity of ...
  12. [12]
    Rayleigh-wave Motion- Incorporated Research Institutions for ... - IRIS
    Rayleigh Waves—surface waves that move in an elliptical motion, producing both a vertical and horizontal component of motion in the direction of wave ...Missing: definition | Show results with:definition
  13. [13]
    [PDF] Surface waves - CSUN
    Rayleigh Wave Particle Motion. P and Sv amplitudes are out of phase by /2. This results in elliptical motion (retrograde) ... Rayleigh Wave Phase Velocities.<|separator|>
  14. [14]
    [PDF] Propagation of Rayleigh waves in thin films - CORE
    A Rayleigh wave is a surface wave that propagates within one wavelength from the specimen surface. Since its energy is concentrated very near to the surface, ...
  15. [15]
    Mathematical foundation of elastic wave propagation - SEG Wiki
    Sep 12, 2020 · The equation of wave propagation in elastic solids are derived by using Hooke's law and Newton's second law of motion.
  16. [16]
    [PDF] Waves in an Isotropic Elastic Solid - Columbia University
    Essentially, it says that particles within a moving deformable body have a rate of change of momentum. (mass × acceleration) that is driven by the applied body ...
  17. [17]
    Rayleigh Wave Propagation | SpringerLink
    Lord Rayleigh: “On Waves Propagating Along the Plane Surface of an Elastic Solid,” Proc. Lond. Math. Soc. 17, 4 (1885). Article MATH Google Scholar. Lord ...Missing: original | Show results with:original
  18. [18]
    What are the different types of seismic waves? - Geometrics
    Like body waves, they are characterized by particle motion. Rayleigh waves are characterized by elliptical motion perpendicular to the surface: p-wave ...
  19. [19]
  20. [20]
    Basics of wave propagation — GPG 0.0.1 documentation - GeoSci.xyz
    Rayleigh waves propagate along a free surface or on the boundary between two materials. Particle motion is a retrograde ellipse, and in the same plane as wave ...
  21. [21]
    Rayleigh Waves - an overview | ScienceDirect Topics
    The particle motion in Rayleigh waves follows a retrograde elliptical orbit. Below the earth's surface, the displacement of the particles decreases to zero ...
  22. [22]
    [PDF] History of Seismology - Institute of Geophysics and Planetary Physics
    First scientific observation of strike-slip faulting. 1889 Apr. 18. Near south ... often associated with Rayleigh waves, remained unexplained. The 1950s ...
  23. [23]
    Rayleigh and Lamb Waves - Google Books
    Igor A. Viktorov. Plenum Press, 1967 - 154 pages. Bibliographic information. Title, Rayleigh and Lamb Waves: Physical Theory and Applications. Transl. from ...
  24. [24]
    A HISTORY OF SURFACE ACOUSTIC WAVE DEVICES
    Aug 6, 2025 · In the 1970s, many of these became established in professional systems, and the SAW bandpass filter became a standard component for domestic TV.
  25. [25]
    Exact formula for the horizontal-to-vertical displacement ratio of ...
    Aug 22, 2019 · In this paper, we derive an explicit exact formula for the Rayleigh wave H/V ratio for a compressible orthotropic elastic half-space coated with ...
  26. [26]
    Surface Wave Methods | US EPA
    Apr 18, 2025 · The effective depth of penetration is approximately one-third to one-half of the wavelength of the Rayleigh wave. Rayleigh wave. Figure 3.Missing: decay | Show results with:decay
  27. [27]
    Interaction of a Rayleigh wave with a disbond in a material ...
    Reflection and transmission of Rayleigh surface waves by an interphase between identical or different materials, has been investigated.
  28. [28]
    Rayleigh wave interaction with surface-breaking cracks
    Mar 23, 2007 · This paper investigates Rayleigh wave interaction with simulated, surface breaking cracks using a finite element method
  29. [29]
    [PDF] PEAT8002 - SEISMOLOGY Lecture 5: Surface waves and dispersion
    The motion of the Rayleigh wave is therefore retrograde elliptical, because at the top of the orbit, it is opposite to the direction of propagation. Page 18 ...
  30. [30]
    Rayleigh Wave Velocity - an overview | ScienceDirect Topics
    The theoretical upper limit of crack speed in an elastic, isotropic, and homogeneous material is set at the Rayleigh wave speed.
  31. [31]
    [PDF] Estimation of near‐surface shear‐wave velocity by inversion of ...
    Rayleigh-wave phase velocities are influenced much less by changes in P-wave velocities than by changes in density. A 25% increase in P-wave velocities (Table ...
  32. [32]
    Rayleigh Wave Modeling in Laterally Inhomogeneous Media with ...
    The inversion of Rayleigh wave dispersion curves can yield the shear wave velocity profile which gives the image of subsurface. In most of the existing ...
  33. [33]
    Estimation of dispersion and attenuation of Rayleigh waves in ...
    Oct 22, 2024 · The dispersion and attenuation characteristics of Rayleigh waves in a non-uniform viscoelastic half-space with a covering layer are studied in this paper.
  34. [34]
    [PDF] Attenuation coefficients of Rayleigh and Lg waves
    As they travel through the Earth, seismic waves attenuate due to three general factors: geometrical spreading (GS), intrinsic (anelastic) energy dissipation, ...
  35. [35]
    [PDF] Near-surface Qs and Qp estimations from Rayleigh waves using ...
    This research also showed that acquired data and obtained results can be also used for Rayleigh-wave analysis for evaluating quality factors Q. We focus on ...
  36. [36]
    Seismic Velocities of Rocks and Various Materials - GeoSci.xyz
    Seismic Velocities of Rocks and Various Materials ; Dolomites. 3500 - 6500. 1900 - 3600 ; Granite. 4500 - 6000. 2500 - 3300 ; Basalt. 5000 - 6000. 2800 - 2400.
  37. [37]
    Identification of the Rayleigh surface waves for estimation of ... - NIH
    The regression analysis, with the 95% confidence interval, estimates that the apparent wave speed is c s = 20.41 ± 0.27 m/s at 200 Hz. The results are expressed ...Missing: granite steel
  38. [38]
    Rayleigh Wave Calibration of Acoustic Emission Sensors and ...
    Jul 16, 2019 · Transmitter output in displacement or velocity can be measured directly by laser interferometry. For all ultrasonic transducers and AE sensors ...Missing: speed | Show results with:speed
  39. [39]
    Surface Crack Monitoring by Rayleigh Waves with a Piezoelectric ...
    Feb 28, 2023 · This paper presents a method for measuring surface cracks based on the analysis of Rayleigh waves in the frequency domain.
  40. [40]
    An improved approximation for the Rayleigh wave equation
    Aug 6, 2025 · The apparent wave velocity was determined to be about 4000 m/s. This is clearly slower than the bulk velocity of steel (about 6100 m/s), but ...
  41. [41]
    Perturbational and nonperturbational inversion of Rayleigh-wave ...
    We describe two different methods of inversion: a perturbational method based on finite elements and a nonperturbational method based on the recently developed ...<|control11|><|separator|>
  42. [42]
    A global model of Love and Rayleigh surface wave dispersion and ...
    The new dispersion model, GDM52, can be used to calculate internally consistent global maps of phase and group velocity, as well as local and path-specific ...
  43. [43]
    Love and Rayleigh waves in nonhorizontally layered media
    Mar 3, 2017 · A new finite element technique has been developed to study the propagation of Love and Rayleigh waves across two-dimensional nonhorizontally ...
  44. [44]
    Eric M. Dunham: Supershear Rupture Dynamics
    In addition to generating shear Mach waves, a surface breaking rupture also excites Rayleigh Mach waves. Rayleigh waves are surface waves comprised of ...
  45. [45]
    [PDF] RADIATION PATTERN OF RAYLEIGH WAVES FROM A ... - DTIC
    It is noted that for one certain orien- tation a shallow double-couple source generates no Rayleigh waves. INTRODUCTION. Ben-Menahem (1961) has recently treated ...Missing: epicentral | Show results with:epicentral
  46. [46]
    Seismic yield estimates from Rayleigh-wave source radiation pattern
    Mar 3, 2017 · The evaluated average attenuation coefficients were used, with the station epicentral distance to adjust the corrected spectral amplitudes to ...
  47. [47]
    Earthquake ground motion amplification for surface waves - Bowden
    Dec 15, 2016 · It is well known, for example, that shallow sediments or soils can give rise to amplification and resonances as seismic waves propagate near ...
  48. [48]
    An Exact Solution for Effects of Topography on Free Rayleigh Waves
    Mar 3, 2017 · An exact solution for the effects of topography on Rayleigh wave amplification is presented. The solution is obtained by incorporating ...
  49. [49]
    What was the magnitude? - Earthquake Hazards Program
    The 1906 earthquake's magnitude is estimated at 7.7 (surface wave), 7.9 (moment magnitude), and 8.3 (traditional Richter).
  50. [50]
    A Rayleigh wave back‐projection method applied to the 2011 ...
    Jan 18, 2012 · We study the rupture process of the 2011 Tohoku megathrust by analyzing 384 regional strong-motion records using a novel back-projection method for Rayleigh ...
  51. [51]
    Determination of Near Surface Shear‐Wave Velocities in the Central ...
    Apr 20, 2021 · In this study, we investigate the shallow shear wave velocity structure of the Los Angeles Basin in southern California, using ambient noise correlations.
  52. [52]
    Evaluation of 2-D shear-wave velocity models and VS30at six strong ...
    Apr 18, 2024 · Rayleigh-wave fundamental mode dispersion curve picks across the entire length of the profile (fig. 2.1A) generally coincide with phase ...
  53. [53]
    ARRAY SEISMOLOGY: METHODS AND APPLICATIONS - Rost - 2002
    Dec 5, 2002 · Most array methods use the ability of seismic arrays to measure the vector velocity of an incident wave front, i.e., slowness and back azimuth.
  54. [54]
    Rayleigh phase velocities in Southern California from beamforming ...
    The beamforming processing requires an array of seismometers that are spread spatially within one or two wavelengths to provide an incidence angle and phase ...Missing: slowness | Show results with:slowness
  55. [55]
    [PDF] Rayleigh and Love wave phase velocity maps - Fan-Chi Lin
    Feb 19, 2008 · SUMMARY. We present the results of Rayleigh wave and Love wave phase velocity tomography in the western United States using ambient seismic ...
  56. [56]
    application of Rayleigh-wave ellipticity extracted from noise cross ...
    Abstract. We present a new 3-D seismic model of the western United States crust derived from a joint inversion of Rayleigh-wave phase velocity and elliptic.
  57. [57]
    Attenuation tomography of the upper mantle - Adenis - 2017
    Jul 18, 2017 · In this study, we present QsADR17, a new 3-D shear wave attenuation model, based on a global data set of attenuation measurements for Rayleigh ...
  58. [58]
    Anelastic tomography: a new perspective on upper mantle thermal ...
    Anelastic tomography provides new insights into the thermal structure of the earth's upper mantle. Lateral variations in Q are largest in the upper 300 km, ...
  59. [59]
    Global centroid moment tensor solutions in a heterogeneous earth
    In this paper, we modify the current procedure for the full-waveform inversion of seismic traces for the six moment-tensor parameters.
  60. [60]
    Centroid Moment Tensor - an overview | ScienceDirect Topics
    The CMT method (Dziewonski et al., 1981) is based on the linear relationship that exists between the six independent elements of the zeroth-order moment tensor.
  61. [61]
    An alternative Rayleigh wave excitation method using an ultrasonic ...
    Aug 6, 2023 · The simplest way to generate Rayleigh waves is to couple a longitudinal or shear wave transducer directly to the surface of a test specimen [5].
  62. [62]
    Measurement of Rayleigh Wave Beams Using Angle Beam Wedge ...
    Jun 20, 2017 · The Rayleigh wave generation process using an angle beam wedge transducer is analyzed, and the actual Rayleigh wave sound source distributions ...
  63. [63]
    Laser-induced ultrasonic measurements for the detection and ...
    As a result, Rayleigh waves offer the potential for the detection of surface and subsurface defects. When Rayleigh waves interact with defects, part of the ...
  64. [64]
    [PDF] Laser generation of Rayleigh and Lamb waves for ultrasonic ...
    The main advantage of the technique is its noncontact nature which alleviates the problem of sensor coupling inherent in conventional techniques. However, laser ...Missing: destructive | Show results with:destructive
  65. [65]
    Estimating the Depths of Normal Surface Notches Using Mode ...
    When the Rayleigh surface wave meets discontinuous interfaces, reflection, transmission, and mode conversion waves are generated, which can be used to determine ...Missing: flaws | Show results with:flaws
  66. [66]
    The Characterization of Surface Defects Using Rayleigh Wave ...
    Flaws can be detected up to depths of one or two wavelengths, corresponding to 0.31–6.1 mm over a frequency range of 0.50–10.0 MHz [1].
  67. [67]
    Laser-generated Rayleigh wave for width gauging of subsurface ...
    Aug 10, 2018 · In this study, the widths of subsurface lateral rectangular defects were successfully gauged using laser-generated Rayleigh waves.
  68. [68]
    An alternative Rayleigh wave excitation method using an ultrasonic ...
    This work demonstrates a method that provides an alternative Rayleigh wave excitation method. In this, a conventional ultrasonic phased array transducer is used ...
  69. [69]
    Standard Guide for Electromagnetic Acoustic Transducers (EMATs)
    Jul 6, 2022 · This guide describes an alternate technique in which electromagnetic energy is used to produce acoustic energy inside an electrically conductive or ...
  70. [70]
    [PDF] Application of ultrasonic methods for manufacture of pipelines and ...
    Application of both TOFD and phased array techniques together for zone discrimination of welds, decrease time testing and increase probability of detection for.
  71. [71]
    Focused Rayleigh wave EMAT for characterisation of surface ...
    Rayleigh waves propagate predominantly within one wavelength of the sample surface [2], hence the sensitivity is dependent on the frequency used for inspection.Missing: flaws | Show results with:flaws
  72. [72]
    [PDF] RAIL INSPECTION USING ULTRASONIC SURFACE WAVES
    This section will only give an overview of selected references relevant for the use of surface waves in non-destructive testing (NDT) applications and more ...
  73. [73]
    Fatigue damage evaluation in A36 steel using nonlinear Rayleigh ...
    This research uses nonlinear Rayleigh waves to characterize the damage due to plastic deformation in A36 steel specimens subjected to quasi-static, ...
  74. [74]
    Air-Coupled Detection of Nonlinear Rayleigh Surface Waves to ...
    Aug 6, 2025 · This research presents a new technique for nonlinear Rayleigh surface wave measurements that uses a non-contact, air-coupled ultrasonic ...
  75. [75]
  76. [76]
    Surface Acoustic Wave (SAW) for Chemical Sensing Applications of ...
    Nov 24, 2017 · Surface acoustic wave (SAW) resonators represent some of the most prominent acoustic devices for chemical sensing applications.
  77. [77]
    Design of ST-cut quartz surface acoustic wave chemical sensors
    ST-cut quartz devices present excellent temperature stability. Previous investigations related to polymers on acoustic sensors focused on the selection of a ...
  78. [78]
    Seismic communication - ElephantVoices
    Seismic energy transmits most efficiently between 10-40 Hz - in the same range as the fundamental frequency and 2nd harmonic of an elephant rumble.
  79. [79]
    Evidence for the use of reflected self-generated seismic waves for ...
    Feb 15, 2005 · The results revealed that the mole-rat is able to detect low frequency seismic waves using only its paws, and in most cases the mole-rats ...
  80. [80]
    Evidence for the use of reflected self-generated seismic waves for ...
    The results revealed that the mole-rat is able to detect low frequency seismic waves using only its paws, and in most cases the mole-rats determined ...Missing: foot pads
  81. [81]
    Did Animals Sense Tsunami Was Coming? - National Geographic
    Jan 3, 2005 · Before the tsunami in Sri Lanka, coastal animals seemed to sense something was coming and fled to safety.
  82. [82]
    Nature | Can Animals Predict Disaster? - Destructive Forces
    Jun 10, 2016 · Hours before the deadly Indian Ocean tsunami of 2004, elephants ran for ... zoo animals rushed into their shelters. This mysterious behavior, coupled ...
  83. [83]
    The use of seismic signals by fossorial southern African mammals
    Both geophone and microphone recordings in the substrate suggest that the golden mole is able to detect termite colonies and other prey items solely using ...
  84. [84]
    The golden mole middle ear: A sensor for airborne and substrate ...
    Functionally, golden moles appear to be low-frequency specialists, and it is likely that they detect prey principally through substrate conduction.
  85. [85]
    Seismic Signal Use by Fossorial Mammals1 - Oxford Academic
    Several species of fossorial rodents use drumming behaviors for conspecific signaling. One example of such a rodent is an African mole-rat in the family ...
  86. [86]
    Vibration-induced Behavioral Responses and Response Threshold ...
    Abstract. Despite documented adverse effects, limits for rodent exposure to vibration in the laboratory animal facility have not been established.Missing: avoidance seismic
  87. [87]
    Vibrating Frequency Thresholds in Mice and Rats - ResearchGate
    Aug 9, 2025 · The goal of this study was to determine the range of frequency values, where animals are unable to attenuate vibrations, and therefore may be most susceptible ...
  88. [88]
    (PDF) Review: Can Animals Predict Earthquakes? - ResearchGate
    Aug 10, 2025 · In this study, we reviewed 180 publications regarding abnormal animal behavior before earthquakes and analyze and discuss them.
  89. [89]
    Animals & Earthquake Prediction | U.S. Geological Survey - USGS.gov
    But many animals with more keen senses are able to feel the P wave seconds before the S wave arrives. As for sensing an impending earthquake days or weeks ...
  90. [90]
    Understanding Animal Detection of Precursor Earthquake Sounds
    Aug 31, 2017 · We use recent research to provide an explanation of how animals might detect earthquakes before they occur. While the intrinsic value of such ...
  91. [91]
    Tuning and sensitivity of the human vestibular system to ... - PubMed
    Oct 17, 2008 · In this article, we demonstrate for the first time that the human vestibular system is also extremely sensitive to low-frequency and infrasound ...
  92. [92]
    Impairment of human proprioception by high-frequency cutaneous ...
    Pacinian corpuscles respond to vibration at high frequencies, 80–450 Hz (Sato, 1961; Talbot et al. 1968; LaMotte & Mountcastle, 1975), whereas rapidly adapting ...
  93. [93]
    How Earthquakes Make the Ground Go Boom | Live Science
    Oct 23, 2012 · A new study finds the Earth's surface acts like a speaker for low-frequency vibrations, transmitting an earthquake's infrasonic tumult into the air.
  94. [94]
    [PDF] Thresholds for the perception of fore-and-aft, lateral and vertical ...
    Absolute thresholds for the perception of whole- body vibration by seated persons have been determined over the frequency range from 0.5 to 300 Hz for vertical.
  95. [95]
    Post-earthquake dizziness syndrome following the 2016 Kumamoto ...
    Aug 5, 2021 · The primary earthquake and repetitive aftershocks caused low frequency vibrations (0.1–3.5 Hz), which may cause motion sickness. We observed ...