Characterising ice slab in firn using seismic full waveform inversion, a sensitivity study

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted s...

Full description

Bibliographic Details
Published in:Journal of Glaciology
Main Authors: Pearce, Emma, Booth, Adam, Rost, Sebastien, Sava, Paul, Konuk, Tuğrul, Brisbourne, Alex, Hubbard, Bryn, Jones, Ian
Format: Article in Journal/Newspaper
Language:English
Published: International Glaciological Society/Cambridge University Press 2023
Subjects:
Online Access:http://nora.nerc.ac.uk/id/eprint/533135/
https://nora.nerc.ac.uk/id/eprint/533135/1/characterising-ice-slabs-in-firn-using-seismic-full-waveform-inversion-a-sensitivity-study.pdf
https://www.cambridge.org/core/journals/journal-of-glaciology/article/characterising-ice-slabs-in-firn-using-seismic-full-waveform-inversion-a-sensitivity-study/B2372F01CC95E8417886227BC7AB9ECC
Description
Summary:The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey.