Lithospheric Waveguides

Physical Mechanism

A waveguide requires two parallel conducting boundaries. In the Earth, this occurs when a resistive layer (low conductivity) is bounded above and below by more conductive layers. The contrast in conductivity creates a potential well for electromagnetic waves, analogous to a dielectric slab waveguide in optics.

The key parameter is the skin depth — the distance an electromagnetic wave penetrates into a conductive medium before its amplitude decays by 1/e:

Skin depth:
  δ = √(2 / (ω μ σ))

Where:
  ω = angular frequency (rad/s)
  μ = magnetic permeability (H/m)
  σ = conductivity (S/m)

For typical crustal rocks at 1 kHz:
  Granite:     σ ≈ 10⁻⁶ S/m  → δ ≈ 500 m
  Sandstone:   σ ≈ 10⁻³ S/m  → δ ≈ 16 m
  Shale:       σ ≈ 10⁻² S/m  → δ ≈ 5 m
  Seawater:    σ ≈ 4 S/m      → δ ≈ 0.25 m
  Salt water:  σ ≈ 0.1 S/m    → δ ≈ 1.6 m

A resistive layer (σ₁) sandwiched between
conductive layers (σ₂, where σ₂ >> σ₁) forms
a waveguide if the layer thickness d satisfies:

  d ≈ δ₂ / π = √(2/(ωμσ₂)) / π

For σ₂ = 10⁻² S/m at 1 kHz:
  d ≈ 5/π ≈ 1.6 m

Geological Structures

Natural lithospheric waveguides form in several geological settings:

• Sedimentary basins: A resistive limestone or sandstone layer between conductive shale layers. Common in petroleum-bearing regions.
• Volcanic formations: Resistive basalt flows between conductive ash or clay layers. Observed in Hawaii and Iceland.
• Permafrost: Frozen ground (very resistive, σ ~ 10⁻⁵ S/m) overlying conductive unfrozen soil. Natural waveguide in Arctic regions.
• Mineralized zones: Conductive ore bodies between resistive host rock. Used in mining geophysics.
• Crustal layering: The Conrad discontinuity (~15 km depth) and Moho discontinuity (~35 km depth) represent conductivity contrasts that could support lithospheric waveguide modes at extremely low frequencies.

Propagation Characteristics

EM waves in lithospheric waveguides exhibit:

• Low attenuation: Waves trapped in the resistive layer experience lower losses than waves propagating through the conductive surroundings. Attenuation constants of 1–10 dB/km at VLF frequencies.
• Modal dispersion: Multiple modes propagate with different phase velocities. Higher-order modes attenuate faster, leaving the fundamental mode dominant at long range.
• Anisotropy: Geological layering is rarely perfectly horizontal. Dip angle affects coupling efficiency and mode structure.
• Frequency selectivity: The waveguide acts as a bandpass filter, favoring frequencies determined by layer thickness and conductivity contrasts.

Experimental Evidence

Several experiments have demonstrated lithospheric waveguide propagation:

• Mine communication (South Africa): VLF signals at 10–30 kHz were observed to propagate along a quartzite layer between shale formations over distances of 10+ km.
• Arctic propagation: ELF signals at 100–300 Hz traveled through permafrost layers with attenuation as low as 0.1 dB/km, far below surface propagation losses.
• Hawaii volcanic rock: VLF signals propagated through basalt layers over 5 km with 15 dB total attenuation, consistent with a basalt waveguide model.
• Saharan experiments: Desert sand layers (very resistive, σ ~ 10⁻⁵ S/m) over conductive bedrock showed waveguide-like propagation at 1–10 kHz.

Applications

Understanding lithospheric waveguides has practical applications:

• Underground communication: Military and mining applications use VLF/ELF waveguide modes for through-rock communication where conventional radio cannot penetrate.
• Geophysical exploration: Measuring the waveguide properties of rock layers helps map subsurface geology for oil/gas/mineral exploration.
• Seismic precursor studies: Stress changes before earthquakes alter rock conductivity, potentially modifying waveguide propagation. Monitoring these changes could contribute to earthquake prediction research.
• Nuclear test detection: Underground nuclear tests generate electromagnetic pulses that propagate through lithospheric waveguides. Monitoring these signals helps verify test ban treaties.