Which Waves Can Travel Through Both Solids And Liquids

You're sitting in a geology lecture, or maybe watching a documentary about earthquakes, and someone says: "P-waves travel through everything. Think about it: " You nod. But then you wonder — wait, why? S-waves stop at the liquid outer core.It sounds right. What is it about a liquid that lets one wave through and blocks the other cold?

That question doesn't just matter for seismology. Day to day, it shows up in medical ultrasound, in industrial testing, in how we explore for oil and gas. The answer comes down to something simple: how the particles in a material push and pull on each other.

What Is a Wave, Really?

Before we sort solids from liquids, let's get clear on what a wave actually is in a solid or liquid. Here's the thing — it's not the water moving across the ocean. Day to day, it's energy moving through* a medium by making particles oscillate around their equilibrium positions. The particles don't travel with the wave — they just bump their neighbors.

Two main ways that bumping happens:

Compression and expansion — particles squeeze together, then spread apart, parallel to the direction the wave travels. Think of a slinky pushed back and forth. This is a longitudinal* wave.

Shear motion — particles move side to side, perpendicular to the wave's travel direction. Shake that slinky up and down. This is a transverse* wave.

Solids support both. Liquids? Only the first.

The Particle-Level Difference

In a solid, atoms or molecules are locked in a lattice. They have shear stiffness* — resistance to being slid past one another. Push the top layer sideways, and the bonds pull it back. That restoring force lets a shear wave propagate.

In a liquid, there's no lattice. On top of that, molecules slide freely. Push a layer sideways, and it just... stays pushed. Also, no restoring force. No shear wave. But squeeze a liquid — compress it — and the molecules push back hard. That's why sound (a compressional wave) travels through water just fine.

Gases work the same way, only softer. They compress easily, so sound travels slower.

Why It Matters: The Earth Tells Its Own Story

This isn't textbook trivia. It's how we know the Earth has a liquid outer core.

When an earthquake ruptures, it radiates P-waves (Primary, compressional) and S-waves (Secondary, shear). Seismometers around the world record their arrival times. P-waves arrive first — they're faster. S-waves follow. But here's the kicker: there's a shadow zone* for S-waves. Between about 103° and 143° from the epicenter, no S-waves show up at all.

P-waves do show up there — but they arrive later than expected, and their paths bend sharply at the core-mantle boundary.

The only explanation: the outer core is liquid. S-waves hit it and die. That's why p-waves transmit through, refracting like light through a lens. Now, the inner core? Solid again — S-waves reappear as converted phases (PKJKP, if you're into the notation).

That discovery — made in the early 1900s by people like Richard Oldham and Beno Gutenberg — rewrote our picture of the planet. All from watching which waves made it through and which didn't.

Beyond Earth: Medical and Industrial Uses

Same physics, different scale. Which means medical ultrasound uses high-frequency P-waves (2–18 MHz typically). Also, they travel through soft tissue — mostly water — and reflect at boundaries where acoustic impedance changes. That's how you see a fetus, a gallstone, a thickened heart wall.

But ultrasound can't* image through bone or air-filled lungs well. Bone is solid — it supports shear waves, and the impedance mismatch reflects almost everything. Air compresses too easily — the wave just dies.

In non-destructive testing, technicians send ultrasonic P-waves into welds, pipelines, turbine blades. Cracks, voids, inclusions — they all scatter or reflect the wave. Some advanced setups even generate shear waves inside* the solid part by angling the transducer (mode conversion at the surface). That lets them inspect at angles a straight beam can't reach.

Oil exploration? But seismic reflection surveys. Vibroseis trucks or air guns generate low-frequency P-waves that penetrate kilometers of rock. The reflections map sedimentary layers, faults, salt domes — the traps where hydrocarbons accumulate. Shear-wave surveys (using specialized sources) add extra detail: they're sensitive to fracture orientation, fluid saturation, things P-waves alone miss.

How It Works: The Physics Without the Jargon

Let's walk through the mechanics. No tensor notation. Just the logic.

Compressional Waves (P-waves)

Imagine a row of balls connected by springs. That's why push the first ball. But it compresses the spring to its neighbor. The disturbance travels. Day to day, that neighbor moves, compresses the next spring, and so on. The balls oscillate along* the line.

In a continuum, this is described by the bulk modulus (K) — resistance to volume change — and density (ρ). Wave speed:

vₚ = √((K + 4/3 μ) / ρ)

μ is the shear modulus — resistance to shape change. In a liquid or gas, μ = 0. So the formula simplifies to:

vₚ = √(K / ρ)

That's why sound speed in water is ~1,480 m/s. In steel, it's ~5,960 m/s — steel is stiffer and denser, but stiffness wins.

Shear Waves (S-waves)

Now shake the first ball sideways*. Consider this: the spring stretches diagonally. It pulls the neighbor sideways. That neighbor pulls the next. The disturbance travels, but the motion is perpendicular* to the travel direction.

This only works if the springs resist sideways stretching — i.e., if there's shear stiffness.

vₛ = √(μ / ρ)

If μ = 0 (liquid, gas), vₛ = 0. Because of that, no wave. Period.

Surface Waves: The Hybrid Case

At a free surface — the ground, the top of a weld — things get interesting. Two main types:

Rayleigh waves — particle motion is elliptical, retrograde (counter to wave travel) at the surface. They're a mix of P and S motion, confined near the surface. Speed is ~0.92 vₛ for a Poisson solid. They carry most of the shaking you feel in an earthquake.

Love waves — purely horizontal shear motion, trapped in a low-velocity layer over a half-space. They need a velocity contrast. They don't exist in a uniform half-space.

Both are guided* waves — they only exist because of a boundary. In a liquid, you don't get Love waves (no shear). You can get a Rayleigh-like wave at a liquid-solid interface (Scholte wave), but it's different — particle motion is elliptical in the vertical plane, and it decays into both media.

Mode Conversion: When Waves Hit a Boundary

This

Mode Conversion: When Waves Hit a Boundary

It's where things get physically rich — and where the simple ball-and-spring picture needs a reality check.

When a P-wave traveling through rock hits a boundary with a different rock type, it doesn't just reflect as a P-wave or transmit as a P-wave. Part of the energy reflects as a P-wave, part as an S-wave. Part transmits as a P-wave, part as an S-wave. Here's the thing — it can split*. Same goes for an incoming S-wave — it can generate reflected and transmitted P-waves and S-waves.

This is mode conversion, and it happens because the boundary has to satisfy both* stress and displacement continuity simultaneously. The incoming wave brings a certain stress state. The outgoing waves — whatever types they are — have to add up to match that stress on one side, and to produce the right displacement on the other. The only way to satisfy both conditions is to generate multiple wave types.

Think of it like this: you're pushing a ball (P-wave) into a wall that's softer on the other side. The wall doesn't just push back the same way. It also shears* — it wants to slide sideways at the interface. That sideways component is the converted S-wave.

Snell's Law — But for Every Wave Type

At every boundary, the horizontal component of wave motion has to be continuous. If it weren't, the two sides of the boundary would tear apart or overlap — and rock doesn't do that easily. This single constraint links the angles of all the generated waves.

For each wave type, the horizontal slowness (the horizontal component of the slowness vector) is the same:

sin(θ₁) / v₁ = sin(θ₂) / v₂

But here's the key: θ₁ and v₁ refer to the incident* wave, while θ₂ and v₂ can refer to any generated wave — P or S, in either medium. Because of that, since P-waves are usually faster than S-waves in the same material, a single P-wave incident on a boundary produces reflected and transmitted P- and S-waves, each at a different* angle. The faster wave goes at a shallower angle.

We're talking about why seismic sections sometimes show "head waves" — critically refracted P-waves that travel along the boundary and re-enter the upper layer. They're a direct consequence of this angular splitting.

Reflection and Transmission Coefficients

How much energy goes where? That's governed by the impedance contrast at the boundary.

Acoustic impedance for P-waves is:

Z = ρ · vₚ

A big impedance jump (say, sandstone to shale, or rock to fluid) means strong reflection. A small jump means most energy passes through. The reflection coefficient for normal incidence simplifies to:

R = (Z₂ − Z₁) / (Z₂ + Z₁)

For non-normal incidence — the real-world case — the math gets more involved. Which means the Zoeppritz equations describe exactly how P- and S-wave energy partitions at an arbitrary angle. They're a set of four coupled equations, and they're the reason AVO (Amplitude Versus Offset) analysis works: by measuring reflected amplitudes at different angles, geophysicists can infer rock properties — lithology, fluid content, porosity — before ever drilling a well.

Why Shear-Wave Data Matters

Here's the practical payoff. P-wave impedance alone can't distinguish a gas-saturated rock from a water-saturated rock of the same porosity and clay content — the P-wave velocities might overlap. That said, a gas-filled sandstone has a lower* shear-wave velocity than a water-filled one of the same structure. But shear waves don't care about the fluid's bulk modulus (remember, μ = 0 in fluids, so S-waves don't propagate through them). That difference — captured in the Vp/Vs ratio — is one of the most reliable indicators of hydrocarbon presence.

This is why marine seismic surveys sometimes use ocean-bottom seismometers (OBS) that record both P- and S-waves. The extra information is worth the cost.

If you found this helpful, you might also enjoy orbitals that have the same energy are called or what are corresponding angles in geometry.

---

Attenuation: The Fading of Seismic Energy

Waves don't travel forever. As they move through rock, they lose energy — and not just by spreading out (geometric spreading, which follows predictable inverse-square or inverse-distance laws). There's an intrinsic* loss, called anelastic attenuation, caused by internal friction and thermoelastic effects at the grain scale.

We quantify this with the quality factor, Q:

Q = 2π × (Energy stored

Attenuation: The Fading of Seismic Energy

The quality factor, Q, is defined as

[ Q = 2\pi \frac{\text{Energy stored per cycle}}{\text{Energy lost per cycle}} ]

In practice, Q is measured by fitting an exponential decay to observed amplitude versus travel‑time data. g., 150–200) indicates a relatively “loss‑free” medium such as crystalline basement, while low Q (e.In real terms, a high Q (e. g., 20–40) signals strong attenuation in unconsolidated sediments or high‑porosity sandstones.

Attenuation is frequency‑dependent: higher‑frequency components lose energy more rapidly, which is why seismic waveforms become progressively lower‑frequency with distance. This dispersion effect can be modeled with the Klebene‑Ahrens or Mavroeidis‑Mavko attenuation equations, which incorporate both intrinsic mechanisms (grain‑boundary friction, squirt‑flow) and scattering from heterogeneities.

Why it matters for interpretation

• Amplitude‑vs‑offset (AVO) inversion must correct for amplitude loss that mimics a fluid‑induced bright spot but is actually just attenuation.
• Time‑lapse (4‑D) monitoring relies on repeatable amplitude behavior; unaccounted attenuation can masquerade as pressure‑induced velocity changes.
• Velocity‑anisotropy studies use attenuation anisotropy to infer grain‑scale preferred orientation, adding a structural dimension beyond pure velocity mapping.

Modern acquisition designs mitigate these issues: broadband sources, multi‑component receivers, and carefully timed repeat surveys all help isolate true amplitude changes from purely dissipative effects.

---

Scattering and Diffraction: Waves Around Obstacles

When the wavelength of a seismic wave is comparable to the size of a geological feature—such as a fault zone, boulder, or thin sand lens—the wavefront is perturbed in ways that cannot be captured by simple ray theory. Scattering redistributes energy into many directions, creating a “hazy” background that can obscure primary reflections.

Diffraction, a subset of scattering, occurs when a wave encounters a sharp edge or a small aperture and generates secondary wavelets that propagate outward. In practice, diffractions are often observed as hyperbolic events on seismic sections, especially when imaging across faults or around salt bodies.

Understanding these phenomena is crucial for two reasons:

1. Imaging challenges – conventional migration algorithms assume smooth wave propagation; unmodeled diffractions can cause migration artifacts or false “bright spots.”
2. Resolution enhancement – by deliberately preserving diffractions, interpreters can extract high‑frequency information that reveals small‑scale structures (e.g., channel edges, thin beds) that would otherwise be invisible.

Advanced processing sequences—such as diffraction separation, frequency‑wavenumber (FK) filtering, and adaptive subtraction—are employed to isolate or suppress these events depending on the geological context.

---

Practical Workflow: From Acquisition to Interpretation

A modern seismic interpretation project typically follows these steps:

1. Velocity model building – using travel‑time tomography that incorporates both P‑ and S‑wave arrivals, anisotropy parameters, and attenuation constraints.
2. Amplitude‑versus‑offset (AVO) and AVO‑AVO‑Angle (AVA) analysis – extracting fluid‑sensitive attributes from angle‑dependent reflection coefficients derived from the Zoeppritz equations.
3. Attribute conditioning – applying statistical normalization, trend removal, and noise‑reduction techniques to isolate true hydrocarbon‑related anomalies.
4. Rock‑physics tying – linking observed attribute patterns to petrophysical variables (porosity, saturation, lithology) through calibrated rock‑physics models.
5. Uncertainty quantification – using Monte‑Carlo simulations or Bayesian inversion to assess the confidence of each interpreted horizon or anomaly.

The integration of P‑wave and S‑wave datasets, along with attenuation and diffraction information, yields a richer, more strong picture of the subsurface than any single attribute could provide.

---

Conclusion

Seismic waves are far more than simple acoustic pulses; they are a suite of coupled mechanical phenomena that reveal the hidden architecture of the Earth. By mastering the fundamentals of wave propagation—reflection, transmission, mode conversion, and the governing impedance contrasts—geophysicists can extract quantitative information about rock composition, fluid content, and structural geometry. Attenuation, scattering, and diffraction add layers of complexity, but they also furnish valuable diagnostics when interpreted correctly.

When these concepts are woven together with modern acquisition techniques, reliable velocity modeling, and sophisticated attribute analysis, seismic data become a

When these concepts are woven together with modern acquisition techniques, solid velocity modeling, and sophisticated attribute analysis, seismic data become a multidimensional diagnostic tool that can be tuned to highlight both the macro‑scale architecture of a basin and the micro‑scale heterogeneity of its constituent units.

A Holistic Interpretation Paradigm

1. Joint P‑S Wave Analysis – By jointly inverting P‑ and S‑wave travel times, interpreters obtain a unified velocity model that respects the elastic coupling between the two modes. This dual‑mode constraint reduces non‑uniqueness and yields more reliable depth estimates for structurally complex regions.

2. Attenuation‑Weighted Migration – Incorporating amplitude‑loss information into the migration aperture allows the imaging of low‑energy reflections that would otherwise be suppressed, revealing subtle stratigraphic interfaces and fault planes that are critical for basin evolution studies.

3. Diffraction Imaging of Small‑Scale Features – Targeted diffraction separation workflows isolate high‑frequency energy associated with thin beds, channel edges, and fracture swarms. When these diffractions are migrated with velocity models calibrated to local anisotropy, they provide “acoustic lenses” that sharpen the resolution of otherwise blurred horizons.

4. Attribute‑Driven Risk Assessment – The calibrated AVO/AVA attributes, when combined with attenuation and diffraction metrics, feed directly into probabilistic risk maps. Bayesian frameworks can then assign posterior probabilities to each prospect, enabling decision‑makers to prioritize drilling locations based on a quantitative measure of confidence.

5. Dynamic Updating Loop – As new wells are drilled and logged, the resulting velocity, density, and attenuation measurements are fed back into the geophysical models. This iterative workflow continuously refines the subsurface description, ensuring that interpretations remain consistent with observed data throughout the field life.

Implications for Exploration and Production

• Improved Success Rates – By extracting a richer set of elastic parameters, operators can better predict the presence and quality of hydrocarbons, reducing the incidence of dry wells.
• Optimized Reservoir Management – High‑resolution imaging of fractures and subtle stratigraphic traps supports more accurate placement of horizontal wells and hydraulic fracturing stages, maximizing sweep efficiency.
• Reduced Exploration Footprint – Accurate, high‑resolution models diminish the need for extensive wildcat drilling, translating into lower environmental impact and lower capital expenditure.
• Cross‑Disciplinary Integration – The workflow described above naturally bridges geophysics, geology, petrophysics, and reservoir engineering, fostering a shared language and set of objectives across multidisciplinary teams.

Future Directions

The next generation of seismic interpretation will likely be shaped by three converging trends:

1. Full‑Waveform Inversion (FWI) at Field Scale – Advances in computational power and data acquisition (e.g., dense 4‑D arrays, ocean‑bottom nodes) are making it feasible to run FWI pipelines that honor both P‑ and S‑wave physics, attenuation, and scattering simultaneously.

2. Machine‑Learning‑Assisted Attribute Extraction – Deep‑learning models trained on synthetic and field datasets can automatically detect subtle diffraction signatures, classify attenuation patterns, and suggest optimal velocity model perturbations, accelerating the interpretation cycle.

3. Integrated Quantum‑Sensing Concepts – Emerging research into quantum‑enhanced seismic sensors promises higher signal‑to‑noise ratios and the ability to measure minute velocity changes, potentially unlocking new levels of resolution for shallow, high‑frequency targets.

Closing Thought

In essence, seismic waves serve as the Earth’s own acoustic fingerprint, encoding a wealth of information about the materials they traverse. By mastering the fundamentals of wave behavior—reflection, transmission, mode conversion, attenuation, scattering, and diffraction—geophysicists can decode that fingerprint with ever‑greater fidelity. When coupled with rigorous modeling, attribute‑driven analytics, and an iterative learning loop, seismic interpretation transforms from a descriptive art into a predictive science, empowering the energy industry to locate resources more responsibly, extract them more efficiently, and ultimately steward the subsurface for the benefit of society.