Gibbs' Phenomenon: Fourier Truncation, Seismic Wavelet Ringing, and Band-Limited Sidelobes

Gibbs' phenomenon is the persistent overshoot and ringing that appears near a sharp discontinuity when a signal is reconstructed from a truncated or band-limited Fourier series. When a square edge or step is rebuilt from a finite sum of sine and cosine components, the reconstruction does not settle smoothly onto the true value at the jump. Instead it overshoots by a fixed proportion of the step height, roughly 8.95 percent, then oscillates with decaying ripples on either side of the edge. The counterintuitive part, first explained by J. Willard Gibbs in 1899, is that adding more Fourier terms does not reduce the height of that overshoot. More terms squeeze the ripples closer to the discontinuity and narrow them, but the peak overshoot stubbornly stays near nine percent no matter how many frequencies are included. The effect matters in reflection seismology because recorded seismic data is inherently band-limited: the earth, the source, and the recording system pass only a finite window of frequencies, typically something like 8 to 80 Hz for conventional Vibroseis surveys. A sharp acoustic-impedance contrast, such as the top of a tightly cemented carbonate or the base of a gas sand, behaves like a near-discontinuity in the reflectivity series. When that reflectivity is convolved with a band-limited seismic wavelet, the missing high and low frequencies produce exactly the Gibbs-style sidelobes: bright trough-peak-trough ringing that flanks the true reflection. Interpreters can mistake these processing sidelobes for real thin beds or extra geological interfaces, so understanding the phenomenon is central to honest interpretation. The same ringing degrades deconvolution results, complicates thin-bed tuning analysis, and forces the use of tapered windows during spectral processing. Managing Gibbs' phenomenon is therefore not an abstract mathematics problem but a daily concern in Western Canadian Sedimentary Basin seismic work, from Montney resource-play mapping to Duvernay amplitude analysis, where a few metres of apparent thickness can change a well's economics.

Key Takeaways

  • Fixed Nine Percent Overshoot: The defining feature is that the overshoot near a jump converges to about 8.95 percent of the step height and never disappears, regardless of how many Fourier terms are summed. Adding frequencies narrows the ripple and pushes it toward the discontinuity but does not lower the peak. This non-uniform convergence is why simply acquiring more bandwidth cannot fully remove the artifact from a band-limited seismic trace.
  • Band-Limiting Is the Root Cause: Seismic data passes only a finite frequency window, often 8 to 80 Hz, set by the source, earth attenuation, and geophones. Any sharp impedance boundary contains frequencies outside that window, and truncating them produces the ringing. The narrower the recorded bandwidth, the more pronounced the sidelobes, which is why low-frequency-rich sources and broadband acquisition reduce, but never eliminate, Gibbs ringing.
  • Sidelobes Mimic Thin Beds: The trough-peak-trough ripples flanking a strong reflection can be misread as genuine geological interfaces. In WCSB thin-bed plays this is dangerous: an interpreter mapping a 4 m Cardium sand can pick a Gibbs sidelobe as a second sand and overstate net pay. Wavelet-phase knowledge and synthetic ties are essential to separate real reflectivity from processing ringing.
  • Tapered Windows Suppress Ripple: Applying a smooth taper such as a Hann, Hamming, or cosine window to the spectrum before inverse transform trades a small loss of resolution for a large reduction in overshoot. Rather than a hard truncation at the band edge, the taper rolls the amplitude gently to zero, which softens the discontinuity in the frequency domain and calms the time-domain ringing.
  • Deconvolution and Tuning Impact: Because deconvolution tries to compress the wavelet and broaden the spectrum, it can amplify Gibbs ringing when the band edges are left sharp. Aggressive whitening spikes the residual sidelobes. Thin-bed tuning analysis, where two close reflections interfere, is also corrupted because the tuning wavelet already carries ringing, biasing net-thickness estimates unless the wavelet is carefully characterized.

How Band-Limited Wavelets Generate Sidelobe Ringing

A zero-phase band-limited wavelet, such as a Ricker or Ormsby wavelet, is the time-domain fingerprint of the recorded frequency band. Because the spectrum is cut off sharply at both ends, the wavelet is not a clean spike but a central peak flanked by symmetric side troughs, the direct expression of Gibbs' phenomenon. An Ormsby wavelet defined by corner frequencies of 5-10-70-80 Hz shows visible first sidelobes about 20 to 25 percent of the main-peak amplitude. When this wavelet convolves with a single sharp reflector, those sidelobes print into the data as false events roughly one half-period above and below the true reflection, spaced by the dominant period of the recorded band.

Windowing and Broadband Acquisition as Mitigations

The practical defence is to shape the amplitude spectrum smoothly rather than clip it. Processors apply cosine tapers of 5 to 10 Hz at each band edge, or full Hann and Hamming windows, so the transition to zero amplitude is gradual. This lowers sidelobe energy by more than half at the cost of slightly reduced vertical resolution. On the acquisition side, low-dwell Vibroseis sweeps that push energy down to 3 Hz and broadband sensors widen the recorded band, moving the ringing period shorter and the sidelobes closer to the main event. Neither approach removes Gibbs overshoot entirely, so interpreters still validate strong reflectors against well synthetics before mapping thin beds.

Fast Facts

Josiah Willard Gibbs, the Yale physicist better known for founding chemical thermodynamics, published his explanation of the overshoot in the journal Nature in 1899, correcting an earlier note by Albert Michelson, who had built a mechanical harmonic analyzer and was baffled when it kept producing spurious spikes at the ends of a square wave. The effect had actually been described decades earlier by the English mathematician Henry Wilbraham in 1848, but it was Gibbs' name that stuck. The overshoot constant, close to 8.9490 percent, ties back to the Wilbraham-Gibbs integral of the sinc function.

Gibbs ringing is inseparable from the seismic wavelet, whose sidelobes are its direct time-domain signature, and from deconvolution, which can sharpen band edges and worsen the overshoot when whitening is aggressive. It shapes seismic resolution because the ringing period sets the practical limit on separating close reflectors, and it interacts with tuning thickness, where interfering sidelobes bias net-pay estimates. Each concept depends on the same band-limited Fourier reconstruction that creates the phenomenon.

Duvernay Amplitude Pitfall Near Fox Creek

On a 3D survey shot over a Duvernay shale target near Fox Creek, Alberta, a processing team delivered a migrated volume with a bright, sharply band-limited response at the top of a Nisku carbonate at roughly 3,300 m. A junior interpreter mapped a strong trough about 12 ms below the carbonate top as a separate porous interval and flagged it as a secondary target, a call that would have justified a 4,000 m horizontal well costing about 8 million CAD. A reprocessing check with a 4 Hz to 90 Hz cosine-tapered spectrum showed the trough collapsing toward the main event.

The apparent second reflector was a Gibbs sidelobe, not rock. Tying the volume to a nearby control well confirmed no porosity at that level, and the operator cancelled the standalone target, saving the drilling capital and avoiding a dry completion. The lesson repeated across the play: validate every bright band-limited event against a well synthetic before committing metres of lateral.